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 451 records · Page 25

Design of High-deflection Foils for Marine Hydrokinetic Applications (Final Report)

This project used model-scale tank testing and fluid-structure-interaction (FSI) simulations to investigate the behavior of foils with large deflections and the effect of these deflections on crossflow turbine performance with the goal of determining the maximum allowable deflections consonant with efficiency and a robust, durable structure. A validated modeling and simulation approach was used in the design of an Ocean Renewable Power Company (ORPC) full-scale turbine. The methodology is applicable to the design of other marine hydrokinetic (MHK) devices, and to hydrokinetic turbines which experience significant deflections during operation.

16 TIDAL AND WAVE POWER↗

Visualizing the SEI formation between lithium metal and solid-state electrolyte

The solid electrolyte interphase (SEI) is regarded as the most important factor affecting the durability of lithium-metal anode in all-solid-state batteries (ASSBs). Despite its significance, the nucleation and growth mechanism of SEI is not yet well understood. Here, we elucidate the thermodynamics and kinetics governing SEI formation at the Li|β-Li 3 PS 4 interface at the atomic scale via thermodynamic phase equilibrium analysis and machine-learning-potential-assisted molecular dynamics (MD) simulations. An accurate moment tensor potential using the machine-learning method is developed for a reactive model of Li|β-Li 3 PS 4 . This potential enabled us to perform large-scale MD simulations with the model size expanded to the experimental dimensions (~40 nm) while maintaining the same level of accuracy as density functional theory calculations. The results reveal a four-stage evolution process at the Li|β-Li 3 PS 4 interface, namely (i) fast ion diffusion, (ii) nucleation, (iii) Li 2 S growth, and (iv) stabilization. Notably, we demonstrate that the SEI can be categorized into crystalline and amorphous regions. The simulated SEI thickness, structure, and composition closely match experimental findings, validating the accuracy of the MD simulations. We further disclose the significant impact of ion diffusion kinetic limitations on the phase formation and crystallization of interfacial products. Furthermore, we shed light on the detailed potential energy (PE) distribution of lithium along the direction perpendicular to the interface. Finally, this information is crucial for better understanding interfacial ion mobility.

25 ENERGY STORAGE↗

BEYONDPLANCK IV. On end-to-end simulations in CMB analysis — Bayesian versus frequentist statistics

End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely, how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data or as a random realization independent from the data. We refer to these as posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before BEYONDPLANCK, most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BEYONDPLANCK represents the first end-to-end CMB simulation framework that is able to generate both types of simulations and these new capabilities have brought this topic to the forefront. The BEYONDPLANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level BEYONDPLANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.

79 ASTRONOMY AND ASTROPHYSICS↗

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗

Minimizing CGYRO HPC Communication Costs in Ensembles with XGYRO by Sharing the Collisional Constant Tensor Structure

First-principles fusion plasma simulations are both compute and memory intensive, and CGYRO is no exception. The use of many HPC nodes to fit the problem in the available memory thus results in significant communication overhead, which is hard to avoid for any single simulation. That said, most fusion studies are composed of ensembles of simulations, so we developed a new tool, named XGYRO, that executes a whole ensemble of CGYRO simulations as a single HPC job. By treating the ensemble as a unit, XGYRO can alter the global buffer distribution logic and apply optimizations that are not feasible on any single simulation, but only on the ensemble as a whole. The main saving comes from the sharing of the collisional constant tensor structure, since its values are typically identical between parameter-sweep simulations. This data structure dominates the memory consumption of CGYRO simulations, so distributing it among the whole ensemble results in drastic memory savings for each simulation, which in turn results in overall lower communication overhead.

CGYRO↗

Structural phase transitions in a MoWSe 2 monolayer: Molecular dynamics simulations and variational autoencoder analysis

Electrical and optoelectronic properties of two-dimensional (2D) transition metal dichalcogenides (TMDCs) can be tuned by exploiting their structural phase transitions. Here semiconducting (2H) to metallic (1T) phase transition is investigated in a strained MoWSe 2 monolayer using molecular dynamics (MD) simulations. Novel intermediate structures called α and β are found between the 2H and 1T phases. These intermediate structures are similar to those observed in a 2D MoS 2 by scanning transmission electron microscopy. Here, a deep generative model, namely the variational autoencoder (VAE) trained by MD data, is used to generate novel heterostructures with α and β interfaces. Quantum simulations based on density functional theory show that these heterostructures are stable and suitable for novel nanoelectronics applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Crystal structure of a thermophilic O 6-alkylguanine-DNA alkyltransferase-derived self-labeling protein-tag in covalent complex with a fluorescent probe

Highlights: • The labeling of proteins is a valid approach for studying their function. • SNAP-tag{sup ®} is a self-labeling protein tag that binds different probes (e.g. SVG). • SsOGT-H{sup 5} is a thermostable alternative to the SNAP-tag{sup ®} technology. • The SsOGT-H{sup 5}-SVG structure reveals new details of protein-probe interaction. • The SsOGT-H{sup 5}-SVG structure highlights new OGTs structure-function relationships. The self-labeling protein tags are robust and versatile tools for studying different molecular aspects of cell biology. In order to be suitable for a wide spectrum of experimental conditions, it is mandatory that these systems are stable after the fluorescent labeling reaction and do not alter the properties of the fusion partner. SsOGT-H{sup 5} is an engineered variant alkylguanine-DNA-alkyl-transferase (OGT) of the hyperthermophilic archaeon Sulfolobus solfataricus, and it represents an alternative solution to the SNAP-tag{sup ®} technology under harsh reaction conditions. Here we present the crystal structure of SsOGT-H{sup 5} in complex with the fluorescent probe SNAP-Vista Green{sup ®} (SsOGT-H{sup 5}-SVG) that reveals the conformation adopted by the protein upon the trans-alkylation reaction with the substrate, which is observed covalently bound to the catalytic cysteine residue. Moreover, we identify the amino acids that contribute to both the overall protein stability in the post-reaction state and the coordination of the fluorescent moiety stretching-out from the protein active site. We gained new insights in the conformational changes possibly occurring to the OGT proteins upon reaction with modified guanine base bearing bulky adducts; indeed, our structural analysis reveals an unprecedented conformation of the active site loop that is likely to trigger protein destabilization and consequent degradation. Interestingly, the SVG moiety plays a key role in restoring the interaction between the N- and C-terminal domains of the protein that is lost following the new conformation adopted by the active site loop in the SsOGT-H{sup 5}-SVG structure. Molecular dynamics simulations provide further information into the dynamics of SsOGT-H{sup 5}-SVG structure, highlighting the role of the fluorescent ligand in keeping the protein stable after the trans-alkylation reaction.

60 APPLIED LIFE SCIENCES↗

Chapter 2: Blade-Resolved Modeling with Fluid-Structure Interaction

Blade -resolved simulations are a high-fidelity modeling tool that helps understand the fundamental flow physics that create loads on wind turbines. Along with fluid - structure interaction (FSI), they provide a foundation to improve the state-of-the-art aerodynamic models of modern wind turbines with large flexible blades and thus the design of the next generation of wind turbines. One of the most important parts of wind-turbine simulation is modeling the conversion of kinetic energy in the wind into mechanical energy using principles of aerodynamics. The hierarchies of aerodynamic models for wind turbines in increasing the order of complexity, cost, and fidelity are blade element momentum (BEM) theory (BEMT), vortex methods, actuator disc/line methods, and blade -resolved modeling.

17 WIND ENERGY↗

Simulating Dopant Diffusion in a Detalied Porous Structure

The National Ignition Facility (NIF) at Lawrence Livermore National Laboratory (LLNL) uses highpowered lasers to explore cutting edge fusion technology. Specialized optics, known as fiber lasers, generate and amplify the NIF’s high-powered lasers that facilitate fusion experiments. Fiber lasers are glass fibers that are doped with heavy elements to stimulate light emission. The fiber laser doping processing uses a predeposited layer of porous silica nanoparticles on the surface of the glass to absorb a dopant rich salt solution. Our project simulated the diffusion of dopants into porous silica with prescribed colloidal stacking geometries using COMSOL Multiphysics. With these simulations we identified the porosity and tortuosity of the porous nanoparticle structures as the primary parameters that will dictate how the dopants will diffusive into the porous surface of the glass. Porosity and Tortuosity are geometric parameters that can easily be input into more detailed simulations to estimate an effective diffusivity within the porous layer. Future work should account for surface reactions to accurately simulate the full-scale diffusion of the solution doping process.

36 MATERIALS SCIENCE↗

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Abstract Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

36 MATERIALS SCIENCE↗

From Electronic Structure to Ion Transport: Photoelectron Spectroscopy and Molecular Dynamics Simulations Reveal the Role of Anions in Lithium Battery Electrolytes

Electrolyte anions are pivotal for lithium battery performance, yet their fundamental electronic structural properties are not well understood. In this work, we employ a combination of negative-ion photoelectron spectroscopy (NIPES), ab initio calculations, and molecular dynamics (MD) simulations to investigate the electronic structures of three representative electrolyte anions. This multiscale approach enables us to elucidate how their intrinsic electronic properties govern anion–solvent interactions in gas-phase clusters, as well as lithium-ion (Li + ) solvation structures and ion transport behavior in the condensed phase. NIPES reveals that difluoro(oxalato)borate (DFOB – ), bis(fluorosulfonyl)imide (FSI – ), and bis(oxalato)borate (BOB – ) all exhibit high electron binding energies, with vertical/adiabatic detachment energies increasing from DFOB – (6.09/5.70 eV) to FSI – (6.80/6.10 eV) to BOB – (6.82/6.40 eV), correlating with enhanced oxidation stability. Ab initio calculations reveal that DFOB – /FSI – –solvent complexes bind Li + ∼ 10 kcal/mol stronger than BOB – series, aligning with the strength of a Li + –anion model. DFOB – exhibits pronounced charge localization on both oxygen and fluorine atoms, enabling their involvement in Li + coordination. In contrast, fluorine atoms in FSI – are largely electron-depleted and remain excluded from direct Li + binding. MD simulations further demonstrate that LiDFOB and LiFSI systems exhibit Li + diffusion coefficients three and five times higher than those of LiBOB across four common solvents. Notably, LiFSI salt in acetonitrile (AN) exhibits the fastest Li + diffusion among 12 electrolyte systems, highlighting the synergistic effect of FSI – and AN in promoting ion mobility. In conclusion, these findings provide a molecular-level understanding of the critical roles of anion and its microsolvation in optimizing Li + diffusion dynamics, once again emphasizing the positioning of FSI – and DFOB – as prime candidates for next-generation electrolytes.

25 ENERGY STORAGE↗

Directed self-assembly of chiral liquid crystals into biomimetic bouligand structures in thin film

The Bouligand structure, renowned for its helicoidal arrangement and enhanced mechanical properties, has attracted significant research interest for its ability to impart enhanced strength to intrinsically soft materials. Biomimetic approaches have centered on fibrous structures in bulk materials, but translating this architecture into thin-film regime for miniaturized-wearable devices with programmable functions remains challenging. Here, we direct the self-assembly of cholesteric liquid crystals (CLCs) into hierarchical helical structures using chemically patterned surfaces. Alternating surface anchoring regions align uniform lying-down helices at the nanoscale, guiding a secondary microscale helical structure exhibiting both left- and right-handed twists. This mimetic Bouligand structure in CLCs enables optical modulation under applied field and strain with enhanced mechanical response. Simulations reveal the structural evolution from initial Bouligand configuration in LC layers to alternating twist helices. This research provides a basis for designing and manufacturing miniaturized or wearable devices with nanometer-scale precision in regulating electro-optical and mechanical properties. Bouligand structures, which offer strength in natural materials, are of interest but difficult to obtain. Here, the authors report the development of such structures by directed self-assembly of cholesteric single crystals into hierarchical helical structures, with the secondary structure having right and left-handed twists.

Bouligand structure↗

BOUT++ nonlinear simulation for a comparative study with the measured 2D ELM structures in the KSTAR H-mode plasma

Here, time-dependent nonlinear simulation using the BOUT++ 3-field model is performed for a study of edge localized mode (ELM) dynamics in KSTAR H-mode plasmas. The simulated results of the ELM evolution are directly compared with the 2D ELM structures, measured by electron cyclotron emission imaging systems, from the initial growth to crash phase. The simulation results are qualitatively in good agreement with the observations. Specifically, it is notable in simulation and experiment that higher-n modes in a quasi-steady state phase evolve into lower-n modes close to the ELM-crash in time. It indicates that lower-n modes may be responsible for the onset of the crash. A postanalysis confirms that nonlinear mode coupling between toroidal harmonics may drive the lower-n mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interfacial Atomistic Mechanisms of Lithium Metal Stripping and Plating in Solid–State Batteries

All-solid-state batteries based on a Li metal anode represent a promising next-generation energy storage system, but are currently limited by low current density and short cycle life. Further research to improve the Li metal anode is impeded by the lack of understanding in its failure mechanisms at lithium–solid interfaces, in particular, the fundamental atomistic processes responsible for interface failure. Here, using large-scale molecular dynamics simulations, the first atomistic modeling study of lithium stripping and plating on a solid electrolyte is performed by explicitly considering key fundamental atomistic processes and interface atomistic structures. In the simulations, the interface failure initiated with the formation of nano-sized pores, and how interface structures, lithium diffusion, adhesion energy, and applied pressure affect interface failure during Li cycling are observed. By systematically varying the parameters of solid-state lithium cells in the simulations, the parameter space of applied pressures and interfacial adhesion energies that inhibit interface failure during cycling are mapped to guide selection of solid-state cells. Furthermore, this study establishes the atomistic modeling for Li stripping and plating, and predicts optimal solid interfaces and new strategies for the future research and development of solid-state Li-metal batteries.

25 ENERGY STORAGE↗

Numerical Simulation of Infrasound Resonance in Underground Tunnel Structures

Remote observation of infrasound resonant signals emanating from underground tunnel structures could potentially allow the remote quantification of the geometry of the underground structures in which the signals were generated. However, the sensitivity of these observations to tunnel geometry and changes in that geometry are unknown. In this report we outline a numerical simulation study with the following three objectives: 1. Can we model infrasound resonance using numerical simulations? 2. What is the sensitivity of the modeled observations to changes in tunnel geometry or boundary conditions? And 3. Can we accurately model resonant observations of explosions occurring in Redmond Salt Mine in Central Utah. In this report we outline affirmative answers to the first two objectives, but we were unable to accurately model the Redmond Explosions due to numerical instability in the model of the complex structure of the mine. Recommendations for future work emphasize the need to acquire additional datasets and to explore a more data-based approach in which changes in data signatures are detected as a first step towards developing a method which inverts resonant infrasound signals for tunnel geometry.

47 OTHER INSTRUMENTATION↗

Decomposition of L a 2 – x S r x Cu O 4 into several L a 2 O 3 phases at elevated temperatures in ultrahigh vacuum inside a transmission electron microscope

Here, we report the decomposition of La 2–x Sr x CuO 4 into La 2 O 3 and Cu nanoparticles in ultrahigh vacuum, observed by in situ heating experiments in a transmission electron microscope. The analysis of electron diffraction data reveals that the phase decomposition process starts at about 150°C and is considerably expedited in the temperature range of 350°C–450°C. Two major resultant solid phases are identified as metallic Cu and La 2 O 3 by electron diffraction, simulation, and electron energy-loss spectroscopy (EELS) analyses. With the aid of calculations, La 2 O 3 phases are further identified to be derivatives of a fluorite structure—fluorite, pyrochlore, and (distorted) bixbyite—characterized by different oxygen-vacancy order. Additionally, the bulk plasmon energy and the fine structures of the O K and LaM 4,5 EELS edges are reported for these structures, along with simulated O K x-ray absorption near-edge structure. The resultant Cu nanoparticles and La 2 O 3 phases remain unchanged after cooling to room temperature.

36 MATERIALS SCIENCE↗

Simulating Sea-Ice Deformation in Viscous-Plastic Sea-Ice Models With CD-Grids

Linear kinematic features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. In viscous-plastic (VP) sea-ice models, the simulation of LKFs depends on several factors such as the grid resolution, the numerical solver convergence, and the placement of the variables on the mesh. In this study, we compare two recently proposed discretization with a CD-grid placement with respect to their ability to reproduce LKFs. The first (CD1) is based on a nonconforming finite element discretization, whereas the second (CD2) uses a conforming subgrid discretization. To analyze their resolution properties, we evaluate runs from different models (e.g., FESOM, MPAS) on a benchmark problem using quadrilateral, hexagonal and triangular meshes. Our findings show that the CD1 setup simulates more deformation structure than the CD2 setup. This highlights the importance of the type of spatial discretization for the simulation of LKFs. Due to the higher number of degrees of freedom, both CD-grids resolve more LKFs than traditional A, B, and C-grids at fixed mesh level. This is an advantage of the CD-grid approach, as high spatial mesh resolution is needed in VP sea-ice models to simulate LKFs.

54 ENVIRONMENTAL SCIENCES↗

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗