Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural modeling”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Structural Propensities in Cs2MBiX6 (M=Na, Ag; X=Cl, Br) Bismuth Halide Double Perovskites

A previously unreported low-temperature phase transition in the bismuth halide double perovskite Cs2AgBiCl6 is reported, thereby establishing trends in the structural ground state across Cs2NaBiCl6, Cs2AgBiCl6, and Cs2AgBiBr6. Using the combined toolkit of variable-temperature synchrotron X-ray and neutron powder diffraction, Raman spectroscopy, and density-functional theory–based electronic structure modeling, we demonstrate a cubic Fm¯3m → tetragonal I4/m transition upon cooling with distinct onset temperatures. Neutron powder diffraction refinements permit the unambiguously assignment of the low-temperature phase of Cs2NaBiCl6 to I4/m, correcting prior reports of an I4/mmm ground state. Cs2AgBiCl6 is also found to transforms to a structure crystallizing in the I4/m space group at low temperatures. Temperaturedependent Raman data and density-functional theory-based modeling capture the softening and freezing of out-of-phase octahedral-tilt modes and quantify relative instabilities. Solid-state nuclear magnetic resonance spectroscopy at room temperature completes the characterization and helps underpin the subtle differences in covalency across the compounds. Trends in the phase transition temperature Ts and tilt magnitudes emerge from coupled effects of halide identity, M(I)–site bonding character, and a mismatch between interatomic distances. These results establish the structure– dynamics–bonding framework for tuning tilt-driven instabilities in halide double perovskites.

Tian, Haowen↗
DOE OSTI2022

A method for automated modeling of plasma shape control response in tokamaks

Here, a simple method commonly used in modeling of tokamak conductor and plasma response for plasma control has been extended to allow more general modeling of tokamak PF conductor circuits. This extension combines a previously existing method for integrating plasma response and tokamak circuits with a more capable circuit modeling technique derived from the theory that underlies present-day automatic circuit simulation codes. This combined method can incorporate any combination of inductors (e.g., coils), resistors, capacitors, voltage sources, and current sources, while retaining the same model structure and ease of interpretation as the simple method. For many tokamaks, use of this extended method would not be the preferred approach, but for certain devices or for specialized models of a device, it is useful and occasionally essential. For example, the method has enabled fast and automated modeling of the complicated and daily-changing circuit topology used in the DIII-D tokamak. This automated modeling is essential to the support of controller design for the frequently changing requirements for DIII-D plasmas. It is likely to also find application in the workflows for new device designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Structural Polymorphism of Chitin and Chitosan in Fungal Cell Walls From Solid-State NMR and Principal Component Analysis

Chitin is a major carbohydrate component of the fungal cell wall and a promising target for novel antifungal agents. However, it is technically challenging to characterize the structure of this polymer in native cell walls. Here, we recorded and compared 13 C chemical shifts of chitin using isotopically enriched cells of six Aspergillus, Rhizopus, and Candida strains, with data interpretation assisted by principal component analysis (PCA) and linear discriminant analysis (LDA) methods. The structure of chitin is found to be intrinsically heterogeneous, with peak multiplicity detected in each sample and distinct fingerprints observed across fungal species. Fungal chitin exhibits partial similarity to the model structures of α- and γ-allomorphs; therefore, chitin structure is not significantly affected by interactions with other cell wall components. Addition of antifungal drugs and salts did not significantly perturb the chemical shifts, revealing the structural resistance of chitin to external stress. In addition, the structure of the deacetylated form, chitosan, was found to resemble a relaxed two-fold helix conformation. This study provides high-resolution information on the structure of chitin and chitosan in their cellular contexts. The method is applicable to the analysis of other complex carbohydrates and polymer composites.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Crystal structure of the $τ_{11}$ Al 4 Fe 1.7 Si phase from neutron diffraction and ab initio calculations

The intermetallic τ 11 Al 4 Fe 1.7 Si phase is of interest for high-temperature structural application due to its combination of low density and high strength. We determine the crystal structure of the τ 11 phase through a combination of powder neutron diffraction and density functional theory calculations. Using Pawley and Rietveld refinements of the neutron diffraction data provides an initial crystal structure model. Since Al and Si have nearly identical neutron scattering lengths, we use density-functional calculations to determine their preferred site occupations. The τ 11 phase exhibits a hexagonal crystal structure with space group P6 3 /mmc and lattice parameters of a = 7.478 Å and c = 7.472 Å. The structure comprises five Wyckoff positions; Al occupies the 6h and 12k sites, Fe the 2a and 6h sites, and Si the 2a sites. Here we observe site disorder and partial occupancies on all sites with a large fraction of 80% Fe vacancies on the 2d sites, indicating an entropic stabilization of the τ 11 phase at high temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Discovery and characterization of two Neptune-mass planets orbiting HD 212729 with TESS

ABSTRACT We report the discovery of two exoplanets orbiting around HD 212729 (TOI 1052, TIC 317060587), a Teff = 6146 K star with V = 9.51 observed by TESS in Sectors 1 and 13. One exoplanet, TOI-1052b, is Neptune-mass and transits the star, and an additional planet TOI-1052c is observed in radial velocities but not seen to transit. We confirm the planetary nature of TOI-1052b using precise radial velocity observations from HARPS and determined its parameters in a joint RV and photometry analysis. TOI-1052b has a radius of $2.87^{+0.29}_{-0.24}$ R⊕, a mass of 16.9 ± 1.7 M⊕, and an orbital period of 9.14 d. TOI-1052c does not show any transits in the TESS data, and has a minimum mass of $34.3^{+4.1}_{-3.7}$ M⊕ and an orbital period of 35.8 d, placing it just interior to the 4:1 mean-motion resonance. Both planets are best fit by relatively high but only marginally significant eccentricities of $0.18^{+0.09}_{-0.07}$ for planet b and $0.24^{+0.09}_{-0.08}$ for planet c. We perform a dynamical analysis and internal structure model of the planets as well as deriving stellar parameters and chemical abundances. The mean density of TOI-1052b is $3.9^{+1.7}_{-1.3}$ g cm−3 consistent with an internal structure similar to Neptune. A nearby star is observed in Gaia DR3 with the same distance and proper motion as TOI-1052, at a sky projected separation of $\scriptstyle \sim$1500 au, making this a potential wide binary star system.

Armstrong, David J. (ORCID:0000000250804117)↗
DOE OSTI2020

Probing the Atomic-Scale Structure of Amorphous Aluminum Oxide Grown by Atomic Layer Deposition

Atomic layer deposition (ALD) is a well-established technique for depositing nanoscale coatings with pristine control of film thickness and composition. The trimethylaluminum (TMA) and water (H 2 O) ALD chemistry is inarguably the most widely used and yet to date, we have little information about the atomic-scale structure of the amorphous aluminum oxide (AlO x ) formed by this chemistry. This lack of understanding hinders our ability to establish process-structure-property relationships and ultimately limits technological advancements employing AlOx made via ALD. In this work, we employ synchrotron high-energy X-ray diffraction (HE-XRD) coupled with pair distribution function (PDF) analysis to characterize the atomic structure of amorphous AlO x ALD coatings. We combine ex situ and in operando HE-XRD measurements on ALD AlO x and fit these experimental data using stochastic structural modeling to reveal variations in the Al-O bond length, Al and O coordination environment, and extent of Al vacancies as a function of growth conditions. In particular, the local atomic structure of ALD AlO x is found to change with the substrate and number of ALD cycles. The observed trends are consistent with the formation of bulk Al 2 O 3 surrounded by an O-rich surface layer. We deconvolute these data to reveal atomic-scale structural information for both the bulk and surface phases. Overall, this work demonstrates the usefulness of HE-XRD and PDF analysis in improving our understanding of the structure of amorphous ALD thin films and provides a pathway to evaluate how process changes impact the structure and properties of ALD films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Temperature Uncertainty Modeling with Proxy Structural Data as Geostatistical Constraints for Well Siting: An Example Applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aim of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations . The data used to support the modeling are measured temperatures and structural proxies from 3D geologic modeling, the latter considered "secondary" data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally likely temperature models (i.e., realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: Vprior and entropy. Vprior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e., linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g., dilation versus distance to fault termination) should be used given the specific geothermal system. In general, Vprior from local varying mean results identify locations that are close to high values for the structural proxies: areas with highe r probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comp arison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

3D temperature modeling↗
DOE OSTI2022

Two interaction surfaces between XPA and RPA organize the preincision complex in nucleotide excision repair

The xeroderma pigmentosum protein A (XPA) and replication protein A (RPA) proteins fulfill essential roles in the assembly of the preincision complex in the nucleotide excision repair (NER) pathway. We have previously characterized the two interaction sites, one between the XPA N-terminal (XPA-N) disordered domain and the RPA32 C-terminal domain (RPA32C), and the other with the XPA DNA binding domain (DBD) and the RPA70AB DBDs. Here, we show that XPA mutations that inhibit the physical interaction in either site reduce NER activity in biochemical and cellular systems. Combining mutations in the two sites leads to an additive inhibition of NER, implying that they fulfill distinct roles. Our data suggest a model in which the interaction between XPA-N and RPA32C is important for the initial association of XPA with NER complexes, while the interaction between XPA DBD and RPA70AB is needed for structural organization of the complex to license the dual incision reaction. Integrative structural models of complexes of XPA and RPA bound to single-stranded/double-stranded DNA (ss/dsDNA) junction substrates that mimic the NER bubble reveal key features of the architecture of XPA and RPA in the preincision complex. Most critical among these is that the shape of the NER bubble is far from colinear as depicted in current models, but rather the two strands of unwound DNA must assume a U-shape with the two ss/dsDNA junctions localized in close proximity. Our data suggest that the interaction between XPA and RPA70 is key for the organization of the NER preincision complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Thermodynamic Modeling of Intrinsic Defects in MnBi₂Te₄

This repository contains the computational data supporting the manuscript titled “The critical role of intrinsic defects and many-body interactions on the stability of MnBi₂Te₄.” It includes: 1. DFT data generated using VASP, used for training and benchmarking electronic structure models. 2. Quantum Monte Carlo (QMC) data produced with QMCPACK, used to apply many-body corrections and validate the electronic and magnetic properties of MnBi₂Te₄. 3. Relevant scripts used to run, analyze, and process the calculations, enabling reproducibility and transparency of the workflows.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Consequences of long-term water exposure for bulk crystal structure and surface composition/chemistry of nickel-rich layered oxide materials for Li-ion batteries

Water exposure of layered nickel-rich transition metal oxide electrodes, widely used in high-energy lithium-ion batteries, has detrimental effects on the electrochemical performance, which complicates electrode handling and prevents implementation of environmentally benign aqueous processing procedures. Elucidating the degradation mechanisms in play may help rationally mitigate/circumvent key challenges. Here, the bulk structural consequences of long-term (>2.5 years) deuterated water (D2O) exposure of intercalation materials with compositions Li x Ni 0.5 Co 0.2 Mn 0.3 O 2 (NCM523) and Li x Ni 0.8 Co 0.1 Mn 0.1 O 2 (NCM811) are studied by neutron powder diffraction (NPD). Detailed inspection of the NPD data reveals gradual formation of a secondary crystalline phase in all exposed samples, not previously reported for this system. This unknown phase forms faster in liquid- compared to vapor-exposed compounds. Structural modelling of the NPD data shows a stable level of Li/Ni anti-site defects and does not indicate any significant changes in lattice parameters or hydrogen-lithium (D + /Li + ) exchange in the structure. Consequently, the secondary phase formation must take place via transformation rather than modification of the parent material. X-ray photoelectron spectroscopy data indicate formation of LiHCO 3 /Li 2 CO 3 at the surface and a Li-deficient oxide in the sub-surface region of the pristine compounds, and the presence of adsorbed water and transition metal hydroxides at the exposed sample surfaces.

25 ENERGY STORAGE↗
DOE OSTI2020

Evaluation of Leading Modes of Climate Variability in the CMIP Archives

The adequate simulation of internal climate variability is key for our understanding of climate as it underpins efforts to attribute historical events, predict on seasonal and decadal time scales, and isolate the effects of climate change. Here the skill of models in reproducing observed modes of climate variability is assessed, both across and within the CMIP3, CMIP5, and CMIP6 archives, in order to document model capabilities, progress across ensembles, and persisting biases. A focus is given to the well-observed tropical and extratropical modes that exhibit small intrinsic variability relative to model structural uncertainty. These include El Niño–Southern Oscillation (ENSO), the Pacific decadal oscillation (PDO), the North Atlantic Oscillation (NAO), and the northern and southern annular modes (NAM and SAM). Significant improvements are identified in models’ representation of many modes. Canonical biases, which involve both amplitudes and patterns, are generally reduced across model generations. For example, biases in ENSO-related equatorial Pacific sea surface temperature, which extend too far westward, and associated atmospheric teleconnections, which are too weak, are reduced. Stronger tropical expression of the PDO in successive CMIP generations has characterized their improvement, with some CMIP6 models generating patterns that lie within the range of observed estimates. For the NAO, NAM, and SAM, pattern correlations with observations are generally higher than for other modes and slight improvements are identified across successive model generations. Finally, for ENSO and PDO spectra and extratropical modes, changes are small compared to internal variability, precluding definitive statements regarding improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluating alternative ebullition models for predicting peatland methane emission and its pathways via data–model fusion

Abstract. Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structural Insights into Mechanisms Underlying Mitochondrial and Bacterial Cytochrome c Synthases

Mitochondrial holocytochrome c synthase (HCCS) is an essential protein in assembling cytochrome c (cyt c) of the electron transport system. HCCS binds heme and covalently attaches the two vinyls of heme to two cysteine thiols of the cyt c CXXCH motif. Human HCCS recognizes both cyt c and cytochrome c1 of complex III (cytochrome bc1). HCCS is mutated in some human diseases and it has been investigated recombinantly by mutational, biochemical, and reconstitution studies in the past decade. Here, we employ structural prediction programs (e.g., AlphaFold 3) on HCCS and its two substrates, heme and cytochrome c. The results, when combined with spectroscopic and functional analyses of HCCS and variants, provide insights into the structural basis for heme binding, apocyt c binding, covalent attachment, and release of the holocyt c product. Results from in vitro reconstitution of purified human HCCS using cyt c and cyt c1 peptides as acceptors are consistent with the structural modeling of substrate binding. Reconstitution of HCCS and cyt c1 provides an approach to studying cyt c1 assembly, which has been refractile to recombinant in vivo reconstitution (unlike HCCS and cyt c). We propose a structural basis for release of the holocyt c product from HCCS based on in vitro studies and on cryoEM structures of the bacterial cyt c synthase (CcsBA) active site. We analyze the kinetoplastid mitochondrial synthase (KCCS), and hypothesize a molecular evolutionary path from mitochondrial endosymbiosis to the current HCCS.

Biochemistry & Molecular Biology↗
DOE OSTI2024

Parameter identification methods for low-order gray box building energy models: A critical review

The body of knowledge on gray box building energy modeling (GBBEM) has been developed over the past few decades and has undergone some important changes recently. Starting with simple methods and simple buildings, the science of GBBEM has grown to encompass complex and more computationally intensive techniques and complex commercial buildings. Numerous works including a recent review have considered model structure and inputs in a nearly systematic way, but no extant work systematically reviews the approaches for GBBEM parameter search initialization and final identification, despite this being arguably the most difficult and impactful part of the modeling process. To this end, we critically review 55 extant works describing advantages, limitations, and domain of applicability of several classes of parameter initialization and optimization techniques specifically for GBBEM. We categorize the classes of methods and analyze the evidence of their applicability for different applications within the field of GBBEM. We find an emerging consensus that initialization of parameter searches for anything other than the simplest building elements is often challenging and sometimes requires a stochastic approach to begin the parameter identification process. After this initial process, faster methods have been used in some cases but often the nature of the problem requires stochastic methods for this portion of the process as well. For less complex systems, more deterministic and efficient methods have been shown to be effective. Finally, we draw conclusions as to the domain of applicability of different classes of initialization and optimization techniques for GBBEEM and offer suggestions for research directions that are likely to prove fruitful.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Modelling Development in Radical (Co)Polymerization of Multivinyl Monomers

Radical polymerization (RP) of multivinyl monomers (MVMs) provides a facile solution for manipulating polymer topology and has received increasing attention due to their industrial and academic significance. Continuous efforts have been made to understand their mechanism, which is the key to regulating materials structure. Modelling techniques have become a powerful tool that can provide detailed information on polymerization kinetics which is inaccessible by experiments. Many publications have reported the combination of experiments and modelling for free radical polymerization (FRP) and reversible-deactivation radical polymerizations (RDRP) of MVMs. Herein, a minireview is presented for the most important modelling techniques and their applications in FRP/RDRP of MVMs. This review hopes to illustrate that the combination of modelling and wet experiments can be a great asset to polymer researchers and inspire new thinking for the future MVMs experiment optimization and product design.

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

Nearly hyperuniform, nonhyperuniform, and antihyperuniform density fluctuations in two-dimensional transition metal dichalcogenides with defects

Hyperuniform many-body systems in d-dimensional Euclidean space R d are characterized by completely suppressed (normalized) infinite-wavelength density fluctuations, and appear to be endowed with novel exotic physical properties. Recently, hyperuniform systems of disordered varieties have been observed in the context of various atomic-scale two-dimensional (2D) materials. Here, we analyze the effects of localized defects on the density fluctuations across length scales and on the hyperuniformity property of experimental samples of two-dimensional transition metal dichalcogenides. In particular, we extract atomic coordinates from time series annular dark field-scanning transmission electron microscopy (ADF-STEM) imaging data of 2D tungsten chalcogenides with the 2H structure (Te-doped 2H-WSe 2 ) showing continuous development and evolution of electron-beam induced defects, and construct the corresponding chemical-bonding informed coordination networks between the atoms. We then compute a variety of pair statistics and bond-orientational statistics to characterize the samples. At low defect concentrations, the corresponding materials are nearly hyperuniform, characterized by significantly suppressed scattering at the zero wave-number limit (omitting forward/ballistic scattering). As more defects are introduced, the (approximate) hyperuniformity of the materials is gradually destroyed, and the system becomes non-hyperuniform even when the material still contains a significant amount of crystalline regions. At high defect concentrations, the structures become antihyperuniform with diverging (normalized) large-scale density fluctuations, mimicking those typically observed at the thermal critical points associated with phase transitions. Overall, the defected materials possess varying degrees of orientation order, and there is apparently no intermediate hexatic phase emerging. To understand the observed nearly hyperuniform density fluctuations in the slightly defected materials, we construct a minimalist structural model and demonstrate that the experimental samples can be essentially viewed as perturbed honeycomb crystals with small correlated displacements and double chalcogen vacancies. Moreover, the small correlated displacements alone can significantly degrade hyperuniformity of the perfect honeycomb structure. Therefore, even a small amount of vacancies, when coupled with correlated displacements, can completely destroy hyperuniformity of the system.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

PORFLOW modeling of Vadose zone flow and transport for the E-area intermediate level vault

In support of the E-Area Performance Assessment, a two-dimensional model of water flow and radionuclide transport through the E-Area Intermediate Level Vault (ILV) and local vadose zone has been developed using the PORFLOWTM software. The purpose of the model is to calculate flux to the water table for radionuclides eluted from the ILV during its operational life, the period of institutional control, and times following site closure. Results of model calculations will be used by a three-dimensional PORFLOW model of transport through the aquifer to determine radionuclide concentrations at a hypothetical 100 meter well and at the site boundary where contaminated groundwater could be assessed by members of the public following site closure. While newly developed, the model structure closely follows that used in the 2008 PA while incorporating a refined computational mesh, updated material properties for the vadose zone soil and vault concrete, revised infiltration rates, and new solid-liquid distribution coefficients. The model also addresses degradation of concrete hydraulic properties by blending soil and concrete water retention curves over a 500 year period. That is, the hydraulic properties of ILV concrete start out as fresh concrete at the time of site closure and degrade to soil properties over a 500 year time period. This approach has not been used previously in E-Area PA’s. The model has been used to calculate water flow through the ILV vadose zone and, on a trial basis, transport of a limited number of radionuclides from the waste region inside the ILV to the water table. Figure ES-1 provides a view of the model and an example of water flow after the concrete has substantially degraded. Figure dimensions are in centimeters. The red area in the left hand side figure is the waste disposal region within the ILV. Colored bands in the figures indicate different soil and concrete regions. Figure ES-2 shows a calculation of flux to the water table (mol/yr) for U-234 and its one-year half-life decay chain (U-234 → Th-230 → Ra-226 → Pb-210). Pb-210 flux is too small to show on the same scale used for Th-230 and Ra-226 in Figure ES-2. The model has also been used to make a preliminary evaluation of the release of tritium from TPBAR disposal containers that will be place in the ILV. A separately calculated tritium release rate as a function of time was used as a source in the PORFLOW model. The TPBAR disposal schedule and tritium release rate differ from the values used in the 2008 PA. Results from these initial trial runs are compared to calculations made for the 2008 PA ILV analysis. In brief, a new PORFLOW model of flow and transport through the ILV vadose zone has been developed and tested. The model is ready to be used for PA calculations.

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

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗