Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Comparison of fluid neutral models for one-dimensional plasma edge modeling with a finite volume solution of the Boltzmann equation

We derive fluid neutral approximations for a simplified 1D edge plasma model, suitable to study the neutral behavior close to the target of a nuclear fusion divertor, and compare its solutions to the solution of the corresponding kinetic Boltzmann equation. The plasma is considered as a fixed background extracted from a detached 2D simulation. We show that the Maxwellian equilibrium distribution is already obtained very close to the target, justifying the use of a fluid approximation. We compare three fluid neutral models: (i) a diffusion model; (ii) a pressure-diffusion model (i.e., a combination of a continuity and momentum equation) assuming equal neutral and ion temperatures; and (iii) the pressure-diffusion model coupled to a neutral energy equation taking into account temperature differences between neutrals and ions. Partial reflection of neutrals reaching the boundaries is included in both the kinetic and fluid models. We propose two methods to obtain an incident neutral flux boundary condition for the fluid models: one based on a diffusion approximation and the other assuming a truncated Chapman-Enskog distribution. The pressure-diffusion model predicts the plasma sources very well. The diffusion boundary condition gives slightly better results overall. Although including an energy equation still improves the results, the assumption of equal ion and neutral temperature already gives a very good approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Resonance coupling in plasmonic nanomatryoshka homo- and heterodimers

Here, we examine the electromagnetic (EM) energy coupling and hybridization of plasmon resonances between closely spaced concentric nanoshells known as “nanomatryoshka” (NM) units in symmetric and antisymmetric compositions using the Finite Difference Time Domain (FDTD) analysis. Utilizing plasmon hybridization model, we calculated the energy level diagrams and verified that, in the symmetric dimer (in-phase mode in a homodimer), plasmonic bonding modes are dominant and tunable within the considered bandwidth. In contrast, in the antisymmetric dimer (out-of-phase mode in a heterodimer), due to the lack of the geometrical symmetry, new antibonding modes appear in the extinction profile, and this condition gives rise to repeal of dipolar field coupling. We also studied the extinction spectra and positions of the antibonding and bonding modes excited due to the energy coupling between silver and gold NM units in a heterodimer structure. Our analysis suggest abnormal shifts in the higher energy modes. We propose a method to analyze the behavior of multilayer concentric nanoshell particles in an antisymmetric orientation employing full dielectric function calculations and the Drude model based on interband transitions in metallic components. This study provides a method to predict the behavior of the higher energy plasmon resonant modes in entirely antisymmetric structures such as compositional heterodimers.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Multi-scale Material Parameter Identification Using LS-DYNA® and LS-OPT®

Ever-tightening regulations on fuel economy and carbon emissions demand continual innovation in finding ways for reducing vehicle mass. Classical methods for computational mass reduction include sizing, shape and topology optimization. One of the few remaining options for weight reduction can be found in materials engineering and material design optimization. Apart from considering different types of materials by adding material diversity, an appealing option in automotive design is to engineer steel alloys for the purpose of reducing thickness while retaining sufficient strength and ductility required for durability and safety. Such a project was proposed and is currently being executed under the auspices of the United States Automotive Materials Partnership (USAMP) funded by the Department of Energy. Under this program, new steel alloys (Third Generation Advanced High Strength Steel or 3GAHSS) are being designed, tested and integrated with the remaining design variables of a benchmark vehicle Finite Element model. In this project the principal phases identified are (i) material identification, (ii) formability optimization and (iii) multi-disciplinary vehicle optimization. This paper serves as an introduction to the LS-OPT methodology and therefore mainly focuses on the first phase, namely an approach to integrate material identification using material models of different length scales. For this purpose, a multi-scale material identification strategy, consisting of a Crystal Plasticity (CP) material model and a Homogenized State Variable (SV) model, is discussed and demonstrated. The paper concludes with proposals for integrating the multi-scale methodology into the overall vehicle design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

CODE-TO-CODE BENCHMARK STUDY FOR THERMAL STRESS MODELING AND PRELIMINARY ANALYSIS OF THE HIGH-TEMPERATURE SINGLE HEAT PIPE EXPERIMENT

Microreactors are very small nuclear reactors with typical thermal-energy output of up to 20 MWth. The microreactor concept is gaining more and more attention for safe, robust, and reliable supply of electricity for remote locations and to meet industrial process heat needs. Heat-pipe-cooled microreactor is one of the microreactors being investigated at Idaho National Laboratory. In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios for the heat-pipe-cooled microreactor designs. This paper discusses finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory, using two commercial software packages, Abaqus and STAR-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the locations where the largest temperature gradients appeared between heat pipe and electric heater. The temperature fields have good agreements between Abaqus and STAR-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The Toughness of Interlocking Metasurfaces

Interlocking metasurfaces (ILMs) are arrays of autogenous latching unit cells patterned across a surface. These create structural joints similar to bioinspired suture joints but patterned over a 2D surface rather than a 1D seam. This enables ILMs to be an alternative to conventional joining technologies such as bolts, welds, and adhesives. However, compared to conventional joining methods, relatively little is known of the engineering considerations for designing structural ILMs. Herein, the interfacial toughness of an archetypal ILM is examined for the first time. Under the conditions studied here, the ILM is substantially tougher than the material from which it is made, in this case, exhibiting up to a 50% increase in interfacial crack initiation energy over the solid base material, a photocured 3D‐printed polymer. Through experimental tests using in‐situ digital image correlation along with complementary computational analyses, the mechanism of toughening in the ILM structure and the origins of toughness anisotropy are revealed. The increase in toughness is associated with cross‐cell interactions, that is, load‐sharing across unit cells, which give rise to a finite process zone length with different effective material properties. In this way, ILM toughening is analogous to crack blunting in ductile materials or fiber bridging in composites; yet here, the ILM is composed of a single‐phase base material and so the architected toughening is geometric in nature and hence amenable to future topological optimization.

fracture↗
DOE OSTI2025

Atwood effects on nonlocality of the scalar transport closure in Rayleigh-Taylor mixing

The importance of nonlocality is assessed in modeling mean scalar transport for turbulent Rayleigh-Taylor (RT) mixing at different Atwood numbers. Building on the two-dimensional incompressible work of Lavacot et al. [J. Fluid Mech. 985, A47 (2024)], the present work extends the macroscopic forcing method to variable density problems in three-dimensional space to measure moments of the generalized eddy diffusivity kernel in RT mixing for increasing Atwood numbers (𝐴 = 0.05, 0.3, 0.5, 0.8). It is found that as 𝐴 increases, (1) the eddy diffusivity moments become asymmetric and (2) the higher-order eddy diffusivity moments become larger relative to the leading-order diffusivity, indicating that nonlocality becomes more important at higher 𝐴. There is a particularly strong temporal nonlocality at higher 𝐴, suggesting stronger history effects. In conclusion, the implications of these findings for closure modeling for finite-Atwood RT are discussed.

general physics↗
DOE OSTI2015

SU-E-T-31: A Fast Finite Size Pencil Beam (FSPB) Convolution Algorithm for a New Co-60 Arc Therapy Machine

Purpose: Present a fast Finite Size Pencil Beam (FSPB) convolution algorithm for a new Co-60 arc therapy machine. The FSPB algorithm accounts for (i) strong angular divergence (short SAD), (ii) heterogeneity effect for primary attenuation, and (iii) source energy spectrum. Methods: The FSPB algorithm is based on a 0.5×0.5-cm2 dose kernel calculated using the GEPTS (Gamma Electron and Positron Transport System) Monte Carlo code. The dose kernel is tabulated using a thin XYZ mesh (0.1 mm steps in lateral directions) for radius less than 1 cm and using an RZ mesh (with varying steps) for larger radial distance. To account for SSD effect, 11 dose kernels with SSDs varying between 30 cm to 80 cm are calculated. Maynord factor and “lateral stretching” are applied to account for differences between closest and actual SSD. Appropriate rotations and second order interpolation are used to calculate the dose from a given beamlet to a point. Results: Accuracy: Dose distributions in water with 80 cm SSD are calculated using the new FSPB convolution algorithm and full Monte Carlo simulation (gold standard). Figs 1–4 show excellent agreements between FSPB and Monte Carlo calculations for different field sizes and at different depths. The dose distribution for a prostate case is calculated using FSPB (Fig.5). Sixty conformal beams with rectum blocking are assumed. Figs 6–8 show the comparison with Monte Carlo simulation based on the same beam apertures. The excellent agreement demonstrates the accuracy of the new algorithm in handling SSD variation, oblique incidence, and scatter contribution.Speed: The FSPB convolution algorithm calculates 28 million dose points per second using a single 2.2-GHz CPU. The present algorithm is seven times faster than a similar algorithm from Gu et al. (Phys. Med. Biol. 54, 2009, 6287–6297). Conclusion: A fast and accurate FSPB convolution algorithm was developed and benchmarked.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

AZTRAN and AZKIND verification with a BWR study case

The AZTLAN Platform project is a national initiative led by the National Institute of Nuclear Research which gather the main public institutions of higher education in Mexico: the National Polytechnic Institute, the National Autonomous University of Mexico and the Autonomous Metropolitan University. It consists in the development of a modeling platform for the analysis and design of nuclear reactors. This project aims to modernize, improve and integrate the neutronic, thermohydraulics and thermomechanical codes, developed in the national institutions, in a platform integrated and developed by Mexican experts. The development of AZTLAN Platform is divided technically into four working groups: neutronics, thermohydraulics, coupling and uncertainty analysis, and users group. The present work is the result of the activities of the Users Group focused on the follow-up, testing and use of the different codes (modules) being developed under the AZTLAN project; with the objective of providing feedback and to suggest development priorities to the team of developers, as well as to promote collaboration and synergy between the participating institutions. In particular, in this work, the AZTRAN and the AZKIND modules were verified and tested against the results of well-known codes of the nuclear industry for boiling water reactor (BWR) models. AZTRAN is a 3D transport code which solves numerically the multi-group time independent neutron transport equation in Cartesian geometry, based on the discrete ordinates (SN) method and the RTN-0 nodal scheme. AZKIND is a 3D diffusion code that solves numerically the time dependent neutron diffusion equations in Cartesian geometry, based on a Galerkin Finite Element Method using the RTN-0 polynomial nodal approximation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Flexosensitive polarization vortices in thin ferroelectric films

In this work, we present a theoretical study on the influence of the flexoelectric coupling on the spatial distribution and temperature behavior of the spontaneous polarization for several types of stable domain structure in thin ferroelectric films, such as stripe domains and vortices. Finite-element modeling (FEM) for BaTiO 3 films and analytical calculations within the Landau-Ginzburg-Devonshire approach reveals that an out-of-plane polarization component can be very sensitive to the flexoelectric coupling for periodic quasi-2D stripe domains and 3D vortex-antivortex structures. However, the influence is rather different for these structures. The flexoelectric coupling increases significantly the amplitude of a small out-of-plane polarization component in the stripe domains, and the “up” or “down” direction of the component is defined by the sign of the flexoelectric coefficients. Concerning the vortex-antivortex pairs, their antivortices with in-plane antiсirculation have smooth wide dipolar cores through the entire film, whose shape and other features are almost insensitive to the coupling. The vortices with in-plane vorticity have spikelike cores with an out-of-plane quadrupolar moment induced by the flexoelectric coupling. The cores are located near the film–dead-layer interfaces. FEM results corroborated by analytical calculations prove that a change of the flexoelectric coefficient sign leads to a reorientation of the core axial polarization, making the flexosensitive 3D vortices similar to the recently introduced “flexons” in cylindrical nanoparticles. The relatively wide temperature range (from 200 to 400 K) of the flexosensitive vortices' existence gives us confidence that they can be observed experimentally in thin ferroelectric films by scanning-probe and nonlinear optical microscopy methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

A NEW METHOD TO DIRECTLY MEASURE THE JEANS SCALE OF THE INTERGALACTIC MEDIUM USING CLOSE QUASAR PAIRS

Although the baryons in the intergalactic medium (IGM) trace dark matter fluctuations on megaparsec scales, on smaller scales ∼100 kpc, fluctuations are suppressed because the finite temperature gas is pressure supported against gravity, analogous to the classical Jeans argument. This Jeans filtering scale, which quantifies the small-scale structure of the IGM, has fundamental cosmological implications. First, it provides a thermal record of heat injected by ultraviolet photons during cosmic reionization events, and thus constrains the thermal and reionization history of the universe. Second, the Jeans scale determines the clumpiness of the IGM, a critical ingredient in models of cosmic reionization. Third, it sets the minimum mass scale for gravitational collapse from the IGM, and hence plays a pivotal role in galaxy formation. Unfortunately, it is extremely challenging to measure the Jeans scale via the standard technique of analyzing purely longitudinal Lyα forest spectra, because the thermal Doppler broadening of absorption lines along the line-of-sight, is highly degenerate with Jeans smoothing. In this work, we show that the Jeans filtering scale can be directly measured by characterizing the coherence of correlated Lyα forest absorption in close quasar pairs, with separations small enough ∼100 kpc to resolve it. We present a novel technique for this purpose, based on the probability density function (PDF) of phase angle differences of homologous longitudinal Fourier modes in close quasar pair spectra. A Bayesian formalism is introduced based on the phase angle PDF, and Markov Chain Monte Carlo techniques are used to characterize the precision of a hypothetical Jeans scale measurement, and explore degeneracies with other thermal parameters governing the IGM. A semi-analytical model of the Lyα forest is used to generate a large grid (500) of thermal models from a dark matter only simulation. Our full parameter study indicates that a realistic sample of only 20 close quasar pair spectra can pinpoint the Jeans scale to ≅ 5% precision, independent of the amplitude T{sub 0} and slope γ of the temperature-density relation of the IGM T=T{sub 0}(ρ/ ρ-bar ){sup γ-1}. This exquisite sensitivity arises because even long-wavelength one-dimensional Fourier modes ∼10 Mpc, i.e., two orders of magnitude larger than the Jeans scale, are nevertheless dominated by projected small-scale three-dimensional (3D) power. Hence phase angle differences between all modes of quasar pair spectra actually probe the shape of the 3D power spectrum on scales comparable to the pair separation. We show that this new method for measuring the Jeans scale is unbiased and is insensitive to a battery of systematics that typically plague Lyα forest measurements, such as continuum fitting errors, imprecise knowledge of the noise level and/or spectral resolution, and metal-line absorption.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Improved multislope MUSCL reconstruction on unstructured grids for shallow water equations

In shallow water flow and transport modeling, the monotonic upstream-centered scheme for conservation laws (MUSCL) is widely used to extend the original Godunov scheme to second-order accuracy. The most important step in MUSCL-type schemes is MUSCL reconstruction, which calculate-extrapolates the values of independent variables from the cell center to the edge. The monotonicity of the scheme is preserved with the help of slope limiters that prevent the occurrence of new extrema during reconstruction. On structured grids, the calculation of the slope is straightforward and usually based on a 2-point stencil that uses the cell centers of the neighbor cell and the so-called far-neighbor cell of the edge under consideration. On unstructured grids, the correct choice for the upwind slope becomes nontrivial. 2 novel total variation diminishing schemes are developed based on different techniques for calculating the upwind slope and the downwind slope. An additional treatment that stabilizes the scheme is discussed. The proposed techniques are compared to 2 existing MUSCL reconstruction techniques, and a detailed discussion of the results is given. Finally, it is shown that the proposed MUSCL reconstruction schemes obtain more accurate results with less numerical diffusion and higher efficiency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Investigation of nonlinear mechanical behavior of two superfine-grained graphites with DIC assisted disc splitting test

Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.

Lin, Lianshan [ORNL] (ORCID:0000000203399219)↗
DOE OSTI2023

Tensile behavior of diffusion bonded AA6061 - AA6061 with variation in cooling method

Hot isostatically pressed AA6061 cladding is an important structural component of the high performance, Zr-laminated U-10Mo monolithic fuel system for the application in research and test reactors. In this study, the mechanical behavior of two diffusion bonded aluminum alloy, AA6061, was examined using tensile testing. Solid-to-solid diffusion bonding between two pieces of AA6061 was performed by isothermal annealing at 560 °C for 1.5 h, and diffusion couples were subsequently cooled via three different cooling methods: furnace cooling, air cooling, and water quenching. Dog-bone shaped tensile specimens, with 10 mm in gauge length (with diffusion bonded interface in the middle), and 1.5 × 1.5 mm 2 gauge cross-sections, were fabricated from the diffusion bonded AA6061 by electro-discharge machining. Yield strength (% EL at failure) of furnace cooled, air cooled and water quenched tensile specimens determined was 82 ~ 89 MPa (10 ~ 30%), 112 ~ 116 MPa (10 ~ 14%), and 149 ~ 164 MPa (10 ~ 17%), respectively. This variation in mechanical behavior was examined with cooling-rate dependent, concentrated precipitation of Mg 2 Si at the diffusion bonded interface, with due respect for mechanical properties of the AA6061 alloy that inherently vary as a function of cooling rate from 560 °C. Finite element analysis using ABAQUS was employed to augment experimental findings with the appropriate microstructural constituents and alloy properties. Results suggest that the strength is dominated by matrix/bulk properties of AA6061, while ductility is strongly influenced by the cooling method dependent presence of Mg 2 Si precipitates at the interface.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Fundamental gaps with approximate density functionals: The derivative discontinuity revealed from ensemble considerations

The fundamental gap is a central quantity in the electronic structure of matter. Unfortunately, the fundamental gap is not generally equal to the Kohn-Sham gap of density functional theory (DFT), even in principle. The two gaps differ precisely by the derivative discontinuity, namely, an abrupt change in slope of the exchange-correlation energy as a function of electron number, expected across an integer-electron point. Popular approximate functionals are thought to be devoid of a derivative discontinuity, strongly compromising their performance for prediction of spectroscopic properties. Here we show that, in fact, all exchange-correlation functionals possess a derivative discontinuity, which arises naturally from the application of ensemble considerations within DFT, without any empiricism. This derivative discontinuity can be expressed in closed form using only quantities obtained in the course of a standard DFT calculation of the neutral system. For small, finite systems, addition of this derivative discontinuity indeed results in a greatly improved prediction for the fundamental gap, even when based on the most simple approximate exchange-correlation density functional – the local density approximation (LDA). For solids, the same scheme is exact in principle, but when applied to LDA it results in a vanishing derivative discontinuity correction. This failure is shown to be directly related to the failure of LDA in predicting fundamental gaps from total energy differences in extended systems.

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

Flaw Tolerance of Octet Truss Lattices with Random Flaws

We investigate the behavior of octet truss lattices in which some proportion of the bars are randomly removed from a perfect lattice. Using a lattice specimen of 9 × 9 × 9 octet truss unit cells as the nominal design, we simulate the response of prefect and flawed lattices with a finite element method implementation of rod theory. For each of two simple load cases, we perform a Monte Carlo analysis where we consider many instantiations of flawed lattices where between 1% and 20% of the bars are missing, and compare the axial loads in the remaining bars to those in the perfect lattice. We also compare the overall stiffness of the flawed lattices to the perfect lattice. Our results show progressive deterioration of lattice strength and stiffness with increasing proportion of missing bars, as expected, with the location and magnitude of the highest-loaded (most failure-prone) bars depending more on the load case and bar orientation than on the location and connectivity of missing bars. As such, the results reported herein may not be directly transferable to different lattice geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A coupled cluster framework for electrons and phonons

Here, we describe a coupled cluster framework for coupled systems of electrons and harmonic phonons. Neutral and charged excitations are accessed via the equation-of-motion version of the theory. Benchmarks on the Hubbard–Holstein model allow us to assess the strengths and weaknesses of different coupled cluster approximations, which generally perform well for weak to moderate coupling. Finally, we report progress toward an implementation for ab initio calculations on solids and present some preliminary results on finite-size models of diamond with a linear electron–phonon coupling. We also report the implementation of electron–phonon coupling matrix elements from crystalline Gaussian type orbitals within the PySCF program package.

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

Enabling topography-resolving structural dynamic contact simulation

Damping of structures and systems is often dominated by frictional dissipation in connections, the prediction of which remains a longstanding scientific challenge. Previous studies have shown that the actual topography of contact interfaces may have a strong effect, especially in the partial slip/liftoff regime. We recently proposed a multi-scale method, which couples finite element and boundary element modeling. The primary benefit of this approach is that it permits to analyze the effect of the actual contact topography on the dynamics of jointed structures. While this multi-scale modeling method was initially developed for quasi-static analysis, we demonstrate herein how it can be used for time step integration and Harmonic Balance analysis. We cross-verify those fully dynamic analysis methods against each other and quasi-static results, for the S4 Beam benchmark. We compare the multi-scale method against state-of-the-art full-FE analysis, in terms of numerical damping and computational performance. Some discrepancy is found to be of physical origin. Depending on the load history, it is shown that the system settles to a slightly different equilibrium. Finally, transient multi-scale simulations enable the prediction of this interesting phenomenon, for the first time, for a structure with bolted joints.

Frictional-unilateral contact↗