Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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 1,027 records · Page 57

Two Strengths of Ordinary Chondritic Meteoroids as Derived from Their Atmospheric Fragmentation Modeling

The internal structure and strength of small asteroids and large meteoroids is poorly known. Observation of bright fireballs in the Earth’s atmosphere can explore meteoroid structure by studying meteoroid fragmentation during the flight. Earlier evaluations showed that the meteoroid’s strength is significantly lower than that of the recovered meteorites. We present a detailed study of atmospheric fragmentation of seven meteorite falls, all ordinary chondrites, and 14 other fireballs, where meteorite fall was predicted but the meteorites, probably also ordinary chondrites, were not recovered. All observations were made by the autonomous observatories of the European Fireball Network and include detailed radiometric light curves. A model, called the semiempirical fragmentation model, was developed to fit the light curves and decelerations. Videos showing individual fragments were available in some cases. The results demonstrated that meteoroids do not fragment randomly but in two distinct phases. The first phase typically corresponds to low strengths of 0.04–0.12 MPa. In two-thirds of cases, the first phase was catastrophic or nearly catastrophic with at least 40% of mass lost. The second phase corresponds to 0.9–5 MPa for confirmed meteorite falls and somewhat lower strengths, from about 0.5 MPa, for smaller meteoroids. All of these strengths are lower than the tensile strengths of ordinary chondritic meteorites cited in the literature, 20–40 MPa. We interpret the second phase as being due to cracks in meteoroids and the first phase as a separation of weakly cemented fragments, which reaccumulated at the surfaces of asteroids after asteroid collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

Introducing a New Spitzer Master BLISS Map to Remove the Instrument Systematic Phase-curve-parameter Degeneracy, as Demonstrated by a Reanalysis of the 4.5 μm WASP-43b Phase Curve

While Spitzer Infrared Array Camera (IRAC) systematics are generally well understood, each data set can provide its own challenges that continue to teach us about the underlying functional form of these systematics. Multiple groups have analyzed the phase curves of WASP-43b with varying detrending techniques, each obtaining different results. In this work, we take another look at WASP-43b while exploring the degenerate relation between Bilinearly Interpolated Subpixel Sensitivity (BLISS) mapping, point response function (PRF)–FWHM detrending, and phase curve parameters. We find that there is a strong correlation between the detrending parameters in the two models, and best-fit phase curve amplitudes vary strongly when the data are temporally binned. To remove this degeneracy, we present a new Gaussian centroided intrapixel sensitivity map (hereafter fixed sensitivity map), generated using 3,712,830 exposures spanning 5 yr, for a variety of aperture sizes at 4.5 μm. We find evidence for time variability in the sensitivity at 3.6 μm and do not generate a visit-independent map for this channel. With the fixed 4.5 μm intrapixel sensitivity map, the best fits for WASP-43b no longer vary strongly with bin size and PRF–FWHM detrending is no longer required to remove correlated noise. For data sets that do not fall completely within the sweet spot, temporal binning should not be used in the analysis of Spitzer phase curves. We confirm nightside emission for WASP-43b with a disk-integrated nightside temperature of 806 ± 48 K at 4.5 μm. The 4.5 μm maps are available at github.com/kevin218/POET.

79 ASTRONOMY AND ASTROPHYSICS↗

Early-time Light Curves of Type Ia Supernovae Observed with TESS

We present the early-time light curves of Type Ia supernovae (SNe Ia) observed in the first six sectors of Transiting Exoplanet Survey Satellite (TESS) data. Ten of these SNe were discovered by ASAS-SN, seven by ATLAS, six by ZTF, and one by Gaia. For nine of these objects with sufficient dynamic range (>3.0 mag from detection to peak), we fit power-law models and searched for signatures of companion stars. We found a diversity of early-time light-curve shapes, although most of our sources are consistent with fireball models where the flux increases as ∝t 2 . Three SNe displayed a flatter rise with flux ∝t. We did not find any obvious evidence for additional structures, such as multiple power-law components, in the early rising light curves. For assumptions about the SN properties and the observer viewing angle (ejecta mass of 1.4 M⊙, expansion velocity of 10 4 km s –1 , opacity of 0.2 cm 2 g –1 , and viewing angle of 45°) and a further assumption that any companion stars would be in Roche lobe overflow, it is possible to place upper limits on the radii of any companion stars. Six of the nine SNe had complete coverage of the early-time light curves, and we placed upper limits on the radii of companion stars of ≲32 R ⊙ for these SNe, ≲20 R ⊙ for five of the six, and ≲4 R ⊙ for two of the six. The small sample size did not allow us to put limits on the occurrence rate of companion stars in the progenitors of SNe Ia. However, we expect that TESS observed enough SNe in its two-year primary mission (26 sectors) to either detect the signature of a large companion (R > 20 R ⊙ ) or constrain the occurrence rate of such systems, at least for the fiducial SN properties adopted here. Here, we also show that TESS is capable of detecting emission from a 1 R ⊙ companion for an SN Ia within 50 Mpc and has a reasonable chance of doing so after about six years.

79 ASTRONOMY AND ASTROPHYSICS↗

ASASSN-14ko is a Periodic Nuclear Transient in ESO 253-G003

We present the discovery that ASASSN-14ko is a periodically flaring active galactic nucleus at the center of the galaxy ESO 253-G003. At the time of its discovery by the All-Sky Automated Survey for Supernovae (ASAS-SN), it was classified as a supernova close to the nucleus. The subsequent 6 yr of V - and g -band ASAS-SN observations revealed that ASASSN-14ko has nuclear flares occurring at regular intervals. The 17 observed outbursts show evidence of a decreasing period over time, with a mean period of P 0 = 114.2 ± 0.4 days and a period derivative of . The most recent outburst in 2020 May, which took place as predicted, exhibited spectroscopic changes during the rise and had a UV bright, blackbody spectral energy distribution similar to tidal disruption events (TDEs). The X-ray flux decreased by a factor of 4 at the beginning of the outburst and then returned to its quiescent flux after ~8 days. The Transiting Exoplanet Survey Satellite observed an outburst during Sectors 4–6, revealing a rise time of 5.60 ± 0.05 days in the optical and a decline that is best fit with an exponential model. We discuss several possible scenarios to explain ASASSN-14ko’s periodic outbursts, but currently favor a repeated partial TDE. The next outbursts should peak in the optical on UT 2020 September 7.4±1.1 and UT 2020 December 26.5±1.4.

Astronomy & Astrophysics↗

Probing Galaxy Evolution in Massive Clusters Using ACT and DES: Splashback as a Cosmic Clock

In this work, we measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev–Zel'dovich effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius is consistent with CDM-only simulations and is located at 2.4$_{-0.4}^{+0.3}$ Mpch -1 . We split the galaxies on color and find significant differences in their profile shapes. Red and green-valley galaxies show a splashback-like minimum in their slope profile consistent with theory, while the bluest galaxies show a weak feature at a smaller radius. We develop a mapping of galaxies to subhalos in simulations and assign colors based on infall time onto their hosts. We find that the shift in location of the steepest slope and different profile shapes can be mapped to the average time of infall of galaxies of different colors. The steepest slope traces a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in clusters over 3.2 Gyr, green galaxies about 2.2 Gyr, while blue galaxies have been accreted most recently and have not reached apocenter. Using the full radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching occurs after a delay of about a gigayear and that galaxies quench rapidly thereafter with an exponential timescale of 0.6 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

Probing the Extragalactic Mid-infrared Background with HAWC

The extragalactic background light (EBL) contains all the radiation emitted by nuclear and accretion processes in stars and compact objects since the epoch of recombination. Measuring the EBL density directly is challenging, especially in the near-to-far-infrared wave band, mainly due to the zodiacal light foreground. Instead, gamma-ray astronomy offers the possibility to indirectly set limits on the EBL by studying the effects of gamma-ray absorption in the very high energy (VHE: >100 GeV) spectra of distant blazars. The High Altitude Water Cherenkov Gamma Ray Observatory (HAWC) is one of the few instruments sensitive to gamma rays with energies above 10 TeV. This offers the opportunity to probe the EBL in the near/mid-IR region: λ = 1–100 μm. In this study, we fit physically motivated emission models to Fermi-LAT gigaelectronvolt data to extrapolate the intrinsic teraelectronvolt spectra of blazars. We then simulate a large number of absorbed spectra for different randomly generated EBL model shapes and calculate Bayesian credible bands in the EBL intensity space by comparing and testing the agreement between the absorbed spectra and HAWC extragalactic observations of two blazars. The resulting bands are in agreement with current EBL lower and upper limits, showing a downward trend toward higher wavelength values λ > 10 μm also observed in previous measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

Bi-Kappa Proton Mirror and Cyclotron Instabilities in the Solar Wind

Abstract The charged particles in the solar wind are often observed to possess a nonthermal tail in the velocity distribution function, a feature that can be fitted with the Kappa model. The anisotropic, or bi-Kappa, model of protons, electrons, and other charged particles is thus adopted in the literature for interpreting the data as well as in the context of the analysis of wave–particle interactions. The present paper develops an approximate but efficient theory of the mirror and cyclotron instabilities excited by the bi-Kappa protons in the solar wind. A velocity moment-based quasi-linear theory of these instabilities is also formulated in order to investigate the saturation behavior. Applications of the formalism are made for instabilities close to the marginally unstable state, which is typical of the solar wind near 1 au.

79 ASTRONOMY AND ASTROPHYSICS↗

Search for Joint Multimessenger Signals from Potential Galactic Cosmic-Ray Accelerators with HAWC and IceCube

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

79 ASTRONOMY AND ASTROPHYSICS↗

Aggressively Dissipative Dark Dwarfs: The Effects of Atomic Dark Matter on the Inner Densities of Isolated Dwarf Galaxies

We present the first suite of cosmological hydrodynamical zoom-in simulations of isolated dwarf galaxies for a dark sector that consists of cold dark matter and a strongly dissipative subcomponent. The simulations are implemented in GIZMO and include standard baryons following the FIRE-2 galaxy formation physics model. The dissipative dark matter is modeled as atomic dark matter (aDM), which forms a dark hydrogen gas that cools in direct analogy to the Standard Model. Our suite includes seven different simulations of ∼10 10 M ⊙ systems that vary over the aDM microphysics and the dwarf’s evolutionary history. We identify a region of aDM parameter space where the cooling rate is aggressive and the resulting halo density profile is universal. In this regime, the aDM gas cools rapidly at high redshifts, and only a small fraction survives in the form of a central dark gas disk; the majority collapses centrally into collisionless dark “clumps,” which are clusters of subresolution dark compact objects. These dark clumps rapidly equilibrate in the inner galaxy, resulting in an approximately isothermal distribution that can be modeled with a simple fitting function. Even when only a small fraction (∼5%) of the total dark matter is strongly dissipative, the central densities of classical dwarf galaxies can be enhanced by over an order of magnitude, providing a sharp prediction for observations.

cold dark matter↗

PARTICLE RESOLVED DNS STUDY OF TURBULENCE EFFECTS ON HYPORHEIC MIXING IN RANDOMLY PACKED SEDIMENT BEDS

Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations.

Karra, Shashank↗

Examination of Digital-Delay Rossi-$\alpha$ for 252 Cf-Driven Highly Enriched Uranium Using Organic Scintillators

Neutron noise techniques constitute several analysis methods applicable to non-destructive assay. One technique is the Rossi-α method to calculate the prompt neutron decay constant (α) or its inverse, the prompt neutron period (1/α), for assemblies of fissionable material. This work evaluates a high-data-throughput measurement at the National Criticality Experiments Research Center (NCERC) with organic scintillators measuring a subcritical assembly of highly enriched uranium (HEU) metal (93% 235 U). The assembly comprises hemi shells stacked together to form fully closed shells and is driven by a Cf source at the center. The assembly is a total of 59.85 kg HEU and a k eff of 0.99, according to MCNP6.2 KCODE simulations. Measurements were acquired with a three-by-four array of 5.08-cm-diameter by 5.08-cm-length trans-stilbene crystals 166 cm from the assembly center. Two types of coincident binning methods are used to build the Rossi-α distribution of coincident neutron detections: A) type 1 binning, also known as any-and-all forward time differences, and B) type 1 binning with a digital-delay technique that is analogous to Orndoff’s use of delay cabling in early Rossi-alpha measurements. Method B) disregards same-detector coincidences and once all nearest time coincidences are collected between all detectors after a trigger, a time delay of 0.75 µs is implemented. The measured prompt neutron decay constants for both techniques are calculated from single exponential fits and the two methods are compared. Method A) shows an apparent timing discontinuity near 500 ns time differences. This timing discontinuity interferes with the fitting method used to calculate the prompt neutron decay constant. Method B) mitigates the artificial timing discontinuity and removes the disagreement of the fit. Future work will model the time-dependent detector response to discern why this timing discontinuity occurs, discern how Method B) reduces the amplitude of this timing discontinuity, and apply this method to critical assembly measurements to ultimately confirm the recommendation to use Method B) for high-data-throughput measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental Considerations for Estimating Degradation in PV Modules

Carefully controlled laboratory experiments and measurements can enable the determination of acceleration factors suitable for extrapolation to durability and performance of a fielded PV module. Ideally, a single mechanism can be identified with appropriate acceleration factors for extrapolation to the field. However, even with a single mechanism, the inherent uncertainty in these factors leads to uncertainty in the extrapolation which is greater the higher the acceleration factor. This course will explain how because of the wide range of acceleration factors for a given degradation mode, utilizing acceleration factors greater than about 10x will typically lead to unacceptable uncertainty in the results. Therefore, if even just a rank ordering of materials is desired, acceleration factors must be minimized which requires a good general understanding of the scale of the different acceleration factors for the degradation mode of interest. In this tutorial we will discuss what the different purposes are for many of the accelerated stress tests used today. E.g., what is a qualification test, a highly accelerated stress test, a rank ordering test, or a service life prediction test. We will discuss how one can understand the relationship between test results and expected field performance. A single accelerated stress test condition cannot duplicate outdoor exposure for all possible degradation pathways; therefore, one must use targeted evaluation of material properties at different stress levels to determine the relevant acceleration factors and fit it to a model. We will also discuss how to interpret the results of experiments understanding what is relevant/not relevant, or not e valuated in a test. There are many common error people make in their test interpretations because they push the stress levels to be too harsh. This creates biases and can mask the relevant failure modes and mechanisms or will erroneously lead one to over design materials against things that aren't relevant. Several case studies will be presented to illustrate appropriate interpretation of accelerated stress testing results.

degradation↗

Electroluminescence Yield Measurements in Xenon Gas with the NEXT-DEMO++ Detector

The NEXT-DEMO++ detector, a high-pressure xenon gas time projection chamber serving as a prototype for the NEXT-100 experiment, was used to measure the electroluminescence (EL) yield as a function of reduced electric field ($E/p$) across pressures from 2.0 to 9.4 bar, utilizing the 41.5 keV de-excitation peak of $^{83m}$Kr. These measurements were made to examine the pressure dependence of the slope of the reduced EL yield $Y/p$, which has shown inconsistencies in the literature. The reduced yield was fitted with a linear model, revealing a modest ($\sim$5%) change in slope, beginning around 5 bar and increasing with pressure up to 9.4 bar.

Renner, J. [Valencia U., IFIC]↗

Benefits of colorectal cancer screening using fecal immunochemical testing with varying positivity thresholds by age and sex

Background Fecal immunochemical test (FIT) performance for colorectal cancer screening varies by age and sex, yet most FIT-based screening programs use uniform positivity thresholds. This study assessed the potential benefits of stratifying FIT thresholds based on age and sex. Methods We conducted a meta-analysis of FIT sensitivity and specificity at various positivity thresholds by age and sex. We then used these estimates in 2 microsimulation models of colorectal cancer and projected lifetime clinical outcomes, incremental costs, and quality-adjusted life-years (QALYs) gained from age- and sex-stratified FIT strategies. FIT thresholds ranged from 10 to 50 µg hemoglobin per gram of feces. Results For current uniform FIT screening (20 µg hemoglobin/gram of feces), models projected 85.67 to 122.15 QALYs gained at incremental costs of ‒$982 to $504 per 1000 individuals compared with no screening. At equivalent costs to current uniform screening, only 1 model found stratified FIT approaches cost-effective, yielding a marginal increase of 1.04 and 1.10 QALYs gained/1000 female and male individuals, respectively. At a willingness-to-pay threshold of $100 000/QALYs gained, both models found stratified FIT cutoffs to be the best strategy, with cutoffs being equal to or higher for males and lowest at older ages (70-75 years). Uniform strategies showed comparable effectiveness, falling within 1 quality-adjusted life-day per person of efficient strategies at up to $112 more per person. Results were sensitive to FIT test performance characteristics and 1-time setup costs. Conclusion Stratifying FIT thresholds by age and sex may be cost-effective compared to current screening. The gain in expected health benefits with stratified FIT screening, however, is likely small.

Harlass, Matthias↗

Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations

Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.

36 MATERIALS SCIENCE↗

Measurements of the E-mode polarization and temperature-E-mode correlation of the CMB from SPT-3G 2018 data

We present measurements of the E -mode ( E E ) polarization power spectrum and temperature- E -mode ( T E ) cross-power spectrum of the cosmic microwave background using data collected by SPT-3G, the latest instrument installed on the South Pole Telescope. This analysis uses observations of a 1500 deg 2 region at 95, 150, and 220 GHz taken over a four-month period in 2018. We report binned values of the E E and T E power spectra over the angular multipole range 300 ≤ ℓ < 3000 , using the multifrequency data to construct six semi-independent estimates of each power spectrum and their minimum-variance combination. These measurements improve upon the previous results of SPTpol across the multipole ranges 300 ≤ ℓ ≤ 1400 for E E and 300 ≤ ℓ ≤ 1700 for T E , resulting in constraints on cosmological parameters comparable to those from other current leading ground-based experiments. We find that the SPT-3G data set is well fit by a Λ CDM cosmological model with parameter constraints consistent with those from Planck and SPTpol data. From SPT-3G data alone, we find H 0 = 68.8 ± 1.5 km s - 1 Mpc - 1 and σ 8 = 0.789 ± 0.016 , with a gravitational lensing amplitude consistent with the Λ CDM prediction ( A L = 0.98 ± 0.12 ). We combine the SPT-3G and the Planck data sets and obtain joint constraints on the Λ CDM model. The volume of the 68% confidence region in six-dimensional Λ CDM parameter space is reduced by a factor of 1.5 compared to Planck-only constraints, with no significant shifts in central values. We note that the results presented here are obtained from data collected during just half of a typical observing season with only part of the focal plane operable, and that the active detector count has since nearly doubled for observations made with SPT-3G after 2018.

79 ASTRONOMY AND ASTROPHYSICS↗

From vertices to vortices in magnetic nanoislands

We report recent studies in magnetic nanoarrays show that a variety of complex magnetic states and textures emerge as a function of a single magnetic nanoisland's aspect ratio. We propose a model that, in addition to fitting experiments, predicts magnetic states with continuous symmetry at particular aspect ratios and reveals a duality between vortex and vertex states. Our model opens new means of engineering novel types of artificial spin systems, and their application to complex magnetic textures in devices and computing.

97 MATHEMATICS AND COMPUTING↗

Code for multi-shape Gaussian process (GP) fitting with uncertainty quantification (UQ)

This is code associated with the publication “Nonparametric Multi-shape Modeling with Uncertainty Quantification,” authored by Hengrui Luo (Lawrence Berkeley National Laboratory) and Justin Strait (Los Alamos National Laboratory). The code is used to fit multiple-output Gaussian process (GP) models of planar closed curves to collections of ordered point sets, allowing for flexible nonlinear prediction of the underlying curve under dense or sparse point set samplings and with or without noise, as well as tractable uncertainty quantification. To do this, we employ use of a periodic kernel to account for the nonlinear input space of closed curves, and combine with coregionalization models to account for dependence both (i) between curve coordinates, and (ii) between pairs of curves. Functions in the code are capable of fitting these models, as well as performing additional tasks with the fitted curves such as (a) shape registration and alignment, (b) shape averaging, and (c) fitting for curve sub-populations / clusters.

Strait, Justin↗