Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Probability density function”

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,009 records · Page 56

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion↗

Supporting information for Few-Shot Learning Enables Population-Scale Analysis of Leaf Traits in Populus trichocarpa

In this work, we use few-shot learning to segment the body and vein architecture of P. trichocarpa leaves from high-resolution scans obtained in the UC Davis common garden. Leaf and vein segmentation are formulated as separate tasks, in which convolutional neural networks (CNNs) are used to iteratively expand partial segmentations until reaching stopping criteria. Our leaf and vein segmentation approaches use just 50 and 8 manually traced images for training, respectively, and are applied to a set of 2,634 top and bottom leaf scans. We show that both methods achieve high segmentation accuracy, in some cases exceeding even human-level segmentation. The leaf and vein segmentations are subsequently used to extract 68 morphological traits using traditional open-source image processing tools, which are validated using real-world physical measurements. For a biological perspective, we perform a genome-wide association study using the vein density trait to discover novel genetic architectures associated with multiple physiological processes relating to leaf development and function. In addition to sharing all of the few-shot learning code (see https://github.com/jlager/few-shot-leaf-segmentation), we are releasing all images, manual segmentations, model predictions, 68 extracted leaf phenotypes, and a new set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes. The data folder includes all images, ground truth segmentations, predicted segmentations, and extracted leaf traits. All images encode the sample ID in the file name by indicating the treatment, block, row, position, and leaf side, respectively. For example, the file, C_1_1_2_bot.jpeg, indicates the control treatment, block 1, row 1, position 2, and the bottom side of the leaf. Tabulated results include position IDs as well as the corresponding genotype IDs. The images folder includes the 2,906 high-resolution leaf scans taken in the field. The leaf_masks folder includes 50 ground truth segmentations used for training the leaf tracing algorithm. The leaf_preds folder includes the 2,906 predicted segmentations from the leaf tracing algorithm. The vein_masks folder includes 8 ground truth segmentations used for training the vein growing algorithm. The vein_preds folder includes the 1,453 predicted segmentations from the vein growing algorithm. The vein_probs folder includes the 1,453 predicted probability maps from the vein growing algorithm before thresholding. The genomes folder includes the set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes with a README file detailing the steps taken. The results folder includes: (i) raw values of the 68 predicted leaf traits in digital_traits.tsv, (ii) manually measured values of petiole length and width in manual_traits.tsv, (iii) thin plate spline (TPS) adjusted values of the vein density trait in vein_density_tps_adj.tsv, (iv) best linear unbiased prediction (BLUP) adjusted values of the vein density trait in vein_density_blups.tsv, and (v) GWAS results for the vein density trait, including chromosome positions and corresponding P values, in gwas_results.csv.

09 BIOMASS FUELS↗

Gas kinematics and ionization along the extended sight line to HD 116852

We present Goddard High-Resolution Spectrograph intermediate observations of the interstellar medium toward HD 116852, a low halo star at a distance of 4.8 kpc (z = -1.3 kpc) in the direction l = 304. deg 9, b = 16.deg 1. The small science aperture observations have signal-to-noise ratios ranging from 30 to 90 and resolutions of 11 to 18 km/s (FWHM). We confirm the optical MK classification of this star through an analysis of its ultraviolet photosperic and stellar wind profiles. We detect interstellar lines of Al III, Si IV, C IV, and N V together with lines of C I, C I*, C I**, Si II, Ge II, P II, and Ni II. We convert the Mg II, P II, S II, Al III, Si IV, C IV, and N V profiles into measure of apparent column density as a function of LSR velocity. Gas scale height, velocity dispersion, and differential Galactic rotation effects govern the profile shapes. A simple computer model of the expected sight line column density profiles for the low and high ion species indicates that the gas velocity dispersions and scale heights increase as the ionization level of the gas increases. We find scale heigts H greater than or = 1 kpc for the high ions, which are comparable to the z-distance of the star, whereas we find H approximatley = 0.6 to 0.7 kpc for A1 III and H approximatley 0.1 kpc for P II and Ge II. An enhancement in the Al II profile near -15 km/s accounts for approximately 25% of the A1 III column along the sight line and probably arises within gas located approximately 500 pc below the Sagittarius-Carina spiral arm link. Portions of the broad underlying A1 III distribution are associated with the higher ionization lines, perhaps in conductive interfaces. The presence of N V and the column density ratios of Si IV, C IV, and N V favor the interpretation that much of the high ion absorption is produced by collisional ionization in gas with T = 1-3 x 10(exp 5) K. An enhancement near -35 km/s in both the Si IV and C IV profiles may be due to an outflow from the Norma spiral arm at a z-distance of about -1 kpc. The smooth decrease of the N(C IV)/N(Si IV) ratio at negative velocities may be due to an ISM ionization structure for hot gas that changes distance from the Galactic plane or to a two phase gas distribution in which the relative contribution from each phase changes with distance from the Galactic plane.

Sembach, Kenneth R.↗

Observable measurement-induced transitions

One of the main postulates of quantum mechanics is that measurements destroy quantum coherence (wave function collapse). Recently it was discovered that in a many-body system dilute local measurements still preserve some coherence across the entire system. As the measurement density is increased, a phase transition occurs that is characterized by the disentanglement of different parts of the system. Unfortunately, this transition is impossible to observe experimentally for macroscopic systems because it requires an exponentially costly full tomography of the many-body wave function or a comparison with the simulation on an oracle classical computer. In this work we report the discovery of another measurement-induced phase transition that can be observed experimentally if quantum dynamics can be reversed. On one side of this phase transition the quantum information encoded in some part of the Hilbert space is fully recovered after the time inversion. On the other side, all quantum information is corrupted. This transition also manifests itself as the change in the behavior of the probability to observe the same measurement outcome in the process that consists of identical blocks repeated many times. In each block the unitary evolution is followed by the measurement. On one side of the transition the probability decreases exponentially with the number of repetitions, on the other it tends to a constant as the number of repetitions is increased. We confirm the existence of the proposed phase transition through numerical simulations of realistic quantum circuits and analytical calculations using an effective random-matrix theory model.

Measurement-induced phase transitions↗

Dark Energy Survey Year 3 results: Cosmological constraints from galaxy clustering and weak lensing

We present the first cosmology results from large-scale structure using the full 5000 deg 2 of imaging data from the Dark Energy Survey (DES) Data Release 1. We perform an analysis of large-scale structure combining three two-point correlation functions ( 3 × 2 pt ): (i) cosmic shear using 100 million source galaxies, (ii) galaxy clustering, and (iii) the cross-correlation of source galaxy shear with lens galaxy positions, galaxy–galaxy lensing. To achieve the cosmological precision enabled by these measurements has required updates to nearly every part of the analysis from DES Year 1, including the use of two independent galaxy clustering samples, modeling advances, and several novel improvements in the calibration of gravitational shear and photometric redshift inference. The analysis was performed under strict conditions to mitigate confirmation or observer bias; we describe specific changes made to the lens galaxy sample following unblinding of the results and tests of the robustness of our results to this decision. We model the data within the flat Λ CDM and w CDM cosmological models, marginalizing over 25 nuisance parameters. We find consistent cosmological results between the three two-point correlation functions; their combination yields clustering amplitude S 8 = 0.77 6 - 0.017 + 0.017 and matter density Ω m = 0.33 9 - 0.031 + 0.032 in Λ CDM , mean with 68% confidence limits; S 8 = 0.77 5 - 0.024 + 0.026 , Ω m = 0.35 2 - 0.041 + 0.035 , and dark energy equation-of-state parameter w = - 0.9 8 - 0.20 + 0.32 in w CDM . These constraints correspond to an improvement in signal-to-noise of the DES Year 3 3 × 2 pt data relative to DES Year 1 by a factor of 2.1, about 20% more than expected from the increase in observing area alone. This combination of DES data is consistent with the prediction of the model favored by the Planck 2018 cosmic microwave background (CMB) primary anisotropy data, which is quantified with a probability-to-exceed p = 0.13 –0.48. We find better agreement between DES 3 × 2 pt and Planck than in DES Y1, despite the significantly improved precision of both. When combining DES 3 × 2 pt data with available baryon acoustic oscillation, redshift-space distortion, and type Ia supernovae data, we find p = 0.34 . Combining all of these datasets with Planck CMB lensing yields joint parameter constraints of S 8 = 0.81 2 - 0.008 + 0.008 , Ω m = 0.30 6 - 0.005 + 0.004 , h = 0.68 0 - 0.003 + 0.004 , and ∑ m ν < 0.13 eV (95% C.L.) in Λ CDM ; S 8 = 0.81 2 - 0.008 + 0.008 , Ω m = 0.30 2 - 0.006 + 0.006 , h = 0.68 7 - 0.007 + 0.006 , and w = - 1.03 1 - 0.027 + 0.030 in w CDM .

79 ASTRONOMY AND ASTROPHYSICS↗

What if chondritic porous interplanetary dust particles are not the real McCoy

To select a target comet for a Comet Nucleus Sample Return Mission (CNSRM) it is necessary to have an experimental data base to evaluate the extent of diversity and similarity of comets. For example, the physical properties (e.g., low density) of chondritic porous (CP) interplanetary dust particles (IDPs) are believed to resemble these properties of cometary dust although it is yet to be demonstrated that the porous structure of CP IDPs is inherent to presolar dust particles stored in comet nuclei. Porous structures of IDPs could conceivably form during sublimation at the surface of active comet nuclei. Porous structures are also obtained during annealing of amorphous Mg-SiO smokes which initially forms porous aggregates of olivine + platey tridymite and which, upon continued annealing, react to fluffy enstatite aggregates. It is therefore uncertain that CP IDPs are entirely composed of unmetamorphosed presolar dust. Conceivably, new minerals and textures may form in situ in nuclei of active comets as a function of their individual thermal history. Unmetamorphosed comet dust is probably structurally amorphous. Thermal annealing of this dust can produce ultra fine-grained minerals and this ultrafine grain size of CP IDPs should be considered in assessments of aqueous alterations that could affect presolar dust in comet nuclei between 200 and 400 K. Devitrification and hydration may occur in situ in ice-dust mixtures and the mantle of active comet nuclei. Devitrification, or uncontrolled crystallization, of amorphous precursor dust can produce a range of chemical compositions of ultrafine-grained minerals and (non-equilibrium) mineral assemblages and textures in dust contained in comet nuclei as a function of period and trajectory of orbit and number of perihelion passages (not considering internal heating). Thus, experimental data on relevant processes and reaction rates between 200 and 400 K are needed in order to evaluate comet selection, penetration depth for sampling device and curation of samples for CNSRM.

Rietmeijer, Frans J. M.↗

High-throughput validation of phase formability and simulation accuracy of Cantor alloys

High-throughput methods enable accelerated discovery of novel materials in complex systems such as high-entropy alloys, which exhibit intricate phase stability across vast compositional spaces. Computational approaches, including Density Functional Theory (DFT) and calculation of phase diagrams (CALPHAD), facilitate screening of phase formability as a function of composition and temperature. However, the integration of computational predictions with experimental validation remains challenging in high-throughput studies. In this work, we introduce a quantitative confidence metric to assess the agreement between predictions and experimental observations, providing a quantitative measure of the confidence of machine learning models trained on either DFT or CALPHAD input in accounting for experimental evidence. The experimental dataset was generated via high-throughput in-situ synchrotron X-ray diffraction on compositionally varied FeNiMnCr alloy libraries, heated from room temperature to ~1000 °C. Agreement between the observed and predicted phases was evaluated using either temperature-independent phase classification or a model that incorporates a temperature-dependent probability of phase formation. This integrated approach demonstrates where strong overall agreement between computation and experiment exists, while also identifying key discrepancies, particularly in FCC/BCC predictions at Mn-rich regions to inform future model refinement.

36 - MATERIALS SCIENCE↗

A collective elastic fluctuation mechanism for decoupling and stretched relaxation in glassy colloidal and molecular liquids

We propose a microscopic theory for the decoupling of self-diffusion and structural relaxation in glass-forming liquids within the Elastically Collective Nonlinear Langevin Equation (ECNLE) activated dynamics framework. Our central hypothesis is that the heterogeneity relevant to this problem is static fluctuations of local density on the scale of 3–4 particle diameters and how this changes local packing correlations. These fluctuations modify the degree of dynamical cage expansion that mechanistically couples intracage large amplitude hopping and longer range collective elasticity in ECNLE theory. Decoupling only emerges in the deeply supercooled regime where the strongly temperature dependent elastic barrier becomes non-negligible relative to its noncooperative local analog. The theory makes predictions for various aspects of the decoupling phenomenon, including apparent fractional power law Stokes-Einstein behavior, that appear to be consistent with experiments and simulations on hard sphere fluids and molecular liquids. Of central importance is a microscopic connection between the barrier fluctuation variance and most probable barrier height. Sensible results are also obtained for the nonexponential stretching of a generic relaxation time correlation function and its temperature evolution. Nonuniversality can arise from the relative importance of the local and collective barriers (related to fragility) and the precise magnitude of the length scale that defines the transition from local cage to elastic physics. Comparison is made with a traplike model based on a Gaussian distribution of barriers.

Xie, Shi-Jie (ORCID:0000000333174720)↗

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Microgravity or Simulated Launch on Testicular Function in Rats

Testes from flight rats on COSMOS 2044 and simulated-launch, vivarium, or caudal-elevation control rats (5/group) were analyzed by subjective and quantitative methods. On the basis of observations of fixed tissue, it was evident that some rats had testicular abnormalities unassociated with treatment and probably existing when they were assigned randomly to the four treatment groups. Considering rats without preexisting abnormalities, diameter of seminiferous tubules and numbers of germ cells per tubule cross section were lower (P less than 0.05) in flight than in simulated-launch or vivarium rats. However, ratios of germ cells to each other or to Sertoli cells and number of homogenization-resistant spermatids did not differ from values for simulated-launch or vivarium controls. Expression of testis-specific gene products was not greatly altered by flight. Furthermore, there was no evidence for production of stress-inducible transcripts of the hsp7O or hsp9O genes. Concentration of receptors for rat luteinizing hormone in testicular tissue and surface density of smooth endoplasmic reticulum in Leydig cells were similar in flight and simulated-launch rats. However, concentrations of testosterone in testicular tissue or peripheral blood plasma were reduced (P less than 0.05) in flight rats to less than 20% of values for simulated-launch or vivarium controls. Thus spermatogenesis was essentially normal in flight rats, but production of testosterone was severely depressed. Exposure to microgravity for more than 2 wk might result in additional changes. Sequelae of reduced androgen production associated with microgravity on turnover of muscle and bone should be considered.

Amann, R. P.↗

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗

Bayesian Inference of the Symmetry Energy of Superdense Neutron-rich Matter from Future Radius Measurements of Massive Neutron Stars

In this work, using as references the posterior probability distribution functions of the equation of state (EOS) parameters inferred from the radii of canonical neutron stars (NSs) reported by the LIGO/VIRGO and NICER Collaborations based on their observations of GW170817 and PSR J0030+0451, we investigate how future radius measurements of more massive NSs will improve our current knowledge about the EOS of superdense neutron-rich nuclear matter, especially its symmetry energy term. Within the Bayesian statistical approach using an explicitly isospin-dependent parametric EOS for the core of NSs, we infer the EOS parameters of superdense neutron-rich nuclear matter from three sets of imagined mass–radius correlation data representing typical predictions by various nuclear many-body theories, that is, the radius stays the same, decreases, or increases with increasing NS mass within ±15% between 1.4 and 2.0 M$_\odot$. The corresponding NS average density increases quickly or slowly or slightly decreases as the NS mass increases from 1.4 to 2.0 M$_\odot$. While the EOSs of symmetric nuclear matter (SNM) inferred from the three data sets are approximately the same, the corresponding symmetry energies above about twice the saturation density of nuclear matter are very different, indicating that the radii of massive NSs carry important information about the high-density behavior of nuclear symmetry energy with little influence from the remaining uncertainties of the SNM EOS at suprasaturation densities.

79 ASTRONOMY AND ASTROPHYSICS↗

Solar Wind Strahl Observations and Their Implication to the Core-Halo Formation due to Scattering

A study of the kinetic properties of the strahl electron velocity distribution functions (VDF?s) in the solar wind is presented. This study focuses on the mechanisms that control and regulate the electron VDF?s and the stability of the strahl electrons in the solar wind; mechanisms that are not yet well understood. Various parameters are investigated such as the strahl-electron density, temperature anisotropy, and electron heat-flux. These parameters are used to investigate the stability of the strahl population. The analysis check for whether the strahl electrons are constrained by some instability (e.g., the whistler or KAW instabilities), or are maintained by other types of processes. The electron heat-flux and temperature anisotropy are determined by modeling of the 3D-VDF?s from which the moments properties of the various populations are obtained. The results of this study have profound implication on the current hypothesis about the probable formation of the solar wind halo electrons produced from the scattering of the strahl population. This hypothesis is strengthened by direct observations of the strahl electrons being scattered into the core-halo in an isolated event. The observation implies that the scattering of the strahl is not a continuous process but occurs in bursts in regions where conditions for wave growth providing the scattering are optimum. Sometimes, observations indicate that the strahl component is anisotropic (Tper/Tpal approx. 2). This provides a possible free energy source for the excitation of whistler waves as a possible scattering mechanism, however this condition is not always observed. The study is based on high time resolution data from the Cluster/PEACE electron spectrometer.

Vinas, Adolfo F.↗

Schematic model for induced fission in a configuration-interaction approach

We model fission at barrier-top energies in a simplified model space that permits comparison of different components of the residual nucleon-nucleon interaction. The model space is built on particle-hole excitations of reference configurations. These are Slater determinants of uniformly spaced orbitals characterized only by their quantum numbers and orbital energies. The residual interaction in the Hamiltonian includes the diabatic interaction connecting similar orbitals at different deformations, the pairing interaction between like nucleons, and a schematic off-diagonal neutron-proton interaction. We find that the fission reaction probability is sensitive to the off-diagonal neutron-proton interaction much more than to the pairing and the diabatic interactions. In particular, the transmission coefficients become insensitive to the strength of the pairing interaction when the neutron-proton interaction is large. Here, we also find that the branching ratio is insensitive to the final-state scission dynamics, as is assumed in the well-known Bohr-Wheeler theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A reexamination of two-fluid solar wind models

The two-fluid solar-wind equations have been solved by a method which is approximately 50 times faster than any previously developed, through the use of asymptotic expansions which are self-consistently iterated upon to find a solution that passes through the critical point. The energy assumptions in two-fluid solar-wind models are reexamined, and the conclusions are as follows: (1) proton thermal conduction may not be neglected, (2) the Coulomb logarithm must be calculated as a function of radius, and (3) the electron and proton temperatures at the base need not be equal, even when the time scale for energy exchange between the species is an order of magnitude smaller than the expansion time at the base. It is possible to reproduce reasonable quiet-time solar-wind parameters at 1 AU, but only if the proton temperature is approximately twice the electron temperature at 1 solar radius. This may indicate that extended proton heating is important in the outer solar corona. Winds with velocities at 1 AU of 450 km/s are generated without nonthermal energy deposition but require high proton temperatures as well as very low densities at the base. Higher-velocity solutions are not possible in a spherically symmetric geometry for reasonable particle fluxes at 1 AU, and it is suggested that these higher-velocity states probably require additional heating, acceleration mechanisms, or nonradial flow.

Nerney, S.↗

Theoretical uncertainty quantification for heavy-ion fusion

Despite recent advances and focus on rigorous uncertainty quantification for microscopic models of quantum many-body systems, the uncertainty on the dynamics of those systems has been underexplored. To address this, we have used time-dependent Hartree-Fock to examine the model uncertainty for a collection of low-energy, heavy-ion fusion reactions. Fusion reactions at near-barrier energies represent a rich test-bed for the dynamics of quantum many-body systems owing to the complex interplay of collective excitation, transfer, and static effects that determine the fusion probability of a given system. The model uncertainty is sizable for many of the systems studied and the primary contribution arises from static properties that are ill-constrained, such as the neutron radius of neutron-rich nuclei. Furthermore, these large uncertainties motivate the use of information from reactions to better constrain existing models and to infer static properties from reaction data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗

Column Grid Array Rework for High Reliability

Due to requirements for reduced size and weight, use of grid array packages in space applications has become common place. To meet the requirement of high reliability and high number of I/Os, ceramic column grid array packages (CCGA) were selected for major electronic components used in next MARS Rover mission (specifically high density Field Programmable Gate Arrays). ABSTRACT The probability of removal and replacement of these devices on the actual flight printed wiring board assemblies is deemed to be very high because of last minute discoveries in final test which will dictate changes in the firmware. The questions and challenges presented to the manufacturing organizations engaged in the production of high reliability electronic assemblies are, Is the reliability of the PWBA adversely affected by rework (removal and replacement) of the CGA package? and How many times can we rework the same board without destroying a pad or degrading the lifetime of the assembly? To answer these questions, the most complex printed wiring board assembly used by the project was chosen to be used as the test vehicle, the PWB was modified to provide a daisy chain pattern, and a number of bare PWB s were acquired to this modified design. Non-functional 624 pin CGA packages with internal daisy chained matching the pattern on the PWB were procured. The combination of the modified PWB and the daisy chained packages enables continuity measurements of every soldered contact during subsequent testing and thermal cycling. Several test vehicles boards were assembled, reworked and then thermal cycled to assess the reliability of the solder joints and board material including pads and traces near the CGA. The details of rework process and results of thermal cycling are presented in this paper.

high reliability↗