Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “virtual cluster”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

THE Pa{alpha} LUMINOSITY FUNCTION OF H II REGIONS IN NEARBY GALAXIES FROM HST/NICMOS

The H II region luminosity function (LF) is an important tool for deriving the birthrates and mass distribution of OB associations and is an excellent tracer of the newly formed massive stars and associations. To date, extensive work (predominantly in H{alpha}) has been done from the ground, which is hindered by dust extinction and the severe blending of adjacent (spatially or in projection) H II regions. Reliably measuring the properties of H II regions requires a linear resolution <40 pc, but analyses satisfying this requirement have been done only in a handful of galaxies, so far. As the first space-based work using a galaxy sample, we have selected 12 galaxies from our HST/NICMOS Pa{alpha} survey and studied the LF and size distribution of H II regions both in individual galaxies and cumulatively, using a virtually extinction-free tracer of the ionizing photon rate. The high angular resolution and low sensitivity to diffuse emission of NICMOS also offer an advantage over ground-based imaging by enabling a higher degree of de-blending of the H II regions. We do not confirm the broken power-law LFs found in ground-based studies. Instead, we find that the LFs, both individual and co-added, follow a single power law dN(L)/dln L{proportional_to}L {sup -1}, are consistent with the mass function of star clusters in nearby galaxies, and are in agreement with the results of the existing analyses with Hubble Space Telescope (HST) data. The individual and co-added size distributions of H II regions are both roughly consistent with dN(D)/dln D{proportional_to}D {sup -3}, but the power-law scaling is probably contaminated by blended regions or complexes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

High precision quantum-chemical treatment of adsorption: Benchmarking physisorption of molecular hydrogen on graphane

A multilevel hierarchical ab initio protocol for calculating adsorption on non-conducting surfaces is presented. It employs fully periodic treatment, which reaches local Møller-Plesset perturbation theory of second order (MP2) with correction for the basis set incompleteness via the local F12 technique. Post-MP2 corrections are calculated using finite clusters. That includes the coupled cluster treatment in the local and canonical frameworks (up to perturbative quadruples) and correlated core (with MP2). Using this protocol, the potential surface of hydrogen molecules adsorbed on graphane was computed. According to the calculations, hydrogen molecules are adsorbed on graphane in a perpendicular to the surface orientation with the minimum of the potential surface of around −3.6 kJ/mol located at the distance of 3.85 Å between the bond center of the hydrogen molecule and the mid-plane of graphane. The adsorption sites along the path from the downward-pointing carbon to the ring center of the graphane are energetically virtually equally preferable, which can enable nearly free translations of hydrogen molecules along these paths. Consequently, the hydrogen molecules on graphane most likely form a non-commensurate monolayer. The analysis of the remaining errors reveals a very high accuracy of the computed potential surface with an error bar of a few tenths of a kJ/mol. The obtained results are a high-precision benchmark for further theoretical and experimental studies of hydrogen molecules interacting with graphane.

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

Development of a water source heat pump hardware-in-the-loop (HIL) testing facility for smart building applications

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while maintaining the thermal comfort of building occupants. As the pivotalpoint between building electricity consumption and indoor thermal comfort, high-efficiency electrical heatpumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on the heatpump operation andperformance under different weather, building load, and grid requests calls for investigation and verification via experimental tests with actual heat pumps integrated with real-time building and grid responses. This study presents a Water-Source Heat Pump Hardware-in-The-Loop (HIL) Test Facility developed with the capability to perform such tests. The hardware configuration for this testfacility introduces a hydronic system that emulates the conditions for the heat pump water-side, and a duct system that emulates conditions for the heat pump airside. Both data acquisition and emulator control are implemented through the National Instruments (NI) LabVIEW software running on an NI PXIplatform. The HIL mechanism based on the hardware-software integration that allows the testbed to communicate with a generic simulation environment is also discussed. Currently, the test facility setup includes a single heatpump and virtual building model in EnergyPlus coupled with an occupant behavioral model in MATLAB. Preliminary test results of the current setup demonstrate the building load emulator's ability to track the simulated gone temperature with a Root Mean Square Deviation (RSME) below 0.12°C (0.216°F). An uncertainty analysis based on sensor accuracies shows that the heat pump coefficient of performance (COP) can be measured with a relative uncertainty of 10.4% in cooling and 3.7% in heating. Apartfrom the current testing on a single heat pump, the test facility also provides the flexibility to include additional heat pumps to form a heat pump cluster, as well as coupling the heat pump with active thermal storage to provide enhanced demandflexibility.

Calfa, Caleb↗
DOE OSTI2019

An Initial Investigation of the Design Challenges Associated with Reliable 100GigE Packet Capture

Network security researchers often rely on Emulytics TM to provide a way to evaluate the safety and security of real world systems. This work involves running a large number of virtual machines on a distributed platform to observe how software and hardware will respond to different types of attacks. While Emulytics TM software such as minimega provide a scalable system for conducting experiments, the sheer volume of network traffic produced in an experiment can easily exceed the rate at which data can be recorded for offline analysis. As such, researchers must perform live analytics, narrow their monitoring scope or accept that they must run an experiment multiple times to capture all the information they require. In support of Sandia's commitment to Emulytics TM , we are developing new storage components for the Carlin cluster that will enable researchers to capture significantly more network traffic from their experiments. This report provides a summary of Haoda Wang's initial investigation of how new AMD Epyc storage nodes can be adapted to perform packet capture at 100Gbps speeds with minimal loss. This work found that the NVMe storage capabilities of the Epyc architecture are suitable for capturing 100Gbps Ethernet traffic. While capturing traffic with existing libraries was surprisingly challenging, we were able to develop a DPDK-based software tool that recorded network traffic to disk with minimal packet loss.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Subcritical Multiplication Ex-Core Calculations with Shift

This report details the development and testing of Virtual Environment for Reactor Applications (VERA) as applied to full core subcritical multiplication ex-core calculations with Shift. The limiting factor in being able to run these detector calculations with Monte Carlo is the amount of computational memory used by a fully-loaded full core problem. The main development features added to VERA in this effort were enabling source biasing and separable weight windows and integrating multithreading in Shift. With these capabilities, a subcritical multiplication detector response calculation on a full core Watts Bar Nuclear Plant Unit 1 (WBN1) ex-core model was successfully run on a moderately sized computing cluster.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2022

Characterizing Student-Driven Research Investigations Contributed to the GLOBE Program Citizen Science Initiative in a Formal Education Context

The Global Learning and Observations to Benefit the Environment (GLOBE) Program offers citizen science opportunities to participants of all ages, with a focus on youth in formal classroom contexts. This study uses student investigation research reports and posters submitted to the 2018 International Virtual Science Symposium (IVSS) and Student Research Symposium (SRS) as testbeds for characterizing student-driven Earth system citizen science investigations. Secondarily, this study aimed to capture GLOBE’s alignment to existing citizen science outcomes frameworks in the literature, which have primarily focused on adults and non-formal settings. Based on a literature review, the evaluation team identified 89 potential characteristics in 27 categories to typify investigations from both formal education and citizen science perspectives. We coded the artifacts from 207 student projects, conducted quantitative analysis of frequencies, and performed a semantic network analysis. By using this networking approach, we conceptually mapped several clusters of co-occurring characteristics, defining a descriptive framework for GLOBE projects. We identified three tiers of citizen science projects, increasing in the sophistication of participants’ demonstrated science practices. The framework includes additional components that reflect student citizen scientists’ thoughtfulness and connection to context as well as their projects’ reflection of their motivation and self efficacy. Through these findings, we have identified areas where student citizen scientists would benefit from further support, and suggest here further research to incorporate the experiences of students into the broader understanding of citizen science outcomes.

network analysis↗
DOE OSTI2014

Optimization of Composite Cloud Service Processing with Virtual Machines

By leveraging virtual machine (VM) technology, we optimize cloud system performance based on refined resource allocation, in processing user requests with composite services. Our contribution is three-fold. (1) We devise a VM resource allocation scheme with a minimized processing overhead for task execution. (2) We comprehensively investigate the best-suited task scheduling policy with different design parameters. (3) We also explore the best-suited resource sharing scheme with adjusted divisible resource fractions on running tasks in terms of Proportional-share model (PSM), which can be split into absolute mode (called AAPSM) and relative mode (RAPSM). We implement a prototype system over a cluster environment deployed with 56 real VM instances, and summarized valuable experience from our evaluation. As the system runs in short supply, lightest workload first (LWF) is mostly recommended because it can minimize the overall response extension ratio (RER) for both sequential-mode tasks and parallel-mode tasks. In a competitive situation with over-commitment of resources, the best one is combining LWF with both AAPSM and RAPSM. It outperforms other solutions in the competitive situation, by 16 + % w.r.t. the worst-case response time and by 7.4 + % w.r.t. the fairness.

cloud resource allocation↗
NASA NTRS2004

Automated Classification of ROSAT Sources Using Heterogeneous Multiwavelength Source Catalogs

We describe an on-line system for automated classification of X-ray sources, ClassX, and present preliminary results of classification of the three major catalogs of ROSAT sources, RASS BSC, RASS FSC, and WGACAT, into six class categories: stars, white dwarfs, X-ray binaries, galaxies, AGNs, and clusters of galaxies. ClassX is based on a machine learning technology. It represents a system of classifiers, each classifier consisting of a considerable number of oblique decision trees. These trees are built as the classifier is 'trained' to recognize various classes of objects using a training sample of sources of known object types. Each source is characterized by a preselected set of parameters, or attributes; the same set is then used as the classifier conducts classification of sources of unknown identity. The ClassX pipeline features an automatic search for X-ray source counterparts among heterogeneous data sets in on-line data archives using Virtual Observatory protocols; it retrieves from those archives all the attributes required by the selected classifier and inputs them to the classifier. The user input to ClassX is typically a file with target coordinates, optionally complemented with target IDs. The output contains the class name, attributes, and class probabilities for all classified targets. We discuss ways to characterize and assess the classifier quality and performance and present the respective validation procedures. Based on both internal and external validation, we conclude that the ClassX classifiers yield reasonable and reliable classifications for ROSAT sources and have the potential to broaden class representation significantly for rare object types.

McGlynn, Thomas↗
DOE OSTI2012

Constraining primordial non-Gaussianity with cosmological weak lensing: shear and flexion

We examine the cosmological constraining power of future large-scale weak lensing surveys on the model of the ESA planned mission Euclid, with particular reference to primordial non-Gaussianity. Our analysis considers several different estimators of the projected matter power spectrum, based on both shear and flexion. We review the covariance and Fisher matrix for cosmic shear and evaluate those for cosmic flexion and for the cross-correlation between the two. The bounds provided by cosmic shear alone are looser than previously estimated, mainly due to the reduced sky coverage and background number density of sources for the latest Euclid specifications. New constraints for the local bispectrum shape, marginalized over σ{sub 8}, are at the level of Δf{sub NL} ∼ 100, with the precise value depending on the exact multipole range that is considered in the analysis. We consider three additional bispectrum shapes, for which the cosmic shear constraints range from Δf{sub NL} ∼ 340 (equilateral shape) up to Δf{sub NL} ∼ 500 (orthogonal shape). Also, constraints on the level of non-Gaussianity and on the amplitude of the matter power spectrum σ{sub 8} are almost perfectly anti-correlated, except for the orthogonal bispectrum shape for which they are correlated. The competitiveness of cosmic flexion constraints against cosmic shear ones depends by and large on the galaxy intrinsic flexion noise, that is still virtually unconstrained. Adopting the very high value that has been occasionally used in the literature results in the flexion contribution being basically negligible with respect to the shear one, and for realistic configurations the former does not improve significantly the constraining power of the latter. Since the shear shot noise is white, while the flexion one decreases with decreasing scale, by considering high enough multipoles the two contributions have to become comparable. Extending the analysis up to l{sub max} = 20,000 cosmic flexion, while being still subdominant, improves the shear constraints by ∼ 10% when added. However on such small scales the highly non-linear clustering of matter, the impact of baryonic physics, and the non-Gaussian part of the covariance matrix make any error estimation uncertain. By considering lower, and possibly more realistic, values of the flexion intrinsic shape noise results in flexion constraining power being a factor of ∼ 2 better than that of shear, and the bounds on σ{sub 8} and f{sub NL} being improved by a factor of ∼ 3 upon their combination.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1996

IRAS Colors of the Pleiades

We present large scale images of the infrared emission of the region around the Pleiades using the ISSA data product from the IRAS mission. Residual Zodiacal background and a discontinuity in the image due to the scanning strategy of the satellite necessitated special background subtraction methods. The 60/100 color image clearly shows the heating of the ambient interstellar medium by the cluster. The 12/100 and 25/100 images peak on the cluster as expected for exposure of small dust grains to an enhanced UV radiation field; however, the 25/100 color declines to below the average interstellar value at the periphery of the cluster. Potential causes of the color deficit are discussed. A new method of identifying dense molecular material through infrared emission properties is presented. The difference between the 100 micron flux density and the 60 micron flux density scaled by the average interstellar 60/100 color ratio (Delta I(sub 100) is a sensitive diagnostic of material with embedded heating sources (Delta I(sub 100) less than 0) and cold, dense cores (Delta I(sub 100) greater than 0). The dense cores of the Taurus cloud complex as well as Lynds 1457 are clearly identified by this method, while the IR bright but diffuse Pleiades molecular cloud is virtually indistinguishable from the nearby infrared cirrus.

Carey, Sean J.↗
DOE OSTI2016

Simulation Model of Microsphere Distribution for Selective Internal Radiation Therapy Agrees With Observations

Purpose: To perform a detailed analysis of microsphere distribution in biopsy material from a patient treated with {sup 90}Y-labeled resin spheres and characterize microsphere distribution in the hepatic artery tree, and to construct a novel dichotomous bifurcation model for microsphere deposits and evaluate its accuracy in simulating the observed microsphere deposits. Methods and Materials: Our virtual model consisted of arteries that successively branched into 2 new generations of arteries at 20 nodes. The artery diameter exponentially decreased from the lowest generation to the highest generation. Three variable parameters were optimized to obtain concordance between simulations and measure microsphere distributions: an artery coefficient of variation (ACV) for the diameter of all artery generations and the microsphere flow distribution at the nodes; a hepatic tree distribution volume (HDV) for the artery tree; and an artery diameter reduction (ADR) parameter. The model was tested against previously measured activity concentrations in 84 biopsies from the liver of 1 patient. In 16 of 84 biopsies, the microsphere distribution regarding cluster size and localization in the artery tree was determined via light microscopy of 30-μm sections (mean concentration, 14 microspheres/mg; distributions divided into 3 groups with mean microsphere concentrations of 4.6, 14, and 28 microspheres/mg). Results: Single spheres and small clusters were observed in terminal arterioles, whereas large clusters, up to 450 microspheres, were observed in larger arterioles. For 14 microspheres/mg, the optimized parameter values were ACV=0.35, HDV = 50 cm{sup 3}, and ADR=6 μm. For 4.6 microspheres/mg, ACV and ADR decreased to 0.26 and 0 μm, respectively, whereas HDV increased to 130 cm{sup 3}. The opposite trend was observed for 28 microspheres/mg: ACV = 0.49, HDV = 20 cm{sup 3}, and ADR = 8 μm. Conclusion: Simulations and measurements reveal that microsphere clusters are larger and more common in volumes with high microsphere concentrations and indicate that the spatial distribution of the artery tree must be considered in estimates of microsphere distributions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1982

Steps toward the Hubble constant. VIII - The global value

The problem of determining the global Hubble constant is approached by using two distance indicators whose dispersion in absolute magnitude is small, and hence which are virtually free of the Malmquist bias. New data are presented on the mean absolute-magnitude of the brightest individual red and blue supergiant stars in very nearby galaxies, calibrated relative to the Cepheids. The results are used to calibrate the maximum absolute blue magnitude for two well-observed Type I supernovae (SNe I). The velocity-relative distance diagram for 16 SNe I with known maximum apparent blue magnitudes is then calibrated using distances for IC 4182 and NGC 4214, with the result that the global Hubble constant equals 50 + or - 7 km per second per megaparsec. The result is compared with the known age of globular clusters in the Galaxy to show that the cosmological and the evolutionary time scales now closely agree.

Sandage, A.↗
DOE OSTI2021

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

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

Dynamic PRA-Based Estimation of PWR Coping Time Using a Surrogate Model for Accident Tolerant Fuel

In this study, we propose an interpolation-based response surface surrogate methodology to manage a large number of scenarios in dynamic probabilistic risk assessment. It adopts the shape Dynamic Time Warping algorithm to cluster the interpolation neighborhood from time series sample data. The interpolation method was adapted from Taylor Kriging to allow a reduced-order model of the Taylor series. In order to demonstrate its applicability to complex issues in risk assessment for nuclear engineering, an example risk response surface to estimate emergency core cooling system (ECCS) criteria for triplex silicon carbide (SiC) accident-tolerant fuel was constructed. The response surface was exploited to estimate the cumulative failure probability of the fuel cladding structure due to the uncertainties in operator actions and safety systems. The functional failures were assessed based on a combination of individual layer failures computed by coupling Risk Analysis Virtual Environment software with a pressurized water reactor 1000-MW(electric) RELAP5 model and the in-house fuel performance assessment module. Results showed that SiC cladding failure probability spiked less than 1 min after a large-break loss-of- coolant accident whenever the current ECCS criteria for Zircaloy-4 (Zr-4) cladding was used. However, it still provides an increased safety margin of three orders of magnitude compared to Zr-4. This positive margin could be utilized to relax active ECCS requirements by allowing deviations of up to 450 s in its actuation time. The proposed surrogate methodology generated a response surface of SiC cladding failure probability reasonably well, with a significant savings of computation time. This methodology is expected to be useful in the analysis of system response with complex uncertainty sources.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

N 5 -Scaling Excited-State-Specific Perturbation Theory

In this paper, we show that by working in a basis similar to that of the natural transition orbitals and using a modified zeroth-order Hamiltonian, the cost of a recently introduced perturbative correction to excited-state mean field theory can be reduced from seventh to fifth order in the system size. The (occupied) 2 (virtual) 3 asymptotic scaling matches that of ground-state second-order Møller-Plesset theory but with a significantly higher prefactor because the bottleneck is iterative: it appears in the Krylov-subspace-based solution of the linear equation that yields the first-order wave function. Here, we discuss the details of the modified zeroth-order Hamiltonian we use to reduce the cost and the automatic code generation process we used to derive and verify the cost scaling of the different terms. Overall, we find that our modifications have little impact on the method's accuracy, which remains competitive with singles and doubles equation-of-motion coupled cluster.

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

Studying the hadron structure with PANDA and CLAS using machine learning techniques

The hadron spectroscopy and structure are currently very active fields of research to study the non-perturbative regime of quantum chronodynamics. The first one studies the complex structure of excited hadrons by looking at their decay products, while the latter uses lepton scattering on nucleons. Both methods require reconstruction algorithms with great efficiency and good particle identification and background rejection rates. This work aims to provide these by either improving the existing methods or developing new ones. The first part of this document presents a feasibility study of a predicted hybrid charmonium state for the PANDA experiment. Lattice QCD calculations predict the ground state hybrid charmonium to be a spin exotic with quantum numbers of JP C = 1?+ at a mass of around 4.3 GeV with a width to be around 20 MeV. A machine learning based data analysis scheme is proposed to further improve the signal efficiency and the background reduction, alongside with improvements of the analysis software (PandaRoot), that are vital for this study. These improvements include a reworked clustering algorithm for the electromagnetic calorimeter (EMC) and an optimized monte carlo matching for neutral particles. The second part of this document is about studying the proton structure. A multidimensional study of the structure function ratio Fsin(?)LU /FUU has been performed for K±, based on the measurement of beam-spin asymmetries. It uses the high statistics data recorded with the CLAS12 spectrometer at Jefferson Laboratory. Fsin(?)LU is a twist-3 quantity that provides information about the quark gluon correlations in the proton. This document will present for the first time a simultaneous analysis of two kaon channels over a large kinematic range of z, xB , PT and Q2 with virtualities Q2 ranging from 1 GeV2 up to 8 GeV2 using machine learning techniques for improved particle identification.

Kripko, Aron↗
DOE OSTI2017

Flexible Spatial Partitions in MPACT Through Module-Based Data Passing

As part of the Consortium for Advanced Simulation of Light Water Reactors (CASL), the Virtual Environment for Reactor Applications (VERA) is being developed to provide high-fidelity multiphysics simulations of nuclear reactor cores. The MPACT code being developed collaboratively by Oak Ridge National Laboratory and the University of Michigan is the primary deterministic neutron transport solver available within VERA. MPACT employs the 2D/1D method to solve 3D problems using 2D method of characteristics (MOC) radial transport solvers with 1D nodal expansion method-based P3 (NEM-P3) axial transport solvers, which are coupled using axial and radial transverse leakages. A 3D coarse mesh finite difference (CMFD) solver is used to provide accelerated convergence and stability. Quarter-core problems in MPACT are typically run on moderately sized clusters ({approx}500-10,000 cores). As such, it relies heavily on spatial decomposition, both radially and axially, using MPI data communication along spatial domain boundaries to pass angular flux boundary condition and neighboring data. The initial spatial decomposition strategy implemented into MPACT overlaid a Cartesian grid over the core to form the domain boundaries. With this approach, a single buffer spanning across the entire boundary between two spatial domains is passed to a domain's neighbor. Additionally, only one neighboring domain index was allowed in each direction. For example, the entire eastern boundary of a domain must be the entire western boundary of the corresponding neighbour's domain. This implementation helped to execute cases on several thousand cores, where substantial radial decomposition was used. The default for several years was to use assembly decomposition, which yielded 73 radial partitions in the Watts Bar Unit 1 (WBN1) models being analyzed. However, as CASL formed stronger connections with industry partners, there was a strong push to allow MPACT to run on fewer cores, preferably in the 500-1,000 core range. Such cases require much less radial decomposition, typically with 8-16 radial domains on each plane. Unfortunately, the previous decomposition strategy is very limited in its ability to yield well-balanced domains, and it does not permit partitions with an arbitrary number of domains, particularly with small numbers of radial domains (3, 5, 7, 11, 13, etc.). With these limitations, a more flexible approach that provided more even-load balancing was needed. The details of this new strategy are outlined here, along with insights to its advantages. The examples provided are for the WBN1 core, which has been the focus of several initial CASL analyses. It is worth noting that spatial decomposition strategies have been employed in many codes with focus on a number of different transport methodologies. The references cited are not exhaustive as this has been very important research topic across many areas.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1988

Far-infrared observations of a luminous dust-shrouded source in the nucleus of NGC 4945

High-resolution far-infrared observations of the galaxy NGC 4945 have been obtained from the Kuiper Airborne Observatory. Using new observational techniques and nonlinear deconvolution routines, it is found that virtually all of the far-infrared luminosity originates from a nuclear source no larger than 12 arcsec x 9 arcsec (225 pc x 170 pc) in extent. This size constraint, coupled with the far-infrared dust temperature, indicates that the source is deeply embedded in dust: the lower limit for the 100 micron optical depth is 0.35, which is by far the largest yet measured in an external galaxy. Published optical spectra of NGC 4945 reveal a heavily obscured nonthermal source which exhibits broad line profiles typical of a Seyfert 2 active nucleus; it is concluded that the far-infrared emission is probably due to thermal radiation from dust grains surrounding the nonthermal nuclear source. A compact cluster of massive young stars may also contribute to the infrared luminosity, but the evidence for such star-forming activity is weak.

Brock, David↗