Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance estimation”

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 883 records · Page 49

A Comparison of Three Random Number Generators for Aircraft Dynamic Modeling Applications

Three random number generators, which produce Gaussian white noise sequences, were compared to assess their suitability in aircraft dynamic modeling applications. The first generator considered was the MATLAB (registered) implementation of the Mersenne-Twister algorithm. The second generator was a website called Random.org, which processes atmospheric noise measured using radios to create the random numbers. The third generator was based on synthesis of the Fourier series, where the random number sequences are constructed from prescribed amplitude and phase spectra. A total of 200 sequences, each having 601 random numbers, for each generator were collected and analyzed in terms of the mean, variance, normality, autocorrelation, and power spectral density. These sequences were then applied to two problems in aircraft dynamic modeling, namely estimating stability and control derivatives from simulated onboard sensor data, and simulating flight in atmospheric turbulence. In general, each random number generator had good performance and is well-suited for aircraft dynamic modeling applications. Specific strengths and weaknesses of each generator are discussed. For Monte Carlo simulation, the Fourier synthesis method is recommended because it most accurately and consistently approximated Gaussian white noise and can be implemented with reasonable computational effort.

Grauer, Jared A.↗

Activation Neutronics for a Swiss PWR: Simulations and Validation - Paper 102

In Switzerland 40% of the electricity generation is produced by nuclear power. With all four NPPs being already beyond their 30. year of operation, Nagra in collaboration with the utilities periodically contributes to the Swiss NPP decommissioning cost estimation studies. The neutron induced activation of structures and components of the NPPs is of particular relevance for the planning of the dismantling actions, the waste-packaging concept, the waste-management considerations and, therefore, the estimation of decommissioning costs. In this context, Nagra has developed a calculation sequence in order to perform state-of-the-art neutron transport and activation calculations. In this paper, aspects relevant to the neutron transport calculations for estimating the residual radioactivity within a Swiss PWR are presented. More specifically, the MCNP5 modeling approach together with the use of the ADVANTG-hybrid code for accelerating the variance reduction efforts is outlined. Additionally, the validation of the neutron transport calculations with in-situ full-cycle foil activation campaigns is discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The financial attractiveness assessment of large waste management projects registered as clean development mechanism

Highlights: • Projects are not financially attractive without registration as CDMs. • WM benchmarks and indicators are converging and reducing in variance. • A sensitivity analysis reveal that revenue has more of an effect on the financial results. • Results indicate that an extensive database would reduce WM project risk and capital costs. • Disclosure standards would make information more comparable worldwide. - Abstract: This study illustrates the financial analyses for demonstration and assessment of additionality presented in the project design (PDD) and enclosed documents of the 431 large Clean Development Mechanisms (CDM) classified as the ‘waste handling and disposal sector’ (13) over the past ten years (2004–2014). The expected certified emissions reductions (CER) of these projects total 63.54 million metric tons of CO{sub 2}eq, where eight countries account for 311 projects and 43.36 million metric tons. All of the projects declare themselves ‘not financially attractive’ without CER with an estimated sum of negative results of approximately a half billion US$. The results indicate that WM benchmarks and indicators are converging and reducing in variance, and the sensitivity analysis reveals that revenues have a greater effect on the financial results. This work concludes that an extensive financial database with simple standards for disclosure would greatly diminish statement problems and make information more comparable, reducing the risk and capital costs of WM projects.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

ADVANTG—An Automated Variance Reduction Parameter Generator (Rev. 1)

The primary objective of ADVANTG is to reduce both the user effort and the computational time required to obtain accurate and precise tally estimates across a broad range of challenging transport applications. ADVANTG has been applied to simulations of real-world radiation shielding, detection, and neutron activation problems. Examples of shielding applications include material damage and dose rate analyses of the Oak Ridge National Laboratory (ORNL) Spallation Neutron Source and High Flux Isotope Reactor (Risner and Blakeman 2013) and the ITER Tokamak (Ibrahim et al. 2011). ADVANTG has been applied to a suite of radiation detection, safeguards, and special nuclear material movement detection test problems (Shaver et al. 2011). ADVANTG has also been used in the prediction of activation rates within light water reactor facilities (Pantelias and Mosher 2013). In these projects, ADVANTG was demonstrated to significantly increase the tally figure of merit (FOM) relative to an analog MCNP simulation. The ADVANTG-generated parameters were also shown to be more effective than manually generated geometry splitting parameters.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Small-Sample Estimator Decisions – Certainly Uncertain in Human Spaceflight at NASA

Within statistics when estimating means, we rely heavily on the Central Limit Theorem (CLT) to aid in our inference, however in cases where the CLT does not apply we begin to be strongly limited in our choices. If a violation to the CLT comes in the form of small samples and unknown population variance, yet we maintain the assumption of normality, the use of t-distributions is perfectly valid. We explore cases where population normality is not assumed, where sample sizes are small, and where some actionable estimate of the central tendency of a distribution is needed such as for spaceflight operational decisions. Here we identify different scenarios based on measures taken from spaceflight, characterizing performance of estimators of central tendency within our set of simulations. We examine bias and variability of standard estimators of central tendency as they apply to varying small sample sizes and varying population distributions. Further work in this area is needed to develop a framework or set of guidelines for individuals set in these situations. Even with limited information, decisions need to be made.

Central Tendency↗

An Ensemble Covariance Framework for Quantifying Forced Climate Variability and Its Time of Emergence

Climate variability and its response to increasing greenhouse gases are important considerations for impacts and adaptation. Modeling studies commonly assess projected changes in variability in terms of changes in the variance of climate variables. Despite the distant and impactful covariations that climate variables can exhibit, the covariance response has received much less attention. Here, a novel ensemble framework is developed that facilitates a unified assessment of the response of the regional variances and covariances of a climate variable to imposed external forcings and their time of emergence from an unforced climate state. Illustrating the framework, the response of variability and covariability of land and ocean temperatures is assessed in the Community Earth System Model Large Ensemble under historical and RCP8.5 forcing. The results reveal that land temperature variance emerges from its preindustrial state in the 1950s and, by the end of the twenty-first century, grows to 1.5 times its preindustrial level. Demonstrating the importance of covariances for variability projections, the covariance between land and ocean temperature is considerably enhanced by 2100, reaching 1.4 times its preindustrial estimate. The framework is also applied to assess changes in monthly temperature variability associated with the Arctic region and the Northern Hemisphere midlatitudes. Consistent with previous studies and coinciding with sea ice loss, Arctic temperature variance decreases in most months, emerging from its preindustrial state in the late twentieth century. Overall, these results demonstrate the utility of the framework in enabling a comprehensive assessment of variability and its response to external climate forcings.

54 ENVIRONMENTAL SCIENCES↗

Doses to Carotid Arteries After Modern Radiation Therapy for Hodgkin Lymphoma: Is Stroke Still a Late Effect of Treatment?

Purpose: Hodgkin lymphoma (HL) survivors are at an increased risk of stroke because of carotid artery irradiation. However, for early-stage HL involved node radiation therapy (INRT) reduces the volume of normal tissue exposed to high doses. Here, we evaluate 3-dimensional conformal radiation therapy (3D-CRT), volumetric-modulated arc therapy (VMAT), and proton therapy (PT) delivered as INRT along with the extensive mantle field (MF) by comparing doses to the carotid arteries and corresponding risk estimates. Methods and Materials: We included a cohort of 46 supradiaphragmatic stage I-II classical HL patients. All patients were initially treated with chemotherapy and INRT delivered as 3D-CRT (30 Gy). For each patient, we simulated MF (36 Gy) and INRT plans using VMAT and PT (30 Gy). Linear dose-response curves for the 20-, 25-, and 30-year risk of stroke were derived from published HL data. Risks of stroke with each technique were calculated for all patients. Statistical analyses were performed with repeated measures analysis of variance. Results: The mean doses to the right and left common carotid artery were significantly lower with modern treatment compared with MF, with substantial patient variability. The estimated excess risk of stroke after 20, 25, and 30 years was 0.6%, 0.86%, and 1.3% for 3D-CRT; 0.67%, 0.96%, and 1.47% for VMAT; 0.61%, 0.96%, and 1.33% for PT; and 1.3%, 1.72%, and 2.61% for MF. Conclusions: INRT reduces the dose delivered to the carotid arteries and corresponding estimated risk of stroke for HL survivors. Even for the subset of patients with lymphoma close to the carotid arteries, the estimated risk is low.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗

Sparse Overcomplete Denoising: Aggregation Versus Global Optimization

Denoising is often addressed via sparse coding with respect to an overcomplete dictionary. There are two main approaches when the dictionary is composed of translates of an orthonormal basis. The first, traditionally employed by techniques such as wavelet cycle spinning, separately seeks sparsity w.r.t. each translate of the orthonormal basis, solving multiple partial optimizations and obtaining a collection of sparse approximations of the noise-free image, which are aggregated together to obtain a final estimate. The second approach, recently employed by convolutional sparse representations, instead seeks sparsity over the entire dictionary via a global optimization. It is tempting to view the former approach as providing a suboptimal solution of the latter. In this letter, we analyze whether global sparsity is a desirable property, and under what conditions the global optimization provides a better solution to the denoising problem. In particular, our experimental analysis shows that the two approaches attain comparable performance in case of natural images and global optimization outperforms the simpler aggregation of partial estimates only when the image admits an extremely sparse representation. Here, we explain this phenomenon by separately studying the bias and variance of these solutions, and by noting that the variance of the global solution increases very rapidly as the original signal becomes less and less sparse.

97 MATHEMATICS AND COMPUTING↗

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↗

The Role of Gravity Waves in the Formation and Organization of Clouds During TWPICE

All convective clouds emit gravity waves. While it is certain that convectively-generated waves play important parts in determining the climate, their precise roles remain uncertain and their effects are not (generally) represented in climate models. The work described here focuses mostly on observations and modeling of convectively-generated gravity waves, using the intensive observations from the DoE-sponsored Tropical Warm Pool International Cloud Experiment (TWP-ICE), which took place in Darwin, from 17 January to 13 February 2006. Among other things, the research has implications the part played by convectively-generated gravity waves in the formation of cirrus, in the initiation and organization of further convection, and in the subgrid-scale momentum transport and associated large-scale stresses imposed on the troposphere and stratosphere. The analysis shows two groups of inertia-gravity waves are detected: group L in the middle stratosphere during the suppressed monsoon period, and group S in the lower stratosphere during the monsoon break period. Waves belonging to group L propagate to the south-east with a mean intrinsic period of 35 h, and have vertical and horizontal wavelengths of about 5-6 km and 3000-6000 km, respectively. Ray tracing calculations indicate that these waves originate from a deep convective region near Indonesia. Waves belonging to group S propagate to the south-south-east with an intrinsic period, vertical wavelength and horizontal wavelength of about 45 h, 2 km and 2000-4000 km, respectively. These waves are shown to be associated with shallow convection in the oceanic area within about 1000 km of Darwin. The intrinsic periods of high-frequency waves are estimated to be between 20-40 minutes. The high-frequency wave activity in the stratosphere, defined by mass-weighted variance of the vertical motion of the sonde, has a maximum following the afternoon local convection indicating that these waves are generated by local convection. The wave activity is strongest in the lower stratosphere below 22 km and, during the suppressed monsoon period, is modulated with a 3-4-day period. The concentration of the wave activity in the lower stratosphere is consistent with the properties of the environment in which these waves propagate, whereas its 3-4-day modulation is explained by the variation of the convection activity in the TWP-ICE domain. At low rainfall intensity the wave activity increases as rainfall intensity increases. At high values of rainfall intensity, however, the wave activity associated with deep convective clouds is independent of the rainfall intensity. The convection and gravity waves observed during TWP-ICE are simulated with the Weather Research and Forecasting (WRF) Model. These simulations are compared with radiosonde observations described above and are used to determine some of the properties of convectively generated gravity waves. The gravity waves appear to be well simulated by the model. The model is used to explore the relationships between the convection, the gravity waves and cirrus.

54 ENVIRONMENTAL SCIENCES↗

Verification of Combined VR Techniques, Derivation of Future Time Equation, and Integration of LLNL Pulsed Sphere V&V Suite [Slides]

To determine whether the combination of forced-collision and DXTRAN variance-reduction (VR) techniques is unbiased, and to gain insight into the operation of these techniques, proof that first-moment estimates from Monte Carlo simulations employing both techniques are unbiased is developed. A general background on the forced-collision and DXTRAN VR techniques and their combination is given. Proof of an unbiased simulation is outlined by showing the equivalence of the history score moment equations of simulations with these techniques in use. A report with detailed proof of this equivalence is available upon request. The derivation of the future time equation using a similar approach, as well as a summary of the addition of the LLNL Pulsed Sphere experiments to the MCNP verification and validation suite, is also briefly discussed.

97 MATHEMATICS AND COMPUTING↗

Inverse sequential procedures for the monitoring of time series

Climate changes traditionally have been detected from long series of observations and long after they happened. The 'inverse sequential' monitoring procedure is designed to detect changes as soon as they occur. Frequency distribution parameters are estimated both from the most recent existing set of observations and from the same set augmented by 1,2,...j new observations. Individual-value probability products ('likelihoods') are then calculated which yield probabilities for erroneously accepting the existing parameter(s) as valid for the augmented data set and vice versa. A parameter change is signaled when these probabilities (or a more convenient and robust compound 'no change' probability) show a progressive decrease. New parameters are then estimated from the new observations alone to restart the procedure. The detailed algebra is developed and tested for Gaussian means and variances, Poisson and chi-square means, and linear or exponential trends; a comprehensive and interactive Fortran program is provided in the appendix.

Radok, Uwe↗

Modeling Oceanic Primary Production: Photoacclimation and Nutrient Effects on Light-saturated Photosynthesis

In this report, we describe a new model (the 'PhotoAcc' model) for estimating changes in the light-saturated rate of chlorophyll-normalized phytoplankton carbon fixation (Pbmax). The model is based on measurements conducted during the Atlantic Meridional Transect studies and the Bermuda Time Series program. The PhotoAcc model explained 64% to 82% of the observed variability in Pbmax for our data set, whereas none of the previously published Pbmax models described over the past 44 years explained any of the variance. The significance of this result is that a primary limiting factor for extracting ocean carbon fixation rates from satellite measurements of near surface chlorophyll has been errors in the estimate of Pbmax. Our new model should thus result in much improved calculations of oceanic photosynthesis and thus the role of the oceans in the global carbon cycle.

Behrenfeld, Michael J.↗

Miniature EPSPs and sensory encoding in the primary afferents of the vestibular lagena of the toadfish, Opsanus tau

The synaptic activity transmitted from vestibular hair cells of the lagena to primary afferent neurons was recorded in vitro using sharp, intracellular microelectrodes. At rest, the activity was composed of miniature excitatory postsynaptic potentials (mEPSPs) at frequencies from 5 to 20/s and action potentials (APs) at frequencies betwen 0 and 10/s. mEPSPs recorded from a single fiber displayed a large variability. For mEPSPs not triggering APs, amplitudes exhibited an average coefficient of variance (CV) of 0.323 and rise times an average CV of 0.516. APs were only triggered by mEPSPs with larger amplitudes (estimated 4-6 mV) and/or steeper maximum rate of rise (10.9 mV/ms, +/- 3.7 SD, n=4 experiments) compared to (3.50 mV/ms, +/-0.07 SD, n=6 experiments) for nontriggering mEPSPs. The smallest mEPSPs showed a fast rise time (0.99 ms between 10% and 90% of peak amplitude) and limited variability across fibers (CV:0.18) confirming that they were not attenuated signals, but rather represented single-transmitter discharges (TDs). The mEPSP amplitude and rise-time relationship suggests that many mEPSPs represented several, rather than a single pulse of secretion of TDs. According to the estimated overall TD frequency, the coincidence of TDs contributing to the same mEPSP were not statistically independent, indicating a positive interaction between TDs that is reminiscent of the way subminiature signals group to form miniature signals at the neuromuscular junction. Depending on the duration and intensity of efferent stimulation, a complete block of AP initiation occurred either immediately or after a delay of a few seconds. Efferent stimulation did not significantly change AP threshold level, but abruptly decreased mEPSP frequency to a near-complete block that followed the block of APs. Maximum mEPSP rate of rise decreased during, and recovered progressively after, efferent stimulation. After termination of efferent stimulation, mEPSP amplitude did not recover instantly and for a few seconds the amplitude distribution of synaptic events showed fewer large-amplitude events than during the control period. This confirms that mEPSP amplitude and rate of rise properties, which are critical for triggering afferent APs, are modified by efferent activity. The depression of afferent AP firing during efferent stimulation corresponded to a decrease in mEPSP frequency and, to a lesser extent, a decrease in mEPSP amplitude and rate of rise, suggesting, a decrease in the level of interaction among TDs contibuting to a mEPSP.

Non-NASA Center↗

Best estimate of the planetary boundary layer height from multiple remote sensing measurements

Remote sensing measurements have been widely used to estimate the planetary boundary layer height (PBLHT). Each remote sensing approach offers unique strengths and faces different limitations. In this study, we use machine learning (ML) methods to produce a best-estimate PBLHT (PBLHT-BE-ML) by integrating four PBLHT estimates derived from remote sensing measurements at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. Three ML models – random forest (RF) classifier, RF regressor, and light gradient-boosting machine (LightGBM) – were trained on a dataset from 2017 to 2023 that included radiosonde, various remote sensing PBLHT estimates, and atmospheric meteorological conditions. Evaluations indicated that PBLHT-BE-ML from all three models improved alignment with the PBLHT derived from radiosonde data (PBLHT-SONDE), with LightGBM demonstrating the highest accuracy under both stable and unstable boundary layer conditions. Feature analysis revealed that the most influential input features at the SGP site were the PBLHT estimates derived from (a) potential temperature profiles retrieved using Raman lidar (RL) and atmospheric emitted radiance interferometer (AERI) measurements (PBLHT-THERMO), (b) vertical velocity variance profiles from Doppler lidar (PBLHT-DL), and (c) aerosol backscatter profiles from micropulse lidar (PBLHT-MPL). The trained models were then used to predict PBLHT-BE-ML at a temporal resolution of 10 min, effectively capturing the diurnal evolution of PBLHT and its significant seasonal variations, with the largest diurnal variation observed over summer at the SGP site. We applied these trained models to data from the ARM Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) field campaign (EPC), where the PBLHT-BE-ML, particularly with the LightGBM model, demonstrated improved accuracy against PBLHT-SONDE. Analyses of model performance at both the SGP and EPC sites suggest that expanding the training dataset to include various surface types, such as ocean and ice-covered areas, could further enhance ML model performance for PBLHT estimation across varied geographic regions.

Zhang, Damao [Pacific Northwest National Laborator↗

Planck 2015 results: XX. Constraints on inflation

In this paper, we present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey, which includes more than twice the integration time of the nominal survey used for the 2013 release papers. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be n s = 0.968 ± 0.006 and tightly constrain its scale dependence to dn s / dlnk = -0.003 ± 0.007 when combined with the Planck lensing likelihood. When the Planck high-ℓ polarization data are included, the results are consistent and uncertainties are further reduced. The upper bound on the tensor-to-scalar ratio is r 0.002 < 0.11 (95% CL). This upper limit is consistent with the B-mode polarization constraint r< 0.12 (95% CL) obtained from a joint analysis of the BICEP2/Keck Array and Planck data. These results imply that V(φ) ∝ φ 2 and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as R 2 inflation. We search for several physically motivated deviations from a simple power-law spectrum of curvature perturbations, including those motivated by a reconstruction of the inflaton potential not relying on the slow-roll approximation. We find that such models are not preferred, either according to a Bayesian model comparison or according to a frequentist simulation-based analysis. Three independent methods reconstructing the primordial power spectrum consistently recover a featureless and smooth P R (k)over the range of scales 0.008 Mpc -1 ≲ k ≲ 0.1 Mpc -1 . At large scales, each method finds deviations from a power law, connected to a deficit at multipoles ℓ ≈ 20-40 in the temperature power spectrum, but at an uncompelling statistical significance owing to the large cosmic variance present at these multipoles. By combining power spectrum and non-Gaussianity bounds, we constrain models with generalized Lagrangians, including Galileon models and axion monodromy models. The Planck data are consistent with adiabatic primordial perturbations, and the estimated values for the parameters of the base Λ cold dark matter (ΛCDM) model are not significantly altered when more general initial conditions are admitted. In correlated mixed adiabatic and isocurvature models, the 95% CL upper bound for the non-adiabatic contribution to the observed CMB temperature variance is | α non - adi | < 1.9%, 4.0%, and 2.9% for CDM, neutrino density, and neutrino velocity isocurvature modes, respectively. We have tested inflationary models producing an anisotropic modulation of the primordial curvature power spectrum finding that the dipolar modulation in the CMB temperature field induced by a CDM isocurvature perturbation is not preferred at a statistically significant level. We also establish tight constraints on a possible quadrupolar modulation of the curvature perturbation. Lastly, these results are consistent with the Planck 2013 analysis based on the nominal mission data and further constrain slow-roll single-field inflationary models, as expected from the increased precision of Planck data using the full set of observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the heliospheric current sheet at Earth's orbit and model comparisons

IMP 8 magnetic field data for the first half of the year 1994, i.e., for about 6 solar rotations, are analyzed around regions of sector boundary crossings with the purpose of obtaining both gross- and fine-scale characteristics of the related heliospheric current sheets separating the observed sectors. For purposes of estimating the attitudes of the normals to the sector boundaries. analysis intervals (sometimes 30 min or more in length) allowing the field to fully complete an excursion of about 180 deg were used in the study, which consisted of variance analyses of the field within those intervals. The resulting boundary normals were analyzed and compared to known (generic) models of projected heliospheric current sheets and to a coronal field model for the same time period. One of the most outstanding features of the resulting ensemble of estimated boundary normals for this period is that they strongly prefer low inclinations, indicating that the observations do not support a 1 AU model that predicts a current sheet whose surface is approximately parallel with the sun's equator, such as the 'sombrero' model. They instead support a model that predicts a relatively high inclination current sheet at 1 AU. Also the normals assume a surprisingly large range of longitudes, somewhat favoring those consistent with a Parker model (45 deg and 225 deg) and/or radial alignment (0 deg and 180 deg). These boundary structures, as defined, are shown typically to be as broad as several hundred proton gyroradii, but having embedded within them very thin structures associated with stronger currents. Such thin structures have normals usually differing markedly from the gross boundary. For some crossings there are indications of a wave-like structure in the current sheet as it passed the spacecraft.

Lepping, R. P.↗