Search NASA⌕ Search

SEARCH · Search NASA

Results for “Manifold sampling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Improving Limit Surface Search Algorithms in RAVEN Using Acceleration Schemes: Level II Milestone

The RAVEN code is becoming a comprehensive tool to perform Probabilistic Risk Assessment (PRA); Uncertainty Quantification (UQ) and Propagation; and Verification and Validation (V&V). The RAVEN code is being developed to support the Risk-Informed Safety Margin Characterization (RISMC) pathway by developing an advanced set of methodologies and algorithms for use in advanced risk analysis. The RISMC approach uses system simulator codes applied to stochastic analysis tools. The fundamental idea behind this coupling approach to perturb (by employing sampling strategies) timing and sequencing of events, internal parameters of the system codes (i.e., uncertain parameters of the physics model) and initial conditions to estimate values ranges and associated probabilities of figures of merit of interest for engineering and safety (e.g. core damage probability, etc.). This approach applied to complex systems such as nuclear power plants requires performing a series of computationally expensive simulation runs. The large computational burden is caused by the large set of (uncertain) parameters characterizing those systems. Consequently, exploring the uncertain/parametric domain, with a good level of confidence, is generally not affordable, considering the limited computational resources that are currently available. In addition, the recent tendency to develop newer tools, characterized by higher accuracy and larger computational resources (if compared with the presently used legacy codes, that have been developed decades ago), has made this issue even more compelling. In order to overcome to these limitations, the strategy for the exploration of the uncertain/parametric space needs to use at best the computational resources focusing the computational effort in those regions of the uncertain/parametric space that are “interesting” (e.g., risk-significant regions of the input space) with respect the targeted Figures Of Merit (FOM): for example, the failure of the system, subject of the analysis. These methodologies are named, in the RAVEN environment, adaptive sampling strategies. These methodologies infer system responses from surrogate models constructed from already existing samples (produced using high fidelity simulations) and suggest the most relevant location (coordinate in the input space) of the next sampling point to be explored in the uncertain/parametric domain. When using those methodologies, it is possible to understand features of the system response with a small number of carefully selected samples. This report focuses on the development and improvement of the limit surface search. The limit surface is an important concept in system reliability analysis. Without going into the details, which will be covered later in the report, the limit surface could be briefly described as an hyper-surface in the system uncertainty/parametric space separating the regions leading to a prescribed system outcome. For example, if the uncertainty/parametric space is the one generated by the reactor power level and the duration of the batteries, the system is a nuclear power plant and the system outcome discriminating variable is the clad failure in a station blackout scenario, then the limit surface separates the combinations of reactor power level and battery duration that lead to clad failure from the ones that do not.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Fractionation of Oxygen Isotopes in Uranium Oxides during Peroxide Precipitation and Dry Air Calcination

The interaction between uranium and oxygen-containing substrates is nearly ubiquitous within the nuclear fuel cycle. Given the well-known and predictable oxygen stable isotope compositions of atmospheric oxygen and meteoric waters around the world, the use of these isotopes as a potential geolocation or processing signature of uranium compounds has been of interest within the nuclear safeguards community. This study focuses on measuring the oxygen-stable isotope composition of synthetically produced uranium oxides to determine the mechanism and extent of fractionation during materials processing relevant to the nuclear fuel cycle. Here, metastudtite [UO 2 O 2 (H 2 O) 2 ] samples were first produced using water and hydrogen peroxide of a known oxygen isotope composition. This starting material was calcined in dry air at temperatures ranging from 300 °C to 1000 °C for 20 h producing oxides ranging in composition from UO 3 +x (0 ≤ x ≤ 0.5) to U 3 O 8 . In addition, calcinations between 500 and 800 °C were performed at different time intervals to evaluate the rate of isotopic equilibration. The δ 18 O values of water bound to metastudtite were measured by thermogravimetric analysis coupled to isotope ratio infrared spectroscopy (TGA–IRIS). A redesigned fluorination manifold was constructed at the University of Utah and was used to extract the bulk oxygen from precipitated and calcined products using bromine pentafluoride (BrF 5 ). Mineralization water on metastudtite (measured by TGA–IRIS) was closely associated with the aqueous environment of precipitation with a fractionation of +3.5 ± 0.6‰. In contrast, bulk oxygen from metastudtite (measured by fluorination and gas chromatography combined with isotope ratio mass spectrometry) was found to have an inconsistent fractionation with precipitation water. Samples calcined in dry air exhibited a wide range of oxygen isotope compositions which increased from δ 18 O = +1.4 ± 0.9‰ at 300 °C to +17.3 ± 0.9‰ at 1000 °C. This variation in δ 18 O values was a result of oxygen exchange with atmospheric O 2 and was found to be rapid at calcination temperatures greater than 600°C—with equilibration occurring in less than 2 h. These results provide insights into the incorporation and exchange of oxygen isotopes during precipitation and calcination processes comparable to those employed within the nuclear fuel cycle and establish a foundation for future investigations to build upon.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1993

Technical Update: Johnson Space Center system using a solid electrolytic cell in a remote location to measure oxygen fugacities in CO/CO2 controlled-atmosphere furnaces

Details are given for the design and application of a (one atmosphere) redox-control system. This system differs from that given in NASA Technical Memorandum 58234 in that it uses a single solid-electrolytic cell in a remote location to measure the oxygen fugacities of multiple CO/CO2 controlled-atmosphere furnaces. This remote measurement extends the range of sample-furnace conditions that can be measured using a solid-electrolytic cell, and cuts costs by extending the life of the sensors and by minimizing the number of sensors in use. The system consists of a reference furnace and an exhaust-gas manifold. The reference furnace is designed according to the redox control system of NASA Technical Memorandum 58234, and any number of CO/CO2 controlled-atmosphere furnaces can be attached to the exhaust-gas manifold. Using the manifold, the exhaust gas from individual CO/CO2 controlled atmosphere furnaces can be diverted through the reference furnace, where a solid-electrolyte cell is used to read the ambient oxygen fugacity. The oxygen fugacity measured in the reference furnace can then be used to calculate the oxygen fugacity in the individual CO/CO2 controlled-atmosphere furnace. A BASIC computer program was developed to expedite this calculation.

Jurewicz, A. J. G.↗
DOE OSTI2020

Machine Learning Calabi–Yau Metrics

We apply machine learning to the problem of finding numerical Calabi–Yau metrics. Building on Donaldson's algorithm for calculating balanced metrics on Kähler manifolds, we combine conventional curve fitting and machine-learning techniques to numerically approximate Ricci-flat metrics. We show that machine learning is able to predict the Calabi–Yau metric and quantities associated with it, such as its determinant, having seen only a small sample of training data. Using this in conjunction with a straightforward curve fitting routine, we demonstrate that it is possible to find highly accurate numerical metrics much more quickly than by using Donaldson's algorithm alone, with our new machine-learning algorithm decreasing the time required by between one and two orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Performance of a CO 2 sorbent for indoor air cleaning applications: Effects of environmental conditions, sorbent aging, and adsorption of co-occurring formaldehyde

Indoor air cleaning systems that incorporate CO 2 sorbent materials enable HVAC load shifting and efficiency improvements. This study developed a bench-scale experimental system to evaluate the performance of a sorbent under controlled operation conditions. A thermostatic holder containing 3.15 g sorbent was connected to a manifold that delivered CO 2 -enriched air at a known temperature and relative humidity (RH). The air stream was also enriched with 0.8-2.1 ppm formaldehyde. The CO 2 concentration was monitored in real-time upstream and downstream of the sorbent, and integrated formaldehyde samples were collected at different times using DNPH-coated silica cartridges. Sorbent regeneration was carried out by circulating clean air in countercurrent. Almost 200 loading/regeneration cycles were performed in the span of 17 months, from which 104 were carried out at reference test conditions defined by loading with air at 25°C, 38% RH, and 1000 ppm CO 2 , and regenerating with air at 80°C, 3% RH and 400 ppm CO 2 . The working capacity decreased slightly from 43-44 mg CO 2 per g sorbent to 39-40 mg per g over the 17 months. The capacity increased with lower loading temperature (in the range 15-35°C) and higher regeneration temperature, between 40 and 80°C. The CO 2 capacity was not sensitive to the moisture content in the range 6-9 g/m 3 , and decreased slightly when dry air was used. Loading isothermal breakthrough curves were fitted to three simple adsorption models, verifying that pseudo-first-order kinetics appropriately describes the adsorption process. The model predicted that equilibrium capacities decreased with increasing temperature from 15 to 35°C, while adsorption rate constants slightly increased. The formaldehyde adsorption efficiency was 80%-99% in different cycles, corresponding to an average capacity of 86 ± 36 µg/g. Formaldehyde was not quantitatively released during regeneration, but its accumulation on the sorbent did not affect CO 2 adsorption.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS2010

Cryogenic Absorption Cells Operating Inside a Bruker IFS-125HR: First Results for 13CH4 at 7 Micrometers

New absorption cells designed specifically to achieve stable temperatures down to 66 K inside the sample compartment of an evacuated Bruker IFS-125HR Fourier transform spectrometer (FTS) were developed at Connecticut College and tested at the Jet Propulsion Laboratory (JPL). The temperature stabilized cryogenic cells with path lengths of 24.29 and 20.38 cm were constructed of oxygen free high conductivity (OFHC) copper and fitted with wedged ZnSe windows using vacuum tight indium seals. In operation, the temperature-controlled cooling by a closed-cycle helium refrigerator achieved stability of 0.01 K. The unwanted absorption features arising from cryodeposits on the cell windows at low temperatures were eliminated by building an internal vacuum shroud box around the cell which significantly minimized the growth of cryodeposits. The effects of vibrations from the closed-cycle helium refrigerator on the FTS spectra were characterized. Using this set up, several high-resolution spectra of methane isotopologues broadened with nitrogen were recorded in the 1200-1800 per centimeter spectral region at various sample temperatures between 79.5 and 296 K. Such data are needed to characterize the temperature dependence of spectral line shapes at low temperatures for remote sensing of outer planets and their moons. Initial analysis of a limited number of spectra in the region of the R(2) manifold of the v4 fundamental band of 13CH4 indicated that an empirical power law used for the temperature dependence of the N2-broadened line widths would fail to fit the observed data in the entire temperature range from 80 to 296 K; instead, it follows a temperature-dependence similar to that reported by Mondelain et al. [17,18]. The initial test was very successful proving that a high precision Fourier transform spectrometer with a completely evacuated optical path can be configured for spectroscopic studies at low temperatures relevant to the planetary atmospheres.

Sung, K.↗
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 OSTI2023

A review on recent machine learning applications for imaging mass spectrometry studies

Imaging mass spectrometry (IMS) is a powerful analytical technique widely used in biology, chemistry, and materials science fields that continue to expand. IMS provides a qualitative compositional analysis and spatial mapping with high chemical specificity. The spatial mapping information can be 2D or 3D depending on the analysis technique employed. Due to the combination of complex mass spectra coupled with spatial information, large high-dimensional datasets (hyperspectral) are often produced. Therefore, the use of automated computational methods for an exploratory analysis is highly beneficial. The fast-paced development of artificial intelligence (AI) and machine learning (ML) tools has received significant attention in recent years. These tools, in principle, can enable the unification of data collection and analysis into a single pipeline to make sampling and analysis decisions on the go. There are various ML approaches that have been applied to IMS data over the last decade. Here, in this review, we discuss recent examples of the common unsupervised (principal component analysis, non-negative matrix factorization, k-means clustering, uniform manifold approximation and projection), supervised (random forest, logistic regression, XGboost, support vector machine), and other methods applied to various IMS datasets in the past five years. The information from this review will be useful for specialists from both IMS and ML fields since it summarizes current and representative studies of computational ML-based exploratory methods for IMS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

Sign problem and Monte Carlo calculations beyond Lefschetz thimbles

We point out that Monte Carlo simulations of theories with severe sign problems can be profitably performed over manifolds in complex space different from the one with fixed imaginary part of the action (“Lefschetz thimble”). We describe a family of such manifolds that interpolate between the tangent space at one critical point (where the sign problem is milder compared to the real plane but in some cases still severe) and the union of relevant thimbles (where the sign problem is mild but a multimodal distribution function complicates the Monte Carlo sampling). As a result, we exemplify this approach using a simple 0+1 dimensional fermion model previously used on sign problem studies and show that it can solve the model for some parameter values where a solution using Lefschetz thimbles was elusive.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Improved surrogates in inertial confinement fusion with manifold and cycle consistencies

Neural networks have become the method of choice in surrogate modeling because of their ability to characterize arbitrary, high-dimensional functions in a data-driven fashion. This paper advocates for the training of surrogates that are 1) consistent with the physical manifold, resulting in physically meaningful predictions, and 2) cyclically consistent with a jointly trained inverse model; i.e., backmapping predictions through the inverse results in the original input parameters. We find that these two consistencies lead to surrogates that are superior in terms of predictive performance, are more resilient to sampling artifacts, and tend to be more data efficient. Using inertial confinement fusion (ICF) as a test-bed problem, we model a one-dimensional semianalytic numerical simulator and demonstrate the effectiveness of our approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Launch Complex 39B, SWMU 009 Air Sparging System Expansion Implementation Work Plan Kennedy Space Center, Florida

This Implementation Work Plan (IWP) for expansion of the Air Sparging (AS) system at Launch Complex 39B (LC39B), presents the design and coordination specifics to expand the existing AS Interim Measure (IM) at LC39B located at Kennedy Space Center (KSC), Florida. The existing AS system at LC39B was installed in 2017 and included 279 AS wells. LC39B has been designated as Solid Waste Management Unit (SWMU) 009 under the KSC Resource Conservation and Recovery Act (RCRA) Corrective Action Program. Expansion of the AS IM will include installation of an additional five AS wells to depths of approximately 13 and 17 feet below land surface (bls) to remediate contaminated groundwater identified northwest of the launch pad, near the liquid oxygen (LOX) tank. The expansion will target concentrations of cis-1,2-dichloroethene (cDCE) and vinyl chloride (VC) greater than Florida Department of Environmental Protection (FDEP) Natural Attenuation Default Concentrations (NADCs) that were identified during direct push technology (DPT) sampling events between 2019-2022 in shallow intervals between 8 and 16 feet bls. The objective for this expansion is to remediate contaminated groundwater within the treatment zone to below NADCs to support transition to monitored natural attenuation (MNA) through a long-term monitoring (LTM) program. This expansion will treat approximately 0.07-acres of residual groundwater contamination. The IWP design for this expansion will utilize the existing AS system trailer that is already onsite and operating. The five new AS expansion wells will be integrated into the existing Zone Z4 distribution manifold to deliver compressed air to the new AS wells. To deliver airflow, three existing lateral lines will be split in the field at a lateral split junction point (using 2-way or 3-way splits). Each split lateral will supply air to AS wells at the same depth interval to minimize preferential pathways. The trenching layout for the new AS wells was developed to minimize ground disturbance due to the amount of utilities in the area. One new monitoring well will also be installed and added to LC39B’s performance monitoring program to evaluate the effectiveness of the AS system expansion. This well will be sampled on a semi-annual basis for volatile organic compounds, with results evaluated to determine if system adjustments may be required. The main components of this IWP were presented to the KSC Remediation Team (KSCRT) at the October 2022 Meeting (Meeting 2210-M05). DPT delineation data from the LOX area were also presented, which included some inferred boundaries due to physical restrictions of the LOX tank and associated concrete pad and infrastructure. Planned monitoring well sampling in this area will help determine if back-diffusion may be occurring from beneath the pad. Additionally, based on potential petroleum odors and staining identified at two DPT locations during the January 2022 DPT event in the 10-16 feet bls range, a KSCRT action item was assigned during the October 2022 meeting to collect soil samples for petroleum compounds if similar odors or staining are identified during IM implementation. Construction completion and performance monitoring will be incorporated into LC39B’s groundwater monitoring reports, which are currently submitted on an annual basis to FDEP to document operation and maintenance activities and results from the AS IM.

Deborah M. Wilson↗
DOE OSTI2022

Illuminating the Photophysics of Oxygen Atom Transfer in Rare Earth Metal-Organic Complexes

Industrial separations of rare earth metals are reliant on selectivity based on ionic radii, which are time, energy, and resource consuming.1 The unique and discontinuous manifold of energy states associated with the f-electron configurations of the rare earth ions2 offer an opportunity to overcome some of these limitations by incorporating light into the separation process. Here, we use transient absorption spectroscopy to understand the differences in mechanism of light-initiated oxygen atom transfer (OAT) in yttrium (Y) and dysprosium (Dy) rare earth complexes. These metal-organic complexes are functionalized with hexfluoroacetyl-acetonate (hfac) and 4-methylmorpholine-N-oxide (NMMO). After light absorption in the presence of triphenyl phosphine (TPP), the NMMO in the Y-NMMO complex is replaced by the TPPO OAT reaction product. Because the Dy f-electron manifold is accessible to the ligand excited states, while there are no accessible states associated with f-electrons for Y, there are significant differences in reactivity. Transient absorption spectroscopy (TAS) was used to monitor changes after 340 nm photoexcitation of Dy/Y-NMMO or the product complex Dy/Y-TPPO. TAS is a pump-probe technique where a white light probe pulse is delayed relative to a visible pump pulse, which photoexcites a sample. Differences between the ground and excited states are monitored over time. With TAS, we access time scales from 100 fs-400 us to determine the photophysics of the complexes of interest. Photoexcitation of Y-NMMO (Fig. 1A) creates a broad photoinduced absorption (PIA) that evolves into a sharp positive feature over 4.2 ps centered at 450 nm. A secondary feature centered at 550 nm grows in with a 200 ps time constant and remains through the 5 ns time window of our experiment. Control experiments for Y-TPPO reveal that the feature at 450 nm is due to internal dynamics of the hfac ligand, while the 550 nm shoulder is due to the presence of the NMMO. Comparing kinetics at 550 nm for Y-NMMO and Dy-NMMO (Fig. 1B) shows that while Y exhibits a distinct growth in the NMMO feature, the Dy decays almost completely within 5 ns. Control experiments with Dy-TPPO prove that the ligand excited states centered on the hfac ligand are deactivated via energy transfer to the Dy manifold within 200 ps. This deactivation leads to significant differences in OAT reactivity between Y and Dy.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2022

Fast Characterization of Inducible Regions of Atrial Fibrillation Models With Multi-Fidelity Gaussian Process Classification

Computational models of atrial fibrillation have successfully been used to predict optimal ablation sites. A critical step to assess the effect of an ablation pattern is to pace the model from different, potentially random, locations to determine whether arrhythmias can be induced in the atria. In this work, we propose to use multi-fidelity Gaussian process classification on Riemannian manifolds to efficiently determine the regions in the atria where arrhythmias are inducible. We build a probabilistic classifier that operates directly on the atrial surface. We take advantage of lower resolution models to explore the atrial surface and combine seamlessly with high-resolution models to identify regions of inducibility. We test our methodology in 9 different cases, with different levels of fibrosis and ablation treatments, totalling 1,800 high resolution and 900 low resolution simulations of atrial fibrillation. When trained with 40 samples, our multi-fidelity classifier that combines low and high resolution models, shows a balanced accuracy that is, on average, 5.7% higher than a nearest neighbor classifier. We hope that this new technique will allow faster and more precise clinical applications of computational models for atrial fibrillation. All data and code accompanying this manuscript will be made publicly available at: https://github.com/fsahli/AtrialMFclass.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Spin skyrmion gaps as signatures of strong-coupling insulators in magic-angle twisted bilayer graphene

Abstract The flat electronic bands in magic-angle twisted bilayer graphene (MATBG) host a variety of correlated insulating ground states, many of which are predicted to support charged excitations with topologically non-trivial spin and/or valley skyrmion textures. However, it has remained challenging to experimentally address their ground state order and excitations, both because some of the proposed states do not couple directly to experimental probes, and because they are highly sensitive to spatial inhomogeneities in real samples. Here, using a scanning single-electron transistor, we observe thermodynamic gaps at even integer moiré filling factors at low magnetic fields. We find evidence of a field-tuned crossover from charged spin skyrmions to bare particle-like excitations, suggesting that the underlying ground state belongs to the manifold of strong-coupling insulators. From the spatial dependence of these states and the chemical potential variation within the flat bands, we infer a link between the stability of the correlated ground states and local twist angle and strain. Our work advances the microscopic understanding of the correlated insulators in MATBG and their unconventional excitations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Doing the impossible: Why neural networks can be trained at all

As deep neural networks grow in size, from thousands to millions to billions of weights, the performance of those networks becomes limited by our ability to accurately train them. A common naive question arises: if we have a system with billions of degrees of freedom, don't we also need billions of samples to train it? Of course, the success of deep learning indicates that reliable models can be learned with reasonable amounts of data. Similar questions arise in protein folding, spin glasses and biological neural networks. With effectively infinite potential folding/spin/wiring configurations, how does the system find the precise arrangement that leads to useful and robust results? Simple sampling of the possible configurations until an optimal one is reached is not a viable option even if one waited for the age of the universe. On the contrary, there appears to be a mechanism in the above phenomena that forces them to achieve configurations that live on a low-dimensional manifold, avoiding the curse of dimensionality. In the current work we use the concept of mutual information between successive layers of a deep neural network to elucidate this mechanism and suggest possible ways of exploiting it to accelerate training.

deep learning, training, curse of dimensionality, ↗
DOE OSTI2022

Photometric redshifts from SDSS images with an interpretable deep capsule network

ABSTRACT Studies of cosmology, galaxy evolution, and astronomical transients with current and next-generation wide-field imaging surveys like the Rubin Observatory Legacy Survey of Space and Time are all critically dependent on estimates of photometric redshifts. Capsule networks are a new type of neural network architecture that is better suited for identifying morphological features of the input images than traditional convolutional neural networks. We use a deep capsule network trained on ugriz images, spectroscopic redshifts, and Galaxy Zoo spiral/elliptical classifications of ∼400 000 Sloan Digital Sky Survey galaxies to do photometric redshift estimation. We achieve a photometric redshift prediction accuracy and a fraction of catastrophic outliers that are comparable to or better than current methods for SDSS main galaxy sample-like data sets (r ≤ 17.8 and zspec ≤ 0.4) while requiring less data and fewer trainable parameters. Furthermore, the decision-making of our capsule network is much more easily interpretable as capsules act as a low-dimensional encoding of the image. When the capsules are projected on a two-dimensional manifold, they form a single redshift sequence with the fraction of spirals in a region exhibiting a gradient roughly perpendicular to the redshift sequence. We perturb encodings of real galaxy images in this low-dimensional space to create synthetic galaxy images that demonstrate the image properties (e.g. size, orientation, and surface brightness) encoded by each dimension. We also measure correlations between galaxy properties (e.g. magnitudes, colours, and stellar mass) and each capsule dimension. We publicly release our code, estimated redshifts, and additional catalogues at https://biprateep.github.io/encapZulate-1.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Studies of the Virgo Cluster. I - Photometry of 109 galaxies near the cluster center to serve as standards

Attention is given to the technical aspects of photometric measurements of 109 galaxies near the center of the Virgo Cluster, noting various types of radii and surface brightness for about 50 E and dE galaxies in the sample that range in absolute magnitude from -20 to -12. These data are combined with data from the literature for giant E and dwarf E galaxies in the Local Group to study the systematic properties of E galaxies over a range of one million luminosities. The radial intensity profiles derived are fitted to the manifold of King (1978) models to derive model-dependent central surface brightness, core radii, and cutoff radii.

Sandage, A.↗
DOE OSTI2024

Finite elements for Matérn-type random fields: Uncertainty in computational mechanics and design optimization

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on prevalent applications in computational mechanics and design optimization. We leverage Matérn-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models as well as the shape and topology of optimized designs. Here we describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications, emphasizing the potential to influence these domains. The flexibility and efficiency of SPDE-based GRF generation empowers us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate design features and topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model and quantify geometric uncertainties on reconstructed submanifolds, such as the interpolated surfaces of cerebral aneurysms provided by postprocessing CT scans. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements.

97 MATHEMATICS AND COMPUTING↗