Results for “Dynamical cluster approximation” — Page 9

NASA NTRS1989

New computing systems and their impact on structural analysis and design

A review is given of the recent advances in computer technology that are likely to impact structural analysis and design. The computational needs for future structures technology are described. The characteristics of new and projected computing systems are summarized. Advances in programming environments, numerical algorithms, and computational strategies for new computing systems are reviewed, and a novel partitioning strategy is outlined for maximizing the degree of parallelism. The strategy is designed for computers with a shared memory and a small number of powerful processors (or a small number of clusters of medium-range processors). It is based on approximating the response of the structure by a combination of symmetric and antisymmetric response vectors, each obtained using a fraction of the degrees of freedom of the original finite element model. The strategy was implemented on the CRAY X-MP/4 and the Alliant FX/8 computers. For nonlinear dynamic problems on the CRAY X-MP with four CPUs, it resulted in an order of magnitude reduction in total analysis time, compared with the direct analysis on a single-CPU CRAY X-MP machine.

Noor, Ahmed K.↗

DOE OSTI2022

Insights into the origin of halo mass profiles from machine learning

ABSTRACT The mass distribution of dark matter haloes is the result of the hierarchical growth of initial density perturbations through mass accretion and mergers. We use an interpretable machine-learning framework to provide physical insights into the origin of the spherically-averaged mass profile of dark matter haloes. We train a gradient-boosted-trees algorithm to predict the final mass profiles of cluster-sized haloes, and measure the importance of the different inputs provided to the algorithm. We find two primary scales in the initial conditions (ICs) that impact the final mass profile: the density at approximately the scale of the haloes’ Lagrangian patch RL ($R\sim 0.7\, R_L$) and that in the large-scale environment (R ∼ 1.7 RL). The model also identifies three primary time-scales in the halo assembly history that affect the final profile: (i) the formation time of the virialized, collapsed material inside the halo, (ii) the dynamical time, which captures the dynamically unrelaxed, infalling component of the halo over its first orbit, (iii) a third, most recent time-scale, which captures the impact on the outer profile of recent massive merger events. While the inner profile retains memory of the ICs, this information alone is insufficient to yield accurate predictions for the outer profile. As we add information about the haloes’ mass accretion history, we find a significant improvement in the predicted profiles at all radii. Our machine-learning framework provides novel insights into the role of the ICs and the mass assembly history in determining the final mass profile of cluster-sized haloes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Qubit-Efficient Quantum Chemistry with the ADAPT Variational Quantum Eigensolver and Double Unitary Downfolding

Here, in this work, we combine the recently developed double unitary coupled cluster (DUCC) theory with the adaptive, problem-tailored variational quantum eigensolver (ADAPT-VQE) to explore the accuracy of unitary downfolded Hamiltonians for quantum simulation of chemistry. We benchmark the ability of DUCC effective Hamiltonians to recover dynamical correlation energy outside of an active space. We consider the effects of strong correlation, commutator truncation, higher-body terms, and approximate external amplitudes on the accuracy of these effective Hamiltonians. When combining these DUCC Hamiltonians with ADAPT-VQE, we observe similar convergence of the ground state as compared with bare active space Hamiltonians, demonstrating that DUCC Hamiltonians provide increased accuracy without increasing the load on the quantum processor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2000

Glass/Jamming Transition in Colloidal Aggregation

We have studied colloidal aggregation in a model colloid plus polymer system with short-range attractive interactions. By varying the colloid concentration and the strength of the attraction, we explored regions where the equilibrium phase is expected to consist of colloidal crystallites in coexistance with colloidal gas (i.e. monomers). This occurs for moderate values of the potential depth, U approximately equal to 2-5 kT. Crystallization was not always observed. Rather, over an extended sub-region two new metastable phases appear, one fluid-like and one solid-like. These were examined in detail with light scattering and microscopy techniques. Both phases consist of a near uniform distribution of small irregular shaped clusters of colloidal particles. The dynamical and structural characteristics of the ergodic-nonergodic transition between the two phases share much in common with the colloidal hard sphere glass transition.

Segre, Philip N.↗

DOE OSTI2018

Lyapunov-Based Anomaly Detection in Highly-Clustered Networks

Network formation models explain the dynamics of the structure of connections using mechanisms that operate under different principles for establishing and removing edges. The Jackson–Rogers model is a generic framework that applies the principle of triadic closure to networks that grow by the addition of new nodes and new edges over time. Past work describes the limit distribution of the in-degree of the nodes based on a continuous-time approximation. Here, we introduce a discrete-time approach of the dynamics of the in- and out-degree distributions of a variation of the model. Furthermore, we characterize the limit distributions and the expected value of the average degree as equilibria, and prove that the equilibria are asymptotically stable. Finally, we use the stability properties of the model to propose a detection criterion for anomalies in the edge formation process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Investigating uncertainties in human adaptation and their impacts on water scarcity in the Colorado river Basin, United States

The Colorado River Basin (CRB) supports the water supply for seven states and forty million people in the Western United States (US) and has been suffering an extensive drought for more than two decades. As climate change continues to reshape water resources distribution in the CRB, its impact can differ in intensity and location, resulting in variations in human adaptation behaviors. The feedback from human systems in response to the environmental changes and the associated uncertainty is critical to water resources management, especially for water-stressed basins. This paper investigates how human adaptation affects water scarcity uncertainty in the CRB and highlights the uncertainties in human behavior modeling. Our focus is on agricultural water consumption, as approximately 80% of the water consumption in the CRB is used in agriculture. We adopted a coupled agent-based and water resources modeling approach for exploring human-water system dynamics, in which an agent is a human behavior model that simulates a farmer’s water consumption decisions. We examined uncertainties at the system, agent, and parameter levels through uncertainty, clustering, and sensitivity analyses. The uncertainty analysis results suggest that the CRB water system may experience 13 to 30 years of water shortage during the 2019–2060 simulation period, depending on the paths of farmers’ adaptation. The clustering analysis identified three decision-making classes: bold, prudent, and forward-looking, and quantified the probabilities of an agent belonging to each class. The sensitivity analysis results indicated agents whose decision-making models require further investigation and the parameters with the higher uncertainty reduction potentials. Here, by conducting numerical experiments with the coupled model, this paper presents quantitative and qualitative information about farmers’ adaptation, water scarcity uncertainties, and future research directions for improving human behavior modeling.

Agent-based modeling↗

DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗

NASA NTRS1996

Dynamic Coupling of Quasi-Electrostatic Thundercloud Fields to the Mesosphere and Lower Ionosphere: Sprites and Jets

Red Sprites and Blue Jets are two different types of recently discovered optical flashes ob- served above large thunderstorm systems. Sprites are luminous glows occurring at altitudes typically ranging from approximately 50 to 90 km. In video they exhibit a red color at their top which gradually changes to blue at lower altitudes. Sprites may occur singly or in clusters of two or more. The lateral extent of "unit" sprites is typically 5-10 km and they endure for several milliseconds. Jets are upward moving (approximately 100 km/s) highly collimated beams of luminosity, emanating from the tops of thunderclouds, extending up to approximately 50 km altitude and exhibiting a primarily blue color. We propose that sprites result from large electric field transients capable of causing electron heating, breakdown ionization and excitation of optical emissions at mesospheric altitudes following the removal of thundercloud charge by a cloud-to-ground discharge. Depending on the history of charge accumulation and removal, and the distribution of ambient atmospheric conductivity, the breakdown region may have the shape of vertically oriented ionization column(s). Results of a two-dimensional and self consistent quasi-electrostatic (QE) model indicate that most of the observed features of sprites can be explained in terms of the formation and self-driven propagation of streamer type channels of breakdown ionization. Comparison of the optical emission intensities of the 1st and 2nd positive bands of N2, Meinel and 1st negative bands of N2(+) and the 1st negative band of O2(+) demonstrates that the 1st positive band of N2 is the dominant optical emission in the altitude range approximately 50-90 km, which accounts for the observed red color of sprites. Optical emissions of the 1st and 2nd positive bands of N2 occur in carrot-like vertical structures with typical transverse dimension approximately 5-10 km which can span an altitude range from approximately 80 km to well below approximately 50 km. The appearance of optical emissions associated with sprites can be delayed in time (approximately 1-20 ms) with respect to the causative cloud to ground discharge. Theoretical model results are found to be in good agreement with recent video, photometric and spectral measurements of sprites.

Pasko, Victor Petrovich↗

DOE OSTI2024

Evaluation of dynamical downscaling in a fully coupled regional earth system model

A set of decadal simulations has been completed and evaluated for gains using the Regional Arctic System Model (RASM) to dynamically downscale data from a global Earth system model and two atmospheric reanalyses. RASM is a fully coupled atmosphere–land–ocean–sea ice regional Earth system model. Nudging to the forcing data is applied to approximately the top half of the atmospheric domain. RASM simulations were also completed with a modification to the atmospheric physics for evaluating changes to the modeling system. The results show that for the top half of the atmosphere, the RASM simulations follow closely to that of the forcing data, regardless of the forcing data. The results for the lower half of the atmosphere, as well as the surface, show a clustering of atmospheric state and surface fluxes based on the modeling system. At all levels of the atmosphere the imprint of the weather from the forcing data is present as indicated in the pattern of the annual means. Biases, in comparison to reanalyses, are evident in the Earth system model forced simulations for the top half of the atmosphere but are not present in the lower atmosphere. This suggests that bias correction is not needed for fully coupled dynamical downscaling simulations. While the RASM simulations tended to go to the same mean state for the lower atmosphere, there are a differences in the variability and changes of weather patterns across the ensemble of simulations. These differences in the weather result in variances in the sea ice and oceanic states.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2008

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively. The Mach Number distribution shows small yet distinct differences between the two cases such as locations of shocks/shock reflections and a slightly different impingement point on the wall of the diffuser from the expansion at the exit of the nozzle. Similarly the temperature distribution indicates different flow recirculation patterns in the test cell. Both cases capture all the essential flow phenomena such as the shock-boundary layer interaction, plume expansion, expansion of the first stage ejectors, mixing between the engine plume and the first stage ejector flow and pressurization due to the first stage ejectors. The final paper will discuss thermal loads on the walls of the diffuser and cooling mechanisms investigated.

Allgood, Daniel C.↗

NASA NTRS2010

Computational Analyses in Support of Sub-scale Diffuser Testing for the A-3 Facility: Steady Predictions - Part 1

Simulation technology can play an important role in rocket engine test facility design and development by assessing risks, providing analysis of dynamic pressure and thermal loads, identifying failure modes and predicting anomalous behavior of critical systems. Advanced numerical tools assume greater significance in supporting testing and design of high altitude testing facilities and plume induced testing environments of high thrust engines because of the greater inter-dependence and synergy in the functioning of the different sub-systems. This is especially true for facilities such as the proposed A-3 facility at NASA SSC because of a challenging operating envelope linked to variable throttle conditions at relatively low chamber pressures. Facility designs in this case will require a complex network of diffuser ducts, steam ejector trains, fast operating valves, cooling water systems and flow diverters that need to be characterized for steady state performance. In this paper, we will demonstrate with the use of CFD analyses s advanced capability to evaluate supersonic diffuser and steam ejector performance in a sub-scale A-3 facility at NASA Stennis Space Center (SSC) where extensive testing was performed. Furthermore, the focus in this paper relates to modeling of critical sub-systems and components used in facilities such as the A-3 facility. The work here will address deficiencies in empirical models and current CFD analyses that are used for design of supersonic diffusers/turning vanes/ejectors as well as analyses for confined plumes and venting processes. The primary areas that will be addressed are: (1) supersonic diffuser performance including analyses of thermal loads (2) accurate shock capturing in the diffuser duct; (3) effect of turning duct on the performance of the facility (4) prediction of mass flow rates and performance classification for steam ejectors (5) comparisons with test data from sub-scale diffuser testing and assessment of confidence levels in CFD based flowpath modeling of the facility. The analyses tools used here expand on the multi-element unstructured CFD which has been tailored and validated for impingement dynamics of dry plumes, complex valve/feed systems, and high pressure propellant delivery systems used in engine and component test stands at NASA SSC. The analyses performed in the evaluation of the sub-scale diffuser facility explored several important factors that influence modeling and understanding of facility operation such as (a) importance of modeling the facility with Real Gas approximation, (b) approximating the cluster of steam ejector nozzles as a single annular nozzle, (c) existence of mixed subsonic/supersonic flow downstream of the turning duct, and (d) inadequacy of two-equation turbulence models in predicting the correct pressurization in the turning duct and expansion of the second stage steam ejectors. The procedure used for modeling the facility was as follows: (i) The engine, test cell and first stage ejectors were simulated with an axisymmetric approximation (ii) the turning duct, second stage ejectors and the piping downstream of the second stage ejectors were analyzed with a three-dimensional simulation utilizing a half-plane symmetry approximation. The solution i.e. primitive variables such as pressure, velocity components, temperature and turbulence quantities were passed from the first computational domain and specified as a supersonic boundary condition for the second simulation. (iii) The third domain comprised of the exit diffuser and the region in the vicinity of the facility (primary included to get the correct shock structure at the exit of the facility and entrainment characteristics). The first set of simulations comprising the engine, test cell and first stage ejectors was carried out both as a turbulent real gas calculation as well as a turbulent perfect gas calculation. A comparison for the two cases (Real Turbulent and Perfect gas turbulent) of the Ma Number distribution and temperature distributions are shown in Figures 1 and 2 respectively.

Allgood, Daniel C.↗

NASA NTRS1993

Gas dynamic simulations of galaxy formation

Results are presented from a simulation modeling the formation of a group of galaxies in a 'standard' cold, dark matter universe with delta = 1, h sub 0 = 50 km/(s(Mpc)), baryon fraction omega sub b = 0.1 and spectrum normalization sigma sub 8 = 0.6 (bias parameter b = 1.7). Initial conditions are generated within a periodic box with comoving length 16 Mpc in a manner constrained to produce a small cluster of total mass approximately 10 exp 14 solar mass. Two sets of 643 particles are used to model the dark matter and baryon fluids. Each gas particle represents 1.08 x 10 exp -8 solar mass, implying an L* galaxy is resolved by approximately 1000 particles. The system is evolved self-consistently in three dimensions using the combined N-body/hydrodynamic scheme P3MSPH up to a final redshift z = 1. Evolving to the present is prohibited by the fact that the mean density in the simulated volume is above critical and the entire volume would be going nonlinear beyond this point, We are currently analyzing another run with somewhat poorer mass resolution which was evolved to the present.

Evrard, August E.↗

DOE OSTI2023

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

Nuclear spin/parity dependent spectroscopy and predissociation dynamics in v OH = 2 ← 0 overtone excited Ne–H 2 O clusters: Theory and experiment

Vibrationally state selective overtone spectroscopy and state- and nuclear spin-dependent predissociation dynamics of weakly bound ortho- and para-Ne–H 2 O complexes (D 0(ortho) = 34.66 cm -1 and D 0(para) = 31.67 cm -1 ) are reported, based on near-infrared excitation of van der Waals cluster bands correlating with v OH = 2 ← 0 overtone transitions (|02 - &rang; and |02 + &rang;) out of the ortho (101) and para (000) internal rotor states of the H 2 O moiety. Quantum theoretical calculations for nuclear motion on a high level potential energy surface [CCSD(T)/VnZf12 (n = 3, 4)], corrected for basis set superposition error and extrapolated to the complete basis set (CBS) limit, are employed to successfully predict and assign Π–Σ, Σ–Σ, and Σ–Π infrared bands in the spectra, where Σ or Π represent approximate projections of the body-fixed H 2 O angular momentum along the Ne–H 2 O internuclear axis. IR-UV pump-probe experimental capabilities permit real-time measurements of the vibrational predissociation dynamics, which indicate facile intramolecular vibrational energy transfer from the H 2 O v OH = 2 overtone vibrations into the VdWs (van der Waals) dissociation coordinate on the τprediss = 15-25 ns time scale. Whereas all predicted strong transitions in the ortho-Ne–H 2 O complexes are readily detected and assigned, vibrationally mediated photolysis spectra for the corresponding para-Ne–H 2 O bands are surprisingly absent despite ab initio predictions of Q-branch intensities with S/N > 20-40. Such behavior signals the presence of highly selective nuclear spin ortho-para predissociation dynamics in the upper state, for which we offer a simple mechanism based on Ne-atom mediated intramolecular vibrational relaxation in the H 2 O subunit (i.e., |02 ± &rang; → {|01 ± &rang; v 2 = 2}), which is confirmed by the ab initio energy level predictions and the nascent OH rotational (N), spin orbit (Π 1/2,3/2 ), and lambda doublet product distributions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2017

A Chandra X-Ray Census of the Interacting Binaries in Old Open Clusters—Collinder 261

We present the first X-ray study of Collinder 261 (Cr 261), which at an age of 7 Gyr is one of the oldest open clusters known in the Galaxy. Our observation with the Chandra X-Ray Observatory is aimed at uncovering the close interacting binaries in Cr 261, and reaches a limiting X-ray luminosity of L{sub X}≈4×10{sup 29} erg s{sup −1} (0.3–7 keV) for stars in the cluster. We detect 107 sources within the cluster half-mass radius r {sub h}, and we estimate that among the sources with L{sub X}≳10{sup 30} erg s{sup −1}, ∼26 are associated with the cluster. We identify a mix of active binaries and candidate active binaries, candidate cataclysmic variables, and stars that have “straggled” from the main locus of Cr 261 in the color–magnitude diagram. Based on a deep optical source catalog of the field, we estimate that Cr 261 has an approximate mass of 6500 M {sub ⊙}, roughly the same as the old open cluster NGC 6791. The X-ray emissivity of Cr 261 is similar to that of other old open clusters, supporting the trend that they are more luminous in X-rays per unit mass than old populations of higher (globular clusters) and lower (the local neighborhood) stellar density. This implies that the dynamical destruction of binaries in the densest environments is not solely responsible for the observed differences in X-ray emissivity.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

Adiabatic invariants in stellar dynamics, 3: Application to globular cluster evolution

The previous two companion papers demonstrate that slowly varying perturbations may not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with R(sub g) less than or approximately 8 kpc: a concentration c = 1.5 cluster with R(sub g) kpc loses up to 95% of its initial mass in 15 Gyr. Only those with concentration c greater than or approximately 1.3 survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will decrease this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, x, which may be considerably larger than their initial values.

Weinberg, Martin D.↗

DOE OSTI2024

VpROM: a novel variational autoencoder-boosted reduced order model for the treatment of parametric dependencies in nonlinear systems

Reduced Order Models (ROMs) are of considerable importance in many areas of engineering in which computational time presents difficulties. Established approaches employ projection-based reduction, such as Proper Orthogonal Decomposition. The limitation of the linear nature of such operators is typically tackled via a library of local reduction subspaces, which requires the assembly of numerous local ROMs to address parametric dependencies. Our work attempts to define a more generalisable mapping between parametric inputs and reduced bases for the purpose of generative modeling. We propose the use of Variational Autoencoders (VAEs) in place of the typically utilised clustering or interpolation operations, for inferring the fundamental vectors, termed as modes, which approximate the manifold of the model response for any and each parametric input state. The derived ROM still relies on projection bases, built on the basis of full-order model simulations, thus retaining the imprinted physical connotation. However, it additionally exploits a matrix of coefficients that relates each local sample response and dynamics to the global phenomena across the parametric input domain. The VAE scheme is utilised for approximating these coefficients for any input state. This coupling leads to a high-precision low-order representation, which is particularly suited for problems where model dependencies or excitation traits cause the dynamic behavior to span multiple response regimes. Moreover, the probabilistic treatment of the VAE representation allows for uncertainty quantification on the reduction bases, which may then be propagated to the ROM response. The performance of the proposed approach is validated on an open-source simulation benchmark featuring hysteresis and multi-parametric dependencies, and on a large-scale wind turbine tower characterised by nonlinear material behavior and model uncertainty.

Conditional VAEs↗

DOE OSTI2019

Temperature dependence of structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 : an ab initio study

Bismuth telluride (Bi 2 Te 3 ) has garnered significant interest in thermoelectric applications and three-dimensional topological insulators due to its unique electronic, transport, and thermal properties. Bi 2 Te 3 and Sb 2 Te 3 chalcogenide compounds have the same crystal structure. While Sb 2 Te 3 has been shown to be a prototypical phase change memory (PCM) compound along the pseudobinary tie-line of Ge-Sb-Te alloys, whether Bi 2 Te 3 can also exhibit PCM functionality is still not well established. In this work, a systematic study on the structural, dynamical, and electronic properties of amorphous Bi 2 Te 3 during the quenching process has been performed by using ab initio molecular dynamics simulations. Pair correlation function, coordination number, bond-angle distribution functions, and a novel atomistic cluster alignment method are used to explore the structural characteristics of Bi 2 Te 3 as a function of temperature. Our study shows that there are many distorted octahedral clusters in amorphous Bi 2 Te 3 . In comparison with the local structures in Sb 2 Te 3 , we found that the degree of distortion of the octahedrons in the Bi 2 Te 3 system is smaller than that in Sb 2 Te 3 system. Moreover, the changes in the dynamical properties of Bi 2 Te 3 from liquid to glassy state are also explored. The approximate range of liquid-to-glass transition temperature is determined to be between 673 and 723 K. The electronic properties of Bi2Te 3 and Sb 2 Te 3 are also analysed by density-of-states and Bader charge calculations, both of them in glass state are semiconductors. Our studies provide useful insights into the local structure and dynamical properties of Bi 2 Te 3 at the atomistic level during the fast cooling process, and suggest that the compound can be a candidate for PCM materials.

36 MATERIALS SCIENCE↗