Electrolytic Capacitor Reliability Analysis using COMSOL MULTIPHYSICS
Explore the source record for details and available documents.
Negoita, Mihai [Sandia National Laboratories (SNL-↗
SEARCH · Search NASA
Results for “COMSOL Multiphysics”
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.
Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”
This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.
ATS↗
Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)
In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.
divertor↗
Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”
This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.
ATS↗
Design and Analysis of a Mutual Inductance Level Sensor for Liquid Metals
Here, this article describes the design and analysis of an electromagnetic level sensor for use in high-temperature liquid metal systems. The mutual inductance level sensor (MILS) described in this work was fabricated using two single-conductor mineral insulated cables wrapped in a bifilar fashion around a stainless steel tube core and was housed in an isolating thimble that preserved the pressure boundary of the test vessel. Two sensor variations were fabricated that differ only in active length, 1016 and 1778 mm. Experimental data were collected using the 1016-mm sensor (MILS-MKII-040) that demonstrated a sensitivity of 9.2 μ V/mm in a room temperature testing stand that used solid aluminum as a surrogate for liquid metal. Experimental data were collected using the 1778-mm sensor (MILS-MKII-070) that demonstrated a sensitivity of 6.9 μ V/mm in the high-temperature (300 ° C) liquid sodium environment at the mechanisms engineering test loop (METL) of Argonne National Laboratory. The sensor performance was found to be repeatable over the course of several months, with roughly ±1% deviation from nominal output. Finite element models were developed in COMSOL Multiphysics that fully describe each test setup, and the models were validated using experimental data. The validated COMSOL models were used to perform an array of analyses that examined the performance of the sensor in differing environments. Maximizing the coil diameter inside the isolating thimble was found to maximize the signal and sensitivity of the sensor. An optimal operating frequency was found near 1000 Hz using both experimental data and COMSOL. The influence of a metallic thimble surrounding the sensor and a metallic sensor core was quantified and found to be negligible at the optimal operating frequency. The sensitivity of the sensor was quantified when monitoring the level of additional liquid metals. These include lead, lead-bismuth eutectic (LBE), sodium-potassium alloy (NaK), and lithium (in addition to sodium). The sensitivities were quantified using liquid metal properties at 350 ° C and 650 ° C. The geometry of the test stand model, all material properties used in the model, and the results are presented in a manner that allows the reader can replicate the model and perform additional analyses.
COMSOL↗
Precise magnetic field mapping of the EMPHATIC Phase 1 magnet with COMSOL
A compact Halbach array magnet is used to measure the momentum of the secondary particles in EMPHATIC(Experiment to Measure the Production ofHadrons At a Test beam InChicagoland). Hall probe data was taken for the central cylindrical bore of the magnet and a field map was constructed. COMSOL Multiphysics® Software is used for modeling the magnet and constructing the corresponding magnetic field map. Here, we present a fitting approach where the hall probe data is used to determine a 1mm-spacing map of the entire volume of the magnet using COMSOL. The new map will allow for linear interpolation within the volume, and expand the map to outside the measurement volume, thus increasing the acceptance and precision of EMPHATIC's tracking system.
47 OTHER INSTRUMENTATION↗
Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field
Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.
Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗
Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)
This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.
08 HYDROGEN↗
Modeling of Vertical Motor-driven Pump for Simulation of a Fault Signature \\ for Condition Monitoring
As part of the ongoing effort to transition from preventive maintenance strategies to condition-based maintenance strategies in nuclear power plants, there is significant reliance on using machine learning techniques. To develop a robust machine learning model that can diagnose all the fault modes of a vertical motor-driven pump, data capturing the unique signature of each fault mode is required. In practice, it is difficult to collect or capture data that captures all the fault modes from a single plant site. So to address this situation, a computational model of a vertical motor-driven pump is developed using the multipurpose finite element software COMSOL Multiphysics. The developed model is used to generate simulated data under normal operation and is compared with the vibration data collected using vibration sensors. Once the simulation model is verified under normal operating condition, simulated data for the fault mode for which minimal or no evidence is available in historical plant process data is developed. This simulated data is used to develop fault signatures to achieve robust predictive models. This paper presents modeling details and verification of the model that can used to generate data for fault modes that are not available at a plant site for condition monitoring purpose.
46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
Enhancing EnergyPlus capabilities to model dynamic building envelopes using python plugin
Nearly half of the energy consumption in the United States is related to buildings, resulting in an urgent need to develop innovative technologies to improve building energy efficiency. Dynamic building envelopes, comprising switchable insulation and thermal energy storage materials, have been proposed recently as a promising solution to reduce buildings' heating and cooling loads by thermally coupling the indoor environment with the ambient environment when beneficial while decoupling them when outdoor conditions are not favorable. Although various related technologies are still underway, the whole-building energy modeling tools, like EnergyPlus, do not have the capability to simulate the transient and dynamic nature of dynamic envelope materials and components to accurately predict their impact on building energy use. The objective of this study is to formulate a method in EnergyPlus simulation engine to model multilayer envelopes, comprising dynamic building materials with variable thermophysical properties, and discuss the changes made to the program using a Python plugin. Furthermore, the thermal performance of the dynamic envelopes using the proposed method is compared and verified with the results from a well-established commercial code, COMSOL Multiphysics. A parametric assessment is also conducted to evaluate the energy efficiency benefits of dynamic envelopes in a single-family residential building, demonstrating total annual energy savings up to 11.6 %, when a dynamic envelope operates alone, and up to 18.2 % when it is combined with a thin layer of phase change material as a thermal storage medium. Finally, a United States wide energy efficiency assessment is presented to showcase the geographical spread of the energy savings. The method designed and implemented in this study provides the researchers with the ability to implement their dynamic insulation methods in EnergyPlus and evaluate the whole building energy impact.
25 ENERGY STORAGE↗
Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether
Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles
In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.
Adsorption-desorption process↗
Coupling thermal energy storage with a thermally anisotropic building envelope for building demand-side management across various US climate conditions
Here, the thermally anisotropic building envelope (TABE) is a novel active building envelope that enhances energy efficiency and thermal comfort in buildings by transferring heat and cold between building envelopes and hydronic loops. When coupled with thermal energy storage (TES) units, the TABE + TES enables the storage of both heat and cold energy captured by the TABE roof or exterior walls. This stored energy can be later released by the TABE floor for indoor heating and cooling, benefiting both the grid and the end user. This paper evaluates the merits of TABE + TES for building demand-side management across various US climate conditions, focusing on peak load shaving, annual energy savings, and cost savings under time-of-use (TOU) electric rate schedules. Simulations were conducted by integrating time-of-day–informed, rule-based control strategies in MATLAB, TABE components and TES units in COMSOL Multiphysics, and whole-building energy analysis in EnergyPlus. A case study using the US Department of Energy’s prototype single-family detached house model in Birmingham, Alabama; Los Angeles, California; Oak Ridge, Tennessee; and Denver, Colorado, showed that the TABE + TES system achieved (1) 70 % peak load shaving in Los Angeles and Denver and 20 % in Birmingham and Oak Ridge; (2) significant peak electricity savings of 351–497 kWh, reducing peak energy consumption by 38 %–78 %; and (3) annual heating cost savings of 0.79 $\$$/m2–1.17 $\$$/m 2 and cooling cost savings of 0.60 $\$$/m 2 –1.17 $\$$/m 2 using a normal utility rate or low-TOU rate. The benefits of employing the TABE + TES system are even more significant under high TOU rates.
25 ENERGY STORAGE↗
Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold
Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Improved internal short circuit models for thermal runaway simulations in lithium-ion batteries
Thermal runaway (TR) modeling is one of the primary tools that can be used to overcome challenges associated with lithium-ion battery (LIB) safety. Among all LIB accidents that have occurred over the past decade, Internal Short Circuit (ISC) remains the most common trigger mechanism. Many available models in the literature either use a simplified approach to simulate ISC or completely ignore its contribution. The aim of this study is to understand the nature of the heat released for different types of ISC scenarios, including aluminum-anode, anode–cathode, and copper-cathode ISC. We study ISC behavior using a coupled electrochemical–thermal model with an integrated TR chemical kinetics solver built in the COMSOL Multiphysics framework. The time duration of heat release and the magnitude of the peak ISC current are studied as functions of parameters such as the size of the penetrating filament and the capacity of the cell. The numerical results are used to build an empirical model validated against the published experimental TR propagation data. Our model can be successfully used as a viable low-cost substitute in lower order (lumped) TR simulations to enable TR prevention and mitigation.
Singh, Bakhshish Preet (ORCID:0000000264751992)↗
Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas
Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .
24 POWER TRANSMISSION AND DISTRIBUTION↗
Novel Full-Ceramic Multi-Tubular Membrane Systems for Pre-Combustion CO2 Capture with Simultaneous H2 Production: Fabrication, Performance Testing, and 3D CFD Modeling
Inorganic membranes show promise for application in pre-combustion CO2 capture with simultaneous H2 production. State-of-the-art systems for use under high temperature and pressure conditions consist of multiple membrane tube bundles prepared in a "candle-filter" configuration, in which the membrane tubes are open at one end and sealed at the other. This configuration is used for practical reasons, specifically the need to minimize potential problems due to thermal expansion mismatch, at high temperatures, between the ceramic tube bundle and the steel housing. The primary technical problem with the candle-filter configuration for use in commercial-scale installations is the inability to purge the permeate side (typically the tube side), a feature that is crucial for high H2 recovery. In this study, we fabricated dual-end open, commercial-size ceramic multiple-tube bundles made of zeolite, palladium (Pd), and carbon molecular sieve (CMS) membranes that enable permeate-side (tube-side) purge for gas separation applications. Experimental gas separation data with these membrane bundles under harsh operating conditions (temperatures up to 350 and pressures up to 800 psig), to be presented at the meeting, manifest excellent performance. Parallel to the membrane bundle construction and testing efforts, we have also developed a detailed 3D CFD modeling package using COMSOL Multiphysics software to gain more insight into the effect on the H2 purity and recovery of the detailed geometry of the multi-tubular membrane system, including the number of tubes used, their dimensions, and placement in the bundle, as well as the number, type, and positioning of internal baffles and other flow-enhancement accessories. The CFD package is validated with experimental data from different systems (1-tube, 3-tube, and 19-tube bundles), and shows high accuracy in predicting the experimental results (<5 % error in all cases). The results of our study show that the detailed internal geometry of these multi-tubular membrane systems has a considerable impact on the performance of the system as the flow maldistribution within the shell-side can substantially decrease (>40%) the H2 recovery.
20 FOSSIL-FUELED POWER PLANTS↗
Increasing Radon Daughter Plate-out onto Superconducting Qubits
Radon daughter plate-out may pose a fundamental limitation to realizing large-scale superconducting quantum computers. Radon daughters plated-out onto a superconducting qubit chip can decay, releasing energy into the substrate, producing phonons. Phonon induced quasiparticles generate bit flip errors and reduce coherence times. We seek to increase plate-out on chips to enable future study of the impact of Radon on coherence times. In this project, we utilize COMSOL Multiphysics simulations to finalize design for a device to increase radon plate-out.
Kay, Dax [Northwestern U.]↗