Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model discovery”

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 793 records · Page 44

Automated and High-Throughput Phase Separation Control for Supramolecular Polymer Blends Enabled by Machine Learning

Supramolecular polymer blends (SPBs) offer tunable morphologies that dictate their macroscopic properties, yet their rational design is limited by the absence of predictive structure−morphology models. Here, we introduce a data-driven highthroughput workflow that integrates modular polymer synthesis, robotic formulation, automated morphology characterization, and machine learning (ML) for accelerated SPB discovery. Using a plug-and-play synthetic strategy, 33 hydrogen-bonding endfunctional homopolymers were prepared and orthogonally combined to generate 260 SPBs in 1 day. A fully automated atomic force microscopy (AFM) pipeline enabled systematic imaging, producing 2340 morphology data sets with minimal human intervention. Domain spacings were extracted through complementary imageprocessing methods and used to train ML models. A support vector regression (SVR) model accurately predicted target phase-separation sizes (50, 100, and 150 nm), which were experimentally validated. This work demonstrates the power of coupling high-throughput experimentation with ML to accelerate morphology discovery and provides one of the first large-scale experimental data sets for supramolecular polymer systems.

ML-guided polymer design↗

Neutrinos from the Galactic Center Hosting a Hypernova Remnant

Similar to star-forming galaxies or starburst galaxies, star-forming regions in our Galaxy can host cosmic-ray (CR) accelerators and rich gas as targets of hadronuclear interaction. By our estimations, the IceCube neutrino observatory might detect muon neutrinos from a CR accelerator associated with a molecular cloud complex in our Galaxy. The associated high-energy gamma-ray emission might be observed by the Cherenkov Telescope Array (CTA), High-Altitude Water Cherenkov Gamma-Ray Observatory (HAWC), and Large High Altitude Air Shower Observatory (LHAASO). Furthermore, taking the Galactic Center (GC) region as an example, we assume that a hypernova exploded in the past in the GC. We simulate the acceleration of CRs in the hypernova remnant (HNR) as well as their confinement and escape. The high-energy protons escape from the HNR, diffuse around the GC, interact with molecular clouds, and then produce gamma-rays and neutrinos. In the optimal cases, the GC would be a promising 100 TeV gamma-ray source for LHAASO’s one-month observation. We propose that neutrino-induced searching for starting track-like and high-energy starting events (HESEs) observed by IceCube, from the GC region with a radius of 1°.8, would help us discover the particle accelerator in the GC or constrain our models. Under the constraint from high-energy gamma-ray observations by the H.E.S.S. telescope, we estimate the exposure time needed to make a significant discovery for the optimal cases. The analysis combining observations of IceCube and ANTARES, starting track-like events and HESEs, future observations by neutrino detectors IceCube-Gen2 and KM3net, and gamma-ray telescopes CTA, HAWC, and LHAASO would help to constrain our models.

79 ASTRONOMY AND ASTROPHYSICS↗

EGRET upper limits to the high-energy gamma-ray emission from the millisecond pulsars in nearby globular clusters

We report upper limits to the high-energy gamma-ray emission from the millisecond pulsars (MSPs) in a number of globular clusters. The observations were done as part of an all-sky survey by the energetic Gamma Ray Experiment Telescope (EGRET) on the Compton Gamma Ray Observatory (CGRO) during Phase I of the CGRO mission (1991 June to 1992 November). Several theoretical models suggest that MSPs may be sources of high-energy gamma radiation emitted either as primary radiation from the pulsar magnetosphere or as secondary radiation generated by conversion into photons of a substantial part of the relativistic e(+/-) pair wind expected to flow from the pulsar. To date, no high-energy emission has been detected from an individual MSP. However, a large number of MSPs are expected in globular cluster cores where the formation rate of accreting binary systems is high. Model predictions of the total number of pulsars range in the hundreds for some clusters. These expectations have been reinforced by recent discoveries of a substantial number of radio MSPs in several clusters; for example, 11 have been found in 47 Tucanae (Manchester et al.). The EGRET observations have been used to obtain upper limits for the efficiency eta of conversion of MSP spin-down power into hard gamma rays. The upper limits are also compared with the gamma-ray fluxes predicted from theoretical models of pulsar wind emission (Tavani). The EGRET limits put significant constraints on either the emission models or the number of pulsars in the globular clusters.

Michelson, P. F.↗

Propane absorption band intensities and band model parameters from 680 to 1580/cm at 296 and 200 K

Band intensities and profiles have been measured for the propane absorption bands from 680 to 1580/cm at 296 and 200 K. This work was stimulated by the discovery of several propane bands in the spectrum of Titan by the Voyager 1 spacecraft. The low temperature laboratory data show that the bands become narrower and the Q branches of the bands somewhat stronger than they are at room temperature. Random band model parameters were determined over the entire region from the 42 spectra obtained at room temperature.

Giver, L. P.↗

Quantitative Studies of the Optical and UV Spectra of Galactic Early B Supergiants

We undertake an optical and ultraviolet spectroscopic analysis of a sample of 20 Galactic B0-B5 supergiants of luminosity classes Ia, Ib, Iab, and II. Fundamental stellar parameters are obtained from optical diagnostics and a critical comparison of the model predictions to observed UV spectral features is made. Methods. Fundamental parameters (e.g., T(sub eff), log L(sub *), mass-loss rates and CNO abundances) are derived for individual stars using CMFGEN, a nLTE, line-blanketed model atmosphere code. The impact of these newly derived parameters on the Galactic B supergiant Ten scale, mass discrepancy, and wind-momentum luminosity relation is examined. Results. The B supergiant temperature scale derived here shows a reduction of about 1000-3000 K compared to previous results using unblanketed codes. Mass-loss rate estimates are in good agreement with predicted theoretical values, and all of the 20 BO-B5 supergiants analysed show evidence of CNO processing. A mass discrepancy still exists between spectroscopic and evolutionary masses, with the largest discrepancy occuring at log (L/(solar)L approx. 5.4. The observed WLR values calculated for B0-B0.7 supergiants are higher than predicted values, whereas the reverse is true for B1-B5 supergiants. This means that the discrepancy between observed and theoretical values cannot be resolved by adopting clumped (i.e., lower) mass-loss rates as for O stars. The most surprising result is that, although CMFGEN succeeds in reproducing the optical stellar spectrum accurately, it fails to precisely reproduce key UV diagnostics, such as the N v and C IV P Cygni profiles. This problem arises because the models are not ionised enough and fail to reproduce the full extent of the observed absorption trough of the P Cygni profiles. Conclusions. Newly-derived fundamental parameters for early B supergiants are in good agreement with similar work in the field. The most significant discovery, however, is the failure of CMFGEN to predict the correct ionisation fraction for some ions. Such findings add further support to revising the current standard model of massive star winds, as our understanding of these winds is incomplete without a precise knowledge of the ionisation structure and distribution of clumping in the wind. Key words. techniques: spectroscopic - stars: mass-loss - stars: supergiants - stars: abundances - stars: atmospheres - stars: fundamental parameters

Searle, S. C.↗

Discovery of the Optical Afterglow and Host Galaxy of Short GRB 181123B at z = 1.754: Implications for Delay Time Distributions

We present the discovery of the optical afterglow and host galaxy of the Swift short-duration gamma-ray burst (SGRB) GRB 181123B. Observations with Gemini-North starting ≈9.1 hr after the burst reveal a faint optical afterglow with i ≈ 25.1 mag at an angular offset of 0 59 ± 0." 16 from its host galaxy. Using grizYJHK observations, we measure a photometric redshift of the host galaxy of z = 1.77+0.30-0.17. From a combination of Gemini and Keck spectroscopy of the host galaxy spanning 4500–18000 Å, we detect a single emission line at 13390 Å, inferred as Hβ at z = 1.754 ± 0.001 and corroborating the photometric redshift. The host galaxy properties of GRB 181123B are typical of those of other SGRB hosts, with an inferred stellar mass of ≈9.1 × 109Msun, a mass-weighted age of ≈0.9 Gyr, and an optical luminosity of ≈0.9L*. At z = 1.754, GRB 181123B is the most distant secure SGRB with an optical afterglow detection and one of only three at z > 1.5. Motivated by a growing number of high-z SGRBs, we explore the effects of a missing z > 1.5 SGRB population among the current Swift sample on delay time distribution (DTD) models. We find that lognormal models with mean delay times of ≈4–6 Gyr are consistent with the observed distribution but can be ruled out to 95% confidence, with an additional ≈one to five Swift SGRBs recovered at z > 1.5. In contrast, power-law models with ∝t−1 are consistent with the redshift distribution and can accommodate up to ≈30 SGRBs at these redshifts. Under this model, we predict that ≈1/3 of the current Swift population of SGRBs is at z > 1. The future discovery or recovery of existing high-z SGRBs will provide significant discriminating power on their DTDs and thus their formation channels.

K Paterson↗

A Scientist-in-the-Loop Data Analytics Framework for Intelligent Simulation Model Tuning and Validation

This project developed a scientist-in-the-loop data analytics framework for intelligent simulation model tuning and validation, targeting the Weather Research and Forecasting (WRF) model and its solar energy variant, WRF-Solar-BNL. Domain experts, such as climate scientists, depend on large-scale numerical simulations for knowledge discovery and decision-making, yet the complexity of parameter tuning and the disconnect between automated optimization and domain expertise pose significant challenges. We extended an interactive visual analytics framework that enables domain experts to observe and intervene in the computational steering process by identifying disagreements between the simulation model, surrogate model, and the expert’s domain knowledge. Using Bayesian Optimization with Gaussian Process Regression as the surrogate model, our system allows users to probe parameter relationships, analyze correlation patterns, and adjust tuning parameters in real time. We developed use cases for solar irradiance forecasting through sustained collaboration with Brookhaven National Laboratory, resolving critical model configuration challenges and achieving meaningful reductions in prediction error. The project supported one PhD student, one MS student, and eight undergraduate students across three Data Science Capstone projects, resulting in one master’s thesis.

Dasgupta, Aritra [New Jersey Institute of Technolo↗

A global biogeocenotical biosphere simulation

This model of the D. Forrester type, constructed in differential equations, predicts the food and mineral supply for the factors biosphere population, depending on two socio-economic factors, until about the year 2500. If contemporary rates of natural resources utilization are maintained and there is no management of the environment, food resources will begin to limit human population growth after 2200, and mineral resources will after 2300. A decrease in the biosphere pollution, increase in effective agricultural production, and discovery of new energy sources may forestall or completely avert the onset of a crisis situation. Conservation measures, according to the model, are to a considerable extent realizable only if carried out simultaneously in both areas.

Moiseyev, N. N.↗

Large-Scale NASA Science Applications on the Columbia Supercluster

Columbia, NASA's newest 61 teraflops supercomputer that became operational late last year, is a highly integrated Altix cluster of 10,240 processors, and was named to honor the crew of the Space Shuttle lost in early 2003. Constructed in just four months, Columbia increased NASA's computing capability ten-fold, and revitalized the Agency's high-end computing efforts. Significant cutting-edge science and engineering simulations in the areas of space and Earth sciences, as well as aeronautics and space operations, are already occurring on this largest operational Linux supercomputer, demonstrating its capacity and capability to accelerate NASA's space exploration vision. The presentation will describe how an integrated environment consisting not only of next-generation systems, but also modeling and simulation, high-speed networking, parallel performance optimization, and advanced data analysis and visualization, is being used to reduce design cycle time, accelerate scientific discovery, conduct parametric analysis of multiple scenarios, and enhance safety during the life cycle of NASA missions. The talk will conclude by discussing how NAS partnered with various NASA centers, other government agencies, computer industry, and academia, to create a national resource in large-scale modeling and simulation.

Brooks, Walter↗

Coding and noncoding variants in EBF3 are involved in HADDS and simplex autism

Abstract Background Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737. Results We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10 −3 ), and combined dataset (p = 1.1 × 10 −4 ). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10 −35 , loss-of-function p = 2.26 × 10 −13 ) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10 −6 , OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia. Conclusions In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs.

59 BASIC BIOLOGICAL SCIENCES↗

The Role of Theory and Modeling in the International Living with a Star Program

Today, theory and modeling play a critical role in our quest to understand the connection between solar eruptive phenomena, and their impacts in interplanetary space and in the near-Earth space environment. This new role is based on two developments, one related to the goal of basic physical understanding, and the other to space weather-related applications. When targeting physical our focus is shifting away from investigations aiming at basic discoveries, to missions and studies that address our basic understanding of processes we know to be important. For these studies, theory and models provide physical explanations that need to be verified or falsified by empirical evidence. Within this paradigm, a much more tight integration between theory modeling, and space flight mission design and execution is not only beneficial, but essential. One of the prime objectives of space weather research, on the other hand, is the prediction of space environmental conditions for the benefit of humans and their assets in near-Earth space and on the ground, as well as on solar system bodies like Mars that are of interest to exploration by humans. By its very nature, prediction requires modeling, which, in turn, requires understanding. We will present an overview of the role of theory and modeling within the International Living With a Star program. Specifically, we will focus on an assessment of present-day and future capabilities, as well as on strategies for tight integration of theory and modeling in space science investigations.

Hesse, M.↗

AI to Automate ModEx for Optimal Predictive Improvement and Scientific Discovery

Focal Areas: Data acquisition and assimilation enabled by machine learning, AI, advanced experimental optimization, unsupervised learning, and hardware-related AI efforts; Predictive modeling through AI techniques and AI-derived model components; Using AI to design a hierarchical model prediction system consisting and model selection; and, Interrogating complex data (observed and simulated) using AI, big data analytics, and other advanced methods such as explainable AI and physics- or knowledge-guided AI.

54 ENVIRONMENTAL SCIENCES↗

Magnitude-squared Coherence: A Powerful Tool for Disentangling Doppler Planet Discoveries from Stellar Activity

Abstract If Doppler searches for Earth-mass, habitable planets are to succeed, observers must be able to identify and model out stellar activity signals. Here we demonstrate how to diagnose activity signals by calculating the magnitude-squared coherence C ˆ xy 2 ( f ) between an activity-indicator time series x t and the radial-velocity (RV) time series y t . Since planets only cause modulation in RV, not in activity indicators, a high value of C ˆ xy 2 ( f ) indicates that the signal at frequency f has a stellar origin. We use Welch’s method to measure coherence between activity indicators and RVs in archival observations of GJ 581, α Cen B, and GJ 3998. High RV-H α coherence at the frequency of GJ 3998 b and high RV-S index coherence at the frequency of GJ 3998c, indicate that the planets may actually be stellar signals. We also replicate previous results showing that GJ 581 d and g are rotation harmonics and demonstrate that α Cen B has activity signals that are not associated with rotation. Welch’s power spectrum estimates have cleaner spectral windows than Lomb–Scargle periodograms, improving our ability to estimate rotation periods. We find that the rotation period of GJ 581 is 132 days, with no evidence of differential rotation. Welch’s method may yield unacceptably large bias for data sets with N < 75 observations and works best on data sets with N > 100. Tapering the time-domain data can reduce the bias of the Welch’s power spectrum estimator, but observers should not apply tapers to data sets with extremely uneven observing cadence. A software package for calculating magnitude-squared coherence and Welch’s power spectrum estimates is available on github.

Dodson-Robinson, Sarah E. (ORCID:0000000287964974)↗

Finding predictive models for singlet fission by machine learning

Singlet fission (SF), the conversion of one singlet exciton into two triplet excitons, could significantly enhance solar cell efficiency. Molecular crystals that undergo SF are scarce. Computational exploration may accelerate the discovery of SF materials. However, many-body perturbation theory (MBPT) calculations of the excitonic properties of molecular crystals are impractical for large-scale materials screening. We use the sure-independence-screening-and-sparsifying-operator (SISSO) machine-learning algorithm to generate computationally efficient models that can predict the MBPT thermodynamic driving force for SF for a dataset of 101 polycyclic aromatic hydrocarbons (PAH101). SISSO generates models by iteratively combining physical primary features. The best models are selected by linear regression with cross-validation. The SISSO models successfully predict the SF driving force with errors below 0.2 eV. Based on the cost, accuracy, and classification performance of SISSO models, we propose a hierarchical materials screening workflow. Three potential SF candidates are found in the PAH101 set.

36 MATERIALS SCIENCE↗

Hydrostatic Stress Effects Incorporated Into the Analysis of the High-Strain-Rate Deformation of Polymer Matrix Composites

Procedures for modeling the effect of high strain rate on composite materials are needed for designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and strain rate dependence of the composite response is primarily due to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. An experimental program has been carried out through a university grant with the Ohio State University to obtain tensile and shear deformation data for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. This information has been used at the NASA Glenn Research Center to develop, characterize, and correlate a material model in which the strain rate dependence and nonlinearity (including hydrostatic stress effects) of the polymer are correctly analyzed. To obtain the material data, Glenn s researchers designed and fabricated test specimens of a representative toughened epoxy resin. Quasi-static tests at low strain rates and split Hopkinson bar tests at high strain rates were then conducted at the Ohio State University. The experimental data confirmed the strong effects of strain rate on both the tensile and shear deformation of the polymer. For the analytical model, Glenn researchers modified state variable constitutive equations previously used for the viscoplastic analysis of metals to allow for the analysis of the nonlinear, strain-rate-dependent polymer deformation. Specifically, we accounted for the effects of hydrostatic stresses. An important discovery in the course of this work was that the hydrostatic stress effects varied during the loading process, which needed to be accounted for within the constitutive equations. The model is characterized primarily by shear data, with tensile data used to characterize the hydrostatic stress effects.

Goldberg, Robert K.↗

Materials representation and transfer learning for multi-property prediction

The adoption of machine learning in materials science has rapidly transformed materials property prediction. Hurdles limiting full capitalization of recent advancements in machine learning include the limited development of methods to learn the underlying interactions of multiple elements as well as the relationships among multiple properties to facilitate property prediction in new composition spaces. To address these issues, we introduce the Hierarchical Correlation Learning for Multi-property Prediction (H-CLMP) framework that seamlessly integrates: (i) prediction using only a material's composition, (ii) learning and exploitation of correlations among target properties in multi-target regression, and (iii) leveraging training data from tangential domains via generative transfer learning. The model is demonstrated for prediction of spectral optical absorption of complex metal oxides spanning 69 three-cation metal oxide composition spaces. H-CLMP accurately predicts non-linear composition-property relationships in composition spaces for which no training data are available, which broadens the purview of machine learning to the discovery of materials with exceptional properties. This achievement results from the principled integration of latent embedding learning, property correlation learning, generative transfer learning, and attention models. The best performance is obtained using H-CLMP with transfer learning [H-CLMP(T)] wherein a generative adversarial network is trained on computational density of states data and deployed in the target domain to augment prediction of optical absorption from composition. H-CLMP(T) aggregates multiple knowledge sources with a framework that is well suited for multi-target regression across the physical sciences.

36 MATERIALS SCIENCE↗

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coalescence dynamics of prismatic dislocation loops due to vacancy supersaturation

Recovery processes in metals with a high defect concentration are at the heart of a structural alloy's usability in engineering applications. However, such processes take place via nonconservative mechanisms mediated by point defect diffusion and are governed by their interactions with the underlying microstructure. Indeed, one of the main challenges in studying recovery processes is their intrinsic multiscale complexity, which often negates the direct use of experimental techniques. Here, in this Research Letter, we study the mechanism of dislocation loop coalescence in metals by simulating the interaction dynamics between vacancies and prismatic loops using a stochastic framework that seamlessly bridges point defect kinetics with dislocation mechanics. Our simulations reveal two fundamental discoveries: first, that climb processes are governed by local fluctuations that are not captured by continuum elastodiffusion models, and, second, that prismatic loops create internal traps for vacancies, leading to the formation of rings around the perimeter of the loops that affect how coalescence takes place. Our results show an excellent agreement with in situ experimental observations of loop coalescence in Fe under irradiation.

36 MATERIALS SCIENCE↗