Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spreading models”

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 847 records · Page 47

Using convolutional neural networks to accelerate three-dimensional coherent synchrotron radiation computations

Calculating the effects of coherent synchrotron radiation (CSR) is one of the most computationally expensive tasks in accelerator physics. Here, we use convolutional neural networks (CNNs), along with a latent conditional diffusion (LCD) model, trained on physics-based simulations to speed up calculations. Specifically, we produce the 3D CSR wakefields generated by electron bunches in circular orbit in the steady-state condition. Two datasets are used for training and testing the models: wakefields generated by three-dimensional Gaussian electron distributions and wakefields from a sum of up to 25 three-dimensional Gaussian distributions. The CNNs are able to accurately produce the 3D wakefields ∼250–1000 times faster than the numerical calculations, while the LCD achieves a gain of a factor of ∼34. We also test the extrapolation and out-of-distribution generalization ability of the models. They generalize well on distributions with larger spreads than what they were trained on but struggle with smaller spreads.

43 PARTICLE ACCELERATORS↗

Deriving Iodine-free Spectra for High-resolution Echelle Spectrographs

We describe a new method to derive clean, iodine-free spectra directly from observations acquired using high-resolution echelle spectrographs equipped with iodine cells. The main motivation to obtain iodine-free spectra is to use portions of the spectrum that are superimposed with the dense forest of iodine absorption lines, in order to retrieve lines that can be used to monitor the magnetic activity of the star, helping to validate candidate planets. In short, we provide a straightforward methodology to clean the spectra using the forward model used to derive radial velocities, the line spread function information plus the stellar spectrum without iodine to reconstruct and subtract the iodine spectrum from the observations. We show our results using observations of the star τ Ceti acquired with the Planet Finder Spectrograph (PFS), High Resolution Echelle Spectrometer (HIRES), and University College London Echelle Spectrograph (UCLES), reaching an iodine-free spectrum correction at the ∼1% rms level. We additionally discuss the limitations and further applications of the method.

79 ASTRONOMY AND ASTROPHYSICS↗

Photometric and Kinematic Study of the Open Clusters SAI 44 and SAI 45

We carry out a detailed photometric and kinematic study of the poorly studied sparse open clusters SAI 44 and SAI 45 using ground-based BVR {sub c} I {sub c} data supplemented by archival data from Gaia eDR3 and Pan-STARRS. The stellar memberships are determined using a statistical method based on Gaia eDR3 kinematic data, and we found 204 members in SAI 44 while only 74 members are identified in SAI 45. The average distances to SAI 44 and SAI 45 are calculated to be 3670 ± 184 and 1668 ± 47 pc. The logarithmic age of the clusters are determined to be 8.82 ± 0.10 and 9.07 ± 0.10 yr for SAI 44 and SAI 45, respectively. The color–magnitude diagram of SAI 45 hosts an extended main-sequence turnoff (eMSTO). The apparent age spread is found to be similar to the apparent age spread predicted on the basis of the age spread and cluster age relation predicted by rotation models. This indicates that eMSTO is a stellar evolution rather than star formation phenomenon in SAI 45. We conclude that eMSTO in SAI 45 is mainly caused by the different rotation rates of stars as the SYCLIST synthetic population with different rotation rates was able to reproduce the observed eMSTO, and stars in the red part of the eMSTO were preferentially concentrated in the inner region, which again hints at different rotations being the reason for the extension in the upper MS. This finding supports the theory attributing the origin of eMSTO to the different rotations of eMSTO stars. The mass function slopes are obtained as -2.24 ± 0.66 and -2.58 ± 3.20 in the mass rages 2.426–0.990 M{sub ⊙} and 2.167–1.202 M{sub ⊙} for SAI 44 and SAI 45, respectively. SAI 44 exhibits the signature of mass segregation while we found weak evidence of mass segregation in SAI 45 possibly due to tidal stripping. The dynamical relaxation times of these clusters indicate that both clusters are in a dynamically relaxed state. Using the AD-diagram method, the apex coordinates are found to be ([[Formula]], -[[Formula]]) for SAI 44 and (-[[Formula]], -[[Formula]]) for SAI 45. The average space velocity components of the clusters SAI 44 and SAI 45 are calculated in units of km s{sup -1} as (-15.14 ± 3.90, -19.43 ± 4.41, -20.85 ± 4.57) and (28.13 ± 5.30, -9.78 ± 3.13, -19.59 ± 4.43), respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dissecting Nearby Galaxies with piXedfit. I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

We study spatially resolved properties (on spatial scales of ~1–2 kpc out to at least 3 effective radii) of the stars, dust, and gas in 10 nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point-spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy-balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX–βrelation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of the star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kiloparsec scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), and gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of the molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center.

79 ASTRONOMY AND ASTROPHYSICS↗

Isochronal Age-mass Discrepancy of Young Stars: SCExAO/CHARIS Integral Field Spectroscopy of the HIP 79124 Triple System

We present SCExAO/CHARIS 1.1-2.4 μm integral field direct spectroscopy of the young HIP 79124 triple system. HIP 79124 is a member of the Scorpius-Centaurus association, consisting of an A0V primary with two low-mass companions at a projected separation of < 1″. Thanks to the high quality wavefront corrections provided by SCExAO, both companions are decisively detected without the employment of any PSF-subtraction algorithm to eliminate quasi-static noise. The spectrum of the outer C object is very well matched by Upper Scorpius M4 ± 0.5 standard spectra, with a T(eff) = 2945 ± 100 K and a mass of ∼350 M(Jup). HIP 79124 B is detected at a separation of only 180 mas in a highly-correlated noise regime, and it falls in the spectral range M6 ± 0.5 with T(eff) = 2840 ± 190 K and ∼100 M(Jup). Previous studies of stellar populations in Sco-Cen have highlighted a discrepancy in isochronal ages between the lower-mass and higher-mass populations. This could be explained either by an age spread in the region, or by conventional isochronal models failing to reproduce the evolution of low-mass stars. The HIP 79124 system should be coeval, and therefore it provides an ideal laboratory to test these scenarios. We place the three components in a color-magnitude diagram and find that the models predict a younger age for the two low-mass companions (∼3 Myr) than for the primary star (∼6 Myr). These results imply that the omission of magnetic effects in conventional isochronal models inhibit them from reproducing early low-mass stellar evolution, which is further supported by the fact that new models that include such effects provide more consistent ages in the HIP 79124 system.

Ruben Asensio-Torres↗

Reappraisal of the Climate Impacts of Ozone‐Depleting Substances

We assess the effective radiative forcing due to ozone‐depleting substances using models participating in the Aerosols and Chemistry and Radiative Forcing Model Intercomparison Projects (AerChemMIP, RFMIP). A large intermodel spread in this globally averaged quantity necessitates an “emergent constraint” approach whereby we link the radiative forcing to ozone declines measured and simulated during 1979–2000, excluding two volcanically perturbed periods. During this period, ozone‐depleting substances were increasing, and several merged satellite‐based climatologies document the ensuing decline of total‐column ozone. Using these analyses, we find an effective radiative forcing of −0.05 to 0.13 W m−2. Our best estimate (0.04 W m−2) is on the edge of the “likely” range given by the Fifth Assessment Report of IPCC of 0.03 to 0.33 W m−2 but is in better agreement with two other literature results.

ozone depletion↗

Coronado Ecological Conservation: Assessing Vegetation Change Due to Border Wall Construction and Shifting Social Trails

Species monitoring is essential in mitigating the impacts of plant invasion, such as radical changes in an area’s ecosystem, degraded soil health, increased wildfire severity, landslides, and increased flooding. NASA DEVELOP partnered with the National Park Service (NPS) to investigate invasive species in disturbed lands: specifically, areas affected by off-trail walking and US-Mexico border construction activities. The team assessed how construction has impacted the distribution of Lehmann’s lovegrass and Russian thistle invasives throughout Coronado National Memorial, AZ from 1986 to 2022. Using data from Landsat 5 and 8, Sentinel-2, the National Agriculture Imagery Program, and PlanetScope, the team computed vegetation indices including the Normalized Difference Vegetation Index, Normalized Difference Moisture Index, Modified Soil Adjusted Vegetation Index 2, Enhanced Vegetation Index, and Tasseled Cap Wetness, Brightness, and Greenness transformations as vegetation health indicators to input into various machine learning algorithms. To minimize noise, the team conducted Principal Component Analysis on the vegetation indices and spectral bands before running k-means++ clustering and random forest classification algorithms. Between all datasets, we found the median area fully overtaken by invasive plants was 5.37% of the park’s total area in 2022. The NPS will use the end products to help increase restoration efforts in disturbed areas with high concentrations of invasive plants. The NPS’s collection of ground data for 2022–2023, in conjunction with future data collection, will notably improve the accuracy of classification models, leading to more precise monitoring of invasive spread over time.

Carson Schuetze↗

Neutron imaging with coded sources: new challenges and the implementation of a simultaneous iterative reconstruction technique

The limitations in neutron flux and resolution (L/D) of current neutron imaging systems can be addressed with a Coded Source Imaging system with magnification (xCSI). More precisely, the multiple sources in an xCSI system can exceed the flux of a single pinhole system for several orders of magnitude, while maintaining a higher L/D with the small sources. Moreover, designing for an xCSI system reduces noise from neutron scattering, because the object is placed away from the detector to achieve magnification. However, xCSI systems are adversely affected by correlated noise such as non-uniform illumination of the neutron source, incorrect sampling of the coded radiograph, misalignment of the coded masks, mask transparency, and the imperfection of the system Point Spread Function (PSF). We argue that a model-based reconstruction algorithm can overcome these problems and describe the implementation of a Simultaneous Iterative Reconstruction Technique algorithm for coded sources. Design pitfalls that preclude a satisfactory reconstruction are documented.

Santos-Villalobos, Hector J↗

A Monte Carlo model for 3D grain evolution during welding

Welding is one of the most wide-spread processes used in metal joining. However, there are currently no open-source software implementations for the simulation of microstructural evolution during a weld pass. Here we describe a Potts Monte Carlo based model implemented in the SPPARKS kinetic Monte Carlo computational framework. The model simulates melting, solidification and solid-state microstructural evolution of material in the fusion and heat-affected zones of a weld. The model does not simulate thermal behavior, but rather utilizes user input parameters to specify weld pool and heat-affect zone properties. Weld pool shapes are specified by Bezier curves, which allow for the specification of a wide range of pool shapes. Pool shapes can range from narrow and deep to wide and shallow representing different fluid flow conditions within the pool. Surrounding temperature gradients are calculated with the aide of a closest point projection algorithm. Furthermore, the model also allows simulation of pulsed power welding through time-dependent variation of the weld pool size. Example simulation results and comparisons with laboratory weld observations demonstrate microstructural variation with weld speed, pool shape, and pulsed-power.

36 MATERIALS SCIENCE↗

Spatially distributed infection increases viral load in a computational model of SARS-CoV-2 lung infection

A key question in SARS-CoV-2 infection is why viral loads and patient outcomes vary dramatically across individuals. Because spatial-temporal dynamics of viral spread and immune response are challenging to study in vivo, we developed Spatial Immune Model of Coronavirus (SIMCoV), a scalable computational model that simulates hundreds of millions of lung cells, including respiratory epithelial cells and T cells. SIMCoV replicates viral growth dynamics observed in patients and shows how spatially dispersed infections can lead to increased viral loads. The model also shows how the timing and strength of the T cell response can affect viral persistence, oscillations, and control. By incorporating spatial interactions, SIMCoV provides a parsimonious explanation for the dramatically different viral load trajectories among patients by varying only the number of initial sites of infection and the magnitude and timing of the T cell immune response. When the branching airway structure of the lung is explicitly represented, we find that virus spreads faster than in a 2D layer of epithelial cells, but much more slowly than in an undifferentiated 3D grid or in a well-mixed differential equation model. These results illustrate how realistic, spatially explicit computational models can improve understanding of within-host dynamics of SARS-CoV-2 infection.

59 BASIC BIOLOGICAL SCIENCES↗

Cold atomic beam ion source for focused ion beam applications

We report measurements and modeling of an ion source that is based on ionization of a laser-cooled atomic beam. We show a high brightness and a low energy spread, suitable for use in next-generation, high-resolution focused ion beam systems. Our measurements of total ion current as a function of ionization conditions support an analytical model that also predicts the cross-sectional current density and spatial distribution of ions created in the source. The model predicts a peak brightness of 2 × 10{sup 7} A m{sup −2} sr{sup −1} eV{sup −1} and an energy spread less than 0.34 eV. The model is also combined with Monte-Carlo simulations of the inter-ion Coulomb forces to show that the source can be operated at several picoamperes with a brightness above 1 × 10{sup 7} A m{sup −2} sr{sup −1} eV{sup −1}. We estimate that when combined with a conventional ion focusing column, an ion source with these properties could focus a 1 pA beam into a spot smaller than 1 nm. A total current greater than 5 nA was measured in a lower-brightness configuration of the ion source, demonstrating the possibility of a high current mode of operation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantifying Uncertainty in Epidemiological Models

Modern epidemiology has made use of a number of mathematical models, including ordinary differential equation (ODE) based models and agent based models (ABMs) to describe the dynamics of how a disease may spread within a population and enable the rational design of strategies for intervention that effectively contain the spread of the disease. Although such predictions are of fundamental importance in preventing the next global pandemic, there is a significant gap in trusting the outcomes/predictions solely based on such models. Hence, there is a need to develop approaches such that mathematical models can be calibrated against historical data. In addition, there is a need to develop rigorous uncertainty quantification approaches that can provide insights into when a model will fail and characterize the confidence in the (possibly multiple) model outcomes/predictions, when such retrospective analysis cannot be performed. In this paper, we outline an approach to develop uncertainty quantification approaches for epidemiological models using formal methods and model checking. By specifying the outcomes expected from a model in a suitable spatio-temporal logic, we use probabilistic model checking methods to quantify the probability with which the epidemiological model satisfies the specification. We argue that statistical model checking methods can solve the uncertainty quantification problem for complex epidemiological models.

model checking↗

Thermal and chemical convection in planetary mantles

Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Dupeyrat, L.↗

Reduced-order modeling for nonlocal diffusion problems

In several settings, diffusive behavior is observed to not follow the rate of spread predicted by parabolic partial differential equations (PDEs) such as the heat equation. Such behaviors, often referred to as anomalous diffusion, can be modeled using nonlocal equations for which points at a finite distance apart can interact. An example of such models is provided by fractional derivative equations. Because of the nonlocal interactions, discretized nonlocal systems have less sparsity, often significantly less, compared with corresponding discretized PDE systems. As such, the need for reduced-order surrogates that can be used to cheaply determine approximate solutions is much more acute for nonlocal models compared with that for PDEs. In this paper, we consider the construction, application, and testing of proper orthogonal decomposition (POD) reduced models for an integral equation model for nonlocal diffusion. For certain modeling parameters, the model we consider allows for discontinuous solutions and includes fractional Laplacian kernels as a special case. Preliminary computational results illustrate the potential of using POD to obtain accurate approximations of solutions of nonlocal diffusion equations at much lower costs compared with, for example, standard finite element methods. Copyright © 2016 John Wiley & Sons, Ltd.

97 MATHEMATICS AND COMPUTING↗

Swift and extensive Omicron outbreak in China after sudden exit from ‘zero-COVID’ policy

In late 2022, China transitioned from a strict ‘zero-COVID’ policy to rapidly abandoning nearly all interventions and data reporting. This raised great concern about the presumably-rapid but unreported spread of the SARS-CoV-2 Omicron variant in a very large population of very low pre-existing immunity. By modeling a combination of case count and survey data, we show that Omicron spread extremely rapidly, at a rate of 0.42/day (95% credibility interval: [0.35, 0.51]/day), translating to an epidemic doubling time of 1.6 days ([1.6, 2.0] days) after the full exit from zero-COVID on Dec. 7, 2022. Consequently, we estimate that the vast majority of the population (97% [95%, 99%], sensitivity analysis lower limit of 90%) was infected during December, with the nation-wide epidemic peaking on Dec. 23. Overall, our results highlight the extremely high transmissibility of the variant and the importance of proper design of intervention exit strategies to avoid large infection waves.

60 APPLIED LIFE SCIENCES↗

Comparison of structurally diverse simulation models for prediction of epidemic outcomes caused by a long-distance dispersed pathogen

Long-distance dispersal (LDD) pathogens pose substantial challenges for epidemic control due to their ability to generate new infection foci at great distances. While various modeling approaches have been developed to understand and manage such outbreaks, little work has compared how models of different structures behave under shared conditions. Here, in this study, we compare four structurally distinct epidemiological models — EPIMUL, GEMF, PoPS, and Warwick — each adapted to simulate the spread of wheat stripe rust (WSR), a wind-dispersed LDD pathogen, under identical epidemiological parameters and dispersal kernel. Using data from a controlled field experiment, we evaluate the ability of each model to replicate disease prevalence under nine intervention scenarios that vary in timing and culling area. While the models differ substantially in design — ranging from spatial grid-based to network-based and raster-based frameworks — the shared dispersal kernel allowed for close alignment in their predictions. All models accurately captured general epidemic trends, particularly the strong effect of early intervention on disease suppression. We qualitatively compared their behavioral responses across scenarios and also evaluated an ensemble prediction by averaging across model outputs. Our findings highlight how integrating shared epidemiological components into distinct modeling frameworks can improve consistency and accuracy, while reinforcing the importance of early culling in managing LDD pathogen outbreaks.

Dispersal kernel↗

Burning Questions in Gravity-Dependent Combustion Science

Building upon a long history of spaceflight and ground based research, NASA's Combustion Science program has accumulated a significant body of accomplishments on the ISS. Historically, NASAs low-gravity combustion research program has sought: to provide a more complete understanding of the fundamental controlling processes in combustion by identifying simpler one-dimensional systems to eliminate the complex interactions between the buoyant flow and the energy feedback to the reaction zone to provide realistic simulation of the fire risk in manned spacecraft and to enable practical simulation of the gravitational environment experienced by reacting systems in future spacecraft. Over the past two decades, low-gravity combustion research has focused primarily on increasing our understanding of fundamental combustion processes (e.g. droplet combustion, soot, flame spread, smoldering, and gas-jet flames). This research program was highly successful and was aided by synergistic programs in Europe and in Japan. Overall improvements were made in our ability to model droplet combustion in spray combustors (e.g. jet engines), predict flame spread, predict soot production, and detect and prevent spacecraft fires. These results provided a unique dataset that supports both an active research discipline and also spacecraft fire safety for current and future spacecraft. These experiments have been conducted using the Combustion Integrated Rack (CIR), the Microgravity Science Glovebox and the Express Rack. In this paper, we provide an overview of the earlier space shuttle experiments, the recent ISS combustion experiments in addition to the studies planned for the future. Experiments in combustion include topics such as droplet combustion, gaseous diffusion flames, solid fuels, premixed flame studies, fire safety, and super critical oxidation processes.

Urban, David↗

A Novel Centrality Measure for Network-wide Cyber Vulnerability Assessment

In this work we propose a novel formulation that models the attack and compromise on a cyber network as a combination of two parts - direct compromise of a host and the compromise occurring through the spread of the attack on the network from a compromised host. The model parameters for the nodes are a concise representation of the host profiles that can include the risky behaviors of the associated human users while the model parameters for the edges are based on the existence of vulnerabilities between each pair of connected hosts. The edge models relate to the summary representations of the corresponding attack-graphs. This results in a formulation based on Random Walk with Restart (RWR) and the resulting centrality metric can be solved for in an efficient manner through the use of sparse linear solvers. Thus the formulation goes beyond mere topological considerations in centrality computations by summarizing the host profiles and the attack graphs into the model parameters. The computational efficiency of the method also allows us to also quantify the uncertainty in the centrality measure through Monte Carlo analysis.

Centrality↗