Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling variability”

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,477 records · Page 82

Phase-field modeling and experiments of dynamic fracture in single crystal quartz

Predicting the onset and characteristics of brittle fracture is important for a wide range of engineering and geological material applications. In this paper, we study important aspects of brittle fracture in α-quartz by phase-field modeling and experiments using a top-down approach. In the modeling framework, the work term in the Griffith energy balance is replaced with internal energy contributions that represent surface energy, thermal energy, and elastic strain energy stored in defects. This allows parametrization of individual energy contributions in terms of internal state variables and keeps track of energy partitioning after the onset of fracture. The path and history dependence of fracture is included in evolution laws for internal state variables, e.g., entropy evolution, while the energy remains a true potential. In the experimental part, dynamic compression experiments coupled with X-ray phase contrast imaging are performed on cube-like samples with a hole. In the top-down analysis, dynamic compression and three point bending experiments from the literature are simulated with the developed phase-field damage model. In conclusion, the fitted model highlights the strain rate, size, and stress state dependence of damage nucleation and evolution in single crystal α-quartz.

36 MATERIALS SCIENCE↗

Structural Conservation of the A 1 Binding Site in Photosystem I across Cyanobacteria and Green Algae

Time-resolved step-scan Fourier transform infrared (FTIR) difference spectroscopy was used to obtain (A 1 − − A 1 ) FTIR difference spectra from photosystem I (PSI) samples isolated from eight phylogenetically diverse cyanobacterial strains and one green alga, totaling 13 PSI preparations. These included samples from cells grown under farred light and PSI in monomeric, dimeric, trimeric, and tetrameric states. Spectral profiles were shown to be independent of oligomeric state. Remarkably, all (A 1 − − A 1 ) FTIR difference spectra exhibited high similarity, underscoring the robustness of the technique and indicating minimal experimental variability. This congruence reveals a highly conserved environment for the phylloquinone cofactor at the A 1 binding site across diverse taxa. Conserved bands associated with the A 0 pigment further suggest structural continuity from A 0 to A 1 . To leverage this consistency, we constructed a composite (A 1 − − A 1 ) FTIR difference spectrum by averaging all 13 spectra. This composite spectrum provides enhanced resolution, enabling unambiguous identification of previously unresolved bands. The fact that a highly resolved composite spectrum can be obtained by averaging demonstrates the similarity in the spectra from the different types of samples. Band assignments were refined using prior studies, yielding an improved spectral framework for future investigations of PSI electron transfer cofactors.

Charge transfer↗

Low Thermal Conductivity in Heteroanionic Materials with Layers of Homoleptic Polyhedra

Although BiAgOSe, an analogue of a well-studied thermoelectric material BiCuOSe, is thermodynamically stable, its synthesis is complicated by the low driving force of formation from the stable binary and ternary intermediates. Here we have developed a “subtraction strategy” to suppress byproducts and produce pure phase BiAgOSe using hydrothermal methods. Electronic structure calculations and optical characterization show that BiAgOSe is an indirect bandgap semiconductor with a bandgap of 0.95 eV. The prepared sample exhibits lower lattice thermal conductivities (0.61 W·m –1 ·K –1 at room temperature and 0.35 W·m –1 ·K –1 at 650 K) than BiCuOSe. Further, lattice dynamical simulations and variable temperature diffraction measurements demonstrate that the low lattice thermal conductivity arises from both the low sound velocity and high phonon–phonon scattering rates in BiAgOSe. These in turn result primarily from the soft Ag–Se bonds in the edge-sharing AgSe 4 tetrahedra and large sublattice mismatch between the quasi-two-dimensional [Bi 2 O 2 ] 2+ and [Ag 2 Se 2 ] 2– layers. These results highlight the advantages of manipulating the chemistry of homoleptic polyhedra in heteroanionic compounds for electronic structure and phonon transport control.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mixed Element Microparticle Characterization by Electron Probe Microanalysis

Abstract Microparticle compositional characterization for nuclear forensics has traditionally been achieved by “gold-standard” destructive analytical techniques such as large geometry-secondary ion mass spectrometry and fission track-thermal ionization mass spectrometry, but long processing times and total sample consumption eliminate the option of subsequent analysis by other methods. Electron probe microanalysis has long been widely employed for rapid, nondestructive, micrometer-scale compositional measurements and imaging in fields such as material science and geology and may therefore need to be reevaluated as an alternative tool for microparticle analysis for the nuclear forensics community. This study presents the use of electron probe microanalysis for imaging and quantitative characterization of homogeneous microparticles utilizing a calibration curve based on high-precision quadrupole-inductively coupled plasma-mass spectrometry analyses. Samples of opportunity synthesized at Savannah River National Laboratory, nickel-doped cerium oxide microparticles, were selected as analogs for plutonium-doped uranium oxide microparticles. Positive detection and accurate quantification of variable amounts of nickel dopant down to trace levels (10s of parts per million) suggest applicability of this technique to other mixed element systems (e.g., actinides). Quantitative single-particle characterization via electron probe microanalyzer may thus provide a high fidelity, nondestructive complement to techniques currently in use.

Riche, Alexis T. (ORCID:0009000121978126)↗

The SAMI Galaxy Survey: stellar population and structural trends across the Fundamental Plane

ABSTRACT We study the Fundamental Plane (FP) for a volume- and luminosity-limited sample of 560 early-type galaxies from the SAMI survey. Using r-band sizes and luminosities from new multi-Gaussian expansion photometric measurements, and treating luminosity as the dependent variable, the FP has coefficients a = 1.294 ± 0.039, b = 0.912 ± 0.025, and zero-point c = 7.067 ± 0.078. We leverage the high signal-to-noise ratio of SAMI integral field spectroscopy, to determine how structural and stellar population observables affect the scatter about the FP. The FP residuals correlate most strongly (8σ significance) with luminosity-weighted simple stellar population (SSP) age. In contrast, the structural observables surface mass density, rotation-to-dispersion ratio, Sérsic index, and projected shape all show little or no significant correlation. We connect the FP residuals to the empirical relation between age (or stellar mass-to-light ratio Υ⋆ ) and surface mass density, the best predictor of SSP age amongst parameters based on FP observables. We show that the FP residuals (anti)correlate with the residuals of the relation between surface density and Υ⋆ . This correlation implies that part of the FP scatter is due to the broad age and Υ⋆ distribution at any given surface mass density. Using virial mass and Υ⋆, we construct a simulated FP and compare it to the observed FP. We find that, while the empirical relations between observed stellar population relations and FP observables are responsible for most (75 per cent) of the FP scatter, on their own they do not explain the observed tilt of the FP away from the virial plane.

D’Eugenio, Francesco↗

Differentiating spatial from environmental effects on foliar fungal communities of Populus trichocarpa

Abstract Aim Foliar fungi – pathogens, endophytes, epiphytes – form taxonomically diverse communities that affect plant health and productivity. The composition of foliar fungal communities is variable at spatial scales both small (e.g. individual plants) and large (e.g. continents), yet few studies have attempted to tease apart spatial from climatic factors influencing these communities. Moreover, few studies have sampled in more than 1 year to gauge interannual variation in community structure. Location The Pacific Northwest of western North America. Taxon Foliar fungi associated with the deciduous tree Populus trichocarpa . Methods In two consecutive years, we used DNA metabarcoding to characterize foliar fungal communities of Populus trichocarpa across its geographic range, which encompasses a sharp climatic transition as it crosses the Cascade Mountain Range. We used multivariate analyses to (a) test for and differentiate spatial and environmental factors affecting community composition and (b) test for temporal variation in community composition across spatial and environmental gradients. Results In both study years, we found that foliar fungal community composition varied among sites and between regions (east vs. west of the Cascades). We found that climate explained more variation in community composition than geographic distance, although the majority of variation explained by each was shared. We also found that interannual variation in community composition depended on environmental context: communities located in the dry, eastern portion of the tree's geographic range varied more between study years than those located in the wet, western portion of the tree's range. Main conclusions Our results suggest that the environment plays a greater role in structuring foliar fungal communities than dispersal limitation.

Barge, Edward G.↗

Impacts of the Additive Manufacturing Process on the Roughness of Engine Scale Vanes and Cooling Channels

Abstract By leveraging the additive manufacturing (AM) platform, development time and costs for turbine component testing can be reduced relative to traditional investment casting. Surface roughness is a key characteristic of the additive manufacturing process that can impact flow, heat transfer, and mechanical integrity of printed components. There are multiple design and build considerations that result in variability in surface roughness, especially when additively fabricating complicated three-dimensional vanes and internal cooling passages. This study characterizes the surface roughness of internal cooling passages, vanes, and flat external surface samples made using additive manufacturing, specifically the direct metal laser sintering process. The samples were manufactured with various wall thicknesses, layer thicknesses, build locations, build directions, and on different AM machines. A combination of computed tomography scanning and optical profilometry was used to evaluate surface roughness levels. The data indicate that the dominant factors in roughness for a given layer thickness are a function of wall thickness, build location, and build direction.

Engineering↗

Personalized, disease-stage specific, rapid identification of immunosuppression in sepsis

Introduction Data overlapping of different biological conditions prevents personalized medical decision-making. For example, when the neutrophil percentages of surviving septic patients overlap with those of non-survivors, no individualized assessment is possible. To ameliorate this problem, an immunological method was explored in the context of sepsis. Methods Blood leukocyte counts and relative percentages as well as the serum concentration of several proteins were investigated with 4072 longitudinal samples collected from 331 hospitalized patients classified as septic (n=286), non-septic (n=43), or not assigned (n=2). Two methodological approaches were evaluated: (i) a reductionist alternative, which analyzed variables in isolation; and (ii) a non-reductionist version, which examined interactions among six (leukocyte-, bacterial-, temporal-, personalized-, population-, and outcome-related) dimensions. Results The reductionist approach did not distinguish outcomes: the leukocyte and serum protein data of survivors and non-survivors overlapped. In contrast, the non-reductionist alternative differentiated several data groups, of which at least one was only composed of survivors (a finding observable since hospitalization day 1). Hence, the non-reductionist approach promoted personalized medical practices: every patient classified within a subset associated with 100% survival subset was likely to survive. The non-reductionist method also revealed five inflammatory or disease-related stages (provisionally named ‘early inflammation, early immunocompetence, intermediary immuno-suppression, late immuno-suppression, or other’). Mortality data validated these labels: both ‘suppression’ subsets revealed 100% mortality, the ‘immunocompetence’ group exhibited 100% survival, while the remaining sets reported two-digit mortality percentages. While the ‘intermediary’ suppression expressed an impaired monocyte-related function, the ‘late’ suppression displayed renal-related dysfunctions, as indicated by high concentrations of urea and creatinine. Discussion The data-driven differentiation of five data groups may foster early and non-overlapping biomedical decision-making, both upon admission and throughout their hospitalization. This approach could evaluate therapies, at personalized level, earlier. To ascertain repeatability and investigate the dynamics of the ‘other’ group, additional studies are recommended.

Immunology↗

Geological and geochemical controls on subsurface microbial life in the Samail Ophiolite, Oman

Microbial abundance and diversity in deep subsurface environments is dependent upon the availability of energy and carbon. However, supplies of oxidants and reductants capable of sustaining life within mafic and ultramafic continental aquifers undergoing low-temperature water-rock reaction are relatively unknown. We conducted an extensive analysis of the geochemistry and microbial communities recovered from fluids sampled from boreholes hosted in peridotite and gabbro in the Tayin block of the Samail Ophiolite in the Sultanate of Oman. The geochemical compositions of subsurface fluids in the ophiolite are highly variable, reflecting differences in host rock composition and the extent of fluid-rock interaction. Principal component analysis of fluid geochemistry and geologic context indicate the presence of at least four fluid types in the Samail Ophiolite (“gabbro,” “alkaline peridotite,” “hyperalkaline peridotite,” and “gabbro/peridotite contact”) that vary strongly in pH and the concentrations of H 2 , CH 4 , Ca 2+ , Mg 2+ , NO 3 $-$ , SO$2-\atop{4}$, trace metals, and DIC. Geochemistry of fluids is strongly correlated with microbial community composition; similar microbial assemblages group according to fluid type. Hyperalkaline fluids exhibit low diversity and are dominated by taxa related to the Deinococcus-Thermus genus Meiothermus , candidate phyla OP1, and the family Thermodesulfovibrionaceae. Gabbro- and alkaline peridotite- aquifers harbor more diverse communities and contain abundant microbial taxa affiliated with Nitrospira , Nitrosospharaceae, OP3, Parvarcheota, and OP1 order Acetothermales. Wells that sit at the contact between gabbro and peridotite host microbial communities distinct from all other fluid types, with an enrichment in betaproteobacterial taxa. Together the taxonomic information and geochemical data suggest that several metabolisms may be operative in subsurface fluids, including methanogenesis, acetogenesis, and fermentation, as well as the oxidation of methane, hydrogen and small molecular weight organic acids utilizing nitrate and sulfate as electron acceptors. Dynamic nitrogen cycling may be especially prevalent in gabbro and alkaline peridotite fluids. As a result, these data suggest water-rock reaction, as controlled by lithology and hydrogeology, constrains the distribution of life in terrestrial ophiolites.

58 GEOSCIENCES↗

Kinetic features of xylan de-polymerization in production of xylose monomer and furfural during acid pretreatment for kenaf, forage sorghums and sunn hemp feedstocks

A kinetic study of acid pretreatment was conducted for sorghum non-brown mid rib (SNBMR) ( Sorghum bicolor L Moench ), sorghum-brown mid rib (SBMR), sunn hemp ( Crotalaria juncea L ) and kenaf ( Gossypiumhirsutum L ), focusing on rates of xylose monomer and furfural formation. The kinetics was investigated using two independent variables, reaction temperature (150 and 160°C) and acid concentration (1 and 2 wt%), with a constant dry biomass loading of 10 wt% and a treatment time up to 20 min while sampling the mixture every 2 min. The experimental data were fitted using a two-step kinetic model based on irreversible pseudo first order kinetics at each step. Varied kinetic orders on the acid concentration, ranging from 0.2 to >3, were observed for both xylose and furfural formation, the values depending on the feedstock. The crystallinity index of raw biomass was shown to be a major factor influencing the rate of both xylose and furfural formation. As a result, a positive correlation was observed between the activation energy and biomass crystallinity index for xylose formation.

09 BIOMASS FUELS↗

Evaluating Entrainment–Mixing Characteristics through Direct Comparisons of Drop Size Distributions Using In Situ Observations from ACE-ENA

Abstract Constraining the impacts of entrainment and associated mixing (i.e., entrainment–mixing) on cloud properties continues to be difficult, partly due to observational uncertainties as well as a lacking consensus of which methodologies for diagnosing entrainment–mixing are most appropriate. This study introduces a novel method to evaluate the presence and degree of inhomogeneous and homogeneous mixing using ∼100 h of in situ observations from a research aircraft over the northeastern Atlantic. Specifically, drop size distributions are compared between regions containing negligible and significant entrainment for select flight legs, making a direct characterization of the degree of homogeneous and inhomogeneous mixing possible. A measure of drop concentration variance is used as a proxy variable to diagnose entrainment–mixing. Results correspond well with entrainment–mixing metrics, showing lower Damköhler numbers where drop size distributions shift toward smaller drop sizes (i.e., inhomogeneous mixing) and greater transition length scales where drop size distributions do not (i.e., homogeneous mixing). Inhomogeneous mixing occurs in most samples from Aerosol and Cloud Experiment in the Eastern North Atlantic (ACE-ENA) (regardless of homogeneous mixing frequencies increasing with increasing spatial resolution from ∼100 to ∼10 m) and is associated with decreased drop size relative dispersion and both greater aerosol and drop concentrations compared with homogeneous mixing. Precipitating clouds have a greater frequency of homogeneous mixing compared with nonprecipitating clouds. The proposed methodology is similarly applied to in situ observations of southeast Pacific stratocumulus, shallow convective clouds over central Oklahoma and low-level clouds over the Southern Ocean. All four locations are primarily dominated by inhomogeneous mixing with minimal variability among each region.

Clouds↗

Temporal Variability of Aerosol Properties during TCAP: Impact on Radiative Forcing

Ground-based remote sensing and in situ observations of aerosol microphysical and optical properties have been collected during summertime (June-August, 2012) as part of the Two-Column Aerosol Project (TCAP; http://campaign.arm.gov/tcap/), which was supported by the U.S. Department of Energy’s (DOE’s) Atmospheric Radiation Measurement (ARM) Program (http://www.arm.gov/). The overall goal of the TCAP field campaign is to study the evolution of optical and microphysical properties of atmospheric aerosol transported from North America to the Atlantic and their impact on the radiation energy budget. During TCAP, the ground-based ARM Mobile Facility (AMF) was deployed on Cape Cod, an arm-shaped peninsula situated on the easternmost portion of Massachusetts (along the east coast of the United States) and that is generally downwind of large metropolitan areas. The AMF site was equipped with numerous instruments for sampling aerosol, cloud and radiative properties, including a Multi-Filter Rotating Shadowband Radiometer (MFRSR), a Scanning Mobility Particle Sizer (SMPS), an Aerodynamic Particle Sizer (APS), and a three-wavelength nephelometer. In this study we present an analysis of diurnal and day-to-day variability of the column and near-surface aerosol properties obtained from remote sensing (MFRSR data) and ground-based in situ measurements (SMPS, APS, and nephelometer data). In particular, we show that the observed diurnal variability of the MFRSR aerosol optical depth is strong and comparable with that obtained previously from the AERONET climatology in Mexico City, which has a larger aerosol loading. Moreover, we illustrate how the variability of aerosol properties impacts the direct aerosol radiative forcing at different time scales.

aerosol microphysical and optical properties↗

Spatial and Temporal Variability of Vertical Velocity under Shallow Cumulus

Vertical velocity distribution below cloud is one of the key determinants of cloud life cycle, but observations of this variable are extremely sparse in space. Doppler lidar retrievals and large-eddy simulations at the U.S. Department of Energy’s Atmospheric Radiation Measurement User Facility Southern Great Plains site are used to determine whether vertical velocity statistics from temporally dense profiles at a single location can be substituted for spatial vertical velocity statistics. We show that even a small number (five) of widely distributed [ O (1°) latitude/longitude spacing] lidars is sufficient sampling to reconstruct domainwide spatial vertical velocity variance, but not higher moments of the vertical velocity distribution. Spatial and temporal vertical velocity variances in the Doppler lidar observations are nearly interchangeable as long as the spatial variance is temporally averaged and the temporal variance is averaged across lidars. This is true even though the dominant spatial scales of vertical velocity variability are ≲ 3 km, more than an order of magnitude smaller than the spacing between the lidars. Further, in the limit where the temporal variance does not vary across a spatial domain (e.g., if the meteorological and surface forcing of the atmospheric turbulence is homogeneous across the domain) and the domain-mean vertical velocity is zero, the commonly available retrieval of temporal vertical velocity variance at one site is equivalent to the spatial variance over the domain. We use an updraft parcel model to show that substituting temporal for spatial vertical velocity statistics will have a relatively minor effect on cloud droplet number concentrations.

54 ENVIRONMENTAL SCIENCES↗

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗

Oxidation of Organics in the Defense Waste Processing Facility Condensate

Defense Waste Processing Facility (DWPF): The state of South Carolina currently holds 35 million gallons of radioactive waste in 51 underground tanks. DWPF was built in an effort to convert the waste into a more stable form. It turns the sludge waste from the tanks into glass, which will keep the waste stable for 10,000 years. Currently, DWPF uses formic acid in its flowsheet and there is one major problem associated with its use: it decomposes into hydrogen, a flammable gas. Glycolic acid, however, has been found as a possible alternative that performs similarly to formic acid but does not directly form hydrogen. One of the main issues with the glycolic acid flowsheet is the potential for there to be glycolate and other organics in DWPF's Recycle Collection Tank (RCT). These organics could lead to hydrogen generation later in the downstream evaporators and must be destroyed in the RCT in order for glycolic acid to be a suitable replacement for formic acid. Oxidants can be added to the RCT to destroy the organics. Two oxidizers, permanganate and Fenton's reagent (iron catalyzed hydrogen peroxide), are possible candidates. Find optimum processing conditions for the oxidation of organic species in a typical RCT simulant. A reactor served as a small-scale version of the RCT. During each run, an RCT simulant with 8 added organics, including glycolate, was created in the reactor. Either permanganate or peroxide was slowly added for each run through a syringe pump. Each experiment lasts 6 hours and a total of 23 runs are scheduled. Variables of Interest: Oxidizer (Permanganate or Peroxide); Temperature (15 deg. C to 50 deg. C); Oxidant addition time (10 mins to 2 hours); Amount of oxidant (1 x to 1.5 x stoichiometric requirements); Solution pH (3 to 13); Glycolate concentration (125 to 250 mg/kg). Offgas Analyzers: Mass spectrometer and FTIR spectrometer. Samples: pulled at 7 different times throughout the run to determine the amount of organic and inorganic carbon and concentrations of various organics through ion chromatography (IC). Only 3 of the runs have been completed at this point, all with permanganate. Based on the first 3 runs, permanganate seems to oxide the majority of the organics and the run at pH 8 and 32.5 deg. C had the most effective oxidation. The rest of the runs will need to be performed to determine the best oxidant and optimum values of the other variables of interest. Once the experiment is completed, DWPF will select one of the oxidants. Testing will continue to determine its performance during foam-over and other non-ideal conditions. Finally, the oxidant will be tested with actual radioactive waste to confirm its effectiveness.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Understanding the Significance of Biochemistry in the Storage, Handling, Purification, and Sampling of Amphiphilic Mycolactone

Mycolactone, the amphiphilic macrolide toxin secreted by Mycobacterium ulcerans , plays a significant role in the pathology and manifestations of Buruli ulcer (BU). Consequently, it follows that the toxin is a suitable target for the development of diagnostics and therapeutics for this disease. Yet, several challenges have deterred such development. For one, the lipophilic nature of the toxin makes it difficult to handle and store and contributes to variability associated with laboratory experimentation and purification yields. In this manuscript, we have attempted to incorporate our understanding of the lipophilicity of mycolactone in order to define the optimal methods for the storage, handling, and purification of this toxin. We present a systematic correlation of variability associated with measurement techniques (thin-layer chromatography (TLC), mass spectrometry (MS), and UV-Vis spectrometry), storage conditions, choice of solvents, as well as the impact of each of these on toxin function as assessed by cellular cytotoxicity. We also compared natural mycolactone extracted from bacterial culture with synthesized toxins in laboratory (solvents, buffers) and physiologically relevant (serum) matrices. Our results point to the greater stability of mycolactone in organic, as well as detergent-containing, solvents, regardless of the container material (plastic, glass, or silanized tubes). They also highlight the presence of toxin in samples that may be undetectable by any one technique, suggesting that each detection approach captures different configurations of the molecule with varying specificity and sensitivity. Most importantly, our results demonstrate for the very first time that amphiphilic mycolactone associates with host lipoproteins in serum, and that this association will likely impact our ability to study, diagnose, and treat Buruli ulcers in patients.

59 BASIC BIOLOGICAL SCIENCES↗

Prospective Validation of a High Dimensional Shape Model for Organ Motion in Intact Cervical Cancer

Purpose: Validated models are needed to justify strategies to define planning target volumes (PTVs) for intact cervical cancer used in clinical practice. Our objective was to independently validate a previously published shape model, using data collected prospectively from clinical trials. Methods and Materials: We analyzed 42 patients with intact cervical cancer treated with daily fractionated pelvic intensity modulated radiation therapy and concurrent chemotherapy in one of 2 prospective clinical trials. We collected online cone beam computed tomography (CBCT) scans before each fraction. Clinical target volume (CTV) structures from the planning computed tomography scan were cast onto each CBCT scan after rigid registration and manually redrawn to account for organ motion and deformation. We applied the 95% isodose cloud from the planning computed tomography scan to each CBCT scan and computed any CTV outside the 95% isodose cloud. The primary aim was to determine the proportion of CTVs that were encompassed within the 95% isodose volume. A 1-sample t test was used to test the hypothesis that the probability of complete coverage was different from 95%. We used mixed-effects logistic regression to assess effects of time and patient variability. Results: The 95% isodose line completely encompassed 92.3% of all CTVs (95% confidence interval, 88.3%-96.4%), not significantly different from the 95% probability anticipated a priori (P=.19). The overall proportion of missed CTVs was small: the grand mean of covered CTVs was 99.9%, and 95.2% of misses were located in the anterior body of the uterus. Time did not affect coverage probability (P=.71). Conclusions: With the clinical implementation of a previously proposed PTV definition strategy based on a shape model for intact cervical cancer, the probability of CTV coverage was high and the volume of CTV missed was low. This PTV expansion strategy is acceptable for clinical trials and practice; however, we recommend daily image guidance to avoid systematic large misses in select patients.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗