Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical modeling”

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 541 records · Page 30

Gravitational lens models based on submillimeter array imaging of Herschel -selected strongly lensed sub-millimeter galaxies at z > 1.5

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub-)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S {sub 500} > 100 mJy, 21 are strongly lensed (i.e., multiply imaged), 4 are moderately lensed (i.e., singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r {sub half}) and far-infrared luminosities (L {sub FIR}) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z {sub lens} > 0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500 μm flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L {sub FIR} (median L {sub FIR} = 7.9 × 10{sup 12} L {sub ☉}) and two decades in FIR luminosity surface density (median Σ{sub FIR} = 6.0 × 10{sup 11} L {sub ☉} kpc{sup –2}). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

High-speed GPU-based finite element simulations for NDT

The finite element method solved with explicit time increments is a general approach which can be applied to many ultrasound problems. It is widely used as a powerful tool within NDE for developing and testing inspection techniques, and can also be used in inversion processes. However, the solution technique is computationally intensive, requiring many calculations to be performed for each simulation, so traditionally speed has been an issue. For maximum speed, an implementation of the method, called Pogo [Huthwaite, J. Comp. Phys. 2014, doi: 10.1016/j.jcp.2013.10.017], has been developed to run on graphics cards, exploiting the highly parallelisable nature of the algorithm. Pogo typically demonstrates speed improvements of 60-90x over commercial CPU alternatives. Pogo is applied to three NDE examples, where the speed improvements are important: guided wave tomography, where a full 3D simulation must be run for each source transducer and every different defect size; scattering from rough cracks, where many simulations need to be run to build up a statistical model of the behaviour; and ultrasound propagation within coarse-grained materials where the mesh must be highly refined and many different cases run.

36 MATERIALS SCIENCE↗

Statistical considerations in design of spacelab experiments

After making an analysis of experimental error sources, statistical models were developed for the design and analysis of potential Space Shuttle experiments. Guidelines for statistical significance and/or confidence limits of expected results were also included. The models were then tested out on the following proposed Space Shuttle biomedical experiments: (1) bone density by computer tomography; (2) basal metabolism; and (3) total body water. Analysis of those results and therefore of the models proved inconclusive due to the lack of previous research data and statistical values. However, the models were seen as possible guides to making some predictions and decisions.

Robinson, J.↗

Explicit formula for the Holevo bound for two-parameter qubit-state estimation problem

The main contribution of this paper is to derive an explicit expression for the fundamental precision bound, the Holevo bound, for estimating any two-parameter family of qubit mixed-states in terms of quantum versions of Fisher information. The obtained formula depends solely on the symmetric logarithmic derivative (SLD), the right logarithmic derivative (RLD) Fisher information, and a given weight matrix. This result immediately provides necessary and sufficient conditions for the following two important classes of quantum statistical models; the Holevo bound coincides with the SLD Cramér-Rao bound and it does with the RLD Cramér-Rao bound. One of the important results of this paper is that a general model other than these two special cases exhibits an unexpected property: the structure of the Holevo bound changes smoothly when the weight matrix varies. In particular, it always coincides with the RLD Cramér-Rao bound for a certain choice of the weight matrix. Several examples illustrate these findings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Framework for the Analysis of Deep Neural Networks in Autonomous Aerospace Applications using Bayesian Statistics

Deep Neural Networks (DNNs) are considered to be key components in many autonomous systems. Applications range from vision-based obstacle avoidance to intelligent/learning control and planning. Safety-critical applications as found in the aerospace domain require that the behavior of the DNN is validated and tested rigorously for safety of the autonomous system (AUS). In this paper, we present a framework to support testing of DNNs and the analysis of the network structure. Our framework employs techniques from statistical modeling and active learning to effectively generate test cases for DNN safety testing and performance analysis. We will present results of a case study on a physics-based Deep recurrent residual neural network (DR-RNN), which has been trained to emulate the aerodynamics behavior of a fixed-wing aircraft.

Deep Neural networks↗

The Initial Mass Function of the Inner Galaxy Measured from OGLE-III Microlensing Timescales

We use the timescale distribution of ∼3000 microlensing events measured by the OGLE-III survey, together with accurate new made-to-measure dynamical models of the Galactic bulge/bar region, to measure the IMF in the inner Milky Way. The timescale of each event depends on the mass of the lensing object, together with the relative distances and velocities of the lens and source. The dynamical model statistically provides these distances and velocities, allowing us to constrain the lens mass function, and thereby infer the IMF. Parameterizing the IMF as a broken power-law, we find slopes in the main-sequence α{sub ms}=1.31±0.10|{sub stat}±0.10|{sub sys}, and brown dwarf region α{sub bd}=−0.7±0.9|{sub stat}±0.8|{sub sys}, where we use a fiducial 50% binary fraction, and the systematic uncertainty covers the range of binary fractions 0%–100%. Similarly, for a log-normal IMF we conclude M{sub c}=(0.17±0.02|{sub stat}±0.01|{sub sys}) M{sub ⊙} and σ{sub m}=0.49±0.07|{sub stat}±0.06|{sub sys}. These values are very similar to a Kroupa or Chabrier IMF, respectively, showing that the IMF in the bulge is indistinguishable from that measured locally, despite the lenses lying in the inner Milky Way where the stars are mostly ∼10 Gyr old and formed on a fast α-element enhanced timescale. This therefore constrains models of IMF variation that depend on the properties of the collapsing gas cloud.

79 ASTRONOMY AND ASTROPHYSICS↗

X-ray and rotational luminosity correlation and magnetic heating of radio pulsars

Previous works have suggested a correlation between the X-ray luminosity L{sub x} and the rotational luminosity L{sub rot} of radio pulsars. However, none of the obtained regression lines is statistically acceptable due to large scatter. We construct a statistical model that has an intrinsic L{sub x}–L{sub rot} relation and reproduces the observed L{sub x} distribution about it by using a Monte Carlo simulator, which takes into account the effects obscuring the intrinsic relation, i.e., the anisotropy of radiation, additional heating, uncertainty in distance, and the detection limit of the instruments. From the ATNF pulsar catalog we collect 57 “ordinary radio pulsars” with significant detection and 42 with upper limits. The sample does not include high-magnetic-field pulsars (>10{sup 13} G), which are analyzed separately. We obtain a statistically acceptable relation L{sub x}(0.5--10keV)=10{sup 31.69}(L{sub rot}/L{sub 0}){sup c{sub 1}} with c {sub 1} = 1.03 ± 0.27 and L {sub 0} = 10{sup 35.38}. The distribution about the obtained L{sub x}–L{sub rot} relation is reproduced well by the simulator. Pulsars with abnormally high L{sub x} fall into two types: one is the soft gamma-ray pulsars, and the other is pulsars that are thermally bright in comparison with the standard cooling curve. On the other hand, pulsars showing low L{sub x} are found to have dim pulsar wind nebulae (PWNs). We argue that there is an unknown mechanism that governs both the magnetospheric emission and the PWNs, and it might involve the production rate of electron–positron pairs. High-field pulsars form a population that is distinct from ordinary pulsars due to their excess luminosities.

79 ASTRONOMY AND ASTROPHYSICS↗

Exposure to pyrethroids insecticides and serum levels of thyroid-related measures in pregnant women

Possible association between environmental exposure to pyrethroid insecticides and serum thyroid-related measures was explored in 231 pregnant women of 10–12 gestational weeks recruited at a university hospital in Tokyo during 2009–2011. Serum levels of free thyroxine (fT4), thyroid stimulating hormone (TSH) and thyroid biding globulin (TBG) and urinary pyrethroid insecticide metabolite (3-phenoxybenzoic acid, 3-PBA) were measured. Obstetrical information was obtained from medical records and dietary and lifestyle information was collected by self-administered questionnaire. Geometric mean concentration of creatinine-adjusted urinary 3-PBA was 0.363 (geometric standard deviation: 3.06) μg/g cre, which was consistent with the previously reported levels for non-exposed Japanese adult females. The range of serum fT4, TSH and TBG level was 0.83–3.41 ng/dL, 0.01–27.4 μIU/mL and 16.4–54.4 μg/mL, respectively. Multiple regression analysis was carried out by using either one of serum levels of thyroid-related measures as a dependent variable and urinary 3-PBA as well as other potential covariates (age, pre-pregnancy BMI, parity, urinary iodine, smoking and drinking status) as independent variables: 3-PBA was not found as a significant predictor of serum level of thyroid-related measures. Lack of association may be due to lower pyrethroid insecticide exposure level of the present subjects. Taking the ability of pyrethroid insecticides and their metabolite to bind to nuclear thyroid hormone (TH) receptor, as well as their ability of placental transfer, into consideration, it is warranted to investigate if pyrethroid pesticides do not have any effect on TH actions in fetus brain even though maternal circulating TH level is not affected. -- Highlights: • Pyrethroid exposure and thyroid hormone status was examined in pregnant women. • Urinary 3-phenoxybenzoic acid was used as a biomarker of exposure. • Iodine nutrition, age and other covariates were included in statistical models. • No association was found between levels of thyroid hormone and pyrethroid exposure. • The result may be ascribed to lower exposure level.

54 ENVIRONMENTAL SCIENCES↗

Numerical simulation of plasma insulator interactions in space. Part 2: Dielectric effects

Any process which may be modeled statistically for a simple plasma particle may be modeled by the particle in cell technique. The success of the calculation is then dependent on having a large enough number of particles that the statistical treatment is meaningful. Thus it is possible to include the effects of secondary emission, backscattering, charge sticking and possibly dielectric breakdown, photoemission, and spallation. The first plasma dielectric interaction included in the computer code for simulating plasma insulator interactions is secondary electron emission. A calculated current density vs. voltage curve is presented and compared to an experimental curve.

Chaky, R. C.↗

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗

Individualized Prediction of Overall Survival After Postoperative Radiation Therapy in Patients With Early-Stage Cervical Cancer: A Korean Radiation Oncology Group Study (KROG 13-03)

Purpose: A nomogram is a predictive statistical model that generates the continuous probability of a clinical event such as death or recurrence. The aim of the study was to construct a nomogram to predict 5-year overall survival after postoperative radiation therapy for stage IB to IIA cervical cancer. Methods and Materials: The clinical data from 1702 patients with early-stage cervical cancer, treated at 10 participating hospitals from 1990 to 2011, were reviewed to develop a prediction nomogram based on the Cox proportional hazards model. Demographic, clinical, and pathologic variables were included and analyzed to formulate the nomogram. The discrimination and calibration power of the model was measured using a concordance index (c-index) and calibration curve. Results: The median follow-up period for surviving patients was 75.6 months, and the 5-year overall survival probability was 87.1%. The final model was constructed using the following variables: age, number of positive pelvic lymph nodes, parametrial invasion, lymphovascular invasion, and the use of concurrent chemotherapy. The nomogram predicted the 5-year overall survival with a c-index of 0.69, which was superior to the predictive power of the International Federation of Gynecology and Obstetrics (FIGO) staging system (c-index of 0.54). Conclusions: A survival-predicting nomogram that offers an accurate level of prediction and discrimination was developed based on a large multi-center study. The model may be more useful than the FIGO staging system for counseling individual patients regarding prognosis.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Characterizing performances of solder paste printing process at flexible manufacturing lines

Solder paste printing (SPP) has been a challenge on printed circuit board (PCB) manufacturing, evident by the proliferation of solder paste inspection equipment, or substituted by rigorous non-value added activity of manual inspections. The objective of this study is to characterize the SPP performance of various products manufactured in flexible production lines with different equipment configurations, and determine areas for process improvement. The study began by collecting information on SPP performance relative to component placement (CP) process, and to the proportion of mixed products. Using a clustering algorithm to group similar elements together, SPP performance across all product-production line pairs are statistically modeled to discover the trend and the influential factors. The main findings are: (a) Ratio of overall dpku for CP and SPP processes are 2:1; (b) logistic regression models of SPP performance indicated that only effects of product-production line and solder paste printer configuration are significant; (c) PCB circuitry design with BGA components and single solder paste printer line configurations generated the highest monthly defects, with the highest variation in the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effect of nano-SiO{sub 2} particles and curing time on development of fiber-matrix bond properties and microstructure of ultra-high strength concrete

Bond properties between fibers and cementitious matrix have significant effect on the mechanical behavior of composite materials. In this study, the development of steel fiber-matrix interfacial bond properties in ultra-high strength concrete (UHSC) proportioned with nano-SiO{sub 2} varying between 0 and 2%, by mass of cementitious materials, was investigated. A statistical model relating either bond strength or pullout energy to curing time and nano-SiO{sub 2} content was proposed by using the response surface methodology. Mercury intrusion porosimetry (MIP) and backscatter scanning electron microscopy (BSEM) were used to characterize the microstructure of the matrix and the fiber-matrix interface, respectively. Micro-hardness around the embedded fiber and hydration products of the matrix were evaluated as well. Test results indicated that the optimal nano-SiO{sub 2} dosage was 1% in terms of the bond properties and the microstructure. The proposed quadratic model efficiently predicted the bond strength and pullout energy with consideration of curing time and nano-SiO{sub 2} content. The improvement in bond properties associated with nano-silica was correlated with denser matrix and/or interface and stronger bond and greater strength of hydration products based on microstructural analysis.

36 MATERIALS SCIENCE↗

Nomograms for Prediction of Outcome With or Without Adjuvant Radiation Therapy for Patients With Endometrial Cancer: A Pooled Analysis of PORTEC-1 and PORTEC-2 Trials

Background: Postoperative radiation therapy for stage I endometrial cancer improves locoregional control but is without survival benefit. To facilitate treatment decision support for individual patients, accurate statistical models to predict locoregional relapse (LRR), distant relapse (DR), overall survival (OS), and disease-free survival (DFS) are required. Methods and Materials: Clinical trial data from the randomized Post Operative Radiation Therapy for Endometrial Cancer (PORTEC-1; N=714 patients) and PORTEC-2 (N=427 patients) trials and registered group (grade 3 and deep invasion, n=99) were pooled for analysis (N=1240). For most patients (86%) pathology review data were available; otherwise original pathology data were used. Trial variables which were clinically relevant and eligible according to data constraints were age, stage, given treatment (pelvic external beam radiation therapy (EBRT), vaginal brachytherapy (VBT), or no adjuvant treatment, FIGO histological grade, depth of invasion, and lymph-vascular invasion (LVSI). Multivariate analyses were based on Cox proportional hazards regression model. Predictors were selected based on a backward elimination scheme. Model results were expressed by the c-index (0.5-1.0; random to perfect prediction). Two validation sets (n=244 and 291 patients) were used. Results: Accuracy of the developed models was good, with training accuracies between 0.71 and 0.78. The nomograms validated well for DR (0.73), DFS (0.69), and OS (0.70), but validation was only fair for LRR (0.59). Ranking of variables as to their predictive power showed that age, tumor grade, and LVSI were highly predictive for all outcomes, and given treatment for LRR and DFS. The nomograms were able to significantly distinguish low- from high-probability patients for these outcomes. Conclusions: The nomograms are internally validated and able to accurately predict long-term outcome for endometrial cancer patients with observation, pelvic EBRT, or VBT after surgery. These models facilitate decision support in daily clinical practice and can be used for patient counseling and shared decision making, selecting patients who benefit most from adjuvant treatment, and generating new hypotheses.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Exact representations of many-body interactions with restricted-Boltzmann-machine neural networks

Restricted Boltzmann machines (RBMs) are simple statistical models defined on a bipartite graph which have been successfully used in studying more complicated many-body systems, both classical and quantum. In this work, we exploit the representation power of RBMs to provide an exact decomposition of many-body contact interactions into one-body operators coupled to discrete auxiliary fields. This construction generalizes the well known Hirsch's transform used for the Hubbard model to more complicated theories such as pionless effective field theory in nuclear physics, which we analyze in detail. Finally, we also discuss possible applications of our mapping for quantum annealing applications and conclude with some implications for RBM parameter optimization through machine learning.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Diabatic paths through the scission point in nuclear fission

An outstanding problem in the theory of nuclear fission is to understand the Hamiltonian dynamics at the scission point. In this work the fissioning nucleus is modeled in self-consistent mean-field theory as a set of generator coordinate method (GCM) configurations passing through the scission point. In contrast to previous methods, the configurations are constructed in the Hartree-Fock approximation with axially symmetric mean fields and using the K -partition numbers as additional constraints. The goal of this work is to find paths through the scission point where the overlaps between neighboring configurations are large. A measure of distance along the path is proposed that is insensitive to the division of the path into short segments. For most of the tested K partitions two shape degrees of freedom are adequate to define smooth paths. However, some of the configurations and candidate paths have sticking points where there are substantial changes in the many-body wave function, especially if quasiparticle excitations are present. The excitation energy deposited in fission fragments arising from thermal excitations in the pre-scission configurations is determined by tracking orbital occupation numbers along the scission paths. This allows us to assess the validity of the well-known scission-point statistical model, in which the scission process is assumed to be fully equilibrated up to the separated fission fragments. The nucleus 236 U is taken as a representative example in the calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Observation of incomplete fusion reactions at l < l {sub crit}

In order to understand the presence of incomplete fusion at low energies i.e. 4-7MeV/nucleon and also to study its dependence on various entrance-channel parameters, the two type of measurements (i) excitation function for {sup 12}C+{sup 159}Tb, and (ii) forward recoil ranges for {sup 12}C+{sup 159}Tb systems have been performed. The experimentally measured excitation functions have been analyzed within the framework of compound nucleus decay using statistical model code PACE4. Analysis of data suggests the production of xn/px)n-channels via complete fusion, as these are found to be well reproduced by PACE4 predictions, while, a significant enhancement in the excitation functions of α-emitting channels has been observed over the theoretical ones, which has been attributed due to the incomplete fusion processes. Further, the incomplete fusion events observed in case of forward recoil range measurements have been explained on the basis of the breakup fusion model, where these events may be attributed to the fusion of {sup 8}Be and/or {sup 4}He from {sup 12}C projectile to the target nucleus. In the present work, the SUMRULE model calculations are found to highly underestimate the observed incomplete fusion cross-sections which indicate that the l-values lower than l {sub crit} (limit of complete fusion) significantly contribute to the incomplete fusion reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Radiological Characterisation of VLLW Produced at High-Energy Particle Accelerators - 19118

In this paper we describe the characterisation method of very-low-level radioactive waste (VLLW) produced at the European Organization for Nuclear Research (CERN). The method is based on the calculation of a preliminary radionuclide inventory using either analytical or Monte Carlo codes, on the measurement of γ-emitters and on the quantification of pure-α, pure-β and low-energy X-ray emitters, called difficult-to-measure (DTM) radionuclides, using the so-called scaling factor (SF) technique. The first stage of the characterisation process is the calculation of the radionuclide inventory, which is a list of radionuclides that can potentially be produced via the interaction of primary and secondary particles with the surrounding media. Once the preliminary radionuclide inventory is obtained, the γ-emitting radionuclides are measured via γ-ray spectrometry on each package of the waste population. The major γemitter, called key nuclide (KN), is also identified. The scaling factor method estimates the activity of DTM radionuclides by checking for a consistent and repeated relationship between the key nuclide and the activity of the difficult to measure radionuclides from samples. Between 2016 and 2018, several hundred samples were collected at CERN and over 1000 radiochemical analyses were performed to estimate robust statistics to quantify scaling factors for the most common DTM radionuclides, such as H-3, Fe-55 and Ni-63. This paper illustrates the sample collection strategy, the tests performed to choose the most appropriate statistical models describing the behaviour of scaling factors and a complete quantification of the uncertainties associated with the characterisation process. Finally, we present the results obtained during the characterisation of over 2000 m{sup 3} radioactive metals and cables disposed of at the French National repository for VLLW. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗