Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling methodology”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Photoluminescence of Cis-Polyacetylene Semiconductor Material

Photoluminescence (PL) is one of the key experimental characterizations of optoelectronic materials, including conjugated polymers (CPs). In this study, a simplified model of an undoped cis-polyacetylene (cis-PA) oligomer was selected and used to explain the mechanism of photoluminescence (PL) of the CPs. Using a combination of the ab initio electronic structure and a time-dependent density matrix methodology, the photo-induced time-dependent excited state dynamics were computed. We explored the phonon-induced relaxation of the photoexcited state for a single oligomer of cis-PA. Here, the dissipative Redfield equation of the motion was used to compute the dissipative excited state dynamics of electronic degrees of freedom. This equation used the nonadiabatic couplings as parameters. The computed excited state dynamics showed that the relaxation rate of the electron is faster than the relaxation rate of the hole. The dissipative excited-state dynamics were combined with radiative recombination channels to predict the PL spectrum. The simulated results showed that the absorption and emission spectra both have a similar transition. The main result is that the computed PL spectrum demonstrates two mechanisms of light emission originating from (i) the inter-band transitions, corresponding to the same range of transition energies as the absorption spectrum and (ii) intra-band transitions not available in the absorption spectra. However, the dissipative Redfield equation of the motion was used to compute the electronic degrees of freedom of the nonadiabatic couplings, which helped to process the time propagation of the excited dynamic state. This excited dynamic state shows that the relaxation rate of the electron is faster than the relaxation rate of the hole, which can be used for improving organic semiconductor materials for photovoltaic and LED applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗
DOE OSTI2023

Variational deep learning of equilibrium transition path ensembles

Here, we present a time-dependent variational method to learn the mechanisms of equilibrium reactive processes and efficiently evaluate their rates within a transition path ensemble. This approach builds off of the variational path sampling methodology by approximating the time-dependent commitment probability within a neural network ansatz. The reaction mechanisms inferred through this approach are elucidated by a novel decomposition of the rate in terms of the components of a stochastic path action conditioned on a transition. This decomposition affords an ability to resolve the typical contribution of each reactive mode and their couplings to the rare event. The associated rate evaluation is variational and systematically improvable through the development of a cumulant expansion. We demonstrate this method in both over- and under-damped stochastic equations of motion, in low-dimensional model systems, and in the isomerization of a solvated alanine dipeptide. In all examples, we find that we can obtain quantitatively accurate estimates of the rates of the reactive events with minimal trajectory statistics and gain unique insights into transitions through the analysis of their commitment probability.

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

Neural simulation-based inference of the neutron star equation of state directly from telescope spectra

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Towards Automated and High-Throughput Quantitative Sizing and Isotopic Analysis of Nanoparticles via Single Particle-ICP-TOF-MS

The work described herein assesses the ability to characterize gold nanoparticles (Au NPs) of 50 and 100 nm, as well as 60 nm silver shelled gold core nanospheres (Au/Ag NPs), for their mass, respective size, and isotopic composition in an automated and unattended fashion. Here, an innovative autosampler was employed to mix and transport the blanks, standards, and samples into a high-efficiency single particle (SP) introduction system for subsequent analysis by inductively coupled plasma–time of flight–mass spectrometry (ICP-TOF-MS). Optimized NP transport efficiency into the ICP-TOF-MS was determined to be >80%. This combination, SP-ICP-TOF-MS, allowed for high-throughput sample analysis. Specifically, 50 total samples (including blanks/standards) were analyzed over 8 h, to provide an accurate characterization of the NPs. This methodology was implemented over the course of 5 days to assess its long-term reproducibility. Impressively, the in-run and day-to-day variation of sample transport is assessed to be 3.54 and 9.52% relative standard deviation (%RSD), respectively. The determination of Au NP size and concentration was of <5% relative difference from the certified values over these time periods. Isotopic characterization of the 107 Ag/ 109 Ag particles (n = 132,630) over the course of the measurements was determined to be 1.0788 ± 0.0030 with high accuracy (0.23% relative difference) when compared to the multi-collector–ICP-MS determination.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A consumer-centric approach to quantify efficiency of receiving goods purchased via online

Virtual participation in shopping activities has increased exponentially in the past four years compared to the last couple of decades. E-tailing or online shopping offers the convenience of goods reaching a consumer instead of a consumer traveling to a store, but it has downsides like geographical service variability and negative social externalities such as increased energy consumption and emissions. This study proposes a novel approach to quantify e-tailing efficiency from the consumers’ viewpoint. The methodology is innovative in its integration of accessibility theory with energy and cost impedance factors and consideration of delivering and picking up goods purchased via online. The methodology is implemented for the San Francisco Bay Area and is subject to scenarios that see enhancements to various facets of online shopping delivery. Results demonstrate that increasing the frequency of e-commerce deliveries helps improve e-tailing efficiency in rural locations, while improvements in energy and cost aspects of delivery modes are seen to improve e-tailing efficiencies in the central parts of the region. The approach proposed can provide valuable insights on where people have limited benefits from online shopping and how emerging delivery mechanisms can change the quality of the e-commerce experience within a city. This research offers a replicable framework for assessing e-commerce systems in diverse geographic contexts, contributing to the development of equitable and environmentally sustainable urban freight systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Risk to Ecological Resources on Burial Waste Sites on the Hanford Site, Washington - 20248

The US and other nations have a large nuclear and chemical waste remediation task remaining from World War II, the Cold War, and industrial activities. The US Department of Energy (DOE) has the largest task in the U.S., with Hanford (Washington State) being the most contaminated site, with the longest and most costly cleanup program projected to extend to the end of the century. One of the greatest remediation problems at Hanford are the burial waste facilities, some of which are inside the 200 Areas on the Central Plateau, while others are outside the 200 Areas. These areas are associated with waste sites, pipes, and contaminated trenches, posing technical cleanup challenges to meet DoE's remediation goals, including protection of human health and the environment commensurate with Tri- Party Agreement and current regulatory requirements. In this paper we select one of the legacy waste sites (the BC Cribs and Trenches) as a case study to explore how risk to ecological resources varies during and as a consequence of remediation on sites primarily considered legacy waste sites. We used the risk methodology developed by the Consortium for Risk Evaluation with Stakeholder Participation (CRESP) during the Hanford Site-wide Risk Review. The CRESP methodology included evaluating the risk to receptors, including human health, the Columbia River, groundwater, and ecological and cultural resources. The risk evaluation process for ecological receptors included defining the boundaries of each evaluation unit (EU) to be remediated, definition of a suitable buffer area that might be affected, and consideration of the known resource value based on field work. This information, coupled with known remediation options or strategies, allowed us to develop a risk rating for ecological resources. Site risks were evaluated on a scale of 'not discernible' to 'very high' risk. Consideration of the risks at the buried waste sites was complicated by uncertainties in remediation methods and the degree of restoration that would occur as part of the remediation process. We use the BC Cribs and Trenches as a case study to explore why the risks are low currently but increase during remediation. The current 'not discernible to low' risk rating is due to the lack of high quality ecological resources on the BC Cribs and Trenches EU, except for edge areas (with some state-listed species). During remediation the risks increase to 'medium', mainly due to high quality resources on adjacent buffer lands, which could be disrupted during active remediation. Medium and high quality resources on the edge of the site are continuous with higher quality native shrub-steppe habitat on the buffer. Ecological resources on some of the areas surrounding the 200 East and West Areas are some of the highest quality shrub-steppe areas on the Hanford Site, and risk to these resources could be prevented by planning before and during remediation to minimize disruption. For example, lay down areas for the BC Cribs and Trenches can be minimized by careful planning. Heavy truck and other vehicular traffic can be minimized. Spread of the noxious weeds now on some of the site can be prevented, and the spread of invasive species into the site (particularly near undisturbed shrub-steppe habitat) can be reduced. Further, the areas around and between the 200 East and West Areas have been largely undisturbed and can serve as a source of seeds and propagules for any revegetation that occurs as part of the remediation process. (authors)

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

Unraveling complex causal processes that affect sustainability requires more integration between empirical and modeling approaches

Scientists seek to understand the causal processes that generate sustainability problems and determine effective solutions. Yet, causal inquiry in nature–society systems is hampered by conceptual and methodological challenges that arise from nature–society interdependencies and the complex dynamics they create. Here, we demonstrate how sustainability scientists can address these challenges and make more robust causal claims through better integration between empirical analyses and process- or agent-based modeling. To illustrate how these different epistemological traditions can be integrated, we present four studies of air pollution regulation, natural resource management, and the spread of COVID-19. The studies show how integration can improve empirical estimates of causal effects, inform future research designs and data collection, enhance understanding of the complex dynamics that underlie observed temporal patterns, and elucidate causal mechanisms and the contexts in which they operate. These advances in causal understanding can help sustainability scientists develop better theories of phenomena where social and ecological processes are dynamically intertwined and prior causal knowledge and data are limited. The improved causal understanding also enhances governance by helping scientists and practitioners choose among potential interventions, decide when and how the timing of an intervention matters, and anticipate unexpected outcomes. Methodological integration, however, requires skills and efforts of all involved to learn how members of the respective other tradition think and analyze nature–society systems.

42 ENGINEERING↗
DOE OSTI2024

A numerical study of process complexity in permafrost dominated regions

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Calibration↗
DOE OSTI2025

Prediction of Redox Potentials for Different Oxidation States of U, Np, Pu, and Am in Alkaline Aqueous Solution

The redox potentials for U, Np, Pu, and Am for oxidation states +III up to +VIII in alkaline aqueous solutions were predicted using density functional theory (DFT) and small-core pseudopotentials and their basis sets, with a hybrid explicit/implicit solvent model using SHE = 4.28 V. For each oxidation state, various oxo/hydroxo complexes were evaluated resulting in a variety of one-electron redox pathways. For An(VIII/VII) couples, the predicted redox potentials for the [An(VIII)O 5 (OH)] -3 /[An(VII)O 4 (OH) 2 ] -3 , or [An(VIII)O 4 (OH) 2 ] -2 /[An(VII)O 4 (OH) 2 ] -3 couples are in good agreement with existing estimates. For An(VII/VI) redox couples, all couples, particularly [An(VII)O 4 (OH) 2 ] -3 /[An(VI)O 2 (OH) 4 ] -2 , were in agreement with experimental values for U, Np, and Pu, but the results for Am showed larger differences from the estimated potentials. The An(VI/V) couples were consistent with experiment for dioxo/tetrahydroxo couples, and the An(V/IV) couples showed acceptable agreement based on actinide-specific couples, with neutral hydroxides often favored in the +IV state. The An(IV/III) couples were consistent with literature values when modeled as soluble neutral hydroxides. The use of our approach yielded calculated redox potentials that were within ±0.2 V from experimental or estimated values consistent with our prior calculations on redox potentials of actinides from Ac to Am in acidic aqueous solutions. Finally, as a result, this supports the robustness of our DFT-based methodology for predicting actinide redox potentials, offering valuable insights for actinide chemistry in aqueous solutions.

Actinides↗
DOE OSTI2021

DOE Energy Frontier Research Centers Center for Direct Catalytic Conversion of Biomass to Biofuels (C3Bio)

New capabilities to predict, design and control the chemistries of carbon could answer a global imperative to transition from fossil-based to sustainable transportation fuels. While the use of inexpensive hydrocarbons has been an unparalleled achievement and enabler of economic prosperity for many nations, singular dependence upon crude oil has given rise to systemic vulnerabilities in climate, energy, economic, and national security. Lignocellulosic biomass, a renewable and carbon-neutral resource, has the potential to displace an estimated annual equivalent of three billion barrels of oil in the U.S. alone (National Research Council 2009, U. S. Department of Energy, 2011). However, biomass has only one-third the energy density of crude oil (Agrawal and Singh 2009, Richard 2010) and lacks petroleum’s versatility as a feedstock for fuels and chemicals. These limitations keep biomass conversion below the efficiency level needed for strategic impact while the scientific challenge of routing carbon from one molecular context to another remains unmet. In 2009, the Center for Direct Catalytic Conversion of Biomass to Biofuels (C3Bio) recognized the potential of chemical catalysis and fast pyrolysis to overcome such limitations by transforming the main components of biomass (cellulose, xylan, and lignin) from grasses and trees directly to liquid hydrocarbons and aromatic co-products. Enabled by the EFRC high-risk, high-reward approach to grand challenge science, C3Bio researchers have been key national players in disrupting the conventional paradigm of the cellulosic biorefinery into a new future of “no carbon left behind”—the full utilization of carbon from plant cell walls in energy-dense fuels (Fig. 1). We identified catalytic and fast-pyrolytic pathways that utilize cellulose, xylan and, most significantly, lignin. We developed catalytic processes that deoxygenate and transform monomers and isolated polymers into useful products and tested their use with intact biomass. We gained control of lignin synthesis within plants and initiated tailoring biomass to its end-use through the tools of plant molecular biology and genetic engineering. C3Bio breakthroughs have increased the energy density of biomass-derived substrates via catalytic and pyrolytic conversions into products such as benzoquinones, furfural and hydroxymethylfurfural, levoglucosan and levulinic acid, methoxypropylphenols and propylbenzene. Such advances in biomass conversion would not have been possible without simultaneous advances in analytical instrumentation and methodologies, and imaging technologies and applications. The legacy science developed by C3Bio enables design and control strategies for achieving a targeted product portfolio of fuel and chemical feedstocks from a diverse range of native and tailored biomass. Our research provides the knowledge base required for a bio-economy with product streams as diverse in functionality as those of the petrochemical industry. Coupling targeted computational modeling with experimentation, we achieved: (1) fundamental understanding of biopolymers and cell wall architecture assembly, (2) discovery of new chemistries that allow the development of highly selective pathways to fuels and desirable chemicals, and (3) an integrated systems-level understanding to control catalytic and pyrolytic pathways.

09 BIOMASS FUELS↗
DOE OSTI2022

Multimode Characterization Approach for Understanding Cell-Level PV Performance and Degradation

Cell-level degradation processes impact the economic viability and large-scale deployment prospects for both established and emerging photovoltaic (PV) technologies. This project addresses the need to develop experimental and device-modeling approaches for studying cell-level degradation processes in photovoltaic (PV) devices using a variety of characterization techniques that provide access to complementary material and device properties. Our results demonstrate that by coupling characterization results with device modeling it is possible to develop comprehensive understanding of processes leading to performance limitations and degradation. This project developed a suite of novel measurement techniques including pulsed-light-bias operando X-ray and photoelectron spectroscopy (popXPS), light-biased scanning microwave impedance microscopy (sMIM), and near-field transport imaging (TI). In addition, operando characterization methodologies and in situ stressing capabilities have been developed and applied for techniques including electron-beam-induced current (EBIC), cathodoluminescence (CL), and Kelvin probe force microscopy (KPFM). Device-physics models were developed and applied to simulate correlative, multi-mode measurements to extract material and device parameters that control performance degradation. These characterization and modeling techniques were applied in a multi-mode approach to probe cell-level degradation mechanisms in Cd(Se,Te) and hybrid perovskite PV devices. Together these efforts contribute to foundational PV degradation science by establishing a framework for understanding PV performance degradation at the cell level and benefit the U.S. PV industry by providing resources in the form of novel experimental capabilities, knowledge gained, and available expertise that can accelerate research and development of improved PV device materials and architectures. The project provided a comprehensive understanding of degradation in baseline Cd(Se,Te) solar cells provided by our collaborators at Colorado State University (CSU). EBIC and CL-based measurements and revealed unusual collection and recombination profiles in these devices, which underwent significant changes with during stressing. KPFM and operando XPS measurements showed that device stressing permanently alters energy-band alignments at the (Mg,Zn)O/Cd(Se,Te) interface, which in turn account for an observed loss in fill factor. Studies on hybrid perovskite devices were hampered to a significant extent by delays related to the pandemic. Nevertheless, a set of hybrid perovskite devices (supplied through an NREL-industry partnership) were stress tested and characterized with techniques including EBIC, sMIM, popXPS/popUPS and optically excited TI. Available results from these measurements informed the device modeling effort and suggest that defects and related band offsets at the C60/LiF/hybrid perovskite interface are the primary source of degradation in these devices.

14 SOLAR ENERGY↗
DOE OSTI2026

Systematic Benchmarking of Climate Models: Methodologies, Applications, and New Directions

As climate models become increasingly complex, there is a growing need to comprehensively and systematically assess model performance with respect to observations. Given the increasing number and diversity of climate model simulations in use, the community has moved beyond simple model intercomparison and toward developing methods capable of benchmarking a large number of simulations against a suite of climate metrics. Here, we present a detailed review of evaluation and benchmarking methods and approaches developed in the last decade, focusing primarily on scientific implications for Coupled Model Intercomparison Project (CMIP) simulations and CMIP6 results that contributed to the Intergovernmental Panel on Climate Change (IPCC) Sixth Assessment Report (AR6). Based on this review, we explain the resulting contemporary philosophy of model benchmarking, and provide clear distinctions and definitions of the terms model verification, process validation, evaluation, and benchmarking. While significant progress has been made in model development based on systematic evaluation and benchmarking efforts, some climate system biases still remain. The development of open‐source community software packages has played a fundamental role in identifying areas of significant model improvement and bias reduction. We review the key features of several software packages that have been commonly used over the past decade to evaluate and benchmark global and regional climate models. Additionally, we discuss best practices for the selection of evaluation and benchmarking metrics and for interpreting the obtained results, the importance of selecting suitable sources of reference data and accurate uncertainty quantification.

Environmental sciences↗
DOE OSTI2025

A Practical Approach to Inertia Assessment: Analyzing New York's Grid Response to Real Events

The non-uniform distribution of inverter-based resources (IBRs) such as solar and wind energy leads to the formation of weakly coupled regions within the broader interconnection. This highlights the urgency of investigating regional dynamics. This paper proposes a practical framework for estimating regional inertia using actual event measurements from the New York (NY) region. The regional frequency is determined using measurement data from frequency disturbance recorders (FDRs) deployed across NY. The rate-of-change of frequency (RoCoF) is estimated using a hybrid methodology that combines the benefits of both moving and non-moving window techniques. By utilizing the power mismatch value and the estimated RoCoF, various inertia metrics such as regional RoCoF, inertia arrival time, and inertia percentage relative to the interconnection inertia are calculated. The identified confirmed-events that took place between 2017–2022 in NY are used in the study. This assessment of inertia using real event data reveals the true behavior of the grid under stress. The results indicate that as the generation mix in NY evolves with the rise of IBRs, both the RoCoF levels and the regional inertia undergo significant changes.

Dulal, Saurav [University of Tennessee, Knoxville ↗
DOE OSTI2022

Segmenting water and shadow regions within WorldView imagery using local binary patterns

Many shadow detection algorithms pertaining to remotely sensed imagery exist. Several of these algorithms exploit the spectral characteristics of shadows within imagery to identify shaded regions. However, these algorithms can have problems when water is also present within the imagery because water shares similar spectral characteristics with shadows. Many of these algorithms are applied to small image subsets instead of the image as a whole and frequently are applied to urban environments that require additional use of the normalized difference water index or other features to affect the removal of water from the shadow mask. This diversity of the image scene content coupled with the complexity and wide variety of environmental conditions that satellite imagery can acquire makes reliably separating shadow and water within larger images a complex problem. Thresholding various spectral indices to produce segmentation maps can be challenging when large, complex, and often imbalanced scenes are captured, which may require manual adjustment of algorithmic parameters for different image areas. Here, we present a methodology that makes use of the near-infrared channel using a watershed segmentation algorithm, local binary pattern measurements, and a support vector machine to classify shadow and water within a full WorldView scene. Results are promising with initial accuracies above 97% for both shadow and water.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Investigating Resiliency of Transportation Network under Targeted and Potential Climate Change Disruptions

Ensuring robustness and resilience in intermodal transportation systems is essential for the continuity and reliability of global logistics. These systems are vulnerable to various disruptions, including natural disasters and technical failures. Despite significant research on freight transportation resilience, investigating the robustness of the system after targeted and climate-change-driven disruption remains a crucial challenge. Drawing on network science methodologies, this study models the interdependencies within the rail and water transport networks and simulates different disruption scenarios to evaluate system responses. Here, we use the data from the U.S. Department of Energy Volpe Center for network topology and tonnage projections. The proposed framework quantifies deliberate, stochastic, and climate-driven infrastructure failure, using higher resolution downscaled multiple Earth System Models’ simulations from Coupled Model Intercomparison Project Phase version 6. We show that the disruptions of a few nodes could have a larger impact on the total tonnage of freight transport than on network topology. For example, the removal of targeted 20 nodes can bring the total tonnage carrying capacity to 30% with about 75% of the rail freight network intact. This research advances the theoretical understanding of transportation resilience and provides practical applications for infrastructure managers and policymakers. By implementing these strategies, stakeholders and policymakers can better prepare for and respond to unexpected disruptions, ensuring sustained operational efficiency in transportation networks.

Climate Change Disruptions↗
DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗