Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model inference”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum optimal control of superconducting qubits based on machine-learning characterization

Implementing fast and high-fidelity quantum operations using open-loop quantum optimal control relies on having an accurate model of the quantum dynamics. Any deviations between this model and the complete dynamics of the device, such as the presence of spurious modes or pulse distortions, can degrade the performance of optimal controls in practice. Here, we propose an experimentally simple approach to realize optimal quantum controls tailored to the device parameters and environment while specifically characterizing this quantum system. Concretely, we use physics-inspired machine learning to infer an accurate model of the dynamics from experimentally available data and then optimize our experimental controls on this trained model. We show the power and feasibility of this approach by optimizing arbitrary single-qubit operations in detailed numerical simulations of a superconducting transmon qubit. Furthermore, we demonstrate that this framework produces an accurate description of the device dynamics under arbitrary controls, together with the precise pulses achieving arbitrary single-qubit gates with a high fidelity of ∼99.99%.

Artificial neural networks↗
DOE OSTI2015

Characterization of swift heavy ion irradiation damage in ceria

We have examined microstructural evolution in irradiated ceria (CeO2) using swift heavy ion irradiation, electron microscopy, and atomistic simulation. CeO2, a UO2 fuel surrogate, was irradiated with gold ions at an energy of 1 GeV to fluences up to 1x1014 ions/cm2. Transmission electron microscopy accompanied by electron energy loss spectroscopy showed that the ion tracks were of similar size at all fluences, and that there was no chemical change in the ion track core. Classical molecular dynamics simulations of thermal spikes in CeO2 with energy deposition of 12 and 36 keV/nm show damage consisting of isolated point defects at the lower energy and defect clusters at 36 keV/nm, with no amorphization at either energy. Inferences are drawn from modeling about density changes in the ion track and the formation of interstitial loops that shed light on features observed by electron microscopy of swift heavy ion irradiated ceria.

Yablinsky, Clarissa↗
DOE OSTI2023

Quantifying the Variation in the Number of Donors in Quantum Dots Created Using Atomic Precision Advanced Manufacturing

Atomic-precision advanced manufacturing enables unique silicon quantum electronics built on quantum dots fabricated from small numbers of phosphorus dopants. The number of dopant atoms comprising a dot plays a central role in determining the behavior of charge and spin confined to the dots and thus overall device performance. Here in this work, we use both theoretical and experimental techniques to explore the combined impact of lithographic variation and stochastic kinetics on the number of P incorporations in quantum dots made using these techniques and how this variation changes as a function of the size of the dot. Using a kinetic model of PH3 dissociation augmented with novel reaction barriers, we demonstrate that for a 2 × 3 silicon dimer window the probability that no donor incorporates goes to zero, allowing for certainty in the placement of at least one donor. However, this still comes with some uncertainty in the precise number of incorporated donors (either one or two), and this variability may still impact certain applications. We also examine the impact of the size of the initial lithographic window, finding that the incorporation fraction saturates to δ-layer-like coverage as the circumference-to-area ratio decreases. We predict that this incorporation fraction depends strongly on the dosage of the precursor and that the standard deviation of the number of incorporations scales as ~√n, as would be expected for a sequence of largely independent incorporation events. Finally, we characterize an array of 36 experimentally prepared multidonor 3 × 3 nm lithographic windows with scanning tunneling microscopy, measuring the fidelity of the lithography to the desired array and the final location of PH x fragments within these lithographic windows. We use our kinetic model to examine the expected variability due to the observed lithographic error, predicting a negligible impact on incorporation statistics. We find good agreement between our model and the inferred incorporation locations in these windows from scanning tunneling microscope measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Implementation of a practical Markov chain Monte Carlo sampling algorithm in PyBioNetFit

Abstract Summary Bayesian inference in biological modeling commonly relies on Markov chain Monte Carlo (MCMC) sampling of a multidimensional and non-Gaussian posterior distribution that is not analytically tractable. Here, we present the implementation of a practical MCMC method in the open-source software package PyBioNetFit (PyBNF), which is designed to support parameterization of mathematical models for biological systems. The new MCMC method, am, incorporates an adaptive move proposal distribution. For warm starts, sampling can be initiated at a specified location in parameter space and with a multivariate Gaussian proposal distribution defined initially by a specified covariance matrix. Multiple chains can be generated in parallel using a computer cluster. We demonstrate that am can be used to successfully solve real-world Bayesian inference problems, including forecasting of new Coronavirus Disease 2019 case detection with Bayesian quantification of forecast uncertainty. Availability and implementation PyBNF version 1.1.9, the first stable release with am, is available at PyPI and can be installed using the pip package-management system on platforms that have a working installation of Python 3. PyBNF relies on libRoadRunner and BioNetGen for simulations (e.g. numerical integration of ordinary differential equations defined in SBML or BNGL files) and Dask.Distributed for task scheduling on Linux computer clusters. The Python source code can be freely downloaded/cloned from GitHub and used and modified under terms of the BSD-3 license (https://github.com/lanl/pybnf). Online documentation covering installation/usage is available (https://pybnf.readthedocs.io/en/latest/). A tutorial video is available on YouTube (https://www.youtube.com/watch?v=2aRqpqFOiS4&t=63s). Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1998

Final Report

This grant was made to fund two accepted proposals to observe with the ASCA satellite: (1) Debunking the myth of two-temperature coronae of active stars, and (2) Dynamic coronae of dMe stars. We obtained the requested observations and have now completed two papers that will be published in the Astrophysical Journal. The RS CVn binary star UX Ari was observed for 14 hours with all four detectors onboard the ASCA satellite. The X-ray emission was at a constant, quiescent level during the first 12 hours, after which time a powerful flare with a peak luminosity of 1.4 x 10(exp 32) ergs/s started. The flare was observed until shortly after its peak. We present a spectral and temporal analysis of the UX Ari observations and analyze the data with a two-ribbon flare model including estimates for cooling losses. A time-dependent reconstruction of the emission measure (EM) distribution shows that two separate plasma components evolve during the flare (one being identified with the quiescent EM). Most of the flare EM reaches temperatures between 50 MK and 100 MK or more. Magnetic confinement requires the loop arcade to be geometrically large, with length scales on the order of one stellar radius. The electron densities inferred from the model decrease from initial values around 10(exp 12)/cc early in the flare to about 10(exp 11)/cc at the flare peak. The best-fit models require surface magnetic field strengths of a few hundred G, compatible with the maximum photospheric fields expected from equipartition. The flare parameters imply a (conductive and radiative) cooling loss time of less than one hour at flare peak. The elemental abundances increase significantly during the flare rise, with the abundances of the low-FIP elements Fe, Mg, Si, and Ni typically increasing to higher levels than the high-FIP elements such as S or Ne. The Fe abundance increases from (17 +/- 4)% of the solar photospheric value during quiescence up to (89 +/- 18)% at flare peak. A fractionation process that occurs during the chromospheric evaporation phase may selectively enrich low-FIP elements as in the solar corona; alternatively, the chromospheric evaporation may itself bring metal-rich plasma into the metal-poor corona.

Linsky, Jeffrey L.↗
NASA NTRS1999

A TRMM Rainfall Estimation Method Applicable to Land Areas

Methods developed to estimate rain rate on a footprint scale over land with the satellite-borne multispectral dual-polarization Special Sensor Microwave Imager (SSM/1) radiometer have met with limited success. Variability of surface emissivity on land and beam filling are commonly cited as the weaknesses of these methods. On the contrary, we contend a more significant reason for this lack of success is that the information content of spectral and polarization measurements of the SSM/I is limited. because of significant redundancy. As a result, the complex nature and vertical distribution C, of frozen and melting ice particles of different densities, sizes, and shapes cannot resolved satisfactorily. Extinction in the microwave region due to these complex particles can mask the extinction due to rain drops. Because of these reasons, theoretical models that attempt to retrieve rain rate do not succeed on a footprint scale. To illustrate the weakness of these models, as an example we can consider the brightness temperature measurement made by the radiometer in the 85 GHz channel (T85). Models indicate that T85 should be inversely related to the rain rate, because of scattering. However, rain rate derived from 15-minute rain gauges on land indicate that this is not true in a majority of footprints. This is also supported by the ship-borne radar observations of rain in the Tropical Oceans and Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA-COARE) region over the ocean. Based on these observations. we infer that theoretical models that attempt to retrieve rain rate do not succeed on a footprint scale. We do not follow the above path of rain retrieval on a footprint scale. Instead, we depend on the limited ability of the microwave radiometer to detect the presence of rain. This capability is useful to determine the rain area in a mesoscale region. We find in a given rain event that this rain area is closely related to the mesoscale-average rain rate. Based on this observation, in this study we have developed a method to estimate the mesoscale-average rain rate over land utilizing microwave radiometer data. Because of the high degree of geographic and seasonal variability in the nature and intensity of rain, this method requires some tuning with 15-minute rain gauge data on land. After tuning the method, it can be applied to an independent set of rain events that are close in time and space. We find that the mesoscale rain rates retrieved over the period of a month on land with this method shows a correlation of about 0.85 with respect to the surface rain-gauge observations. This mesoscale-average rain rate estimation method can be useful to extend the spatial and temporal coverage of the rainfall data provided by the Precipitation Radar on board the Tropical Rainfall Measuring Mission (TRMM) satellite.

Prabhakara, C.↗
DOE OSTI2016

G64-12 and G64-37 are carbon-enhanced metal-poor stars

We present new high-resolution chemical-abundance analyses for the well-known high proper-motion subdwarfs G64-12 and G64-37, based on very high signal-to-noise ratio spectra (S/N∼700/1) with resolving power R ∼ 95,000. These high-quality data enable the first reliable determination of the carbon abundances for these two stars; we classify them as carbon-enhanced metal-poor (CEMP) stars based on their carboni cities, which both exceed [C/Fe] = +1.0. They are sub-classified as CEMP-no Group-II stars, based on their location in the Yoon–Beers diagram of absolute carbon abundance, A(C) versus [Fe/H], as well as on the conventional diagnostic [Ba/Fe]. The relatively low absolute carbon abundances of CEMP-no stars, in combination with the high effective temperatures of these two stars (T{sub eff}∼6500 K), weakens their CH molecular features to the point that accurate carbon abundances can only be estimated from spectra with very high S/N. A comparison of the observed abundance patterns with the predicted yields from massive, metal-free supernova models reduces the inferred progenitor masses by factors of ∼2–3, and explosion energies by factors of ∼10–15, compared to those derived using previously claimed carbon-abundance estimates. There are certainly many more warm CEMP-no stars near the halo main-sequence turnoff that have been overlooked in past studies, directly impacting the derived frequencies of CEMP-no stars as a function of metallicity, a probe that provides important constraints on Galactic chemical evolution models, the initial mass function in the early universe, and first-star nucleosynthesis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2014

Turbulent Concentration of MM-Size Particles in the Protoplanetary Nebula: Scaled-Dependent Multiplier Functions

The initial accretion of primitive bodies (asteroids and TNOs) from freely-floating nebula particles remains problematic. Here we focus on the asteroids where constituent particle (read "chondrule") sizes are observationally known; similar arguments will hold for TNOs, but the constituent particles in those regions will be smaller, or will be fluffy aggregates, and are unobserved. Traditional growth-bysticking models encounter a formidable "meter-size barrier" [1] (or even a mm-cm-size barrier [2]) in turbulent nebulae, while nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids [3]. Even if growth by sticking could somehow breach the meter size barrier, other obstacles are encountered through the 1-10km size range [4]. Another clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids [5]; scenarios leading directly from independent nebula particulates to this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios [6-8]. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. The typical sizes of planetesimals and the rate of their formation [7,8] are determined by a statistical model with properties inferred from large numerical simulations of turbulence [9]. Nebula turbulence can be described by its Reynolds number Re = L/eta sup(4/3), where L = ETA alpha sup (1/2) the largest eddy scale, H is the nebula gas vertical scale height, and α the nebula turbulent viscosity parameter, and η is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time tη. In the nebula, Re is far larger than any numerical simulation can handle, so some physical model is needed to extend the results of numerical simulations to nebula conditions.

protoplanetary nebula↗
DOE OSTI2017

Deuterium charge exchange recombination spectroscopy from the top of the pedestal to the scrape off layer in H-mode plasmas

Recent completion of the thirty two channel main-ion (deuterium) charge exchange recombination spectroscopy (CER) diagnostic on DIII-D enables detailed comparisons between impurity and main-ion temperature, density, and toroidal rotation. Sixteen sightlines cover the core of the plasma and another sixteen are densely packed towards the edge, providing high resolution measurements of the pedestal and steep gradient edge region of H-mode plasmas. The complexities of the D α spectrum require fitting with a comprehensive model, as well as using iterative collisional radiative modeling to determine the underlying thermal deuterium ion properties. Large differences in the structure and magnitude of impurity (C 6+ ) and main-ion (D + ) toroidal rotation profiles are seen in the H-mode pedestal. Additionally the D + temperature can be half the value of the C 6+ temperature at the separatrix and shows more of a pedestal structure. Typically only the impurity properties are measured and the main-ion properties are either assumed to be the same, or inferred using neoclassical models, which require validation in the steep gradient region. Furthermore, these measured differences have implications for transport model validation, intrinsic rotation studies, pedestal stability, and the boundary conditions for scrape off layer and plasma material interactions studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analytical methods for superresolution dislocation identification in dark-field X-ray microscopy

In this work, we develop several inference methods to estimate the position of dislocations from images generated using dark-field X-ray microscopy (DFXM)—achieving superresolution accuracy and principled uncertainty quantification. Using the framework of Bayesian inference, we incorporate models of the DFXM contrast mechanism and detector measurement noise, along with initial position estimates, into a statistical model coupling DFXM images with the dislocation position of interest. We motivate several position estimation and uncertainty quantification algorithms based on this model. We then demonstrate the accuracy of our primary estimation algorithm on synthetic realistic DFXM images of edge dislocations in single-crystal aluminum. We conclude with a discussion of our methods’ impact on future dislocation studies and possible future research avenues.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D

Abstract The plasma and neutral density dynamics after an edge localized mode are investigated and utilized to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization and neutral density and shows significant poloidal asymmetries in both. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport. Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s −1 and diffusion coefficient of 0.05 m 2 s −1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

PLASMA INSTABILITIES IN THE CONTEXT OF CURRENT HELIUM SEDIMENTATION MODELS: DYNAMICAL IMPLICATIONS FOR THE ICM IN GALAXY CLUSTERS

Understanding whether Helium can sediment to the core of galaxy clusters is important for a number of problems in cosmology and astrophysics. All current models addressing this question are one-dimensional and do not account for the fact that magnetic fields can effectively channel ions and electrons, leading to anisotropic transport of momentum, heat, and particle diffusion in the weakly collisional intracluster medium (ICM). This anisotropy can lead to a wide variety of instabilities, which could be relevant for understanding the dynamics of heterogeneous media. In this paper, we consider the radial temperature and composition profiles as obtained from a state-of-the-art Helium sedimentation model and analyze its stability properties. We find that the associated radial profiles are unstable to different kinds of instabilities depending on the magnetic field orientation at all radii. The fastest growing modes are usually related to generalizations of the magnetothermal instability (MTI) and the heat-flux-driven buoyancy instability which operate in heterogeneous media. We find that the effect of sedimentation is to increase (decrease) the predicted growth rates in the inner (outer) cluster region. The unstable modes grow quickly compared to the sedimentation timescale. This suggests that the composition gradients as inferred from sedimentation models, which do not fully account for the anisotropic character of the weakly collisional environment, might not be very robust. Our results emphasize the subtleties involved in understanding the gas dynamics of the ICM and argue for the need of a comprehensive approach to address the issue of Helium sedimentation beyond current models.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Formation of the infrared emission lines of Mg I in the solar atmosphere

A non-LTE radiative transfer investigation of the emission lines is conducted at 7 and 12 microns using a realistic atomic model for neutral magnesium. An average quiet sun atmospheric model is used to calculate emission-line profiles that resemble the observed ones, i.e., broad absorption troughs with narrow central emission, and significant limb brightening. The charge exchange rates are found to be significant, but the effects of high-n coupling between Mg and Mg(+) together with radiative low-n transitions are of greater importance. It is confirmed that the emission cores are formed no higher than the temperature minimum region, and that the emission is caused by non-LTE effects rather than by the chromospheric temperature rise. It is inferred from the model calculations that the line core is sensitive to magnetic fields located almost 400 km above those measured in ordinary magnetograms; the gas pressure decreases 20-fold between these two heights.

Chang, E. S.↗
DOE OSTI2013

Atmospheric retrieval analysis of the directly imaged exoplanet HR 8799b

Directly imaged exoplanets are unexplored laboratories for the application of the spectral and temperature retrieval method, where the chemistry and composition of their atmospheres are inferred from inverse modeling of the available data. As a pilot study, we focus on the extrasolar gas giant HR 8799b, for which more than 50 data points are available. We upgrade our non-linear optimal estimation retrieval method to include a phenomenological model of clouds that requires the cloud optical depth and monodisperse particle size to be specified. Previous studies have focused on forward models with assumed values of the exoplanetary properties; there is no consensus on the best-fit values of the radius, mass, surface gravity, and effective temperature of HR 8799b. We show that cloud-free models produce reasonable fits to the data if the atmosphere is of super-solar metallicity and non-solar elemental abundances. Intermediate cloudy models with moderate values of the cloud optical depth and micron-sized particles provide an equally reasonable fit to the data and require a lower mean molecular weight. We report our best-fit values for the radius, mass, surface gravity, and effective temperature of HR 8799b. The mean molecular weight is about 3.8, while the carbon-to-oxygen ratio is about unity due to the prevalence of carbon monoxide. Our study emphasizes the need for robust claims about the nature of an exoplanetary atmosphere to be based on analyses involving both photometry and spectroscopy and inferred from beyond a few photometric data points, such as are typically reported for hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Reaction rates for mesoscopic reaction-diffusion kinetics

The mesoscopic reaction-diffusion master equation (RDME) is a popular modeling framework frequently applied to stochastic reaction-diffusion kinetics in systems biology. The RDME is derived from assumptions about the underlying physical properties of the system, and it may produce unphysical results for models where those assumptions fail. In that case, other more comprehensive models are better suited, such as hard-sphere Brownian dynamics (BD). Although the RDME is a model in its own right, and not inferred from any specific microscale model, it proves useful to attempt to approximate a microscale model by a specific choice of mesoscopic reaction rates. In this paper we derive mesoscopic scale-dependent reaction rates by matching certain statistics of the RDME solution to statistics of the solution of a widely used microscopic BD model: the Smoluchowski model with a Robin boundary condition at the reaction radius of two molecules. We also establish fundamental limits on the range of mesh resolutions for which this approach yields accurate results and show both theoretically and in numerical examples that as we approach the lower fundamental limit, the mesoscopic dynamics approach the microscopic dynamics. Finally, we show that for mesh sizes below the fundamental lower limit, results are less accurate. Thus, the lower limit determines the mesh size for which we obtain the most accurate results.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2015

Turbulent Concentration of mm-Size Particles in the Protoplanetary Nebula: Scale-Dependent Cascades

The initial accretion of primitive bodies (here, asteroids in particular) from freely-floating nebula particles remains problematic. Traditional growth-by-sticking models encounter a formidable "meter-size barrier" (or even a mm-to-cm-size barrier) in turbulent nebulae, making the preconditions for so-called "streaming instabilities" difficult to achieve even for so-called "lucky" particles. Even if growth by sticking could somehow breach the meter size barrier, turbulent nebulae present further obstacles through the 1-10km size range. On the other hand, nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids. Theoretical understanding of nebula turbulence is itself in flux; recent models of MRI (magnetically-driven) turbulence favor low-or- no-turbulence environments, but purely hydrodynamic turbulence is making a comeback, with two recently discovered mechanisms generating robust turbulence which do not rely on magnetic fields at all. An important clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids; scenarios leading directly from independent nebula particulates to large objects of this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. There is evidence that at least the ordinary chondrite parent bodies were initially composed entirely of a homogeneous mix of such particles. Thus, while they are arcane, turbulent concentration models acting directly on chondrule size particles are worthy of deeper study. The typical sizes of planetesimals and the rate of their formation can be estimated using a statistical model with properties inferred from large numerical simulations of turbulence. Nebula turbulence is described by its Reynolds number Re = (L/eta)(exp 4/3), where L = H alpha(exp 1/2) is the largest eddy scale, H is the nebula gas vertical scale height, alpha the turbulent viscosity parameter, and eta is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time t(sub eta). In the nebula, Re is far larger than any numerical simulation can handle, so some physical arguments are needed to extend the results of numerical simulations to nebula conditions. In this paper, we report new physics to be incorporated into our statistical models.

mm-size particles↗
DOE OSTI2022

Multimodal Freight Energy Model for Emerging Freight Technology Analysis

While freight movement represents a small portion of the volume in the transportation sector, they are a critical contributor for energy consumed in that sector. To reduce logistics costs, energy consumption and negative environmental impact, emerging technologies, such as digitalization of logistics and alternative powertrain, have been developed and extended applications in freight. These trends are expected to grow and provide opportunities for greater efficiencies in freight movement and corresponding energy use. Still, the complexity of freight systems present challenges in evaluating these benefits, especially in the multimodal inter-city freight. Addressing the research need, this paper develops a multimodal freight energy modeling framework for the analysis of emerging freight technology scenarios. The framework is a bi-level optimization problem: network cost minimization problem (lower-level) and energy minimization problem (upper-level). The lower-level problem is a mode-path assignment problem in multimodal inter-city freight networks, where commodity-specific congestion effects on trans-shipment links are considered. For this model, an inverse modeling approach is applied to infer parameters of the lower-level model. The upper-level problem is designed to search for an optimal scenario that has the lowest energy consumption among different levels of technology applications. The proposed model is empirically tested to analyze truck load-pooling and multimodal load-pooling scenarios through stand-alone and mixed applications using freight shipments originating from or destined to the Chicago region. This framework can be used to explore the impact of emerging freight technologies on mode-path freight flow and energy consumption in the national multimodal freight network.

ADVANCED PROPULSION SYSTEMS↗
NASA NTRSDate not supplied

ARCTRON: A Rapid Experimental Proving Ground for TPS Experiments and Arcjet Technology Development

Innovation in high-enthalpy facilities is fundamentally limited by the cost and risk of experimentation. New concepts for plasma control, diagnostics, facility components, and plasma-material interaction often require repeated iterations that are impractical to perform in production arcjets. As a result, promising ideas may remain unexplored or reach operational facilities only after significant development effort. ARCTRON is being developed as a rapid experimental proving ground where new ideas in plasma science, arcjet engineering, diagnostics, and material response can be conceived, tested, and quantitatively evaluated before transition to large-scale facilities. The platform combines radio-frequency (RF) and DC arc plasma generation, externally applied magnetic fields, configurable gas composition, reduced-pressure operation, laser heating, electrical biasing, and modular diagnostic access. These capabilities permit the plasma source, applied forcing, test article, and measurement configuration to be modified independently, allowing individual physical mechanisms to be isolated more readily than in a traditional test environment. One class of investigations addresses fundamental plasma-surface interaction physics. Conventional material tests often expose a specimen simultaneously to convective heating, reactive species, pressure, shear, radiation, and surface-current effects. The resulting material response may be measured accurately, while the contribution of each mechanism remains difficult to identify. ARCTRON is designed to vary these effects selectively. Plasma chemistry can be changed independently through configurable gas mixtures; magnetic fields and electrical biasing can modify charged-particle transport; laser heating can provide a non-plasma thermal input; and pressure, flow, and discharge mode can be varied over a broad operating space. This enables controlled tests of hypotheses involving surface catalycity, reactive-species transport, plasma-assisted oxidation, electromagnetic effects, shear, and the relative contributions of thermal and chemical loading. A second class of investigations enabled by this approach concerns the engineering of high-enthalpy facilities themselves. Arc-heated facilities are limited by electrode erosion, unstable arc attachment, localized heating, and damage to nozzles and other plasma-facing components. ARCTRON provides a lower-cost environment for testing concepts intended to mitigate these limitations. Candidate investigations include the use of applied magnetic fields to alter current paths and reduce plasma interaction with nozzle walls, ExB forcing to introduce controlled plasma rotation, magnetic or geometric approaches for distributing arc attachment, and alternative electrode or discharge configurations intended to reduce erosion and improve stability. Because the platform is reconfigurable, these concepts can be evaluated through repeated design--build--test cycles before they are considered for implementation in operational facilities. The platform also supports the development and validation of diagnostics that may be difficult to introduce initially into a large arcjet. Current and planned measurements include spatially resolved optical emission spectroscopy, electrostatic probes, fast imaging, pyrometry, calorimetry, laser-induced fluorescence, and absorption spectroscopy. These diagnostics are intended not merely to document a nominal operating condition, but to constrain the local plasma state and its relationship to component or material response. The modular facility geometry allows diagnostic concepts to be tested, calibrated, and compared under repeatable conditions before deployment in more demanding environments. ARCTRON is also supported by an integrated software suite. Automated control and data acquisition allow discharge parameters, gas composition, magnetic fields, diagnostic timing, and test configuration to be recorded as part of each experiment (STARDAC - Software for Testing, Analysis, Research Data, and Control). The Backend for Experiment Analysis, Storage, and Traceability (BEAST) is a database that provides the infrastructure needed to associate heterogeneous measurements with facility configuration, specimen identity, calibration state, geometry, and analysis provenance. This backend is particularly important for exploratory campaigns, in which many related configurations may be tested, and the value of an individual experiment depends on its connection to earlier and subsequent iterations. Complementary analysis capabilities, including computer-vision-based transient response measurements (arcjetCV), three-dimensional surface reconstruction (STARSCAN), and model-based Bayesian inference (SHIELD), and tomography data analysis (TOMATO, PuMA) can be incorporated when required by a specific hypothesis without becoming the focus of every campaign. The central objective of ARCTRON is therefore not to maximize heat flux or reproduce a complete flight environment. Its purpose is to reduce the cost and time required to ask consequential questions about plasma behavior, plasma-facing materials, diagnostics, and arcjet technology. By providing a controlled environment for rapid reconfiguration, mechanism isolation, quantitative measurement, and iterative engineering, ARCTRON can help mature concepts that would otherwise remain too speculative or too risky for evaluation in production facilities. The resulting knowledge can then guide the design of material models, focus test objectives in larger arcjets, reduce facility-development risk, and improve the physical basis of high-enthalpy ground testing. This work will present the ARCTRON architecture, operating modes, diagnostic suite, and digital experimental workflow. Initial experimental results from the first integrated operation of the facility will be presented, including flow characterization, power limitations, and deployment of the initial diagnostic suite. Ongoing development efforts aimed at catalycity characterization, magnetic plasma control, and advanced optical diagnostics will also be discussed, illustrating how the platform supports rapid iteration from concept to experiment.

experimental diagnostics↗