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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermo-Mechanical Modeling and Evaluation of the Cracking Response of Additively Manufactured Monolithic SiC Lattice Structures Subjected to Laser Heating

Increasing operating temperatures of solar receivers is paramount to the efficiency of concentrated solar thermal (CST) and solar power (CSP) systems. Owing to its high temperature stability combined with excellent thermal and optical properties, SiC has been the material of choice for application in high-temperature solar receivers. The state-of-the-art SiC volumetric concentrating solar air receivers such as honeycomb design have been demonstrated in field tests to achieve exit air temperatures approaching 800oC. However, successful application of CST systems for decarbonization requires significant increase in temperature capability of SiC receiver technology. Supported by an award from the Solar Technology Office (SETO), US Department of Energy (DOE), GE Research in collaboration with Heliogen Holdings Inc, is engaged in the development of ultra-High Operating Temperature SiC-matrix Solar Thermal Air Receiver (HOTSSTAR) enabled by additive manufacturing. The program objective is to design and to demonstrate a techno-economically viable SiC air receiver technology to achieve exit air temperature of 1100oC for CST applications. This report summarizes the learnings of a computational study on additively manufactured SiC lattice structures subjected to laser heating. The results from this study aim to help guide the design of a SiC receiver through a better understanding of lattice structure geometric parameters and their implications on the cracking response of the structure. In the study, a heat flux was applied to the surface of a 2”-diameter cylindrical lattice structure to simulate a 4kW CO2 laser. The resulting temperature distribution was applied to a structural model to approximate the stress distribution within the lattice structure and a Weibull analysis was performed to gain insight into the probability of failure and to evaluate the cracking response of the structure. With this approach, the influence of lattice density on temperature, stress, and probability of failure was explored for two different lattice beam spacings. A discussion on modeling assumptions, a comparison with experimental results, and an evaluation of the lattice cracking response is provided.

14 SOLAR ENERGY↗
DOE OSTI2022

Improved Axisymmetric and High Temperature Material Structural Modeling in MOOSE and NEML

This report describes improvements made to the solid mechanics formulation in the MOOSE open source finite element simulation environment and the open source Nuclear Engineering Material model Library (NEML) for mechanical constitutive models. The focus of these improvements is to improve the usability and performance of simulations involving one or both pieces of software. Specifically, this work completes a new system for solid mechanics simulations in the MOOSE ecosystem providing exact linearizations and optimal (quadratic) convergence, for a variety of coordinate systems and material types, including large deformation simulations. This work then provides users a framework to build highly efficient mechanical simulations of structures or materials or to couple in additional MOOSE physics modules to build complex, scalable multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Formation of a constructed microbial community in a nutrient-rich environment indicates bacterial interspecific competition

ABSTRACT Understanding the organizational principles of microbial communities is essential for interpreting ecosystem stability. Previous studies have investigated the formation of bacterial communities under nutrient-poor conditions or obligate relationships to observe cooperative interactions among different species. How microorganisms form stabilized communities in nutrient-rich environments, without obligate metabolic interdependency for growth, is still not fully disclosed. In this study, three bacterial strains isolated from the Populus deltoides rhizosphere were co-cultured in complex medium, and their growth behavior was tracked. These strains co-exist in mixed culture over serial transfer for multiple growth-dilution cycles. Competition is proposed as an emergent interaction relationship among the three bacteria based on their significantly decreased growth levels. The effects of different initial inoculum ratios, up to three orders of magnitude, on community structure were investigated, and the final compositions of the mixed communities with various starting composition indicate that community structure is not dependent on the initial inoculum ratio. Furthermore, the competitive relationships within the community were not altered by different initial inoculum ratios. The community structure was simulated by generalized Lotka-Volterra and dynamic flux balance analysis to provide mechanistic predictions into emergence of community structure under a nutrient-rich environment. Metaproteomic analyses provide support for the metabolite exchanges predicted by computational modeling and for highly altered physiologies when microbes are grown in co-culture. These findings broaden our understanding of bacterial community dynamics and metabolic diversity in higher-order interactions and could be significant in the management of rhizospheric bacterial communities. IMPORTANCE Bacteria naturally co-exist in multispecies consortia, and the ability to engineer such systems can be useful in biotechnology. Despite this, few studies have been performed to understand how bacteria form a stable community and interact with each other under nutrient-rich conditions. In this study, we investigated the effects of initial inoculum ratios on bacterial community structure using a complex medium and found that the initial inoculum ratio has no significant impact on resultant community structure or on interaction patterns between community members. The microbial population profiles were simulated using computational tools in order to understand intermicrobial relationships and to identify potential metabolic exchanges that occur during stabilization of the bacterial community. Studying microbial community assembly processes is essential for understanding fundamental ecological principles in microbial ecosystems and can be critical in predicting microbial community structure and function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Incorporating bubble growth volume feedback to improve simulation of the response of a structure containing liquid and gas to sudden energy input

The SNS target module is a stainless-steel structure that contains and directs mercury. The mercury is struck with short, intense proton pulses to create neutrons. The pulses deposit energy and cause loads on the target structure and mercury cavitation. Helium bubbles are introduced into the mercury to reduce loads and cavitation. This gas complicates the prediction of the target module’s physical response. The current method used to simulate the structural response to the proton pulse incorporates a simple but effective approach to approximate mercury cavitation effects. This method cannot account for the changes from injected non-condensable gas. The ability to model the target vessel response is a prerequisite for estimating its life. Known approaches for simulating a bubbly fluid mixture’s response are computationally expensive and impractical for detailed engineering models. A constitutive model is proposed for target mercury with injected gas bubbles, with the intent of providing more accurate mercury–vessel response behavior without the computing costs required to track individual bubbles. The model’s assumptions and a computational model embodiment are described. The bubbly mercury computational model is subjected to test problems to assess its potential utility. Here, the model was able to better predict the effect of gas on wave propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

EM counterparts of structured jets from 3D GRMHD simulations

ABSTRACT GW170817/GRB170817A has offered unprecedented insight into binary neutron star post-merger systems. Its Prompt and afterglow emission imply the presence of a tightly collimated relativistic jet with a smooth transverse structure. However, it remains unclear whether and how the central engine can produce such structured jets. Here, we utilize 3D general relativistic magnetohydrodynamic simulations starting with a black hole surrounded by a magnetized torus with properties typically expected of a post-merger system. We follow the jet, as it is self-consistently launched, from the scale of the compact object out to more than three orders of magnitude in distance. We find that this naturally results in a structured jet, which is collimated by the disc wind into a half-opening angle of roughly 10°; its emission can explain features of both the prompt and afterglow emission of GRB170817A for a 30° observing angle. Our work is the first to compute the afterglow, in the context of a binary merger, from a relativistic magnetized jet self-consistently generated by an accreting black hole, with the jet’s transverse structure determined by the accretion physics and not prescribed at any point.

Kathirgamaraju, Adithan (ORCID:000000028560692X)↗
DOE OSTI2021

Structural stability and artificial buckling modes in topology optimization

Abstract This paper demonstrates how a strain energy transition approach can be used to remove artificial buckling modes that often occur in stability constrained topology optimization problems. To simulate the structural response, a nonlinear large deformation hyperelastic simulation is performed, wherein the fundamental load path is traversed using Newton’s method and the critical buckling load levels are estimated by an eigenvalue analysis. The goal of the optimization is to minimize displacement, subject to constraints on the lowest critical buckling loads and maximum volume. The topology optimization problem is regularized via the Helmholtz PDE-filter and the method of moving asymptotes is used to update the design. The stability and sensitivity analyses are outlined in detail. The effectiveness of the energy transition scheme is demonstrated in numerical examples.

Dalklint, Anna (ORCID:0000000346195205)↗
DOE OSTI2023

Understanding the dynamic evolution of atomically dispersed Cu catalyst for CO 2 electrochemical conversion using integrated XANES analysis and mechanistic studies

Direct electrochemical conversion of CO 2 to ethanol (CH 3 CH 2 OH) offers a promising strategy to lower CO 2 emission while storing energy from renewable electricity. Our recent study reported a carbon-supported atomically dispersed Cu catalyst that achieved the highest reported selectivity for CH 3 CH 2 OH formation (91%) at a relatively low potential (-0.6 V), however, the active site structure that is responsible for such high activity and selectivity has yet to be understood. Here, in this paper, we demonstrate a computational investigation combining X-ray absorption near edge structure (XANES) simulations and a mechanistic study via density functional theory (DFT) to understand the catalyst structures of this Cu catalyst during electrocatalysis and the corresponding reaction mechanisms of the key products. An integrated computational and experimental XANES analysis depicted the dynamic evolution of the catalytic site during electrocatalysis. The as-prepared, atomically dispersed Cu catalyst aggregates and forms metallic clusters/nanoparticles under electrochemical condition, which then break down to smaller oxidized clusters after electrocatalysis. The formed Cu clusters/nanoparticles showed distinct catalytic activity and selectivity as a function of particle size based on the mechanistic investigation using DFT, which is consistent with experimental observations for catalyst samples with different Cu loadings. This comprehensive study which combines experimental and computational XANES investigation, mechanistic study via DFT calculations, and experimental performance of the catalysts, provides unprecedented dynamic and mechanistic insights into the supported atomically dispersed metal catalysts for CO 2 reduction. Such strategy and details gained can further guide discovery of novel catalyst materials for CO 2 electrochemical reduction.

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

Preferential Paths of Air-water Two-phase Flow in Porous Structures with Special Consideration of Channel Thickness Effects

Accurate understanding and predicting the flow paths of immiscible two-phase flow in rocky porous structures are of critical importance for the evaluation of oil or gas recovery and prediction of rock slides caused by gas-liquid flow. A 2D phase field model was established for compressible air-water two-phase flow in heterogenous porous structures. The dynamic characteristics of air-water two-phase interface and preferential paths in porous structures were simulated. The factors affecting the path selection of two-phase flow in porous structures were analyzed. Transparent physical models of complex porous structures were prepared using 3D printing technology. Tracer dye was used to visually observe the flow characteristics and path selection in air-water two-phase displacement experiments. The experimental observations agree with the numerical results used to validate the accuracy of phase field model. The effects of channel thickness on the air-water two-phase flow behavior and paths in porous structures were also analyzed. The results indicate that thick channels can induce secondary air flow paths due to the increase in flow resistance; consequently, the flow distribution is different from that in narrow channels. This study provides a new reference for quantitatively analyzing multi-phase flow and predicting the preferential paths of immiscible fluids in porous structures.

42 ENGINEERING↗
DOE OSTI2018

ATDM/ECP Milestone Memo WBS 2.3.4.04 / SNL ATDM Data and Visualization Projects STDV04-21 - [MS1/YR2] Q3: Prototype Catalyst/ParaView in-situ viz for unsteady RV flow on ATS-1

ParaView Catalyst is an API for accessing the scalable visualization infrastructure of ParaView in an in-situ context. In-situ visualization allows simulation codes to access data post-processing operations while the simulation is running. In-situ techniques can reduce data post-processing time, allow computational steering, and increase the resolution and frequency of data output. For a simulation code to use ParaView Catalyst, adapter code needs to be created that interfaces the simulations data structures to ParaView/VTK data structures. Under ATDM, Catalyst is to be integrated with SPARC, a code used for simulation of unsteady reentry vehicle flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simulating Hydrodynamics in Cosmology with CRK-HACC

Abstract We introduce CRK-HACC, an extension of the Hardware/Hybrid Accelerated Cosmology Code (HACC), to resolve gas hydrodynamics in large-scale structure formation simulations of the universe. The new framework couples the HACC gravitational N -body solver with a modern smoothed-particle hydrodynamics (SPH) approach called conservative reproducing kernel SPH (CRKSPH). CRKSPH utilizes smoothing functions that exactly interpolate linear fields while manifestly preserving conservation laws (momentum, mass, and energy). The CRKSPH method has been incorporated to accurately model baryonic effects in cosmology simulations—an important addition targeting the generation of precise synthetic sky predictions for upcoming observational surveys. CRK-HACC inherits the codesign strategies of the HACC solver and is built to run on modern GPU-accelerated supercomputers. In this work, we summarize the primary solver components and present a number of standard validation tests to demonstrate code accuracy, including idealized hydrodynamic and cosmological setups, as well as self-similarity measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Theoretical Investigation of the Hydrogenation of Cyclohexene Catalyzed by Supported Single-Atom Sites on Redox Noninnocent LiMn 2 O 4 and Li 2 Mn 2 O 4 Surfaces

Here, the tuning of catalyst activity via stereoelectronic modulation of the active-site structure remains a grand challenge in heterogeneous catalysis. In homogeneous catalysis, the redox noninnocent ligands can be introduced to organometallic fragments to donate electrons to the metal center and fine-tune the catalytic activity. Analogously, lithium-ion battery materials, such as lithium manganese oxide (LMO), lithium titanium oxides (LTO), etc., can serve as redox noninnocent catalyst support to modulate the electronic structures of the active site via lithiation, hence tuning the catalytic activity of supported active sites. Experimentally, the Ni single-atom site was supported on LiMn 2 O 4 via oxidative grafting and exhibited no catalytic activity toward the hydrogenation of cyclohexene. After introducing additional Li into the catalyst support forming Li 2 Mn 2 O 4 , the Ni single-atom site becomes active, with the catalytic rates increasing as a function of the lithiation for Li/Mn ratios >0.9. In this paper, density functional theory (DFT) calculations are performed to study the Ni site structure via X-ray absorption near edge structure (XANES) simulations and investigate the electronic properties of Ni single-atom site before and after the addition of intercalated lithium in the LMO spinel structure. Furthermore, the study of the reaction mechanism is also carried out to understand the thermodynamically and kinetically favored pathways. XANES simulation suggests that the Ni single-atom site is likely to stay in the Li channel of the spinel support structure and form an octahedral structure. After Li intercalation, the Ni site becomes less positively charged, indicating the partial reduction of the Ni site. The simulated reaction energy profile over the LiMn 2 O 4 support exhibits high-energy barriers (1.07 eV) for the hydrogenation of cyclohexene; however, the Li 2 Mn 2 O 4 support is able to better stabilize low-coordinated ion sites and improve ion mobility, leading to lower overall energy barriers (0.64 eV). The reduced and low-coordinated ion site, thus, can better stabilize the reaction intermediates and promote hydrogenation reaction. Similarly, other transition metal ions (Fe, Co, and Cu) are also considered over the LiMn 2 O 4 and Li 2 Mn 2 O 4 catalyst supports for hydrogenation reaction.

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

The impact of thermoelectric leg geometries on thermal resistance and power output

Thermoelectric devices enable direct, solid-state conversion of heat to electricity and vice versa. Rather than designing the shape of thermoelectric units or legs to maximize this energy conversion, the cuboid shape of these legs has instead remained unchanged in large part because of limitations in the standard manufacturing process. However, the advent of additive manufacturing (a technique in which freeform geometries are built up layer-by-layer) offers the potential to create unique thermoelectric leg geometries designed to optimize device performance. This work explores this new realm of novel leg geometry by simulating the thermal and electrical performance of various leg geometries such as prismatic, hollow, and layered structures. The simulations are performed for two materials, a standard bismuth telluride material found in current commercial modules and a higher manganese silicide material proposed for low cost energy conversion in high-temperature applications. The results include the temperature gradient and electrical potential developed across individual thermoelectric legs as well as thermoelectric modules with 16 legs. Even simple hollow and layered leg geometries result in larger temperature gradients and higher output powers than the traditional cuboid structure. As a result, the clear dependence of thermal resistance and power output on leg geometry provides compelling motivation to explore additive manufacturing of thermoelectric devices.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Precipitation Hardening in Ferroelectric Ceramics

Domain wall motion in ferroics, similar to dislocation motion in metals, can be tuned by well-concepted microstructural elements. In demanding high-power applications of piezoelectric materials, the domain wall motion is considered as a lossy hysteretic mechanism that should be restricted. Current applications for so-called hard piezoelectrics are abundant and hinge on the use of an acceptor-doping scheme. However, this mechanism features severe limitations due to enhanced mobility of oxygen vacancies at moderate temperatures. By analogy with metal technology, the authors present here a new solution for electroceramics, where precipitates are utilized to pin domain walls and improve piezoelectric properties. Through a sequence of sintering, nucleation, and precipitate growth, intragranular precipitates leading to a fine domain structure are developed as shown by transmission electron microscopy, piezoresponse force microscopy, and phase-field simulation. This structure impedes the domain wall motion as elucidated by electromechanical characterization. As a result, the mechanical quality factor is increased by ≈50% and the hysteresis in electrostrain is suppressed considerably. This is even achieved with slightly increased piezoelectric coefficient and electromechanical coupling factor. This novel process can be smoothly implemented in industrial production processes and is accessible to simple laboratory experimentation for microstructure optimization and implementation in various ferroelectric systems.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Halogen Bond Structure and Dynamics from Molecular Simulations

Halogen bonding has emerged as an important noncovalent interaction in a myriad of applications, including drug design, supramolecular assembly, and catalysis. Current understanding of the halogen bond is informed by electronic structure calculations on isolated molecules and/or crystal structures that are not readily trans- ferable to liquids and disordered phases. To address this issue, we present a first-principles simulation-based approach for quantifying halogen bonds in molecular systems rooted in an understanding of nuclei-nuclei and electron-nuclei spatial correlations. We then demonstrate how this approach can be used to quantify the structure and dynamics of halogen bonds in condensed phases, using solid and liquid molecular chlorine as prototypical examples with high concentrations of halogen bonds. We close with a discussion of how the knowledge generated by our first-principles approach may inform the development of classical empirical models, with a consistent representation of halogen bonding.

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

Simulation Study of the Influence of Experimental Variations on the Structure and Quality of Plasma Liners

Simulation studies of a section of a spherically imploding plasma liner, formed by the merger of six hypersonic plasma jets, have been performed at conditions relevant to the Plasma Liner Experiment (PLX) [S. C. Hsu et al., IEEE Trans. Plasma Sci. 46, 1951 (2018)]. The main aim of simulations was the sensitivity study of the detailed structure of plasma liners and their global properties to experimental mass variations and timing jitter across the six plasma jets. Experimentally observable synthetic quantities have been computed using simulation data and compared with the available experimental data. Simulations predicted that the primary oblique shock wave structure is preserved at small experimental variations. At later phases of the liner implosion, primary shocks and, especially, secondary shocks are more sensitive to experimental variations. These conclusions follow from the simulation data as well as the experimental CCD camera images. Small displacements of the shock wave structures may also cause significant changes in the synthetic interferometer data at early time. Our studies showed that the global properties of the plasma liners (averaged Mach number and averaged ram pressure along leading edges of plasma liners) are not very sensitive to experimental variations. Simulation data of the liner structure were largely confirmed by the PLX experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Structure and dynamics of water confined in nanoporous carbon

We study the structure and dynamics of water confined in a nanoporous graphitic carbon structure using molecular-dynamics (MD) simulations. The carbon structure is generated by a reactive MD simulation of oxidation of a silicon carbide nanoparticle. We embed water molecules in the nanopores and study structural and dynamical properties of nanoconfined water as a function of temperature. MD simulation results indicate the presence of high-density water (HDW) and low-density water. Radial distribution functions and spatial density functions indicate that the second solvation shell of the HDW is broken. We calculate the self-diffusion coefficient of confined water molecules as a function of temperature and find a significant decrease in the diffusion of water molecules around T = 190K. The cage correlation function c(t) of confined water molecules at T = 200K exhibits stretched exponential decay, c(t) = exp[–(t/τ) β ], with β = 0.43, which matches exactly with the theoretical prediction β = 3/7. Finally, the self-intermediate scattering function at T = 200K indicates differences in small-scale and large-scale dynamics of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2019

High-throughput simulations for insight into grain boundary structure-property relationships and other complex microstructural phenomena

High-throughput simulations can be a powerful tool in the discovery of new materials and behaviors. As part of a special issue on Rising Stars in Computational Materials Science, this paper uses the work of the author to show how high-throughput simulations have had an impact in grain boundary structure-property relationships and other complex microstructural phenomena. The work demonstrates how new tools designed to analyze large datasets produced by high-throughput simulations are enabling comprehensive grain boundary structure-property relationships to be obtained. High-throughput simulations are also used to demonstrate the impact descriptive and inferential statistics has had in extracting key aspects of deformation in metallic glasses. Finally, several different examples are used to show a balanced approach between simulations designed to survey and simulations designed for detailed analysis. Together, the two approaches provide a comprehensive picture of the variety of behaviors that exist, while ensuring that the physics underlying the behaviors are thoroughly understood.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simulation of double buffer layer on CIGS solar cell with SCAPS software

In a CIGS thin film solar cell, the buffer layer is interposed between the absorber layer and the window layer, which plays an important role in interface electricity. CdS is often used as a buffer layer in CIGS cells. Beside the toxicity of Cd, the main drawback of CdS is its relatively narrow band gap (2.4 eV), which can lead to the current loss due to parasitic absorption. In order to reduce the use of CdS, we try other buffer layers instead of CdS and optimize buffer layer structure. Numerical simulation has been used to explore the effect of different buffer layer materials (CdS and ZnSe) on the performance of CuInGaSe{sub 2} solar cells with SCAPS-1D software. The main photovoltaic parameters of the analog device: open circuit voltage, short circuit current density, fill factor and conversion efficiency (η) were compared and analyzed. The results show that different buffer materials and thickness have great influence on the performance of the cell. In this simulation, CdS buffer layer was replaced by ZnSe material, which has a wider band gap. However, the efficiency of this kind of cell is lower than that of CIGS device with CdS buffer layer because of the bad energy band alignment. In order to improve the short-wave absorption and the efficiency of CIGS cell, we simulated a novel CdS/ZnSe double buffer layer structure. The efficiency of the optimized cell with double buffer layer increases from 18.47% (conventional cell with the same CdS buffer layer thickness) to 19.01%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗