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 199 records · Page 11

Synthetic Access to a Framework-Stabilized and Fully Sulfided Analogue of an Anderson Polyoxometalate that is Catalytically Competent for Reduction Reactions

Polyoxometalates (POMs) featuring 7, 12, 18, or more redox-accessible transition metal ions are ubiquitous as selective catalysts, especially for oxidation reactions. The corresponding synthetic and catalytic chemistry of stable, discrete, capping-ligand-free polythiometalates (PTMs), which could be especially attractive for reduction reactions, is much less well developed. Among the challenges are the propensity of PTMs to agglomerate and the tendency for agglomeration to block reactant access of catalyst active sites. Nevertheless, the pervasive presence of transition metal sulfur clusters metalloenzymes or cofactors that catalyze reduction reactions and the justifiable proliferation of studies of two-dimensional (2D) metal-chalcogenides as reduction catalysts point to the promise of well-defined and controllable PTMs as reduction catalysts. Here, we report the fabrication of agglomeration-immune, reactant-accessible, capping-ligand-free Co II Mo 6 IV S 24 n– clusters as periodic arrays in a water-stable, hierarchically porous Zr-metal–organic framework (MOF; NU1K) by first installing a disk-like Anderson polyoxometalate, Co III Mo 6 VI O 24 m– , in size-matched micropores where the siting is established via difference electron density (DED) X-ray diffraction (XRD) experiments. Flowing H 2 S, while heating, reduces molybdenum(VI) ions to Mo(IV) and quantitatively replaces oxygen anions with sulfur anions (S 2– , HS – , S 2 2– ). DED maps show that MOF-templated POM-to-PTM conversion leaves clusters individually isolated in open-channel-connected micropores. Importantly, the structure of the immobilized cluster as determined, in part, by X-ray photoelectron spectroscopy (XPS), X-ray absorption fine structure (XAFS) analysis, and pair distribution function (PDF) analysis of total X-ray scattering agrees well with the theoretically simulated structure. PTM@MOF displays both electrocatalytic and photocatalytic competency for hydrogen evolution. Nevertheless, the initially installed PTM appears to be a precatalyst, gaining competency only after the loss of ~3 to 6 sulfurs and exposure to hydride-forming metal ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Establishing model credibility for process-microstructure-property relationships in additive manufacturing using exascale computing

Additive Manufacturing (AM) of alloys holds significant promise as a disruptive technology in various industries, yet its adoption is often hindered by challenges in achieving consistent part quality. These issues are primarily due to the complex process-microstructure-property (PSP) relationships inherent to AM. Computational models can greatly aid in understanding these relationships, but their widespread impact and adoption has been limited by a lack of validated, open-source, and computationally efficient PSP modeling frameworks and hardware limitations. Here, this study leverages the ExaAM software suite and data from the AMBench-2018 series of laser powder bed fusion (LPBF) benchmark experiments to perform a comprehensive model assessment, including verification, validation, sensitivity analysis, and uncertainty quantification. The RADICAL-EnTK workflow manager was used to perform an ensemble of heat transport, solidification, and mechanical response simulations on the exascale computer Frontier, considering uncertainties in critical model inputs such as laser spot size and nucleation parameters, and consisting of 125 explicit grain structure simulations and 7875 crystal plasticity simulations. For a selected location within the Inconel 625 AMBench-2018 test artifact, sensitivity analysis and uncertainty quantification were performed using the predicted distributions of grain structure and mechanical properties. Qualitative agreement was found between the predicted grain size and texture and the observed AMBench-2018 microstructure, the mean predicted yield stress was within 5% of the experimental measurement mean, and the mean predicted engineering stress at 5% strain was within 10% of the experimental measurement mean. The insights gained from development and validation of the ExaAM PSP modeling framework will help guide future directions for enhancing the credibility and reliability of PSP models in AM, thereby accelerating the adoption of AM technologies in various industries.

Additive manufacturing↗

MPEX AI Digital Twins

All magnetically confined plasma fusion power plant concepts (Tokamak, Spherical Tokamak, Stellarator, Mirror, ...) must exhaust the heat and plasma from the core confinement region to the material walls. The primary channel for this exhaust is through a plasma divertor which directs plasma along open magnetic field lines to a material target. The Material Plasma Exposure eXperiment (MPEX) illustrated in Figure 1, is a high-power, steady-state linear plasma device designed to produce the plasma material interaction (PMI) conditions of the divertor of future magnetic confinement fusion power plants: energy flux 20MW/m 2 , ion fluence 1031/m 2 , pulse duration 106 sec. These goals of plasma exposure in MPEX are well beyond those achieved in magnetic fusion experimental devices. Successfully achieving these high power steady state conditions for long pulses requires operational control of the heating and particle sources and the plasma flux to the walls and target. The MPEX AI Hot Spot Controller, proposed in this project, will help achieve the operational milestones of MPEX. The MPEX device will begin commissioning at the end of FY26. A smaller proto-MPEX was operated for 14,666 plasma discharges and will resume operation in September of 2025 as proto-MPEX-lite, with reduced capability, to test a new window for the Helicon plasma source. The proto-MPEX data has undergone surrogate modeling with machine learning methods (R. Archibald, 2022 IEEE International Conference on Big Data). This proto-MPEX data will be used to begin development of the AI digital twins described in this white paper. The scientific mission of MPEX is to qualify materials of different composition for use in the high energy and plasma flux conditions of a fusion power plant. The materials exposed in MPEX will in some cases be exposed to high neutron fluxes at other ORNL facilities to measure the changes to their PMI properties. The targets exposed in MPEX will be transported under vacuum to a Surface Analysis Station (SAS). The SAS will be equipped with the following diagnostics: Focused Ion Beam (FIB) for trench milling, 100-400 angstrom resolution scanning electron microscope (SEM), surface mapping x-ray spectrometer, high resolution camera, and a future upgrade to a laser induced breakdown spectroscopy quadruple mass spectrometer (LIBS-QMS). The MPEX experiments will generate diverse pre- and post-exposure measurement data of detailed material properties down to the crystal grain level in 3D for post-exposure assessment of PMI damage (e.g. cracking, melting, erosion and redeposition of the material). Physics models for the PMI, and how the material composition and manufacturing impact its performance under high energy plasma exposure, need to be validated with MPEX data to guide the selection of new candidate materials. Our vision for the MPEX AI Digital Twins project is to supply experimental and physics model simulation data to train Artificial Intelligence (AI) models for data processing, analysis, operational control, PMI and materials simulation to maximize the scientific output of the MPEX device. Ultimately, an AI digital twin of MPEX material assessment metrics for tested and synthetic material types with simulated PMI will be trained by the AI Modeling Teams on the experimental and physics simulation data submitted to the American Science Cloud by this project. A purely empirical search for the best material is inefficient given the finite number of samples that can be tested on MPEX. In order to expand the material properties database for training the MPEX Material Assessment AI Digital Twin, and to gain physics understanding of the PMI processes, physics models of the material properties and PMI processes are required. The physics simulations provide detailed simulation data, like impact angles for plasma ions, sputtering yields, transport of the ionized sputtered target material in the plasma, and redeposition locations. This simulation data expands the measurement data for deeper physics understanding. The experimental data is essential to validate the PMI and material structure simulation models. The validated models can then be used to generate new simulation data of MPEX material assessments for synthetic material compositions that have not been exposed in MPEX. These predictive simulations, plus the whole experimental dataset, will be used to train the MPEX Material Assessment AI Digital Twin allowing a rapid generative AI search for new materials with reduced PMI damage by interpolating the domain of the training set. These new optimum materials can be simulated with the physics codes and/or tested in MPEX. The ability of AI neural networks to interpolate multi-dimensional parameter spaces and generate virtual data is exploited for a more efficient search for optimum materials. The advent of the Transformational AI Models Consortium (TAIMC) is an opportunity to engage with state of the art private and public AI developers to achieve the goals of the AI digital twins and AI accelerated physics models proposed in this project. Our partners at ORNL from the Advance Scientific Computing Research (ASCR) organization will collaborate in accelerating the integrated plasma material interaction simulation framework. This simulation framework will provide a platform for generating simulation data across a range of physical fidelities, including hybrid methods that produce multi-fidelity results. This data will be leveraged for AI model development, both for generation of surrogates and the automation of simulation campaigns. A part of the research below will include collaborative efforts with the TAIMC to (i) adapt data storage approaches to ensure AI-readiness, (ii) provide a protypical exemplar to inform and exercise constructed workflows, and (iii) generate and share data, using the TAIMC unified AI data standard, for foundational models that will be trained from multiple sources across the DOE complex. We will also collaborate with the TAIMC, as well as the planned AI modeling teams, to develop approaches for reducing the cost of data generation. These include tailored multi-fidelity approaches as well as fine-tuning strategies to augment general, large-scale foundational models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗

Numerical simulation of actively controlled space structures

The numerical simulation of large, actively controlled space structures is a problem that has received considerable attention. A new computational method is presented in this paper, which simulates the vibration suppression of a large space structure by a linear compensator. The structure is assumed to be modeled by a large linear finite element model. Nested circularly linked lists and the trapezoidal rule are key aspects of this algorithm. There is an emphasis on generality, while establishing an equivalence between numerical stability and system stability.

Quan, Ralph↗

OC7 phase I: CFD investigation of viscous forces on rectangular members of semisubmersibles

In Phase I of the Offshore Code Comparison Collaboration 7 (OC7) project, we aim to address key questions related to the modeling of viscous hydrodynamic loads on floating platforms. Here, this paper presents results from a collaborative effort within Phase I of the OC7 project. We investigate viscous forces acting on the rectangular members of semisubmersible platforms using high-fidelity computational fluid dynamics (CFD) simulations. The parametric CFD analyses include two-and-a-half-dimensional (2.5D) simulations of pontoon sections of two different semisubmersible platform designs (the VolturnUS-S and the INO WINDMOOR) and full three-dimensional (3D) structure simulations of the VolturnUS-S platform. We extract both the added mass and drag coefficients that are typically used in engineering-level hydrodynamic models from the 2.5D forced oscillation simulations with various Keulegan-Carpenter numbers and compare them to the coefficients from the 3D forced oscillation simulations. We observe consistent variation of the coefficients along the pontoon legs due to 3D platform and free-surface effects. With the breadth of the simulations and component-level detailed analyses, this work establishes a useful database and several empirical formulas for selecting hydrodynamic coefficients when modeling semisubmersible platforms in engineering tools.

13 HYDRO ENERGY↗

A-SST Initial Specification

The U.S. Army Research Office (ARO), in partnership with IARPA, are investigating innovative, efficient, and scalable computer architectures that are capable of executing next-generation large scale data-analytic applications. These applications are increasingly sparse, unstructured, non-local, and heterogeneous. Under the Advanced Graphic Intelligence Logical computing Environment (AGILE) program, Performer teams will be asked to design computer architectures to meet the future needs of the DoD and the Intelligence Community (IC). This design effort will require flexible, scalable, and detailed simulation to assess the performance, efficiency, and validity of their designs. To support AGILE, Sandia National Labs will be providing the AGILE-enhanced Structural Simulation Toolkit (A-SST). This toolkit is a computer architecture simulation framework designed to support fast, parallel, and multi-scale simulation of novel architectures. This document describes the A-SST framework, some of its library of simulation models, and how it may be used by AGILE Performers.

97 MATHEMATICS AND COMPUTING↗

3D radiated power analysis of JET SPI discharges using the Emis3D forward modeling tool

Abstract Precise values for radiated energy in tokamak disruption experiments are needed to validate disruption mitigation techniques for burning plasma tokamaks like ITER and SPARC. Control room analysis of radiated power ( P rad ) on JET assumes axisymmetry, since fitting 3D radiation structures with limited bolometry coverage is an under-determined problem. In mitigated disruptions, radiation is toroidally asymmetric and 3D, due to fast-growing 3D MHD modes and localized impurity sources. To address this problem, Emis3D adopts a physics motivated forward modeling (‘guess and check’) approach, comparing experimental bolometry data to synthetic data from user-defined radiation structures. Synthetic structures are observed with the Cherab modeling framework and a best fit chosen using a reduced χ 2 statistic. 2D tomographic inversion models are tested, as well as helical flux tubes and 3D MHD simulated structures from JOREK. Two nominally identical pure neon shattered pellet injection (SPI) mitigated discharges in JET are analyzed. 2D tomographic inversions with added toroidal freedom are the best fits in the thermal quench (TQ) and current quench (CQ). In the pre-TQ, 2D reconstructions are statistically the best fits, but are likely over-optimized and do not capture the 3D radiation structure seen in fast camera images. The next-best pre-TQ fits are helical structures that extend towards the high-field side, consistent with an impurity flow under the magnetic nozzle effect also observed in JOREK simulations. Whole-disruption radiated fractions of 0.98 + 0.03 / − 0.29 and 1.01 + 0.02 / − 0.17 are found, suggesting that the stored energy may have been fully mitigated by each SPI, although mitigation efficiencies well below ITER and SPARC requirements for high energy pulses are still within the large uncertainties. Emis3D is also used to validate JOREK SPI simulations, and confirms improvements in matching experiment from changes to impurity modeling. Time-dependent toroidal peaking factors are calculated and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Application of Composite Mechanics to Composites Enhanced Concrete Structures

A new and effective method is described to design composites to repair damage or enhance the overload strength of concrete infrastructures. The method is based on composite mechanics which is available in computer codes. It is used to simulate structural sections made from reinforced concrete which are typical in infrastructure as well as select reinforced concrete structures. The structural sections are represented by a number of layers through the thickness where different layers are used in concrete, and for the composite. The reinforced concrete structures are represented with finite elements where the element stiffness parameters are from the structural sections which are represented by composite mechanics. The load carrying capability of the structure is determined by progressive structural fracture. Results show up to 40 percent improvements for damage and for overload enhancement with relatively small laminate thickness for the structural sections and up to three times for the composite enhanced select structures (arches and domes).

Chamis, Christos C.↗

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗

Difference in load predictions obtained with effective turbulence vs. a dynamic wake meandering modeling approach

Abstract. According to the international standard for wind turbine design, the effects of wind turbine wakes on structural loads can be considered in two ways: (1) by augmenting the ambient turbulence levels with the effective turbulence model (EFF) and then calculating the resulting loads and (2) by performing dynamic wake meandering (DWM) simulations, which compute wake effects and loads for all turbines on a farm at once. There is no definitive answer in scientific literature as to the consequences of choosing one model over the other, but the two approaches are unarguably very different. The work presented here expounds on these differences and investigates to what extent they affect the simulated structural loads. We consider an idealized 4×4 rectangular array of National Renewable Energy Laboratory 5 MW wind turbines with a spacing of 5 by 8 rotor diameters and three wind speed scenarios at high ambient turbulence. Load simulations are performed in OpenFAST with EFF and in FAST.Farm with the DWM implementation. We compare ambient turbulence, wind farm turbulence, and loads between both approaches. When omnidirectional results are compared, EFF wind farm turbulence intensity is consistently higher by 0.2 % (above-rated wind speed) to 2.7 % (below-rated wind speed). However, for certain wind directions, the EFF turbulence can be lower than FAST.Farm by almost 9 %. Wind speeds within the farm were found to differ by up to 3 m s−1 due to the lack of wake deficits in the EFF approach, leading to longer tails toward low values in the FAST.Farm mean load distributions. Consistent with the turbulence results, the median EFF load standard deviations are also consistently higher, by a maximum of 20 % and 17 % for blade-root out-of-plane and tower-base fore-aft moments, respectively.

17 WIND ENERGY↗

Dynamic analysis of an unrestrained, rotating structure through nonlinear simulation

An alternate derivation based on Hamilton's principle of equations of motion of a flexible body for use in nonlinear simulation modeling is presented. The equations include terms representing variation in the inertial properties of the body resulting from flexibility as well as terms coupling rigid body and flexible momenta that arise from satisfaction of only the linear part of the mean axis conditions. An approximate method for including the centrifugal stiffening in a way that permits a variable spin vector is given. The equations are exercised in a simulation program on an example structure initially spinning about its unstable principal axis in a gravity-free vacuum.

Zeiler, Thomas A.↗

Airblast High Explosive Sensitivity Study

The ability to accurately simulate the flow produced inside the LANL blast tube is an important capability that is currently under development in W-13. This capability will be used to produce boundary conditions for structural simulations of blast tube tests. This capability could also be used to improve future blast tube test campaigns and inform future blast tube upgrades.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Numerical Simulations of Two-Phase Reacting Flow in a Single-Element Lean Direct Injection (LDI) Combustor Using NCC

A series of numerical simulations of Jet-A spray reacting flow in a single-element lean direct injection (LDI) combustor have been conducted by using the National Combustion Code (NCC). The simulations have been carried out using the time filtered Navier-Stokes (TFNS) approach ranging from the steady Reynolds-averaged Navier-Stokes (RANS), unsteady RANS (URANS), to the dynamic flow structure simulation (DFS). The sub-grid model employed for turbulent mixing and combustion includes the well-mixed model, the linear eddy mixing (LEM) model, and the filtered mass density function (FDF/PDF) model. The starting condition of the injected liquid spray is specified via empirical droplet size correlation, and a five-species single-step global reduced mechanism is employed for fuel chemistry. All the calculations use the same grid whose resolution is of the RANS type. Comparisons of results from various models are presented.

Liu, Nan-Suey↗

SPARC: Structural properties associated with residue constraints

SPARC facilitates the generation of plausible hypotheses regarding underlying biochemical mechanisms by structurally characterizing protein sequence constraints. Such constraints appear as residues co-conserved in functionally related subgroups, as subtle pairwise correlations (i.e., direct couplings), and as correlations among these sequence features or with structural features. SPARC performs three types of analyses. First, based on pairwise sequence correlations, it estimates the biological relevance of alternative conformations and of homomeric contacts, as illustrated here for death domains. Second, it estimates the statistical significance of the correspondence between directly coupled residue pairs and interactions at heterodimeric interfaces. Third, given molecular dynamics simulated structures, it characterizes interactions among constrained residues or between such residues and ligands that: (a) are stably maintained during the simulation; (b) undergo correlated formation and/or disruption of interactions with other constrained residues; or (c) switch between alternative interactions. We illustrate this for two homohexameric complexes: the bacterial enhancer binding protein (bEBP) NtrC1, which activates transcription by remodeling RNA polymerase (RNAP) containing σ 54 , and for DnaB helicase, which opens DNA at the bacterial replication fork. Based on the NtrC1 analysis, we hypothesize possible mechanisms for inhibiting ATP hydrolysis until ADP is released from an adjacent subunit and for coupling ATP hydrolysis to restructuring of σ 54 binding loops. Based on the DnaB analysis, we hypothesize that DnaB ‘grabs’ ssDNA by flipping every fourth base and inserting it into cavities between subunits and that flipping of a DnaB-specific glutamine residue triggers ATP hydrolysis.

97 MATHEMATICS AND COMPUTING↗

Dynamic test/analysis correlation using reduced analytical models

Test/analysis correlation is an important aspect of the verification of analysis models which are used to predict on-orbit response characteristics of large space structures. This paper presents results of a study using reduced analysis models for performing dynamic test/analysis correlation. The reduced test-analysis model (TAM) has the same number and orientation of DOF as the test measurements. Two reduction methods, static (Guyan) reduction and the Improved Reduced System reduction, are applied to the test/analysis correlation of a laboratory truss structure. Simulated test results and modal test data are used to examine the performance of each method. It is shown that selection of DOF to be retained in the TAM is critical when large structural masses are involved. In addition, the use of modal test results may provide difficulties in TAM accuracy even if a large number of DOF are retained in the TAM.

Mcgowan, Paul E.↗