Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution”

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 919 records · Page 51

A distributed cellulosic biorefinery system in the US Midwest based on corn stover

Corn stover supply chains in a distributed biorefinery system are explored here. The distributed cellulosic biorefinery uses pre-processed and densified cellulosic feedstock from a geographically separated facility (a depot) as raw material. A network of small-scale depot facilities supplies pre-processed feedstock to a distributed biorefinery. Depot facilities are assumed to be located at existing grain elevators, while distributed biorefineries are located adjacent to coal-fired power plants in areas with high gasoline consumption (urban areas) in the Midwest. The county level corn stover projections in 2022 by the US Billion-Ton Update report (2011) are used to estimate ethanol selling price and greenhouse gas (GHG) emissions of the ethanol fuel. The supply chain for each distributed biorefinery is determined by minimizing the ethanol selling price. Approximately ten distributed biorefineries based on corn stover could be established in the Midwest. Over 700 individual depot facilities participate in supplying the distributed biorefinery systems which collectively can produce greater than 12 hm3 of ethanol (3.3 billion gallons) per year. Ethanol selling price in the distributed system ranges from US$\$ $0.66 to US$\$ $1.03 per liter. Some distributed biorefineries are economically competitive with a centralized biorefinery. However, not every region can support a distributed biorefinery system due to inadequate corn stover availability. Cradle-to-gate GHG emissions of ethanol in the distributed systems are 22.1–46.6 g CO 2 per MJ. The external energy consumption in the depot facilities is the major GHG source. Optimizing process energy use in the depot facility is required to reduce both operation costs and GHG emissions.

09 BIOMASS FUELS↗

Compressive deformation behavior and strain rate sensitivity of Al-cenosphere hybrid foam with mono-modal, bi-modal and tri-modal cenosphere size distribution

Highlights: • Cenospheres of different size ranges and distributions are used in hybrid aluminium foam. • Microstructure and deformation response of hybrid foam cenosphere size distribution. • Hybrid foam with bi or tri-modal cenosphere distribution exhibits higher strength and energy absorption. • The plateau stresses and energy absorption increased marginally with strain rate. • In case of bi- or tri-modal size distribution, strain rate sensitivity and strain rate sensitivity parameter is marginally lower. - Abstract: AlSi12Cu1Mg1 aluminium alloy-cenosphere closed cell hybrid foams with different kinds of cenosphere size distribution (mono-modal, bi-modal and tri-modal) are used to see the effect of cenosphere size distribution on their deformation behavior under quasi static compressive loading conditions at different strain rates. The stress strain curves are almost similar type irrespective of the cenosphere size distribution and strain rate. However, the plateau stress and energy absorption increased significantly when one used bi-modal and or tri-modal cenosphere size distribution instead of mono-modal cenosphere size distribution. These plateau stresses and energy absorption also increased marginally with strain rate. The densification strain remained almost invariant with strain rate and cenosphere size distribution. The strain rate sensitivity and strain rate sensitivity parameter is relatively lower in case of bi-modal and tri-modal cenosphere size distribution particularly at higher relative density.

36 MATERIALS SCIENCE↗

Image quality of microcalcifications in digital breast tomosynthesis: Effects of projection-view distributions

Purpose: To analyze the effects of projection-view (PV) distribution on the contrast and spatial blurring of microcalcifications on the tomosynthesized slices (X-Y plane) and along the depth (Z) direction for the same radiation dose in digital breast tomosynthesis (DBT). Methods: A GE GEN2 prototype DBT system was used for acquisition of DBT scans. The system acquires PV images from 21 angles in 3 deg. increments over a {+-}30 deg. range. From these acquired PV images, the authors selected six subsets of PV images to simulate DBT of different angular ranges and angular increments. The number of PV images in each subset was fixed at 11 to simulate a constant total dose. These different PV distributions were subjectively divided into three categories: uniform group, nonuniform central group, and nonuniform extreme group with different angular ranges and angular increments. The simultaneous algebraic reconstruction technique (SART) was applied to each subset to reconstruct the DBT slices. A selective diffusion regularization method was employed to suppress noise. The image quality of microcalcifications in the reconstructed DBTs with different PV distributions was compared using the DBT scans of an American College of Radiology phantom and three human subjects. The contrast-to-noise ratio (CNR) and the full width at half maximum (FWHM) of the line profiles of microcalcifications within their in-focus DBT slices (parallel to detector plane) and the FWHMs of the interplane artifact spread function (ASF) in the Z-direction (perpendicular to detector plane) were used as image quality measures. Results: The results indicate that DBT acquired with a large angular range or, for an equal angular range,with a large fraction of PVs at large angles yielded superior ASF with smaller FWHM in the Z-direction. PV distributions with a narrow angular range or a large fraction of PVs at small angles had stronger interplane artifacts. In the X-Y focal planes, the effect of PV distributions on spatial blurring depended on the directions. In the X-direction (perpendicular to the chestwall), the normalized line profiles of the calcifications reconstructed with the different PV distributions were similar in terms of FWHM; the differences in the FWHMs between the different PV distributions were less than half a pixel. In the Y-direction (x-ray source motion), the normalized line profiles of the calcifications reconstructed with PVs acquired with a narrow angular range or a large fraction of PVs at small angles had smaller FWHMs and thus less blurring of the line profiles. In addition, PV distributions with a narrow angular range or a large fraction of PVs at small angles yielded slightly higher CNR than those with a wide angular range for small, subtle microcalcifications; however, PV distributions had no obvious effect on CNR for relatively large microcalcifications. Conclusions: PV distributions affect the image quality of DBT. The relative importance of the impact depends on the characteristics of the signal and the direction (perpendicular or parallel) relative to the direction of x-ray source motion. For a given number of PVs, the angular range and the distribution of the PVs affect the degree of in-plane and interplane blurring in opposite ways. The design of the scan parameters of tomosynthesis systems would require proper consideration of the characteristics of the signals of interest and the potential trade-off of the image quality of different types of signals.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Hydroacoustic Evaluation of Juvenile Salmonid Passage and Distribution at Lookout Point Dam, 2010

This report presents the results of an evaluation of juvenile salmonid passage and distribution at Lookout Point Dam (LOP) on the Middle Fork Willamette River. The study was conducted by the Pacific Northwest National Laboratory for the U.S. Army Corps of Engineers, Portland District (USACE). The goal of the study was to provide fish passage and distribution data to support decisions on long-term measures to enhance downstream passage at LOP and others dams in USACE’s Willamette Valley Project in response to the listing of Upper Willamette River Spring Chinook salmon (Oncorhynchus tshawytscha) and Upper Willamette River steelhead (O. mykiss) as threatened under the Endangered Species Act. During the year-long study period - February 1, 2010 to January 31, 2011the objectives of the hydroacoustic evaluation of fish passage and distribution at LOP were to: 1. Estimate passage rates, run timing, horizontal distribution, and diel distribution at turbine penstock intakes for smolt-size fish. 2. Estimate passage rates, run timing and diel distribution at turbine penstock intakes for small-size fish. 3. Estimate passage rates and run timing at the regulating outlets for smolt-size fish. 4. Estimate vertical distribution of smolt-size fish in the forebay near the upstream face of the dam. The fixed-location hydroacoustic technique was used to accomplish the objectives of this study. Transducers (420 kHz) were deployed in each penstock intake, above each RO entrance, and on the dam face; a total of nine transducers (2 single-beam and 7 split-beam) were used. We summarize the findings from the hydroacoustic evaluation of juvenile salmonid passage and distribution at LOP during February 2010 through January 2011 as follows. • Fish passage rates for smolt-size fish (> ~90 mm) were highest during December-January and lowest in mid-summer through early fall. • During the entire study period, an estimated total of 142,463 fish ± 4,444 (95% confidence interval) smolt-size fish passed through turbine penstock intakes. • Diel periodicity of smolt-size fish showing crepuscular peaks was evident in fish passage into turbine penstock intakes. • Run timing for small-size fish (~65-90 mm) peaked (702 fish) on December 18. Downstream passage of small-size juvenile fish was variable, occurring on two days in the spring, eight days in the summer, and at times throughout late fall and winter. A total of 7,017 ± 690 small-size fish passed through the turbine penstock intakes during the study period. • Relatively few fish passed into the ROs when they were open in summer (2 fish/d) and winter (8 fish/d). • Fish were surface-oriented with 62-80% above 10 m deep. The highest percentage of fish (30-60%) was in the 5-10 m depth bin. We draw the following conclusions from the study. • The non-obtrusive hydroacoustic data from this study are reliable because passage estimates and patterns were similar with those observed in the direct capture data from the tailrace screw trap and were consistent with distribution patterns observed in other studies of juvenile salmonid passage at dams. • Fish passage at LOP was apparently affected but not dominated by dam operations and reservoir elevation. • The surface-oriented vertical distribution of fish we observed supports development of surface passage or collector devices. In summary, the high-resolution spatially and temporally data reported herein provide detailed estimates of vertical, horizontal, diel, daily, and seasonal passage and distributions at LOP during March 2010 through January 2011. This information is applicable to management decisions on design and development of surface passage and collections devices to help restore Chinook salmon populations in the Middle Fork Willamette River watershed above Lookout Point Dam.

13 HYDRO ENERGY↗

A hybrid architecture for volt-var control in active distribution grids

Modern active distribution grids are characterized by the increasing penetration of distributed energy resources (DERs). The proper coordination and scheduling of a large numbers of these small-scale and spatially distributed DERs is necessary, and warrants the use of novel distributed approaches. In this paper, we propose a hybrid volt-var control architecture for the distribution grid, which leverages existing centralized and local approaches to planning, decision making, and control, and augments it with distributed optimization and distributed control for DER management. First, we propose a convex model to describe the power physics of distribution grids of meshed topology and unbalanced structure, based on current injection and McCormick Envelopes. Second, we employ the distributed proximal atomic coordination (PAC) algorithm to coordinate DERs to provide voltage support. We implement volt-var optimization by optimally coordinating DERs including PV smart inverters and demand response. We present results using the IEEE-34 bus network, using real data from a distribution feeder in Hawaii, to model load and PV generation. Different levels of DER penetration and objective functions are simulated. Finally, our results show the need for the coordination of DERs to improve voltage profiles, even in networks with existing voltage control devices. Further, we show the need for flexible reactive power capabilities to achieve desired grid performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Equity-driven Planning of Distributed Solar PV using Optimal Transport

Typically, distribution system planning processes do not explicitly incorporate energy equity considerations, such as identifying consumers most affected by energy costs and determining how investments in the distribution system can address existing energy burden imbalances. This paper proposes a novel optimal transport (OT)-based method to improve the energy burden distribution of consumers. The approach involves the strategic siting and sizing of solar PV in order to assist customers with high energy burden and improve the overall energy burden distribution of the community. The desired energy burden distribution is defined using the equal distribution equivalent (EDE) concept. The OT-based method is then used to estimate the distributed solar PV capacity to be installed at various locations and the tariffs to be adjusted, all while improving the energy burden distribution and providing valuable insights into distributed generation (DG) planning. The results on IEEE 37 bus test system demonstrate how DG planning, considering EDE and OT, can help reduce the energy burden of low-income consumers. Additionally, the approach also reveals optimal tariff adjustments needed to ensure revenue neutrality for distribution utilities.

Optimal transport, equal distribution equivalent, ↗

Experimental investigation of the effect of pore size distribution on nano-particle capture efficiency within ceramic particulate filters

The effect of the pore size distribution on size-resolved filtration efficiency was investigated for two ceramic particulate filters using particulate matter (PM) generated by a spark-ignition direct-injection engine fueled with gasoline. The cordierite filter tested had a porosity of 43%, a median pore diameter (d_50) of 12 µm, and a wide pore size distribution with a lognormal standard deviation (s’) of 0.3. The aluminum titanate filter had very similar porosity, d_50, and thickness, but significantly narrower pore size distribution (s’ = 0.1). The testing of two filters under identical experimental conditions enabled the impact of the pore size distribution on filtration performance to be evaluated. Filtration experiments were performed focusing on just the filter wall, starting from a clean filter until the transition to cake filtration (filtration efficiency > 99%). Time-resolved particle size distribution measurements were used to evaluate the progression of filtration performance and estimate trapped mass within the filter. The aluminum titanate filter, with a narrow pore size distribution, exhibited significantly better diffusion capture efficiency. The negative impact of higher flow velocity on diffusion capture efficiency was more pronounced for a narrower pore size distribution. Flow distributions measured using capillary-flow porometry were used to develop a cylindrical pore flow model to understand the impact of the differences in pore size distribution on observed trends in diffusion capture efficiency within a clean filter. The model helps demonstrate how large pores on the high end of a wider pore size distribution can carry a large fraction of the flow, with less efficient capture of small particles by diffusion. The experimental results and data demonstrate that the bubble point diameter and width of the pore size distribution significantly influence diffusion capture efficiencies for filters with very similar median pore diameter, porosity, and thickness.

aerosol filtration, Aftertreatment, ceramic exhaus↗

Joint scheduling of energy, fast and primary frequency response reserves in integrated transmission–distribution networks

Inverter-based distributed energy resources (DERs) connected to distribution networks (DNs) can provide fast frequency support, but their reserve deliverability depends on feeder constraints and differs from synchronous primary frequency response (PFR). Existing transmission–distribution coordination studies usually treat reserve generically or neglect feeder-level feasibility, while frequency-security scheduling studies rarely represent distribution feeders explicitly. This paper develops a bi-level day-ahead scheduling framework for integrated transmission–distribution networks that jointly clears energy, transmission-side PFR, and distribution-side fast frequency response (FFR) under exogenous hourly inertia and largest-loss inputs from an external unit commitment (UC) schedule. The transmission problem is modeled with DC-optimal power flow (OPF) and closed-form second-order cone (SOC) frequency-security constraints, whereas each DN is represented by a reserve-aware branch-flow AC-OPF so that scheduled fast reserves remain deliverable during activation. The bi-level problem is reformulated through Karush–Kuhn–Tucker (KKT) conditions into a mixed-integer SOC program, and a penalty term is used to tighten the distribution-network relaxation. In the reduced test system, lower exogenous inertia increased the required primary response from 179.64 MW to 191.08 MW, distribution-side fast response reduced total frequency-response procurement by up to 4.9%, and neglecting distribution constraints overstated the combined distribution-side energy and reserve award by up to 18%. In the expanded study, the largest case was solved in 2.02 s with a 0.00% optimality gap. Time-domain simulations kept the frequency nadir above 59.0 Hz in all tested hours. These results demonstrate the value of fast-response modeling and distribution-feasible reserve delivery in coordinated market clearing.

Noh, Seung-Gil↗

Establishing metrics to quantify spatial similarity in spherical and red blood cell distributions

As computational power increases and systems with millions of red blood cells can be simulated, it is important to note that varying spatial distributions of cells may affect simulation outcomes. Since a single simulation may not represent the ensemble behavior, many different configurations may need to be sampled to adequately assess the entire collection of potential cell arrangements. In order to determine both the number of distributions needed and which ones to run, we must first establish methods to identify well-generated, randomly placed cell distributions and to quantify distinct cell configurations. We utilize metrics to assess (1) the presence of any underlying structure to the initial cell distribution and (2) similarity between cell configurations. We propose the use of the radial distribution function to identify long-range structure in a cell configuration and apply it to a randomly distributed and structured set of red blood cells. To quantify spatial similarity between two configurations, we make use of the Jaccard index, and characterize sets of red blood cell and sphere initializations. As an extension to our work submitted to the International Conference on Computational Science, we significantly increase our data set size from 72 to 1048 cells, include a similar set of studies using spheres, compare the effects of varying sphere size, and utilize the Jaccard index distribution to probe sets of extremely similar configurations. Our results show that the radial distribution function can be used as a metric to determine long-range structure in both distributions of spheres and RBCs. We determine that the ideal case of spheres within a cube versus bi-concave shaped cells within a cylinder affects the shape of the Jaccard index distributions, as well as the range of Jaccard values, showing that both the shape of particle and the domain may play a role. Furthermore, we also find that the distribution is able to capture very similar configurations through Jaccard index values greater than 95% when appending several nearly identical configurations into the data set.

59 BASIC BIOLOGICAL SCIENCES↗

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydroacoustic Evaluation of Juvenile Salmonid Passage and Distribution at Detroit Dam, 2011

Pacific Northwest National Laboratory evaluated juvenile salmonid passage and distribution at Detroit Dam (DET) on the North Santiam River, Oregon for the U.S. Army Corps of Engineers (USACE) to provide data to support decisions on long-term measures to enhance downstream passage at DET and others dams in USACE’s Willamette Valley Project. This study was conducted in response to regulatory requirements necessitated by the listing of Upper Willamette River Spring Chinook salmon (Oncorhynchus tshawytscha) and Upper Willamette River steelhead (O. mykiss) as threatened under the Endangered Species Act. The goal of the study was to provide information of juvenile salmonid passage and distribution at DET from February 2011 through February 2012. The results of the hydroacoustic study provide new and, in some cases, first-ever data on passage estimates, run timing, distributions, and relationships between fish passage and environmental variables at the dam. This information will inform management decisions on the design and development of surface passage and collection devices to help restore Chinook salmon populations in the North Santiam River watershed above DET. During the entire study period, an estimated total of 182,526 smolt-size fish (±4,660 fish, 95% CI) passed through turbine penstock intakes. Run timing peaked in winter and early spring months. Passage rates were highest during late fall, winter and early spring months and low during summer. Horizontal distribution for hours when both turbine units were operated simultaneously indicated Unit 2 passed almost twice as much fish as Unit 1. Diel distribution for smolt-size fish during the study period was fairly uniform, indicating fish were passing the turbines at all times of the day. A total of 5,083 smolt-size fish (± 312 fish, 95% CI) were estimated passed via the spillway when it was open between June 23 and September 27, 2011. Daily passage was low at the spillway during the June-August period, and increased somewhat in September 2011. When the spillway was operated simultaneously with the turbines, spillway efficiency (efficiency is estimated as spillway passage divided by total project passage) was 0.72 and effectiveness (fish:flow ratio—proportion fish passage at a route (e.g., spillway) divided by proportion water through that route out of the total project) was 2.69. That is, when the spillway was open, 72% of the fish passing the dam used the spillway and 28% passed into the turbine penstocks. Diel distribution for smolt-size fish at the spillway shows a distinct peak in passage between mid-morning and mid-afternoon and low passage at night. We estimated that 23,339 smolt-size fish (± 572 fish, 95% CI) passed via the Regulating Outlet (RO) when it was open from October 29 through November 12, 2011, January 2-6, and January 20 through February 3, 2012. During the October–November period, RO passage peaked at 1,086 fish on November 5, with a second peak on November 7 (1,075 fish). When the RO was operated simultaneously with the turbines, RO efficiency was 0.33 and effectiveness was 0.89. In multiple regression analyses, a relatively parsimonious model was selected that predicted the observed fish passage data well. The best model included forebay temperature at depth, forebay elevation, total discharge, hours of daylight, and the operation period. The vertical distribution of fish in the forebay near the face of the dam where the transducers sampled showed fish were generally distributed throughout the water column during all four operational periods. During the refill and full pool periods, vertical distribution was bi-modal with surface-layer and mid-water modes. Patterns for day and night distributions were variable. Fish were distributed above and below the thermocline when it was present (full pool and drawdown periods).

Hydroacoustics↗

Software-Based Challenges of Developing the Future Distribution Grid

The software that the utility industry currently uses may be insufficient to analyze the distribution grid as it rapidly modernizes to include active resources such as distributed generation, switch and voltage control, automation, and increasingly complex loads. Although planners and operators have traditionally viewed the distribution grid as a passive load, utilities and consultants increasingly need enhanced analysis that incorporates active distribution grid loads in order to ensure grid reliability. Numerous commercial and open-source tools are available for analyzing distribution grid systems. These tools vary in complexity from providing basic load-flow and capacity analysis under steady-state conditions to time-series analysis and even geographical representations of dynamic and transient events. The need for each type of analysis is not well understood in the industry, nor are the reasons that distribution analysis requires different techniques and tools both from those now available and from those used for transmission analysis. In addition, there is limited understanding of basic capability of the tools and how they should be practically applied to the evolving distribution system. The study reviews the features and state of the art capability of current tools, including usability and visualization, basic analysis functionality, advanced analysis including inverters, and renewable generation and load modeling. We also discuss the need for each type of distribution grid system analysis. In addition to reviewing basic functionality current models, we discuss dynamics and transient simulation in detail and draw conclusions about existing software?s ability to address the needs of the future distribution grid as well as the barriers to modernization of the distribution grid that are posed by the current state of software and model development. Among our conclusions are that accuracy, data transfer, and data processing abilities are key to future distribution grid modeling, and measured data sources are a key missing element . Modeling tools need to be calibrated based on measured grid data to validate their output in varied conditions such as high renewables penetration and rapidly changing topology. In addition, establishing a standardized data modeling format would enable users to transfer data among tools to take advantage of different analysis features. ?

99 GENERAL AND MISCELLANEOUS↗

The Value of Distributed Wind: A Valuation Framework

While there has been an increase in distributed energy resource (DER) adoption, particularly for solar photovoltaic (PV) DERs, distributed wind deployment is still relatively low. However, distributed wind has significant growth potential for behind-the-meter applications and can provide improved resource diversity and resilience in high-DER grid systems. In order to make optimized decisions about distributed wind deployment, it is important to identify and characterize the value streams of distributed wind in a way that will help policy makers, microgrid developers, utilities, project developers, and others understand the potential value of distributed wind. In our paper, we develop a framework to provide a consistent and comprehensive approach to finding the net benefits of distributed wind to various stakeholders in different scenarios. We do this by considering the main use cases for distributed wind, and then developing a framework which is flexible for these various use cases and the relevant stakeholders while maintaining the same overall structure for analysis. Ongoing work includes applying this framework to valuing two distributed wind power systems, one in a remote village in Alaska and one in a front-of-the-meter deployment in Iowa. Our analysis aims to understand the full value that distributed wind is providing in these systems, including the resilience benefits it is contributing. As a better understanding of distributed wind valuation is gained, it will be easier for stakeholders to make informed decisions in order to have resilient renewable energy systems.

17 WIND ENERGY↗

A Tool for Automatic Data Distribution for CFD Applications on Structured Grids

Development of HPF versions of NPB and ARC3D has shown that HPF provides an efficient, concise way to express parallelism and to organize data traffic. The use of HPF, as noted in the papers, requires an intimate knowledge of the applications and a detailed analysis of data affinity, data movement, and data granularity. To simplify and accelerate the task of developing HPF versions of existing CFD applications we have designed and implemented ADAPT (Automatic Data Alignment and Placement Tool). ADAPT analyzes a CFD application working on a single structured grid and generates HPF TEMPLATE, (RE)DISTRIBUTION, ALIGNMENT, and INDEPENDENT directives. The directives can be generated on the nest level, subroutine level, application level, or on the application interface level. ADAPT annotates an existing CFD FORTRAN application, performing computations on single or multiple grids. On each grid the application is considered as a sequence of operators, each applied to a set of variables defined in a particular grid domain. ADAPT automatically detects implicit operators (i.e., having data dependences) and explicit operators (without data dependences). For parallelization of an explicit operator ADAPT creates a template for the operator domain, aligns arrays used in the operator with the template, distributes the template, and declares the loops over the distributed dimensions as INDEPENDENT. For parallelization of an implicit operator, the distribution of the operator's domain should be consistent with the operator's dependences. Any dependence between sections distributed on different processors would preclude parallelization if the compiler does not have an ability to pipeline computations. If a data distribution is "orthogonal" to the dependences of an implicit operator, then the loop which implements the operator can be declared as INDEPENDENT. ADAPT starts with an analysis of array index expressions of the loop nests. For each pair of arrays referenced in an assignment statement, it generates an arc in the alignment graph and annotates it with an affinity relation. The template, alignment, and distribution directives for a particular loop nest are then derived from a transitive closure of the affinity relation. A compromise of data distributions in different nests and subroutines is achieved by merging annotated alignment graphs for adjacent nests/stibroutine calls in the nest/call graph of the application in the process called distribution lifting. ADAPT has been implemented as a C++ program running in conjunction with a parallelization tool called CAPTools. ADAPT uses the parse tree, interprocedural analysis and application database generated by CAPTools. It also uses the Directed Graph class, initially implemented in p2d2 (parallel debugger oi distributed programs), and some other classes supporting symbolic computations. ADAPT uses data distribution techniques described. ADAPT was tested with ARC3D and the FT benchmark and has demonstrated a code performance within a factor of 1.5 of handwritten versions.

Frumkin, Michael↗

Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study

Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.

Convery, P. D.↗

Implementing Access to Data Distributed on Many Processors

A reference architecture is defined for an object-oriented implementation of domains, arrays, and distributions written in the programming language Chapel. This technology primarily addresses domains that contain arrays that have regular index sets with the low-level implementation details being beyond the scope of this discussion. What is defined is a complete set of object-oriented operators that allows one to perform data distributions for domain arrays involving regular arithmetic index sets. What is unique is that these operators allow for the arbitrary regions of the arrays to be fragmented and distributed across multiple processors with a single point of access giving the programmer the illusion that all the elements are collocated on a single processor. Today's massively parallel High Productivity Computing Systems (HPCS) are characterized by a modular structure, with a large number of processing and memory units connected by a high-speed network. Locality of access as well as load balancing are primary concerns in these systems that are typically used for high-performance scientific computation. Data distributions address these issues by providing a range of methods for spreading large data sets across the components of a system. Over the past two decades, many languages, systems, tools, and libraries have been developed for the support of distributions. Since the performance of data parallel applications is directly influenced by the distribution strategy, users often resort to low-level programming models that allow fine-tuning of the distribution aspects affecting performance, but, at the same time, are tedious and error-prone. This technology presents a reusable design of a data-distribution framework for data parallel high-performance applications. Distributions are a means to express locality in systems composed of large numbers of processor and memory components connected by a network. Since distributions have a great effect on the performance of applications, it is important that the distribution strategy is flexible, so its behavior can change depending on the needs of the application. At the same time, high productivity concerns require that the user be shielded from error-prone, tedious details such as communication and synchronization.

James, Mark↗

Conceptual Design of a Counter-Rotating Fan System for Distributed Boundary Layer Ingesting Propulsion

The present paper details the design of the counter rotating fans for a Turboelectric Distributed Propulsion (TeDP) system. Sixteen propulsors installed in mail-slot-shape nacelles are embedded on an aerodynamically optimized hybrid wing-body configuration. The hybrid-wing/body (HWB) configuration which was previously designed to satisfy the conditions of trim, longitudinally static stability and specific cargo space is employed as the baseline configuration in pursuing an optimal distributed propulsion system. A set of distributed propulsors is conceptually designed and the collective performance is evaluated against the target thrust mandated by the mission requirements. The concept of the distributed propulsion allows the fan pressure ratio to be around 1.27~1.32 for the target thrust. In addition, further splitting of the fan pressure ratio by using the counter-rotating fans for each slot realizes the target pressure ratio with low tip speed. In the distributed propulsion system, the nature of the flow conditions and/or the thickness of the ingested boundary layer may differ and result in different propulsive reaction of each individual propulsor. The optimization is, thus, approached from both the propulsion system and individual propulsor perspectives. An optimal distribution of the thrust and power output is determined by how the system utilizes each passage's propulsive characteristics and its interaction with the airframe. These system level analysis and optimization are conducted using an actuator disk model to account for the propulsion-airframe integration numerically. With respect to the propulsor level, aerodynamic shape optimizations of the fan blades are performed in a sequential multi-objective optimization process for various design objectives, such as mass flow rate condition, fan pressure ratio, efficiency and the exit flow angle of the fan stage by using a genetic algorithm, NSGA-II. The radial chord distribution, and meanline distribution of the rotors are designed on the circumferentially averaged axi-symmetric inlet profiles and tested on the six inlet profiles from six divided sectors to reckon flow distortion. The performances of the counter rotating fans are, thus, evaluated accordingly for obtaining distortion tolerant fan. The performance of the distributed propulsion system is evaluated by two CFD tools, i.e., a multi-stage turbo-machinery CFD code and one propulsion-airframe integration flow solver coupled with a body-force model. The optimized boundary layer ingestion propulsion system of 16 distributed slots not only reaches the system target thrust, but also delivers a close to 20% fuel saving benefit against its counterpart 12 distributed clean inlet propulsion system.

Boundary-Layer-Ingestion Propulsion↗

Radionuclide Risk Coefficient Uncertainty Distributions for Use in Performance Assessments - 19156

Traditionally, radiation dose has been used as the regulatory criterion for decision making in performance assessments (PAs), despite the concept of cancer risk arguably being a more easily understood endpoint than dose for the public, and a preferable metric for decision analysis. However, the United States Environmental Protection Agency (EPA) has developed risk parameters for use in risk assessments and PAs in the form of radionuclide-specific risk coefficients, or RCs, which are published in various Federal Guidance Reports (FGRs). RCs are employed as scalars to convert environmental radionuclide exposure (via ingestion, inhalation, and external exposure) estimates to incremental lifetime cancer risk (ILCR). ILCR is the additional cancer risk averaged over a lifetime, over-and-above all-cause cancer risk, associated with carcinogen exposure. Uncertainties associated with RCs are discussed in FGR 13 and in a joint report published by EPA and Oak Ridge National Laboratory, but statistical distributions for the RCs were not directly developed. The uncertainty assessments in the EPA/ONRL report were based on dosimetric modeling and elicitation of nine experts. From these sources, five uncertainty categories were defined based on ratios of the 5. and 95. quantiles of uncertainty. Thus, the representations of uncertainty defined by EPA/ORNL are not distributions in a traditional sense. Furthermore, the authors only addressed uncertainty in RCs for inhalation and ingestion pathways. We develop a methodology for assigning statistical distributions to EPA's published mortality RC values (i.e. to estimate the probability of dying from radiation-induced cancers) for the inhalation (slow, medium, and fast absorption types), ingestion (water, food, and soil), and external exposure pathways associated with all cancer types. We employed the same dosimetric concepts and assumptions as EPA/ORNL, including cancer risk (dose-response) distributional information from a previous expert elicitation, which was also employed by EPA/ORNL in their work. We assumed that uncertainty in RCs can be described by lognormal distributions based on the modeling conducted previously by EPA/ORNL, and that point estimates tabulated in FGR 13 represent median values of the distributions. We also extended EPA/ORNL's uncertainty ratios to statistical distributions for external exposure RCs. The results are a set of radionuclide-specific distributions that are applicable in PAs that focus on cancer mortality. These distributions are currently and will subsequently be employed in PAs at various sites, in addition to the required regulatory dose limits. Quantification of ILCRs and associated uncertainty in PAs will provide regulators and the public with additional information for good decision making. Future directions include further refinement of the distributions and development of statistical distributions for morbidity RCs, which are used to estimate the probability of developing, but not necessarily dying as a result of, radiation-induced cancer. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗