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 1,477 records · Page 82

Fluctuating Hydrodynamics Methods for Dynamic Coarse-Grained Implicit-Solvent Simulations in LAMMPS

We present a software package integrated with the molecular dynamics software LAMMPS for fluctuating hydrodynamics simulations of fluid-structure interactions subject to thermal fluctuations. The package is motivated to provide dynamic thermostats to extend implicit-solvent coarse-grained (IS-CG) models by incorporating kinetic contributions from the solvent to facilitate their use in a wider range of applications. To capture the thermal and hydrodynamic contributions of the solvent to dynamics, we introduce momentum conserving thermostats and computational methods based on fluctuating hydrodynamics and the stochastic Eulerian Lagrangian method (SELM). SELM couples the coarse-grained microstructure degrees of freedom to continuum stochastic fields to capture both the relaxation of hydrodynamic modes and thermal fluctuations. Features of the SELM software include (i) numerical time-step integrators for SELM fluctuating hydrodynamics in inertial and quasi-steady regimes, (ii) Lees--Edwards-style methods for imposing shear, (iii) a Java-based graphical user interface (GUI) for setting up models and simulations, (iv) standardized XML formats for parametrization and data output, and (v) standardized formats VTK for continuum fields and microstructures. The SELM software package facilitates for pre-established models in LAMMPS easy adoption of the SELM fluctuating hydrodynamics thermostats. We offer here an overview of the SELM software package, computational methods, and applications.

97 MATHEMATICS AND COMPUTING↗

Fingerprints of anomalous primordial Universe on the abundance of large scale structures

We study the predictions of anomalous inflationary models on the abundance of structures in large scale structure observations. The anomalous features encoded in primordial curvature perturbation power spectrum are (a): localized feature in momentum space, (b): hemispherical asymmetry and (c): statistical anisotropies. We present a model-independent expression relating the number density of structures to the changes in the matter density variance. Models with localized feature can alleviate the tension between observations and numerical simulations of cold dark matter structures on galactic scales as a possible solution to the missing satellite problem. In models with hemispherical asymmetry we show that the abundance of structures becomes asymmetric depending on the direction of observation to sky. In addition, we study the effects of scale-dependent dipole amplitude on the abundance of structures. Using the quasars data and adopting the power-law scaling k{sup n{sub A}-1} for the amplitude of dipole we find the upper bound n{sub A} < 0.6 for the spectral index of the dipole asymmetry. In all cases there is a critical mass scale M{sub c} in which for M M{sub c}) the enhancement in variance induced from anomalous feature decreases (increases) the abundance of dark matter structures in Universe.

79 ASTRONOMY AND ASTROPHYSICS↗

Phase coexistence and domain configuration in Pb(Mg 1/3 Nb 2/3 )O 3 -0.34PbTiO 3 single crystal revealed by synchrotron-based X-ray diffractive three-dimensional reciprocal space mapping and piezoresponse force microscopy

The crystalline phases and domain configuration in the morphotropic phase boundary composition Pb(Mg 1/3 Nb 2/3 )O 3 -0.34PbTiO 3 (PMN-0.34PT) single crystal have been investigated by synchrotronbased X-ray 3D Reciprocal Space Mapping (3D-RSM) and Piezoresponse Force Microscopy. The coexistence of tetragonal (T) and monoclinic MC phases in this PMN-0.34PT single crystal is confirmed. The affiliation of each diffraction spot in the 3D-RSM was identified with the assistance of qualitative simulation. Most importantly, the twinning structure between different domains in such a mixed phase PMN-PT crystal is firmly clarified, and the spatial distribution of different twin domains is demonstrated. Additionally, the lattice parameters of T and MC phases in PMN-0.34PT single crystal as well as the tilting angles of crystal lattices caused by the interfacial lattice mismatch are determined.

36 MATERIALS SCIENCE↗

Short- and medium-range orders in Al90Tb10 glass and their relation to the structures of competing crystalline phases

Molecular dynamics simulations using an interatomic potential developed by artificial neural network deep machine learning are performed to study the local structural order in Al90Tb10 metallic glass. We show that more than 80% of the Tb-centered clusters in Al90Tb10 glass have short-range order (SRO) with their 17 first coordination shell atoms stacked in a ‘3661’ or ‘15551’ sequence. Medium-range order (MRO) in Bergman-type packing extended out to the second and third coordination shells is also clearly observed. Analysis of the network formed by the ‘3661’ and ‘15551’ clusters show that ~82% of such SRO units share their faces or vertexes, while only ~6% of neighboring SRO pairs are interpenetrating. Such a network topology is consistent with the Bergman-type MRO around the Tb-centers. Moreover, crystal structure searches using genetic algorithm and the neural network interatomic potential reveal several low-energy metastable crystalline structures in the composition range close to Al 90 Tb 10 . Some of these crystalline structures have the ‘3661’ SRO while others have the ‘15551’ SRO. While the crystalline structures with the ‘3661’ SRO also exhibit the MRO very similar to that observed in the glass, the ones with the ‘15551’ SRO have very different atomic packing in the second and third shells around the Tb centers from that of the Bergman-type MRO observed in the glassy phase.

36 MATERIALS SCIENCE↗

SPONTANEOUS FORMATION OF SURFACE MAGNETIC STRUCTURE FROM LARGE-SCALE DYNAMO IN STRONGLY STRATIFIED CONVECTION

We report the first successful simulation of spontaneous formation of surface magnetic structures from a large-scale dynamo by strongly stratified thermal convection in Cartesian geometry. The large-scale dynamo observed in our strongly stratified model has physical properties similar to those in earlier weakly stratified convective dynamo simulations, indicating that the α {sup 2}-type mechanism is responsible for the dynamo. In addition to the large-scale dynamo, we find that large-scale structures of the vertical magnetic field are spontaneously formed in the convection zone (CZ) surface only in cases with a strongly stratified atmosphere. The organization of the vertical magnetic field proceeds in the upper CZ within tens of convective turnover time and band-like bipolar structures recurrently appear in the dynamo-saturated stage. We consider several candidates to be possibly be the origin of the surface magnetic structure formation, and then suggest the existence of an as-yet-unknown mechanism for the self-organization of the large-scale magnetic structure, which should be inherent in the strongly stratified convective atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermodynamics of concentrated solid solution alloys

This study reviews the three main approaches for predicting the formation of concentrated solid solution alloys (CSSA) and for modeling their thermodynamic properties, in particular, utilizing the methodologies of empirical thermo-physical parameters, CALPHAD method, and first-principles calculations combined with hybrid Monte Carlo/Molecular Dynamics (MC/MD) simulations. In order to speed up CSSA development, a variety of empirical parameters based on Hume-Rothery rules have been developed. Herein, these parameters have been systematically and critically evaluated for their efficiency in predicting solid solution formation. The phase stability of representative CSSA systems is then illustrated from the perspectives of phase diagrams and nucleation driving force plots of the σ phase using CALPHAD method. The temperature-dependent total entropies of the FCC, BCC, HCP, and σ phases in equimolar compositions of various systems are presented next, followed by the thermodynamic properties of mixing of the BCC phase in Al-containing and Ti-containing refractory metal systems. First-principles calculations on model FCC, BCC and HCP CSSA reveal the presence of both positive and negative vibrational entropies of mixing, while the calculated electronic entropies of mixing are negligible. Temperature dependent configurational entropy is determined from the atomic structures obtained from MC/MD simulations. Current status and challenges in using these methodologies as they pertain to thermodynamic property analysis and CSSA design are discussed.

36 MATERIALS SCIENCE↗

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

58 GEOSCIENCES↗

Predissociation dynamics of lithium iodide

The predissociation dynamics of lithium iodide (LiI) in the first excited A-state is investigated for molecules in the gas phase and embedded in helium nanodroplets, using femtosecond pump-probe photoionization spectroscopy. In the gas phase, the transient Li{sup +} and LiI{sup +} ion signals feature damped oscillations due to the excitation and decay of a vibrational wave packet. Based on high-level ab initio calculations of the electronic structure of LiI and simulations of the wave packet dynamics, the exponential signal decay is found to result from predissociation predominantly at the lowest avoided X-A potential curve crossing, for which we infer a coupling constant V{sub XA} = 650(20) cm{sup −1}. The lack of a pump-probe delay dependence for the case of LiI embedded in helium nanodroplets indicates fast droplet-induced relaxation of the vibrational excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collective modes in two-dimensional one-component-plasma with logarithmic interaction

The collective modes of a familiar two-dimensional one-component-plasma with the repulsive logarithmic interaction between the particles are analysed using the quasi-crystalline approximation (QCA) combined with the molecular dynamic simulation of the equilibrium structural properties. It is found that the dispersion curves in the strongly coupled regime are virtually independent of the coupling strength. Arguments based on the excluded volume consideration for the radial distribution function allow us to derive very simple expressions for the dispersion relations, which show excellent agreement with the exact QCA dispersion over the entire domain of wavelengths. Comparison with the results of the conventional fluid analysis is performed, and the difference is explained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temperature structure of the intracluster medium from smoothed-particle hydrodynamics and adaptive-mesh refinement simulations

Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10%-30%. The largest bias originates from both violation of hydrostatic equilibrium (HE) and an additional temperature bias caused by inhomogeneities in the X-ray-emitting intracluster medium (ICM). To elucidate this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle hydrodynamics (SPH) or adaptive-mesh refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in nonradiative simulations, whereas it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with that observed in the central regions of clusters. Around R {sub 500} the temperature inhomogeneities of the SPH simulations can generate twice the typical HE mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high-sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization and structural analysis of a thermophilic GH11 xylanase from compost metatranscriptome

Xylanase is efficient for xylan degradation and widely applied in industries. We found a GH11 family xylanase (Xyn11A) with high thermostability and catalytic activity from compost metatranscriptome. This xylanase has the optimal reaction temperature at 80 °C with the activity of 2907.3 U/mg. The X-ray crystallographic structure shows a typical “right hand” architecture, which is the characteristics of the GH11 family enzymes. Comparing it with the mesophilic XYN II, a well-studied GH11 xylanase from Trichoderma reesei, Xyn11A is more compact with more H-bonds. Our mutagenic results show that the electrostatic interactions in the thumb and palm region of Xyn11A could result in its high thermostability and activity. Introducing a disulfide bond at the N-terminus further increased its optimal reaction temperature to 90 °C with augmented activity.

59 BASIC BIOLOGICAL SCIENCES↗

Detecting underscreening and generalized Kirkwood transitions in aqueous electrolytes

We establish the connection between the measured small angle x-ray scattering signal and the charge–charge correlations underlying Kirkwood transitions (KTs) in 1:1, 2:1, and 3:1 aqueous electrolytes. These measurements allow us to obtain underscreening lengths for bulk electrolytes independently verified by theory and simulations. Furthermore, we generalize the concept of KTs beyond those theoretically predicted for 1:1 electrolytes, which involves the inverse screening length, a 0 , and the inverse periodicity length, Q 0 . Above the KTs, we find a universal scaling of a 0 ∝ c – $\sqrt{ζ/3}$ and Q 0 ∝ c 1/3 for the studied electrolyte solutions, where ζ is the ionic strength factor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THERMAL/STRUCTURAL ANALYSIS OF THE AXION QUANTUM METROLOGY CAVITY AND ITS COMPONENTS

This research was centered around maximizing the capability to cool dielectric material within a containment unit, or Photonic Band Gap (PBG) cavity, designed for detecting axion dark matter and identifying the unit’s thermal properties. There are multiple types of PBG cavities, but the latest version that axion researchers wish to use has been theorized to contain possible issues related to its thermal properties. Thermal conductivity is an issue with the dielectric material because it is made from alumina which is highly insulative. This is important since the research is being done in a cryogenic environment and the thermal noise affects the quantum bit used for detecting the axion to photon conversion process. Therefore, any improvements to this unit should be justified and implemented but are not entirely limited to thermal contact related aspects of PBG cavities. A prospect of using sapphire in place of alumina also exists, but this is a more expensive, less tested, and more elusive material to justify the creation of a full dielectric structure out of sapphire. Thermal aspects of the cavity were analyzed using finite element method (FEM) and an experiment designed to test different thermal joint materials. FEM was used to check the contraction of the cavity during cooling, the contact quality between the dielectric material and cavity wall, theheat flow rate through the assembly, and helps visualize the cavity’s reaction to different design changes. A simple comparison between thermal conductivity curves justified the usage of sapphire over alumina from a thermal property point of view. Thus, the method of experimentation had an additional dielectric material, sapphire, to test with and compare to alumina. The thermal test identified which material is best to use for a thermal joint but also simultaneously found the conductance of the joint and dielectric material as well as an estimate for what the temperature is inside the larger scale cavity since a thermometer cannot be placed inside the PBG when it cools. Initially, the cavity simulation was tested for structural deformation properties and stress distribution since it is made of copper which contracts heavily in comparison to most other materials. The reduction in volume of the cavity was less than anticipated which gave more room for a possible usage of a modified thermal joint for additional contact area and higher conductivity.

43 PARTICLE ACCELERATORS↗

Light Water Reactor Sustainability Program Advanced Seismic Soil Structure Modeling

Risk calculations should focus on providing best estimate results, and associated insights, for evaluation and decision-making. Specifically, seismic probabilistic risk assessments (SPRAs) are intended to provide best estimates of the various combinations of structural and equipment failures that can lead to a seismic induced core damage event. However, in some instances the current SPRA approach has large uncertainties, and potentially masks other important events (for instance, it was not the seismic motions that caused the Fukushima core melt events, but the tsunami ingress into the facility). SPRA’s are performed by convolving the seismic hazard (this is the estimate of all likely damaging earthquakes at the site of interest) with the seismic fragility (the conditional probability of failure of a structure, system, or component given the occurrence of earthquake ground motion). In this calculation, there are three main pieces to seismic risk quantification, 1) seismic hazard and nuclear power plants (NPPs) response to the hazard, 2) fragility or capacity of structures, systems and components (SSC), and 3) systems analysis. Two areas where NLSSI effects may be important in SPRA calculations are, 1) when calculating in-structure response at the area of interest, and 2) calculation of seismic fragilities (current fragility calculations assume a lognormal distribution for probability of failure of components). Some important effects when using NLSSI in the SPRA calculation process include, 1) gapping and sliding, 2) inclined seismic waves coupled with gapping and sliding of foundations atop soil, 3) inclined seismic waves coupled with gapping and sliding of deeply embedded structures, 4) soil dilatancy, 5) soil liquefaction, 6) surface waves, 7) buoyancy, 8) concrete cracking and 9) seismic isolation The focus of the research task presented here-in is on implementation of NLSSI into the SPRA calculation process when calculating in-structure response at the area of interest. The specific nonlinear soil behavior included in the NLSSI calculation presented in this report is gapping and sliding. Other NLSSI effects are not included in the calculation. The results presented in this report document initial model runs in the linear and nonlinear analysis process. Final comparisons between traditional and advanced SPRA will be presented in the September 30th deliverable.

54 ENVIRONMENTAL SCIENCES↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

HACC: Simulating sky surveys on state-of-the-art supercomputing architectures

Current and future surveys of large-scale cosmic structure are associated with a massive and complex datastream to study, characterize, and ultimately understand the physics behind the two major components of the ‘Dark Universe’, dark energy and dark matter. In addition, the surveys also probe primordial perturbations and carry out fundamental measurements, such as determining the sum of neutrino masses. Large-scale simulations of structure formation in the Universe play a critical role in the interpretation of the data and extraction of the physics of interest. Just as survey instruments continue to grow in size and complexity, so do the supercomputers that enable these simulations. In this paper, we report on HACC (Hardware/Hybrid Accelerated Cosmology Code), a recently developed and evolving cosmology N-body code framework, designed to run efficiently on diverse computing architectures and to scale to millions of cores and beyond. HACC can run on all current supercomputer architectures and supports a variety of programming models and algorithms. It has been demonstrated at scale on Cell- and GPU-accelerated systems, standard multi-core node clusters, and Blue Gene systems. HACC’s design allows for ease of portability, and at the same time, high levels of sustained performance on the fastest supercomputers available. We present a description of the design philosophy of HACC, the underlying algorithms and code structure, and outline implementation details for several specific architectures. We show selected accuracy and performance results from some of the largest high resolution cosmological simulations so far performed, including benchmarks evolving more than 3.6 trillion particles.

79 ASTRONOMY AND ASTROPHYSICS↗

Studies on structural and electrical properties of nanostructured GdMnO{sub 3}

In this communication, we report the results of the studies on structural and electrical properties of cost effective sol-gel grown nanostructured GdMnO{sub 3} manganites. Samples were sintered at different temperatures. Structural phases have been analyzed by performing the X-ray diffraction (XRD) measurements while crystallite size has been calculated using Scherer’s formula. AC conductivity (σ{sub AC}) has been measured for all the samples and fitted theoretically by using power law equation. Variation in σ{sub AC} has been understood on the basis of crystallite boundaries and oxygen ordering in all the samples. Nature of conduction in the samples has been discussed in the light of correlated barrier hopping mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗