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

Structural and biochemical analyses of selectivity determinants in chimeric Streptococcus Class A sortase enzymes

Abstract Sequence variation in related proteins is an important characteristic that modulates activity and selectivity. An example of a protein family with a large degree of sequence variation is that of bacterial sortases, which are cysteine transpeptidases on the surface of gram‐positive bacteria. Class A sortases are responsible for attachment of diverse proteins to the cell wall to facilitate environmental adaption and interaction. These enzymes are also used in protein engineering applications for sortase‐mediated ligations (SML) or sortagging of protein targets. We previously investigated SrtA from Streptococcus pneumoniae , identifying a number of putative β7–β8 loop‐mediated interactions that affected in vitro enzyme function. We identified residues that contributed to the ability of S. pneumoniae SrtA to recognize several amino acids at the P1′ position of the substrate motif, underlined in LPXT G , in contrast to the strict P1′ Gly recognition of SrtA from Staphylococcus aureus . However, motivated by the lack of a structural model for the active, monomeric form of S. pneumoniae SrtA, here, we expanded our studies to other Streptococcus SrtA proteins. We solved the first monomeric structure of S. agalactiae SrtA which includes the C‐terminus, and three others of β7–β8 loop chimeras from S. pyogenes and S. agalactiae SrtA. These structures and accompanying biochemical data support our previously identified β7–β8 loop‐mediated interactions and provide additional insight into their role in Class A sortase substrate selectivity. A greater understanding of individual SrtA sequence and structural determinants of target selectivity may also facilitate the design or discovery of improved sortagging tools.

Gao, Melody↗
DOE OSTI2021

Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication

Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lithium-Ion Battery Life Model with Electrode Cracking and Early-Life Break-in Processes

This paper develops a physically justified reduced-order capacity fade model from accelerated calendar- and cycle-aging data for 32 lithium-ion (Li-ion) graphite/nickel-manganese-cobalt (NMC) cells. The large data set reveals temperature-, charge C-rate-, depth-of-discharge-, and state of charge (SOC)-dependent degradation patterns that would be unobserved in a smaller test matrix. Model structure is informed by incremental capacity analysis that shows loss of lithium inventory and cathode-material loss as the dominant capacity fade mechanisms. The model includes terms attributable to solid-electrolyte interface (SEI) growth, electrode cracking, cycling-driven acceleration of SEI growth, and "break-in" mechanisms that slightly decrease or increase available Li inventory early in life. The study explores what mathematical couplings of these mechanisms best describe calendar aging, cycle aging, and mixed calendar/cycle aging. Various approaches are discussed for extracting relevant stress factors from complex cycling profiles to predict lifetime during real-world battery loads using models trained on constant-current laboratory test results. The complexity of the present human-driven model identification process motivates future work in machine learning to more widely search and statistically discern the optimal model that correctly extrapolates capacity fade based on physical knowledge.

25 ENERGY STORAGE↗
DOE OSTI2020

Deep Learning Prediction of Interspecies Interactions from Self-organized Spatiotemporal Patterns of Co-evolving Organisms

Microorganisms colonizing natural habits such as soils co-evolve to form specific spatial patterns through interspecies interactions. These self-organized patterns are a key ecological phenotype, which provides critical information on their interaction mechanisms. However, conventional network inference techniques that analyze species population data in bulk have yet to be extended to account for such spatial heterogeneity. Here we proposed supervised deep learning as a new network inference tool for predicting interspecies interactions from spatiotemporal patterns of microbial evolution. Due to lack of biological imaging data that can be used for training deep learning networks, we used in silico data generated from high-fidelity agent-based models to determine model structure and parameters. Even though networks were trained under simple configurations where interaction coefficients are assumed to be spatially invariant, we demonstrated that the resulting model can be utilized to successfully predict spatial variation of interactions in more complex domains (i.e., configured with a context-dependent mixture of interaction coefficients) as well as in simple domains without further training. In the further test against real biological data obtained through imaging experiments of a binary consortium (Pseudomonas fluorescens and a mutant of Escherichia coli), our model also predicted the dramatic shifts in interactions of the two organisms across different environmental contexts. Through various successful demonstrations in this work, the combined use of the agent-based model and machine learning algorithm provides a means to use new type of data - microscopic images - for extracting microbial interactions, therefore presenting itself as a useful tool for the analysis of more complex microbial community interactions.

Lee, Joon-Yong↗
DOE OSTI2026

Multipoint Aerostructural Optimization of Wind Turbine Rotors Using a Coupled Blade‐Resolved Aerostructural Solver

Physics‐based design optimization workflows thread the needle between computational cost limitations and simulation complexity, often compromising between modeling detail and the range of operating design conditions. Multipoint aerostructural optimization of wind turbine rotors has so far been confined to low‐fidelity analyses or to high‐fidelity studies with simplified structural models, leaving the most complex design trade‐offs unexplored. We close this gap by performing the first tightly coupled gradient‐based multipoint aerostructural rotor optimization using 3D aerodynamic and structural solvers with discrete coupled adjoints. The optimizer simultaneously varies blade planform, airfoil shapes, and structural thickness through more than 270 design variables, minimizing a weighted combination of rotor mass and power across multiple wind speeds. Applied to a modified DTU 10‐MW benchmark under conservative structural and aerodynamic constraints, our multipoint optimization reduces rotor mass by up to 36% and increases power by 12%–15% across the main operating conditions; biasing the objective toward power yields power gains up to 18% and a 17% mass reduction. For a nominal wind distribution, 3‐point rotor designs accounting for low RPM and high thrust conditions capture dominant trade‐offs and outperform single‐point designs. Adding two off‐design points changes individual‐condition power by less than 3% but leaves the weighted average within 0.5%, and the mass‐power bias has a stronger effect on the final design than the operating‐point weighting itself. Our framework extends naturally to richer load cases and site‐specific wind distributions, providing a basis for high‐fidelity multipoint design earlier in industrial workflows.

17 WIND ENERGY↗
DOE OSTI2021

Quantifying Turbine-Level Risk to Golden Eagles Using a High-Fidelity Updraft Model and a Stochastic Behavioral Model

To minimize the effects of wind farms on Golden Eagle (Aquila chrysaetos) populations while enabling sustainable development of renewable energy resources, it is important to understand how eagles interact with atmospheric flows, terrain features, and anthropogenic structures. Models that predict migratory flight paths provide one tool that helps us grasp how the location of wind farms may influence interactions and impacts on migrating Golden Eagles. The current state-of-the-art in predicting migratory flight paths uses a deterministic fluid-flow analogy to predict eagle trajectory using only an orographic updraft potential computed from topographical features. This model does not take into account variables, such as thermal updrafts and time varying atmospheric conditions that are known to influence migratory behavior. In this work, we improve on the model with the objective of developing tools that advance our understanding of how atmospheric flows and terrain features affect migratory eagle behavior and their interactions with wind farms. Specifically, we 1) incorporate both orographic and thermal updraft information in simulating eagle flight paths; 2) incorporate stochasticity into eagle travel patterns to better capture the influence of exogenous factors on, and the inherent stochasticity of eagle behavior; 3) consider spatio-temporal atmospheric data at wind-farm-scale when computing updraft potential; and 4) account for how atmospheric conditions and the direction of migration change seasonally and how these changes affect eagle migratory flight behavior. We tested the model using a 50km by 50km region with 50 m resolution in the western United States. We simulated 900 independent, probabilistic eagle tracks during southerly and northerly migration, assuming eagles solely rely on orographic updrafts. The preliminary results indicate that the inclusion of finer resolution atmospheric data allows for the inclusion of realistic conditions that an eagle experiences. The stochasticity in eagle tracks provides a platform to include uncertainty in eagle decision making and help produce robust eagle presence maps. We will deploy updraft and downdraft velocities computed using a high-fidelity, wind farm scale, computational fluid dynamics solver under development at National Renewable Energy Laboratory. This work is a first step in the development of a predictive and generalizable eagle behavior model at the wind farm scale that does not rely on empirical data collection. Although the current model is intended for migratory eagles, we will extend and refine this model to inform the development of additional behavioral modes, including resident eagle behavior. This modeling approach improves our ability to understand eagle use of the landscape at a fine scale, and it is our hope that this work will ultimately help advance strategies that minimize the impact of wind development on Golden Eagle populations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

On the Representation of Hyporheic Exchange in Models for Reactive Transport in Stream and River Corridors

Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Deformation, dislocation evolution and the non-Schmid effect in body-centered-cubic single- and polycrystal tantalum

A physically-informed continuum crystal plasticity model is presented to elucidate deformation mechanisms, dislocation evolution and the non-Schmid effect in body-centered-cubic (bcc) tantalum widely used as a key structural material for mechanical and thermal extremes. We show the unified structural modeling framework informed by mesoscopic dislocation dynamics simulations is capable of capturing salient features of the large inelastic behavior of tantalum at quasi-static (10 -3 s -1 ) to extreme strain rates (5000 s -1 ) and at low (77 K) to high temperatures (873 K) at both single- and polycrystal levels. Here we also present predictive capabilities of the model for microstructural evolution in the material. To this end, we investigate the effects of dislocation interactions on slip activities, instability and the non-Schmid behavior at the single crystal level. Furthermore, ex situ measurements on crystallographic texture evolution and dislocation density growth are carried out for polycrystal tantalum specimens at increasing strains. Numerical simulation results also support that the modeling framework is capable of capturing the main features of the polycrystal behavior over a wide range of strains, strain rates and temperatures. The theoretical, experimental and numerical results at both single- and polycrystal levels provide critical insight into the underlying physical pictures for micro- and macroscopic responses and their relations in this important class of refractory bcc materials undergoing large inelastic deformations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Octahedral oxide glass network in ambient pressure neodymium titanate

Abstract Rare-earth titanates form very fragile liquids that can be made into glasses with useful optical properties. We investigate the atomic structure of 83TiO 2 -17Nd 2 O 3 glass using pair distribution function (PDF) analysis of X-ray and neutron diffraction with double isotope substitutions for both Ti and Nd. Six total structure factors are analyzed (5 neutron + 1 X-ray) to obtain complementary sensitivities to O and Ti/Nd scattering, and an empirical potential structure refinement (EPSR) provides a structural model consistent with the experimental measurements. Glass density is estimated as 4.72(13) g cm −3 , consistent with direct measurements. The EPSR model indicates nearest neighbor interactions for Ti-O at $$\overline{r}_{TiO}$$ r ¯ TiO = 1.984(11) Å with coordination of $$n_{TiO}$$ n TiO = 5.72(6) and for Nd-O at $$\overline{r}_{NdO}$$ r ¯ NdO = 2.598(22) Å with coordination of $$n_{NdO}$$ n NdO = 7.70(26), in reasonable agreement with neutron first order difference functions for Ti and Nd. The titanate glass network comprises a mixture of distorted Ti-O 5 and Ti-O 6 polyhedra connected via 71% corner-sharing and 23% edge-sharing. The O-Ti coordination environments include 15% nonbridging O-Ti 1 , 51% bridging O-Ti 2 , and 32% tricluster O-Ti 3 . This structure is highly unusual for oxide glasses melt-quenched at ambient pressure, as it consists of Ti-O x predominantly in octahedral (with nearly no tetrahedral) coordination.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Characterizing Performance of Freshwater Wetland Methane Models Across Time Scales at FLUXNET-CH 4 Sites Using Wavelet Analyses

Process-based land surface models are important tools for estimating global wetland methane (CH 4 ) emissions and projecting their behavior across space and time. So far there are no performance assessments of model responses to drivers at multiple time scales. In this study, we apply wavelet analysis to identify the dominant time scales contributing to model uncertainty in the frequency domain. We evaluate seven wetland models at 23 eddy covariance tower sites. Our study first characterizes site-level patterns of freshwater wetland CH 4 fluxes (FCH 4 ) at different time scales. A Monte Carlo approach was developed to incorporate flux observation error to avoid misidentification of the time scales that dominate model error. Our results suggest that (a) significant model-observation disagreements are mainly at multi-day time scales (<15 days); (b) most of the models can capture the CH 4 variability at monthly and seasonal time scales (>32 days) for the boreal and Arctic tundra wetland sites but have significant bias in variability at seasonal time scales for temperate and tropical/subtropical sites; (c) model errors exhibit increasing power spectrum as time scale increases, indicating that biases at time scales <5 days could contribute to persistent systematic biases on longer time scales; and (d) differences in error pattern are related to model structure (e.g., proxy of CH 4 production). Our evaluation suggests the need to accurately replicate FCH 4 variability, especially at short time scales, in future wetland CH 4 model developments.

58 GEOSCIENCES↗
DOE OSTI2023

Desmearing small-angle scattering data by central moment expansions

This article proposes a robust method for desmearing 1D small-angle scattering (SAS) intensity profiles. Using the central-moment-expansion technique, smearing in SAS intensity profiles is achieved in a model-free manner, eliminating the need for presumptive input during the desmearing. This innovative approach enables data analysis of SAS data without the resolution being convoluted with the structural model, thereby facilitating faster implementation and performance of fitting functions directly from the scattering models.

42 ENGINEERING↗
DOE OSTI2021

Physics-Informed Machine Learning from Observations for Clouds, Convection, and Precipitation Parameterizations and Analysis

Fusing learning from Earth observations from space and from the ground (e.g., ARM) with newly developed interpretable and generalizable physics-informed model structures to improve the parameterization of clouds and convection and advance the simulation, understanding and analysis of hydrological extreme events. This covers the prescribed foci of (2) predictive modeling using a hierarchy of models and (3) insight gleaned from complex data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Intrusive Uncertainty Quantification and Optimal Experiment Design in the Open-Source Pyomo Ecosystem

This contribution describes ParmEst and Pyomo.DoE, two pillars of the open-source Python-based Pyomo ecosystem for computational optimization with (partial differential) algebraic equation mathematical models. Specifically, ParmEst facilitates intrusive frequentist parameter estimation (PE) and uncertainty quantification (UQ) through built-in features, such as covariance matrix estimation, bootstrapping, and likelihood ratio tests. Complementary, Pyomo.DoE enables optimal experiment design by maximizing various metrics of the Fisher information matrix, such as A-optimality (trace), D-optimality (determinant), E-optimality (minimum eigenvalue), and ME-optimality (condition number). ParmEst and Pyomo.DoE can solve high-dimensional optimization problems by leveraging the model structure and exact derivative information. Finally, we will discuss future opportunities to integrate PE and UQ capabilities with optimization under uncertainty, including robust optimization with non-convex models via PyROS.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Vaterite Optical Petrography in Lake Sturgeon Otoliths

Vaterite occurring in fish otoliths exhibits plumose and spherulitic habits, the latter being like those grown from aqueous solutions. Vaterite spherulites are known to have sheaf-like texture and can be up to 100 μm in size. In thin section, the spherulites typically show uniaxial positive interference figures between cross-polarizers using a polarizing light microscope but without conoscopic or Bertrand lenses because the spherulites mimic the effect of these lenses. Plumose vaterite can show both uniaxial (+) and biaxial (+) interference figures with a small 2V, which suggests that crystal structural models for vaterite should have the plane of the carbonate group parallel to the optic axis or acute bisectrix. Finally, vaterite spherulites grown from aqueous solution exhibit similar structure and growth habits.

58 GEOSCIENCES↗
DOE OSTI2020

Are we modeling the properties of unconventional shales correctly?

Technological advancements have led to impressive growth in hydrocarbon (HC) production from unconventional shale reservoirs. However, major limiting factors in shale gas development are i) low recovery of gas in place (<20%), (ii) a rapid decline in well productivity, iii) release of contaminants, despite the use of advanced hydraulic fracturing fluids and multiple fracturing stages. This is mainly due to a lack of understanding of the nanoscale physicochemical properties of shales, especially that of kerogen, the macromolecule that is not only the source of the majority of the HC’s in shales but also holds most of the HC’s in adsorbed state. Over the last few years, a large number of studies have been published on simulating the physicochemical properties of shale utilizing molecular models of kerogen. However, the molecular models of kerogen input into these simulations are based on the kerogen “type” derived from a very limited number of shale samples. In this paper, we examine the variations that can exist in kerogen structure within a particular kerogen “type” across different shale basins and single shale basin at similar thermal maturity levels. We propose that using kerogen “type” based structural models for molecular simulations could lead to inaccurate estimation of HC reserves, HC recovery, HC production, frackability, and quantity and quality of produced water. Our study highlights the need for developing a better classification of kerogen based on its “molecular structure” instead of “type” for a more accurate prediction of physicochemical properties of shales.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2022

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

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

Convergent close-coupling calculations of electron scattering on LiH

Cross sections for electron scattering on the ground state of LiH are calculated using the molecular convergent close-coupling method. The fixed-nuclei approximation is utilized, and calculations are performed at the mean internuclear separation of 3.06 𝑎 0 for projectile energies ranging from 0.1 to 500 eV. Here, a model-potential approach is utilized to treat the LiH molecule as a two-electron system, and good agreement is found with the literature for various properties of the structure model. Cross sections are presented for elastic scattering, excitation to several electronic states, and ionization. Comparison is made with existing results, and qualitative agreement is found with small close-coupling calculations at low projectile energies. Good agreement is found with first-order results at high projectile energies.

74 ATOMIC AND MOLECULAR PHYSICS↗