SEARCH · Engineering Papers
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.
The poster will be presented at the American Geophysical Union Fall Meeting, during the session on Developments in Earthquake Ground Motion Modeling and Advancements toward Fault-to-Structure Simulations. https://www.agu.org/Fall-Meeting
Higher plants synthesize cellulose using membrane-bound, six-lobed cellulose synthase complexes, each lobe containing trimeric cellulose synthases (CESAs). Although molecular biology reports support heteromeric trimers composed of different isoforms, a homomeric trimer was reported for in vitro studies of the catalytic domain of CESA1 of Arabidopsis (AtCESA1CatD) and confirmed in cryoEM structures of full-length CESA8 and CESA7 of poplar and cotton, respectively. In both structures, a small portion of the plant-conserved region (P-CR) forms the only contacts between catalytic domains of the monomers. Here, we report inter-subunit lysine-crosslinks that localize to the small P-CR, negative-stain EM structure, and modeling data for homotrimers of AtCESA1CatD. Molecular dynamics simulations for AtCESA1CatD trimers based on the CESA8 cryoEM structure were stable and dependent upon a small set of residue contacts. The results suggest that homomeric CESA trimers may be important for the synthesis of primary and secondary cell walls and identify key residues for future mutagenic studies.
Quantitative structure–odor relationships are critically important for studies related to the function of olfaction. Current literature data sets contain expert-labeled molecules but lack feature data. This paper introduces QuantumScents, a quantum mechanics augmented derivative of the Leffingwell data set. QuantumScents contains 3.5k structurally and chemically diverse molecules ranging from 2 to 30 heavy atoms (CNOS) and their corresponding 3D coordinates, total PBE0 energy, molecular dipole moment, and per-atom Hirshfeld charges, dipoles, and ratios. The authors demonstrate that Hirshfeld charges and ratios contain sufficient information to perform molecular classification by training a Message Passing Neural Network with chemprop (Heid, E.; et al. ChemRxiv, 2023, DOI: 10.26434/chemrxiv-2023-3zcfl) to predict scent labels. Finally, the QuantumScents data set is freely available on Zenodo along with the authors’ code, example models, and data set generation workflow (https://zenodo.org/doi/10.5281/zenodo.8239853).
The goal of this project is to model and control a novel unidirectional cascaded multilevel bridgeless rectifier as an active front end in medium and high voltage applications. This topology has many advantages over a conventional cascaded H-bridge rectifier, such as lower implementation cost, higher reliability, and greater flexibility with similar power quality. The complete design process of the proposed converter is developed step by step in order to meet all the desired objectives. The steady-state mathematical model is used to develop a method for the voltage balancing of dc cells. Power factor analysis is discussed to mathematically derive requirements for the number of partially controlled and fully controlled H-bridges in the proposed H-bridge converter. Power loss, efficiency, and cost comparison studies between the traditional cascaded H-Bridge converter and the proposed bridgeless converter demonstrate the advantages. After exploring various well-established control methods, a novel control strategy is proposed to achieve dc voltage balancing, fast and robust grid synchronization, power factor correction, and elimination of zero crossing current distortion under both balanced and unbalanced load conditions. The converter can also be used for reactive power compensation in a grid tied power system if a sufficient number of fully controlled H-bridge modules are included. Processor-In-the-Loop (PIL) simulation has been the utilized to validate the performance of discrete control structure. Simulation and experimental results validate the models and control method.
Abstract This paper presents the development of a multi‐fidelity digital twin structural model (virtual model) of an as‐built wind turbine blade. The goal is to develop and demonstrate an approach to produce an accurate and detailed model of the as‐built blade for use in verifying the performance of the operating two‐bladed, downwind rotor. The digital twin model development methodology, presented herein, involves a novel calibration process to integrate a wide range of information including design specifications, manufacturing information, and structural testing data (modal and static) to produce a multi‐fidelity digital twin structural model: a detailed high‐fidelity model (i.e., 3D finite element analysis [FEA]) and consistent beam‐type models for aeroelastic simulation. A key element is that the multi‐fidelity structural digital twin method follows the rotor from the stages of design, to manufacturing, then to the ground testing and field operation. The result of this comprehensive approach is an accurate multi‐fidelity digital twin structural model for the geometric, structural, and structural dynamic properties of the as‐built blade within a 1% match in mass properties, 3.2% in blade frequencies, and 6% in deflection. The different stages of processing this information within the methodology are discussed. The rotor examined is the SUMR‐Demonstrator (SUMR‐D), which was installed on the Controls Advanced Research Testbed (CART‐2) wind turbine at the National Wind Technology Center. The digital twin model developed here was utilized to design controllers to safely operate SUMR‐D in field tests, which are providing additional data for further evaluation and development of the multi‐fidelity digital twin structural model.
The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.
The vertebrate inner ear, responsible for hearing and balance, is able to sense minute mechanical stimuli originating from an extraordinarily broad range of sound frequencies and intensities or from head movements. Integral to these processes is the tip-link protein complex, which conveys force to open the inner-ear transduction channels that mediate sensory perception. Protocadherin-15 and cadherin-23, two atypically large cadherins with 11 and 27 extracellular cadherin (EC) repeats, are involved in deafness and balance disorders and assemble as parallel homodimers that interact to form the tip link. Here we report the X-ray crystal structure of a protocadherin-15 + cadherin-23 heterotetrameric complex at 2.9-Å resolution, depicting a parallel homodimer of protocadherin-15 EC1-3 molecules forming an antiparallel complex with two cadherin-23 EC1-2 molecules. In addition, we report structures for 10 protocadherin-15 fragments used to build complete high-resolution models of the monomeric protocadherin-15 ectodomain. Molecular dynamics simulations and validated crystal contacts are used to propose models for the complete extracellular protocadherin-15 parallel homodimer and the tip-link bond. Steered molecular dynamics simulations of these models suggest conditions in which a structurally diverse and multimodal protocadherin-15 ectodomain can act as a stiff or soft gating spring. These results reveal the structural determinants of tip-link–mediated inner-ear sensory perception and elucidate protocadherin-15’s structural and adhesive properties relevant in disease.
Metal oxide/water interfaces play a crucial role in many electrochemical and photocatalytic processes, such as photoelectrochemical water splitting, the creation of fuel from sunlight, and electrochemical CO 2 reduction. First-principles electronic structure calculations can reveal unique insights into these processes, such as the role of the alignment of the oxide electronic energy levels with those of liquid water. An essential prerequisite for the success of such calculations is the ability to predict accurate structural models of these interfaces, which in turn requires careful experimental validation. Here we report a general, quantitative validation protocol for first-principles molecular dynamics simulations of oxide/aqueous interfaces. The approach makes direct comparisons of interfacial x-ray reflectivity (XR) signals from experimental measurements and those obtained from ab initio simulations with semilocal and van der Waals functionals. The protocol is demonstrated here for the case of the Al 2 O 3 (001)/water interface, one of the simplest oxide/water interfaces. We discuss the technical requirements needed for validation, including the choice of the density functional, the simulation cell size, and the optimal choice of the thermodynamic ensemble. Our results establish a general paradigm for the validation of structural models and interactions at solid/water interfaces derived from first-principles simulations. Furthermore, while there is qualitative agreement between the simulated structures and the experimental best-fit structure, direct comparisons of simulated and measured XR intensities show quantitative discrepancies that derive from both bulk regions (i.e., alumina and water) as well as the interfacial region, highlighting the need for accurate density functionals to properly describe interfacial interactions. Our results show that XR data are sensitive not only to the atomic structure (i.e., the atom locations) but also to the electron-density distributions in both the substrate and at the interface.
This project will ultimately combine observations from Raman lidar, Doppler lidar, AERI, eddy correlation measurements, and airborne and ground-based high spectral resolution lidar systems with many large eddy simulations (LES) of realistic boundary layer cases to study: The diurnal evolution of the boundary layer thermodynamic structure; The turbulent structure of the CBL and the entrainment zone; The ability of LES models to simulate the thermodynamic and turbulent structure of the CBL; How perturbations in the surface properties and CBL properties influence boundary layer clouds; The similarities and differences in the boundary layer and cloud interaction in the marine vs. continental boundary layer.
Additive manufacturing (AM) offers opportunities to produce complex part geometries not possible with conventional processing and in some cases even improve part performance. However, adoption has been slowed by difficulties assessing microstructure variability and there is no straightforward approach to relate processing to grain structure characteristics. In this study, datasets from AdditiveFOAM heat transport simulations of laser powder bed fusion (LPBF) are used to drive ExaCA simulations of grain structure. The GPU utilization of ExaCA and an algorithmic update for modeling melt pool overlap region solidification enabled rapid and parallel simulation across a wider range of process conditions than previously explored with cellular automata-based solidification models. A texture selection angle $θ_s$ is defined based on melt pool overlap geometry, and the range of $θ_s$ over which a commonly observed texture transition occurs in characterized AM builds was well-reproduced by ExaCA simulations over a wide range of melt pool shape, hatch spacing, and layer height. ExaCA simulations with 90 degree rotation of the scan direction on every other layer reproduced a number of trends from the AM literature including grain refinement, the dominance of layers with larger melt pools on the final grain structure, and the weakening or strengthening of texture depending on odd and even layer melt pool overlap geometry. EBSD data from a benchmark AM part is used to validate the simulated mechanism of a layer rotation-induced texture strengthening effect. Importantly, these results expand the understanding of the mechanisms for texture selection in alloys with cubic crystal symmetry and offer an approach to easily evaluate processing conditions. With this new understanding, these modeling tools will enable anticipation of previously unexpected variations in grain structure and target specific microstructures and properties.
Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing.
Microstructure largely dictates macroscopic material properties and is strongly affected by processing. Therefore, the simulation of microstructure evolution in response to thermal fields during processing is of significant interest within the computational materials science community. Additive manufacturing (AM) has emerged as a technique for producing complex geometries and unique microstructures. Yet, complex and rapid thermal cycles in AM pose computational challenges for existing microstructure models. This work proposes a discrete event inspired cellular automata (CA) approach, titled DECA, to accelerate simulation of grain structure evolution in AM. In contrast to conventional time-stepped CA models, this model directly solves the times capture events would take place allowing for stepping in events rather than time (a technique also found in the field of discrete-event simulation). In comparison to purely serial discrete-event models, DECA allows for temporary violation of the causality constraint, but detects and corrects these violations, leading to an emergent phenomenon dubbed causality rippling, in which previously calculated capture events are overwritten. The amount of repeated calculations, defined by the capture ratio, is taken as a measure of computational inefficiency, and the model parameters that affect this ratio are evaluated. The new DECA approach was found to be more computationally efficient than conventional time-stepped CA models while guaranteeing an accurate solution, which can only be achieved in the conventional models for vanishingly small time steps. Finally, opportunities for parallelization and scaling of the new approach are discussed.
This report describes a collaborative project to develop a validated, predictive model for the high temperature mechanical properties of a titanium-matrix, tungsten-carbide/cobalt-reinforced metal matrix composite. The modeling approach was to first develop a detailed, microstructural model linking the material structure and the interfacial debonding properties to the effective properties of the material. The project then completed a throughput simulation campaign to generate a large number of simulations for discrete microstructures and different debonding parameters. Finally, the project trained a fast, Gaussian process surrogate model against this simulation database to provide a quick model linking the material compositions, structure, and processing parameters to the resulting material properties. This model was validated against high temperature tensile test data on a few particular composite compositions. The tests validate the model predictions for ultimate tensile strength and uniform elongation/ductility, meaning the final surrogate model can now be used to tune the material composition and processing parameters to identify optimal composite compositions for particular applications.
Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa–Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa–Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same 'equation of state', which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs.
Many nations responded to the COVID-19 pandemic by restricting travel and other activities during 2020, resulting in temporarily reduced emissions of CO 2 , other greenhouse gases and ozone and aerosol precursors. We perform a coordinated Intercomparison, CovidMIP, of Earth System model simulations to assess the impact on climate of these emissions reductions. Eleven models performed multiple initial-condition ensembles to produce over 280 simulations spanning both initial condition and model structural uncertainty. We find model consensus on reduced aerosol amounts (particularly over East Asia) and associated increases in surface shortwave radiation levels. However, any impact on near-surface temperature or rainfall during 2020-2024 is extremely small and is not detectable in this initial analysis. Regional analyses on a finer scale, and closer attention to extremes (especially linked to changes in atmospheric composition and air quality) are required to test the impact of COVID-19-related emission reductions on near-term climate.
Abstract Floodplain lakes are widespread and ecologically important throughout tropical river systems, however data are rare that describe how temporal variations in hydrological, meteorological and optical conditions moderate stratification and mixing in these shallow lakes. Using time series measurements of meteorology and water‐column temperatures from 17 several day campaigns spanning two hydrological years in a representative Amazon floodplain lake, we calculated surface energy fluxes and thermal stratification, and applied and evaluated a 3‐dimensional hydrodynamic model. The model successfully simulated diel cycles in thermal structure characterized by buoyancy frequency, depth of the actively mixing layer, and other terms associated with the surface energy budget. Diurnal heating with strong stratification and nocturnal mixing were common; despite considerable heat loss at night, the strong stratification during the day meant that mixing only infrequently extended to the bottom at night. Simulations indicated that the diurnal thermocline up and downwelled creating lake‐wide differences in near‐surface temperatures and mixing depths. Infrequent full mixing creates conditions conducive to anoxia in these shallow lakes given their warm temperatures.
Protein structure prediction algorithms, such as AlphaFold, have accelerated protein design and advanced the understanding of the relationship between amino acid sequence and protein structure. However, these algorithms are limited in their ability to predict the structures of conformationally dynamic, intrinsically disordered, and stimuli-responsive proteins. To evaluate sequence-to-structure predictions of such challenging proteins, we explored a class of conformationally dynamic, repeats-in-toxin (RTX) proteins. RTX proteins adopt intrinsically disordered conformations in the absence of calcium and undergo reversible folding into β-roll structures upon binding to calcium. RTX proteins are characterized by tandem repeats of the sequence GGXGXDXUX, in which X can be any amino acid and U is an aliphatic amino acid. We designed RTX sequence variants with global substitutions of nonconserved amino acids, tandem repeats of consensus sequences GGAGXDTLY, and tandem repeats of scrambled sequences GGAGXDTYL. AlphaFold2 and AlphaFold3 predicted that all of these RTX variants adopt β-roll structures, characteristic of wild-type RTX bound to calcium. However, modeling the predicted structures with molecular dynamics simulations and characterizing the protein variants with circular dichroism spectroscopy, small-angle x-ray scattering, and x-ray crystallography revealed that variants adopt diverse, sequence-dependent structures in the absence and presence of calcium. To better design proteins for applications in biotechnology and sustainability, it is critical to build predictive tools that consider intrinsically disordered protein states and validate these tools with multi-mode, multi-scale experimental data.
Improper lengths of actin-thin filaments are associated with altered contractile activity and lethal myopathies. Leiomodin, a member of the tropomodulin family of proteins, is critical in thin filament assembly and maintenance; however, its role is under dispute. Using nuclear magnetic resonance data and molecular dynamics simulations, we generated the first atomic structural model of the binding interface between the tropomyosin-binding site of cardiac leiomodin and the N-terminus of striated muscle tropomyosin. Our structural data indicate that the leiomodin/tropomyosin complex only forms at the pointed end of thin filaments, where the tropomyosin N-terminus is not blocked by an adjacent tropomyosin protomer. This discovery provides evidence supporting the debated mechanism where leiomodin and tropomodulin regulate thin filament lengths by competing for thin filament binding. Data from experiments performed in cardiomyocytes provide additional support for the competition model; specifically, expression of a leiomodin mutant that is unable to interact with tropomyosin fails to displace tropomodulin at thin filament pointed ends and fails to elongate thin filaments. Together with previous structural and biochemical data, we now propose a molecular mechanism of actin polymerization at the pointed end in the presence of bound leiomodin. In the proposed model, the N-terminal actin-binding site of leiomodin can act as a “swinging gate” allowing limited actin polymerization, thus making leiomodin a leaky pointed-end cap. Results presented in this work answer long-standing questions about the role of leiomodin in thin filament length regulation and maintenance.