SEARCH · Search NASA
Results for “source models”
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.
Matrix Approach to Land Carbon Cycle Modeling
Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchanges among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.
Open Source High Fidelity Modeling Of A Type-5 Wind Turbine Drivetrain
The purpose of this code is to provide a high-fidelity model of Type-5 wind turbine drivetrain for grid integration and transient stability studies. The increasing integration of renewable energy resources in evolving bulk power system (BPS) is impacting the system’s inertia. Type-5 wind turbine generation has the potential to behave like a traditional synchronous generator and can help mitigate the aforementioned impact. A hydraulic torque converter (TC) and gearbox with torque limiting feature are integral parts of a Type-5 wind turbine unit. A generic and high-fidelity model of a Type-5 wind turbine drivetrain is not openly and widely available for grid integration and transient stability studies. This hinders appropriate assessment of Type-5 wind power plant's contribution to the BPS resilience. This work develops and validates a TC model based on those generally used in automobile's transmission system.
Comparison of Open-Source Photovoltaic Performance Models Against Multi-Year Field Data
Explore the source record for details and available documents.
Comparison of Open-Source Photovoltaic Performance Models Against Multi-Year Field Data
Explore the source record for details and available documents.
Review of the 1st EUV Light Sources Code Comparison Workshop
We review the results of the 1 st Extreme Ultraviolet (EUV) Light Sources Code Comparison Workshop. The goal of this workshop was to provide a platform for specialists in EUV light source plasma modeling to benchmark and validate their numerical codes using well-defined case studies. Detailed consideration of a plethora of atomic collisional and radiative processes is required for modeling EUV light source plasmas. Eight institutions spanning four countries contributed data to the workshop. Two topics were addressed, namely (i) the atomic kinetics and radiative properties of tin plasmas under EUV-generating conditions and (ii) laser absorption in a fully ionized, one-dimensional hydrogen plasma. In this paper, we summarize the key findings of the workshop and outline plans for future iterations of the code comparison activity.
One Planet or Two Planets? The Ultra-sensitive Extreme-magnification Microlensing Event KMT-2019-BLG-1953
We present the analysis of a very high-magnification (A ∼ 900) microlensing event KMT-2019-BLG-1953. A single-lens single-source (1L1S) model appears to approximately delineate the observed light curve, but the residuals from the model exhibit small but obvious deviations in the peak region. A binary-lens (2L1S) model with a mass ratio of q ∼ 2 × 10{sup −3} improves the fits by Δχ {sup 2} = 181.8, indicating that the lens possesses a planetary companion. From additional modeling by introducing an extra planetary lens component (3L1S model) and an extra source companion (2L2S model), it is found that the residuals from the 2L1S model further diminish, but claiming these interpretations is difficult due to the weak signals with Δχ {sup 2} = 16.0 and 13.5 for the 3L1S and 2L2L models, respectively. From a Bayesian analysis, we estimate that the host of the planets has a mass of M{sub host}=0.31{sub −0.17}{sup +0.37} M{sub ⊙} and that the planetary system is located at a distance of D{sub L}=7.04{sub −1.33}{sup +1.10} kpc toward the Galactic center. The mass of the securely detected planet is M{sub p}=0.64{sub −0.35}{sup +0.76} M{sub J}. The signal of the potential second planet could have been confirmed if the peak of the light curve had been more densely observed by follow-up observations, and thus the event illustrates the need for intensive follow-up observations for very high-magnification events even in the current generation of high-cadence surveys.
A Systematic Framework for Tuning Open-Source Multifunctional IBR Models To Emulate OEM Black-Box Fault Dynamics
This paper presents a systematic framework to tune a generic IBR EMT model to match with an OEM provided balckbox inverter model based on the fault current responses. The key learnings and findings are summarized as follows: The tunable key parameters include inner control loops and current limiters to align the fault current magnitude, sequence content, and phase trajectories with the OEM models across diverse fault type and locations. The tuned model's fidelity is validated through comparative analysis with an OEM blackbox model, assessing both the fault current response and the responses of multiple relay elements. The results demonstrate the tuned generic model can trigger relay decision logic that is identical or near identical to that of the OEM model, thus generating very good match model for fault studies.
NRAP Open-Source Integrated Assessment Model and Relevant Application
NRAP Phase II Tools Workshop held in conjunction with GWPC 2021 Annual Forum, Virtual, September 29, 2021
NRAP-Open-IAM: Open-Source Integrated Assessment Model
NRAP Open IAM is the latest evolution in a series of quantitative models developed under the U.S. Department of Energy’s National Risk Assessment Partnership (NRAP) to address risk assessment, risk management and containment assurance questions of geological carbon storage (GCS). NRAP Open IAM allows a user to define a site specific GCS scenario, characterize important site features, events, and processes (considering uncertainty), couple together reduced order models of various system components to rapidly forecast response to CO 2 injection, and perform stochastic analysis to quantify containment effectiveness, leakage risk, and impact to receptors of concern.
Studying the Impact of IBR Modeling on the Commonly Applied Transmission Line Protective Elements
It is a well-accepted fact that detailed inverter-based resource (IBR) models are needed for electromagnetic transient (EMT) studies, especially for protection studies and protection designs. However, it is not clear which modeling aspects should be included in the EMT model to accurately represent the IBR fault responses. Therefore, this paper addresses this gap by evaluating the impact of key modeling aspects (DC source, inverter model, power control, current/voltage control, and current limiting) on the response of protective relay elements both qualitatively and quantitatively. Both grid-following (GFL) and grid-forming (GFM) IBR models are comprehensively evaluated through a simplified real-world transmission network by varying the IBR modeling aspects, fault locations, and types. The insights and findings from this paper will benefit the power system protection community because this work provides modeling requirements and guidelines for IBR modeling so that correct IBR models will be used for protection scheme settings design.
Numerical modelling of impact seismic sources using the stress glut theory
SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.
Multiscale modeling for high burnup structure formation and associated pulverization
This report summarizes the lower length scale modeling work performed in the fiscal year 2023 under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to capture the microstrutural evolution and associated pulverization criteria in the high burnup regions. This year, the resolution mechanisms within the cluster dynamics code has been updated and it influence on the bubble growth in HBS regions at the mesoscale is studied. To improve the accuracy of phase-field models of fission gas bubble growth, a new Helmholtz free energy for high-density Xe gas is derived from a virial equation of state. Two previously used strategies for representing net vacancy production in phase-field models of fission gas bubble growth are compared with each other and with analytical models of bubble growth. The vacancy source only model is found to be more convenient to parameterize realistically compared with the vacancy source+sink model. Furthermore, the phase-field-fracture simulation have been performed using MD- informed failure stress values and realistic HBS structure obtained from the phase-field simulations. We also present the uncertainty bands on prediction of the critical stress to account for the effect of the lower length scale variabilities on the failure criteria at the mesoscale. It is observed that pulverization may occur in partially restructured regions with restructuring fraction as low as 17%. In light of this, an update to the BISON’s pulverization criteria is recommended.
SALT3-NIR: Taking the Open-source Type Ia Supernova Model to Longer Wavelengths for Next-generation Cosmological Measurements
Abstract A large fraction of Type Ia supernova (SN Ia) observations over the next decade will be in the near-infrared (NIR), at wavelengths beyond the reach of the current standard light-curve model for SN Ia cosmology, SALT3 (∼2800–8700 Å central filter wavelength). To harness this new SN Ia sample and reduce future light-curve standardization systematic uncertainties, we train SALT3 at NIR wavelengths (SALT3-NIR) up to 2 μ m with the open-source model-training software SALTshaker , which can easily accommodate future observations. Using simulated data, we show that the training process constrains the NIR model to ∼2%–3% across the phase range (−20 to 50 days). We find that Hubble residual (HR) scatter is smaller using the NIR alone or optical+NIR compared to optical alone, by up to ∼30% depending on filter choice (95% confidence). There is significant correlation between NIR light-curve stretch measurements and luminosity, with stretch and color corrections often improving HR scatter by up to ∼20%. For SN Ia observations expected from the Roman Space Telescope, SALT3-NIR increases the amount of usable data in the SALT framework by ∼20% at redshift z ≲ 0.4 and by ∼50% at z ≲ 0.15. The SALT3-NIR model is part of the open-source SNCosmo and SNANA SN Ia cosmology packages.
Modeling glass degradation and release of radionuclides from vitrified waste for performance assessment simulations
The release of radionuclides initially encapsulated in a slowly degrading solid waste form and contained in an eventually corroding canister defines the source term for numerical simulations for the assessment of a geologic repository for high-level radioactive waste. While the details of waste degradation, canister corrosion, and dissolution and mobilization of the radionuclides in pore water include complex chemical reaction and transport processes that are coupled to the thermal, hydrological, microbiological, and mechanical conditions in the repository, the source-term model suitable for use in a numerical performance assessment model should be a defensible abstraction of these mechanisms. We developed a radiological source-term model and implemented it into a non-isothermal flow and transport simulator. While the proposed source-term model is applicable to various waste forms, canister systems, and disposal concepts, we specifically considered radionuclide releases from vitrified high-level waste placed in a cylindrical canister disposed in a deep vertical borehole repository. In this model, waste degradation is a function of temperature, and it can be adjusted to evaluate the influence of and propagate uncertainties in pH, passivation reactions, and chemical conditions as well as geometrical factors. The time-dependent, congruent release of safety-relevant radionuclides present in the decaying inventory is then calculated. Finally, the radionuclides are mobilized by diffusive and advective transport according to the thermo-hydraulic conditions prevailing in the near field of the repository, from where they migrate through the geosphere to the accessible environment. We examine the influence of the source-term model’s parameters on performance assessment calculations through sensitivity and uncertainty propagation analyses, identifying influential factors and confirming the upper bound of their impact. These considerations align with the overarching goal of repository design, which is to demonstrate that engineered and natural barriers can collectively delay radionuclide migration for timescales far exceeding human planning, thereby providing multiple, redundant barriers against environmental contamination.
Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures
Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.
An Open-Source Frequency-Domain Model for Floating Wind Turbine Design Optimization
A new frequency-domain dynamics model has been developed that uses open-source components to efficiently represent a complete floating wind turbine system. The model, called RAFT (Response Amplitudes of Floating Turbines), incorporates quasi-static mooring reactions, strip-theory and potential-flow hydrodynamics, blade-element-momentum aerodynamics, and linear turbine control. The formulation is compatible with a wide variety of support structure configurations and no manual or time-domain preprocessing steps are required, making RAFT very practical in design and optimization workflows. The model is applied to three reference floating wind turbine designs and its predictions are compared with results from time-domain OpenFAST simulations. There is good agreement in mean offsets as well the statistics and spectra of the dynamic response, verifying RAFT’s general suitability for floating wind analysis. Follow-on work will include verification of potential-flow and turbine-control features and application to optimization problems.
Vegetation biogeography is a main source of uncertainty in modelling the land carbon cycle
The terrestrial biosphere exchanges a large amount of CO 2 with the atmosphere through photosynthesis and respiration, determining the magnitude of land carbon sink and consequently influencing the rate of global warming. The magnitudes of global photosynthesis and respiration, however, vary widely across models (100-200 PgC/year), constituting a key and persistent source of uncertainty in carbon cycle and climate modelling. Here, we argue that the uncertainty in the land carbon cycle modelling is largely attributable to the uncertainty in biogeography – the distribution of plant functional types (PFTs). Using an ensemble of dynamic global vegetation models (DGVMs), we find a strong dependence of total photosynthesis on total area for each PFT. The dependence allows us to reduce the spread of land carbon cycle estimates by ~75% using remote sensing-based PFT maps. We further find that 56 ± 21% of climate-driven changes in global photosynthesis modelled by DGVMs are caused by changes in PFT distribution in the last two decades. Our study identifies vegetation biogeography as a main controlling factor of uncertainty in land carbon cycle modelling and highlights the importance of biogeography-climate interactions in carbon cycle and climate studies.