Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

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

Hurricane Eyewall Winds and Structural Response of Wind Turbines

This paper describes the analysis of a wind turbine and support structure subject to simulated hurricane wind fields. The hurricane wind fields, which result from a Large Eddy Simulation of a hurricane, exhibit features such as very high gust factors (> 1.7), rapid direction changes (30 degrees in 30 sec) and substantial veer. Wind fields including these features have not previously been used in an analysis of a wind turbine and their effect on structural loads may be an important driver of enhanced design considerations. With a focus on blade root loads and tower base loads, the simulations show that these features of hurricane wind fields can lead to loads that are substantially in excess of those that would be predicted were wind fields with equally high mean wind speeds but without associated direction change and veer used in the analysis. This result, if further verified for a range of hurricane and tropical storm simulations, should provide an impetus for revisiting design standards.

17 WIND ENERGY↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part II: linesweep ALE rezoning for nonconformal block-structured AMR meshes

The simulation of flows presenting contact discontinuities, vorticity, and large variations in spatial scales can be performed in a framework coupling Arbitrary Lagrangian Eulerian (ALE) algorithms and Adaptive Mesh Refinement (AMR). This coupling requires adaptation of ALE rezoning techniques to meshes containing nonconformal nodes arising from both the AMR topology and the junction of mesh blocks. Here, in this paper, we present an ALE rezoning strategy that is compatible with such meshes, and that can also act as a disentangling algorithm. Emphasis is put on an algorithm that respects intrinsic Lagrangian mesh properties in order to preserve accuracy around discontinuities. To that end, we adapt the weighted linesweep algorithm to nonconformal block-structured AMR meshes. Then, we present control parameters introduced in the method for it to be applicable in practical situations. Notably, the method is coupled to a specific metric optimization in order to palliate some shortcomings of the linesweep method. Finally, numerical test cases are presented that feature the capabilities of the ALE-AMR algorithm for flows that present discontinuities, vorticity, and a variety of scales. Notably, we show that our ALE-AMR algorithm gives results at least similar to Euler-AMR, but provides better accuracy in cases where discontinuities are involved, thanks to a method that respects the Lagrangian features of the mesh. Additionally, it enables Euler-AMR-like computations on domains with temporally varying domain boundaries.

Adaptive mesh refinement↗
DOE OSTI2021

PUMIPic: A mesh-based approach to unstructured mesh Particle-In-Cell on GPUs

Unstructured mesh particle-in-cell, PIC, simulations executing on the current and next generation of massively parallel systems require new methods for both the mesh and particles to achieve performance and scalability on GPUs. The traditional approach to implementing PIC simulations defines data structures and algorithms in terms of particles with a full copy of the unstructured mesh on every process. To effectively scale the unstructured mesh and particles, mesh-based PIC uses the unstructured mesh as the predominant data structure with the particles stored in terms of the mesh entities. Here, this paper details the PUMIPic library, a framework for developing efficient and performance-portable mesh-based PIC simulations on GPU systems. A pseudo physics simulation based on a five-dimensional gyro-kinetic code for modeling plasma physics is used to examine the performance of PUMIPic. Scaling studies of the unstructured mesh partition and number of particles are performed up to 4096 nodes of the Summit system at Oak Ridge National Laboratory. The studies show that mesh-based PIC can utilize a partitioned mesh and maintain scaling up to system limitations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Assessment of MOOSE-Based Tools for Calculating Radial Core Expansion

Radial core expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. An inherently safe reactor design takes advantage of negative reactivity feedback in accident scenarios by utilizing a core restraint system which produces a bowed shape that allows for radial expansion of the fuel regions. Detailed modeling and simulation of core radial expansion itself as well as subsequent reactivity feedback is a challenging task involving contact of many fuel assembly elements and physics feedback from neutronics, thermal hydraulics, and thermal mechanical response. A variety of physics codes have been developed to model aspects of radial core expansion but in general invoke geometrical or physics approximations. No code system currently exists which tightly and robustly couples these physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A high-fidelity code will also be used to benchmark existing lower fidelity, faster running models to understand their benefits, limitations, and range of applications. A code development path using MOOSE-based tools is proposed in order to leverage the detailed geometry capabilities and natural tight coupling and robustness of MOOSE-based applications for modeling this complex phenomena. While simulation of the full phenomenon involves several physics, an assessment has been initiated on the capabilities and readiness of the currently available Tensor Mechanics module within MOOSE for calculation of the structural mechanical responses which occur within the reactor core. This report focuses on modeling the force-deformation response which mimics the physics of a duct contact deformation, as well as differential thermal expansion which produces a thermal bowed shaped for the fuel assemblies in a core. Simple examples were initially performed such as simple supported beam bending under load. The complexity of examples was progressively increased to better mimic the duct behavior by including a differential thermal example and inclusion of hexagonal cross-sections in the geometry. Further assessment of the structural mechanical response simulation capability is still required for modeling duct contact interactions and irradiation creep and swelling. Companion thermal hydraulic and neutronics assessments will also be required; these activities are planned for future years. Finally, integration of the multiple physics components through MOOSE is required to predict the core radial expansion and subsequent feedback effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Real-space inversion and super-resolution of ultrafast scattering

Ultrafast scattering using x-rays or electrons is an emerging method to obtain structure dynamics at the atomic lengthscales and timescales. However, directly resolving in real-space atomic motions is inherently limited by the finite detector range and the probe energy. So, as a result, the time-resolved signal interpretation is mostly done in reciprocal space and relies on modeling and simulations of specific structures and processes. Here, we introduce a model-free approach to directly resolve scattering signals in real space, surpassing the diffraction limit, using scattering kernels and signal priors that naturally arise from the measurement constraints. We demonstrate the approach on simulated and experimental data, recover multiple atomic motions at sub-angstrom resolutions, and discuss the recovery accuracy and resolution limits versus signal fidelity. The approach offers a robust path to obtain high-resolution real-space information of atomic-scale structure dynamics using current time-resolved x-ray or electron scattering sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Advancing coupled fire-atmosphere research with HIGRAD/FIRETEC (w20_firetec) (CY 2021 Institutional Computing Annual Progress Report) [Slides]

LANL’s HIGRAD and FIRETEC models were used in a multi-scale simulation framework to investigate impacts of small-scale fuel heterogeneities on large-scale Pyro-Cb formation. Forest structural diversity is found to increase canopy fire spread, and thus risk of tree mortality in FIRETEC simulations. Structural diversity is associated with an increase in ladder fuels, which facilitate fire spread into the canopy. Low-intensity fire burning in marginal conditions is highly sensitive to small perturbation in atmospheric turbulence. These results highlight predictive limitations of individual fire simulations, and emphasize the need for ensemble simulations, especially for low-intensity, prescribed fire.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Zinc binding alters the conformational dynamics and drives the transport cycle of the cation diffusion facilitator YiiP

YiiP is a secondary transporter that couples Zn 2+ transport to the proton motive force. Structural studies of YiiP from prokaryotes and Znt8 from humans have revealed three different Zn 2+ sites and a conserved homodimeric architecture. These structures define the inward-facing and outward-facing states that characterize the archetypal alternating access mechanism of transport. To study the effects of Zn 2+ binding on the conformational transition, we use cryo-EM together with molecular dynamics simulation to compare structures of YiiP from Shewanella oneidensis in the presence and absence of Zn 2+ . To enable single-particle cryo-EM, we used a phage-display library to develop a Fab antibody fragment with high affinity for YiiP, thus producing a YiiP/Fab complex. To perform MD simulations, we developed a nonbonded dummy model for Zn 2+ and validated its performance with known Zn 2+ -binding proteins. Using these tools, we find that, in the presence of Zn 2+ , YiiP adopts an inward-facing conformation consistent with that previously seen in tubular crystals. After removal of Zn 2+ with high-affinity chelators, YiiP exhibits enhanced flexibility and adopts a novel conformation that appears to be intermediate between inward-facing and outward-facing states. This conformation involves closure of a hydrophobic gate that has been postulated to control access to the primary transport site. Comparison of several independent cryo-EM maps suggests that the transition from the inward-facing state is controlled by occupancy of a secondary Zn 2+ site at the cytoplasmic membrane interface. This work enhances our understanding of individual Zn 2+ binding sites and their role in the conformational dynamics that govern the transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Modeling of Sagging for 3D Printed Layers During the Curing Process

Here, a printed layer of silica-enforced poly(dimethylsiloxane)-co-(diphenylsiloxane) is modeled as a two-phase system consisting of air and polymer with an interface set up during the printing process. The structural geometry changes mostly due to the action of surface tension, while all material properties are strongly temperature dependent. Polymer flow is described using equations of the extended Herschel–Bulkley model, with parameters strongly dependent on temperature and degree of curing. Parameters of the model are determined using flow sweep measurements and separate experiments with vertical structure sagging at different temperatures. The curing process is modeled using dependencies between the curing rate, degree of curing, and temperature obtained in studies by differential scanning calorimetry. The developed model is used for simulations of printed structure deformations with different initial and boundary conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Deflection predictions of involute-shaped fuel plates using a fully-coupled numerical approach

This paper describes the modeling and simulation of fluid structure interactions (FSI) of involute-shaped fuel plates used in nuclear research reactors. We believe this to be the first time that this type of application is described in the literature using a fully-coupled, and monolithic, finite element approach. The simulations are validated against plate deflection data for the conceptual design of the Advanced Neutron Source Reactor (ANSR), which was envisioned to be the world’s most powerful nuclear research reactor for neutron scattering and other applications, but was ultimately never completed. The high performance of the ANSR creates a bounding envelope for involute-shaped research reactors such as that used in the High Flux Isotope Reactor (HFIR) at the Oak Ridge National Laboratory (ORNL) which is undergoing research for the conversion from highly-enriched uranium (HEU) to low-enriched uranium (LEU) fuel. As such, the findings from the present FSI analyses carried out herein for the ANSR plates provide good guidelines and inform designers what should be expected for the next generation of plates in the HFIR. Furthermore it is shown herein that the current approach can accurately capture the leading-edge deflections of the involute-shaped plates and simulations can predict the ‘S-shaped’ deflection of the first mode instilling confidence in the methodology.

42 ENGINEERING↗
DOE OSTI2022

Hybrid particle-in-cell simulations of electromagnetic coupling and waves from streaming burst debris

Various systems can be modeled as a point-like explosion of ionized debris into a magnetized, collisionless background plasma-including astrophysical examples, active experiments in space, and laser-driven laboratory experiments. Debris streaming from the explosion parallel to the magnetic field may drive multiple resonant and non-resonant ion-ion beam instabilities, some of which can efficiently couple the debris energy to the background and may even support the formation of shocks. We present a large-scale hybrid (kinetic ions + fluid electrons) particle-in-cell simulation, extending hundreds of ion inertial lengths from a 3D explosion, that resolves these instabilities. We show that the character of these instabilities differs notably from the 1D equivalent by the presence of unique transverse structure. Additional 2D simulations explore how the debris beam length, width, density, and speed affect debris-background coupling, with implications for the generation of quasi-parallel shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Correlated Anion Disorder in Heteroanionic Cubic TiOF 2

Resolving anion configurations in heteroanionic materials is crucial for understanding and controlling their properties. For anion-disordered oxyfluorides, conventional Bragg diffraction cannot fully resolve the anionic structure, necessitating alternative structure determination methods. We have investigated the anionic structure of anion-disordered cubic (ReO 3 -type) TiOF 2 using X-ray pair distribution function (PDF), 19 F MAS NMR analysis, density functional theory (DFT), cluster expansion modeling, and genetic-algorithm structure prediction. Our computational data predict short-range anion ordering in TiOF 2 , characterized by predominant cis-[O 2 F 4 ] titanium coordination, resulting in correlated anion disorder at longer ranges. To validate our predictions, we generated partially disordered supercells using genetic-algorithm structure prediction and computed simulated X-ray PDF data and 19 F MAS NMR spectra, which we compared directly to experimental data. To construct our simulated 19 F NMR spectra, we derived new transformation functions for mapping calculated magnetic shieldings to predicted magnetic chemical shifts in titanium (oxy)fluorides, obtained by fitting DFT-calculated magnetic shieldings to previously published experimental chemical shift data for TiF 4 . We find good agreement between our simulated and experimental data, which supports our computationally predicted structural model and demonstrates the effectiveness of complementary experimental and computational techniques in resolving anionic structure in anion-disordered oxyfluorides. From additional DFT calculations, we predict that increasing anion disorder makes lithium intercalation more favorable by, on average, up to 2 eV, highlighting the significant effect of variations in short-range order on the intercalation properties of anion-disordered materials.

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

Small-molecule modulation of β-arrestins

β-Arrestins are multifunctional regulators of G-protein-coupled receptor (GPCR) signalling and orchestrate diverse downstream signalling events and physiological responses across the GPCR superfamily. Although GPCR pharmacology has advanced to target orthosteric and allosteric sites, as well as G proteins and GPCR kinases, direct chemical tools to modulate β-arrestin activities have remained conspicuously absent. Here we report the identification of small-molecule inhibitors that selectively target β-arrestins and delineate their mechanism of action through integrated pharmacological, biochemical, biophysical and structural analyses. These inhibitors disrupt β-arrestin engagement with agonist-activated GPCRs, impairing desensitization, internalization and β-arrestin-dependent physiological functions while sparing G protein–receptor coupling. Cryo-electron microscopy, molecular dynamics simulations and structure-guided mutagenesis reveal that one modulator, Cmpd-5, engages a pocket within the central crest of β-arrestin1 formed by the middle, C and lariat loops, a critical receptor-binding interface, stabilizing a distinct conformation that is incompatible with full β-arrestin–receptor engagement. Together, these findings establish a mechanistic framework for β-arrestin modulation, reveal a novel allosteric site for structure-based drug design, and open new avenues for transducer-targeted, pathway-specific GPCR therapeutic agents.

Kahsai, Alem W. [Duke University, Durham, NC (Unit↗
DOE OSTI2021

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗
DOE OSTI2021

Simulation of enhanced light extraction from periodic, disordered, and quasi-periodic OLED structures

Here we have developed a rigorous scattering matrix theory of light emission from periodically structured media using a Green’s function approach. We computationally simulate the spectral power inside the structure, incorporating Purcell factor enhancements, and find internal waveguiding modes and plasmonic losses. We simulate the light-outcoupling factor η out and describe how corrugations and structured media can enhance η out . We have extended our framework to describe non-periodic disordered arrays as well as aperiodic quasi-crystalline arrays. Flat organic light-emitting diodes (OLEDs) have low outcoupling with η out ~20% since most of the light is trapped in waveguided modes in higher index layers or lost to plasmonic excitations at the metal cathode. Periodically corrugated OLEDs can achieve highly enhanced η out ~65%-70% for pitch values between 1000 and 2000 nm, representing an enhancement factor >3 over planar structures. Periodic corrugations strongly diffract trapped waveguided and plasmonic modes to the emissive air cone. Disordered templates also can cause significant enhanced outcoupling with η out ~50%-55% for smaller nearest-neighbor separations of 300–400 nm. Quasi-crystalline tilings can also lead to enhancements of η out ~50%-55%. This framework can be utilized to design novel structured media that can generate high light extraction.

42 ENGINEERING↗
DOE OSTI2019

Inter-subunit coupling enables fast CO 2 -fixation by reductive carboxylases

Enoyl-CoA carboxylases/reductases (ECRs) belong to the most efficient CO 2 -fixing enzymes described to date. However, the molecular mechanisms underlying ECR’s extraordinary catalytic activity on the level of the protein assembly remain elusive. Here we used a combination of ambient temperature X-ray Free Electron Laser (XFEL) and cryogenic synchrotron experiments to study the structural organization of the ECR from Kitasatospora setae. K. setae ECR is a homo-tetramer that differentiates into a dimer of dimers of open- and closed-form subunits in the catalytically active state. Using molecular dynamics simulations and structure-based mutagenesis, we show that catalysis is synchronized in K. setae ECR across the pair of two dimers. This conformational coupling of catalytic domains is conferred by individual amino acids to achieve high CO 2 -fixation rates. Our results provide unprecedented insights into the dynamic organization and synchronized inter- and intra-subunit communications of this remarkably efficient CO 2 -fixing enzyme during catalysis.

DeMirci, Hasan↗
DOE OSTI2020

A fluid-structure interaction solver for investigating torsional galloping in solar-tracking photovoltaic panel arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. Furthermore, these simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects. These results are supplemented by a discussion of stability trends noted during testing and possible implications when considering multi-row array interactions.

14 SOLAR ENERGY↗
DOE OSTI2020

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗