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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

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

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗
DOE OSTI2015

Automated contouring error detection based on supervised geometric attribute distribution models for radiation therapy: A general strategy

Purpose: One of the most critical steps in radiation therapy treatment is accurate tumor and critical organ-at-risk (OAR) contouring. Both manual and automated contouring processes are prone to errors and to a large degree of inter- and intraobserver variability. These are often due to the limitations of imaging techniques in visualizing human anatomy as well as to inherent anatomical variability among individuals. Physicians/physicists have to reverify all the radiation therapy contours of every patient before using them for treatment planning, which is tedious, laborious, and still not an error-free process. In this study, the authors developed a general strategy based on novel geometric attribute distribution (GAD) models to automatically detect radiation therapy OAR contouring errors and facilitate the current clinical workflow. Methods: Considering the radiation therapy structures’ geometric attributes (centroid, volume, and shape), the spatial relationship of neighboring structures, as well as anatomical similarity of individual contours among patients, the authors established GAD models to characterize the interstructural centroid and volume variations, and the intrastructural shape variations of each individual structure. The GAD models are scalable and deformable, and constrained by their respective principal attribute variations calculated from training sets with verified OAR contours. A new iterative weighted GAD model-fitting algorithm was developed for contouring error detection. Receiver operating characteristic (ROC) analysis was employed in a unique way to optimize the model parameters to satisfy clinical requirements. A total of forty-four head-and-neck patient cases, each of which includes nine critical OAR contours, were utilized to demonstrate the proposed strategy. Twenty-nine out of these forty-four patient cases were utilized to train the inter- and intrastructural GAD models. These training data and the remaining fifteen testing data sets were separately employed to test the effectiveness of the proposed contouring error detection strategy. Results: An evaluation tool was implemented to illustrate how the proposed strategy automatically detects the radiation therapy contouring errors for a given patient and provides 3D graphical visualization of error detection results as well. The contouring error detection results were achieved with an average sensitivity of 0.954/0.906 and an average specificity of 0.901/0.909 on the centroid/volume related contouring errors of all the tested samples. As for the detection results on structural shape related contouring errors, an average sensitivity of 0.816 and an average specificity of 0.94 on all the tested samples were obtained. The promising results indicated the feasibility of the proposed strategy for the detection of contouring errors with low false detection rate. Conclusions: The proposed strategy can reliably identify contouring errors based upon inter- and intrastructural constraints derived from clinically approved contours. It holds great potential for improving the radiation therapy workflow. ROC and box plot analyses allow for analytically tuning of the system parameters to satisfy clinical requirements. Future work will focus on the improvement of strategy reliability by utilizing more training sets and additional geometric attribute constraints.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

THE SPECTRAL ENERGY DISTRIBUTION OF THE COLDEST KNOWN BROWN DWARF

WISE J085510.83–071442.5 (hereafter WISE 0855–0714) is the coldest known brown dwarf (∼250 K) and the fourth-closest known system to the Sun (2.2 pc). It has been previously detected only in the J band and two mid-IR bands. To better measure its spectral energy distribution (SED), we have performed deep imaging of WISE 0855–0714 in six optical and near-IR bands with Gemini Observatory, the Very Large Telescope, and the Hubble Space Telescope . Five of the bands show detections, although one detection is marginal (S/N ∼ 3). We also have obtained two epochs of images with the Spitzer Space Telescope for use in refining the parallax of the brown dwarf. By combining astrometry from this work and previous studies, we have derived a parallax of 0.449 ± 0.008″ (2.23 ± 0.04 pc). We have compared our photometry for WISE 0855–0714 to data for known Y dwarfs and to the predictions of three suites of models by Saumon et al. and Morley et al. that are defined by the presence or absence of clouds and nonequilibrium chemistry. Our estimates of Y − J and J − H for WISE 0855–0714 are redder than colors of other Y dwarfs, confirming a predicted reversal of near-IR colors to redder values at temperatures below 300–400 K. In color–magnitude diagrams, no single suite of models provides a clearly superior match to the sequence formed by WISE 0855–0714 and other Y dwarfs. Instead, the best-fitting model changes from one diagram to the next. Similarly, all of the models have substantial differences from the SED of WISE 0855–0714. As a result, we are currently unable to constrain the presence of clouds or nonequilibrium chemistry in its atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyα Forest

The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyα forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit–accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible IGM thermal histories by varying the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen at redshifts 2.2 ≤ z ≤ 5.0 and with the helium Lyα opacity at redshifts 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyα forest data to break model degeneracies, we find that the IGM experienced two main reheating events over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z R ≈ 6.0 with a first IGM temperature peak of T 0 ≃ 1.3 × 10 4 K and is followed by the reionization of He ii that completes by z R ≈ 3.0 and yields a second temperature peak of T 0 ≃ 1.4 × 10 4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

A precessing relativistic jet model for 3C 449

It is shown that the radio structure of 3C 449 can be matched with a model in which the jets are precessing and have relativistic (beta greater-than or equal to 0.4) velocities. The best-fit model implies a precession period of about 100,000 yr and a cone angle which increases with time. A similar model may be relevant for the radio structure of 3C 31. A brief discussion of the implications for 3C 449 is given.

Gower, A. C.↗
NASA NTRS2018

Seeing Double: ASASSN-18bt Exhibits a Two-component Rise in the Early-time K2 Light Curve

On 2018 February 4.41, the All-Sky Automated Survey for SuperNovae (ASAS-SN) discovered ASASSN-18bt in the K2 Campaign 16 field. With a redshift of z=0.01098 and a peak apparent magnitude of B(max)=14.31, ASASSN-18bt is the nearest and brightest SNe Ia yet observed by the Kepler spacecraft. Here we present the discovery of ASASSN-18bt, the K2 light curve, and prediscovery data from ASAS-SN and the Asteroid Terrestrial-impact Last Alert System. The K2 early-time light curve has an unprecedented 30-minute cadence and photometric precision for an SN Ia light curve, and it unambiguously shows a ∼4 day nearly linear phase followed by a steeper rise. Thus, ASASSN-18bt joins a growing list of SNe Ia whose early light curves are not well described by a single power law. We show that a double-power-law model fits the data reasonably well, hinting that two physical processes must be responsible for the observed rise. However, we find that current models of the interaction with a nondegenerate companion predict an abrupt rise and cannot adequately explain the initial, slower linear phase. Instead, we find that existing published models with shallow 56Ni are able to span the observed behavior and, with tuning, may be able to reproduce the ASASSN-18bt light curve. Regardless, more theoretical work is needed to satisfactorily model this and other early-time SNe Ia light curves. Finally, we use Swift X-ray nondetections to constrain the presence of circumstellar material (CSM) at much larger distances and lower densities than possible with the optical light curve. For a constant-density CSM, these nondetections constrain ρ<4.5×10(exp 5)per cu.cm at a radius of 4×10(exp 15) cm from the progenitor star. Assuming a wind-like environment, we place mass loss limits of M˙ < 8 x 10(exp -6) M(ʘ)per yr for v(w)=100 km/s, ruling out some symbiotic progenitor systems. This work highlights the power of well-sampled early-time data and the need for immediate multiband, high-cadence follow-up for progress in understanding SNe Ia.

B. J. Shappee↗
NASA NTRSDate not supplied

A Misperception in Reliability Growth Modelling

The Duane reliability growth model is n(t)/t = k t^-alpha (1) The reliability growth rate is alpha, the downward slope of n(t)/t versus t. It usually varies from 0.2 to 0.6. k is a constant. Crow used a 56-failure data set to illustrate reliability growth.1 A graphical Duane model fit to this data gives n(t)/t = 0.640 t^-0.283 (2) A problem in using the Duane-Crow reliability growth model is that it assumes that reliability growth continues and the failure rate decreases throughout the test period. It is more usual that reliability growth stops when the failure rated is low enough. Growth testing is often followed by testing with a low constant failure rate due to rare or uncorrectable failure modes. As more and more low constant rate acceptable failures accumulate after the period of reliability growth, the reliability growth time exponent alpha decreases toward zero. This occurs if constant rate failures are treated as occurring during the reliability growth period. It is more accurate to model a period of initial reliability growth followed by testing without repair to more accurately determine the final constant failure rate. This is done in the abcd model. n(t)/t = a t^-b + c from t = 0 to td (3) = c + d after td, where d = a td^-b (4) The term a t^-b describes the continuous reliability growth that continues out to time td and c is the constant uncorrected failure rate. The parameter d represents an additional constant failure rate due to correctable but uncorrected failure modes. After the reliability growth process is terminated, the failure rate n(t)/t = c + d.

Harry W Jones↗
NASA NTRS1997

Spectral and Temporal Properties of Galactic Black Hole Systems

Kusunose, Mineshige & Yamada (1996; hereafter KMY) extended the model of Kusunose & Mineshige (1995) to the Galactic black hole candidates by considering nonthermal electron injection with gamma(EQ\0(,\s\up2(less than),\s\do-l(_))) 10. The effects of pair escape and advection on the disk structure and general relativistic effects on the emission spectrum were also examined. They found that the energy spectral index (alpha)(sub x) of the power law X-rays is about-0.8 and-2.0 when 1(sub soft)/1 = 0.2 and 2, respectively, where 1(sub soft)/1 is the ratio of the compactness of the injected soft photons to that of the gravitational energy. The power law index was found to be nearly independent of the mass accretion which is consistent with the observed luminosity independence. The model with small 1(sub soft)/1 (less than 1) shows promise for explaining the low state observed in Galactic black hole candidates. Model fits were provided for GX339-4 and Cyg X- 1 data from COMPTEL and OSSE on the Compton Gamma Ray Observatory. The difference in emission spectra between thermal disks and the model of KMY appears only in the energy range greater than 100 keV. Li, Kusunose and Liang (1996) studied stochastic particle acceleration to produce nonthermal particle distributions which then were used in the model of Kusunose & Mineshige (1995) to model the spectrum above 1 Mev from GBHC's. Under certain conditions, stochastic electron acceleration overcame Coulomb and Compton losses resulting in a suprathermal electron population. Good fits were obtained by COMPTEL and OSSE observations of Cyg X-1 and GRO J0422+23. Kusunose & Mineshige (1996a) examined the role of electron-positron pairs in advection-dominated disks. They found that the results for advection-dominated disks without pairs are not qualitatively changed by including pairs. Summaries of work sponsored by this grant are given in Wheeler, Kim, Moscoso, Kusunose & Mineshige (1996) and Kusunose (1996) Work was also done on developing a model for an e(+-) pair wind from the inner disk region of a black hole. The model consists of three zones: a pair production/annihilation zone at the base of the wind, a pair annihilation zone slightly further out from the disk and a wind zone extending to infinity where no annihilation occurs. The model assumes an input X-ray / gamma-ray power-law spectrum as a function of photon energy and angular distribution. Pairs are created in the pair production/annihilation zone via photon-photon, photonparticle and particle-particle collisions. The bulk velocity of the pairs is obtained from the momentum component perpendicular to the disk taking into account the radiation pressure acceleration of the pairs. Energy balance in the pair production/annihilation zone is calculated by assuming that the momentum component parallel to the disk is thermalized and by taking into account heating/cooling via anisotropic thermal Compton scattering along with the cooling processes of bremsstrahlung and synchrotron radiation. The pair density is also calculated. The outer boundary of the pair production/annihilation zone is where the optical depth for photon-photon pair produciton is unity. Beyond this point only pair annihilation will occur in the pair annihilation zone. When the pair density becomes small, pair annihilation will be negligible and the pairs will flow freely to infinity. In this model we have found that the X-ray power-law spectral index, alpha(sub x) is the primary parameter which determines the density, temperature and velocity of the pair production/annihilation zone (Moscoso, Kusunose & Wheeler 1996).

Wheeler, J. Craig↗
DOE OSTI2025

Q -score as a reliability measure for protein, nucleic acid and small-molecule atomic coordinate models derived from 3DEM maps

Atomic coordinate models are important for the interpretation of 3D maps produced with cryoEM and cryoET (3D electron microscopy; 3DEM). In addition to visual inspection of such maps and models, quantitative metrics can inform about the reliability of the atomic coordinates, in particular how well the model is supported by the experimentally determined 3DEM map. A recently introduced metric, Q-score, was shown to correlate well with the reported resolution of the map for well fitted models. Here, we present new statistical analyses of Q-score based on its application to ∼10 000 maps and models archived in the EMDB (Electron Microscopy Data Bank) and PDB (Protein Data Bank). Further, we introduce two new metrics based on Q-score to represent each map and model relative to all entries in the EMDB and those with similar resolution. We explore through illustrative examples of proteins, nucleic acids and small molecules how Q-scores can indicate whether the atomic coordinates are well fitted to 3DEM maps and also whether some parts of a map may be poorly resolved due to factors such as molecular flexibility, radiation damage and/or conformational heterogeneity. These examples and statistical analyses provide a basis for how Q-scores can be interpreted effectively in order to evaluate 3DEM maps and atomic coordinate models prior to publication and archiving.

B factors↗
DOE OSTI2016

Partitioning controls on Amazon forest photosynthesis between environmental and biotic factors at hourly to interannual timescales

Gross ecosystem productivity (GEP) in tropical forests varies both with the environment and with biotic changes in photosynthetic infrastructure, but our understanding of the relative effects of these factors across timescales is limited. Here, we used a statistical model to partition the variability of seven years of eddy covariance-derived GEP in a central Amazon evergreen forest into two main causes: variation in environmental drivers (solar radiation, diffuse light fraction, and vapor pressure deficit) that interact with model parameters that govern photosynthesis and biotic variation in canopy photosynthetic light-use efficiency associated with changes in the parameters themselves. Our fitted model was able to explain most of the variability in GEP at hourly (R 2 = 0.77) to interannual (R 2 = 0.80) timescales. At hourly timescales, we found that 75% of observed GEP variability could be attributed to environmental variability. When aggregating GEP to the longer timescales (daily, monthly, and yearly), however, environmental variation explained progressively less GEP variability: At monthly timescales, it explained only 3%, much less than biotic variation in canopy photosynthetic light-use efficiency, which accounted for 63%. These results challenge modeling approaches that assume GEP is primarily controlled by the environment at both short and long timescales. Our approach distinguishing biotic from environmental variability can help to resolve debates about environmental limitations to tropical forest photosynthesis. For example, we found that biotically regulated canopy photosynthetic light-use efficiency (associated with leaf phenology) increased with sunlight during dry seasons (consistent with light but not water limitation of canopy development) but that realized GEP was nonetheless lower relative to its potential efficiency during dry than wet seasons (consistent with water limitation of photosynthesis in given assemblages of leaves). Lastly, this work highlights the importance of accounting for differential regulation of GEP at different timescales and of identifying the underlying feedbacks and adaptive mechanisms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

The Y-type Brown Dwarfs: Estimates of Mass and Age from New Astrometry, Homogenized Photometry, and Near-infrared Spectroscopy

The survey of the mid-infrared sky by the Wide-field Infrared Survey Explorer ( WISE ) led to the discovery of extremely cold, low-mass brown dwarfs, classified as Y dwarfs, which extend the T class to lower temperatures. Twenty-four Y dwarfs are known at the time of writing. Here we present improved parallaxes for four of these, determined using Spitzer images. We give new photometry for four late-type T and three Y dwarfs and new spectra of three Y dwarfs, obtained at Gemini Observatory. We also present previously unpublished photometry taken from HST , ESO, Spitzer , and WISE archives of 11 late-type T and 9 Y dwarfs. The near-infrared data are put onto the same photometric system, forming a homogeneous data set for the coolest brown dwarfs. We compare recent models to our photometric and spectroscopic data set. We confirm that nonequilibrium atmospheric chemistry is important for these objects. Nonequilibrium cloud-free models reproduce well the near-infrared spectra and mid-infrared photometry for the warmer Y dwarfs with 425 ≤ T {sub eff} (K) ≤ 450. A small amount of cloud cover may improve the model fits in the near-infrared for the Y dwarfs with 325 ≤ T {sub eff} (K) ≤ 375. Neither cloudy nor cloud-free models reproduce the near-infrared photometry for the T {sub eff} = 250 K Y dwarf W0855. We use the mid-infrared region, where most of the flux originates, to constrain our models of W0855. We find that W0855 likely has a mass of 1.5–8 Jupiter masses and an age of 0.3–6 Gyr. The Y dwarfs with measured parallaxes are within 20 pc of the Sun and have tangential velocities typical of the thin disk. The metallicities and ages we derive for the sample are generally solar-like. We estimate that the known Y dwarfs are 3 to 20 Jupiter-mass objects with ages of 0.6–8.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

High resolution, low temperature photoabsorption cross-section of C2H2 with application to Saturn's atmosphere

New laboratory observations of the VUV absorption cross-section of C2H2, obtained under physical conditions approximating stratospheres of the giant planets, were combined with IUE observations of the albedo of Saturn, for which improved data reduction techniques have been used, to produce new models for that atmosphere. When the effects of C2H2 absorption are accounted for, additional absorption by other molecules is required. The best-fitting model also includes absorption by PH3, H2O, C2H6 and CH4. A small residual disagreement near 1600 A suggests that an additional trace species may be required to complete the model.

Caldwell, John↗
NASA NTRS2004

The Distribution of Cosmic-Ray Sources in the Galaxy, Gamma-Rays and the Gradient in the CO-to-H2 Relation

We present a solution to the apparent discrepancy between the radial gradient in the diffuse Galactic gamma-ray emissivity and the distribution of supernova remnants, believed to be the sources of cosmic rays. Recent determinations of the pulsar distribution have made the discrepancy even more apparent. The problem is shown to be plausibly solved by a variation in the Wco-to-N(H2) scaling factor. If this factor increases by a factor of 5-10 from the inner to the outer Galaxy, as expected from the Galactic metallicity gradient and supported by other evidence, we show that the source distribution required to match the radial gradient of gamma-rays can be reconciled with the distribution of supernova remnants as traced by current studies of pulsars. The resulting model fits the EGRET gamma-ray profiles extremely well in longitude, and reproduces the mid-latitude inner Galaxy intensities better than previous models.

Strong, A. W.↗
DOE OSTI2017

Integrated Nondestructive Assay Systems to Estimate Plutonium in Spent Fuel Assemblies

An integrated nondestructive assay (NDA) system combining active (neutron generator) and passive neutron detection and passive gamma (PG) detection is being analyzed in order to estimate the amount of plutonium, verify initial enrichment, burnup, and cooling time, and detect partial defects in a spent fuel assembly (SFA). Active signals are measured using the differential die-away (DDA), delayed neutron (DN), and delayed gamma (DG) techniques. Passive signals are measured using total neutron (TN) counts and both gross and spectral resolved gamma counts. To quantify how a system of several NDA techniques is expected to perform, all of the relevant NDA techniques listed above were simulated as a function of various reactor conditions such as initial enrichment, burnup, cooling time, assembly shuffling pattern, reactor operating conditions (including temperature, pressure, and the presence of burnable poisons) by simulating the NDA response for five sets of light water reactor assemblies. In this study, we compare the performance of several exploratory model-fitting options (including neural networks, adaptive regression with splines, iterative bias reduction smoothing, projection pursuit regression, and regression with quadratic terms and interaction terms) to relate data simulated with measurement and model error effects from various subsets of the NDA techniques to the total Pu mass. Isotope masses for SFAs and expected detector responses (DRs) for several NDA techniques are simulated using MCNP, and the DRs become inputs to the fitting process. Such responses include eight signals from DDA, one from DN, one from TN, and up to seven from PG; the DG signal will be examined separately. Results are summarized using the root-mean-squared estimation error for plutonium mass in held-out subsets of the data for a range of model and measurement error variances. Different simulation assumptions lead to different spent fuel libraries relating DRs to Pu mass. Some results for training with one library and testing with another library are also given.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2014

Modeling the Dynamics and Migratory Pathways of Virus-Specific Antibody-Secreting Cell Populations in Primary Influenza Infection

The B cell response to influenza infection of the respiratory tract contributes to viral clearance and establishes profound resistance to reinfection by related viruses. Numerous studies have measured virus-specific antibody-secreting cell (ASC) frequencies in different anatomical compartments after influenza infection and provided a general picture of the kinetics of ASC formation and dispersion. However, the dynamics of ASC populations are difficult to determine experimentally and have received little attention. Here, we applied mathematical modeling to investigate the dynamics of ASC growth, death, and migration over the 2-week period following primary influenza infection in mice. Experimental data for model fitting came from high frequency measurements of virus-specific IgM, IgG, and IgA ASCs in the mediastinal lymph node (MLN), spleen, and lung. Model construction was based on a set of assumptions about ASC gain and loss from the sampled sites, and also on the directionality of ASC trafficking pathways. Most notably, modeling results suggest that differences in ASC fate and trafficking patterns reflect the site of formation and the expressed antibody class. Essentially all early IgA ASCs in the MLN migrated to spleen or lung, whereas cell death was likely the major reason for IgM and IgG ASC loss from the MLN. In contrast, the spleen contributed most of the IgM and IgG ASCs that migrated to the lung, but essentially none of the IgA ASCs. This finding points to a critical role for regional lymph nodes such as the MLN in the rapid generation of IgA ASCs that seed the lung. Results for the MLN also suggest that ASC death is a significant early feature of the B cell response. Overall, our analysis is consistent with accepted concepts in many regards, but it also indicates novel features of the B cell response to influenza that warrant further investigation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

PANCHROMATIC IMAGING OF A TRANSITIONAL DISK: THE DISK OF GM AUR IN OPTICAL AND FUV SCATTERED LIGHT

We have imaged GM Aurigae with the Hubble Space Telescope , detected its disk in scattered light at 1400 and 1650 Å, and compared these with observations at 3300 Å, 5550 Å, 1.1 μ m, and 1.6 μ m. The scattered light increases at shorter wavelengths. The radial surface brightness profile at 3300 Å shows no evidence of the 24 au radius cavity that has been previously observed in submillimeter observations. Comparison with dust grain opacity models indicates that the surface of the entire disk is populated with submicron grains. We have compiled a spectral energy distribution from 0.1 μ m to 1 mm and used it to constrain a model of the star + disk system that includes the submillimeter cavity using the Monte Carlo radiative transfer code by Barbara Whitney. The best-fit model image indicates that the cavity should be detectable in the F330W bandpass if the cavity has been cleared of both large and small dust grains, but we do not detect it. The lack of an observed cavity can be explained by the presence of submicron grains interior to the submillimeter cavity wall. We suggest one explanation for this that could be due to a planet of mass <9 M {sub J} interior to 24 au. A unique cylindrical structure is detected in the far-UV data from the Advanced Camera for Surveys/Solar Blind Channel. It is aligned along the system semiminor axis, but does not resemble an accretion-driven jet. The structure is limb brightened and extends 190 ± 35 au above the disk midplane. The inner radius of the limb brightening is 40 ± 10 au, just beyond the submillimeter cavity wall.

79 ASTRONOMY AND ASTROPHYSICS↗