Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional 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 631 records · Page 35

Coupled CFD-PBE Predictions of Renal Stone Size Distributions in the Nephron in Microgravity

In this paper, a deterministic model is developed to assess the risk of critical renal stone formation for astronauts during space travel. A Population Balance Equation (PBE) model is used to compute the size distribution of a population of nucleating, growing and agglomerating renal calculi as they are transported through different sections of the nephron. The PBE model is coupled to a Computational Fluid Dynamics (CFD) model that solves for steady state flow of urine and transport of renal calculi along with the concentrations of ionic species, calcium and oxalate, in the nephron using an Eulerian two-phase mathematical framework. Parametric simulation are performed to study stone size enhancement and steady state volume fraction distributions in the four main sections of the nephron under weightlessness conditions. Contribution of agglomeration to the stone size distribution and effect of wall friction on the stone volume fraction distributions are carefully examined. Case studies using measured astronaut urinary calcium and oxalate concentrations in microgravity as input indicate that under nominal conditions the largest stone sizes developed in Space will be still considerably below the critical range for problematic stone development. However, results also indicate that the highest stone volume fraction occurs next to the tubule and duct walls. This suggests that there is an increased potential for wall adhesion with the possibility of evolution towards critical stone sizes.

Crystal Growth↗

Watershed Modeling with Remotely Sensed Big Data: MODIS Leaf Area Index Improves Hydrology and Water Quality Predictions

Traditional watershed modeling often overlooks the role of vegetation dynamics. There is also little quantitative evidence to suggest that increased physical realism of vegetation dynamics in process-based models improves hydrology and water quality predictions simultaneously. In this study, we applied a modified Soil and Water Assessment Tool (SWAT) to quantify the extent of improvements that the assimilation of remotely sensed Leaf Area Index (LAI) would convey to streamflow, soil moisture, and nitrate load simulations across a 16,860 km2 agricultural watershedin the midwestern United States. We modified the SWAT source code to automatically override the model’s built-in semiempirical LAI with spatially distributed and temporally continuous estimates from Moderate Resolution Imaging Spectroradiometer (MODIS). Compared to a “basic” traditional model with limited spatial information, our LAI assimilation model (i) significantly improved daily streamflow simulations during medium-to-low flow conditions, (ii) provided realistic spatial distributions of growing season soil moisture, and (iii) substantially reproduced the long-term observed variability of daily nitrate loads. Further analysis revealed that the overestimation or underestimation of LAI imparted a proportional cascading effect on how the model partitions hydrologic fluxes and nutrient pools. As such, assimilation of MODIS LAI data corrected the model’sLAI overestimation tendency, which led to a proportionally increased rootzone soil moisture and decreased plant nitrogen uptake. With these new findings, our study fills the existing knowledge gap regarding vegetation dynamics in watershed modeling and confirms that assimilation of MODIS LAI data in watershed models can effectively improve both hydrology and water quality predictions.

Adnan Rajib↗

Local and utility-wide cost allocations for a more equitable wildfire-resilient distribution grid

Climate-induced extreme weather conditions make electricity infrastructure more vulnerable. They increase the risk of power-line-ignited wildfires which can, in turn, jeopardize electric power delivery. Here, leveraging machine learning, we show that lower-income communities in California not only have lower fractions of power distribution lines undergrounded, but overhead lines and poles in their neighbourhoods are also more vulnerable to wildfires. Should they bear the cost of undergrounding fire-prone lines themselves, they would have to pay a disproportionately higher cost per household. We propose a cost allocation scheme with an income threshold below which the cost is borne by utility-wide ratepayers and above which the cost is borne locally. This scheme can not only minimize the average of undergrounding costs per household as a share of income, but also homogenize such cost–income ratios across communities. Furthermore, our research demonstrates the opportunity to appropriately integrate existing policies to make electricity infrastructure affordable, equitable and reliable amidst climate change.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The influence of submarine canyons-related processes on recent benthic foraminiferal distribution, Espírito Santo Basin, Southeastern Brazil

The complex topography of submarine canyons may result in different composition of benthic foraminifera assemblages. To understand how trophic, hydrological and sedimentological conditions in submarine canyons can influence the distribution of benthic foraminifera, and to use this information to corroborate paleoenvironmental interpretations for the Holocene, we investigated recent benthic foraminiferal assemblages (total fauna >63 μm) and sedimentological data in two canyons (CANWN and CAND) in the Espírito Santo Basin (ESB) between 18°20' and 21°20' S. Surface sediment samples (0–2 cm) were collected inside the canyons (150 to 1300 m water depth) and in shelf-slope adjacent transects (50 to 1300 m water depth). The density, taxonomic diversity, and assemblage composition of benthic foraminifera change with depth and location. The distinct ecological preferences of the most abundant taxa allowed us to recognize five benthic foraminiferal groups. Three groups (V, III, and I) are present in different bathymetric sectors; Group V: outer shelf (50 m), Group III: upper, and middle – lower slope (150 to 400 m), and Group I: middle – lower slope (1000 to 1300 m). Groups II and IV show no characteristic bathymetric distribution and are present only in CAND and in CANWN, respectively. Group V consists of Hanzawaia boueana, Peneroplis planatus, and Quinqueloculina lamarckiana; Group III is dominated by Globocassidulina rossensis and Trifarina spp.; Group I consists of Globocassidulina crassa, Bolivina lowmani, Gavelinopsis versiformis, Alabaminella weddellensis, and Epistominella exigua. The main species in Group II (CAND, 150, 1000 to 1300 m) are Trifarina angulosa, Globocassidulina subglobosa, and Discorbis vilardeboanus. Group IV (middle – lower CANWN, 1000 to 1300 m), consists mainly of agglutinated species Glomospira charoides, Rhabdammina abyssorum, and Psammosphaera fusca. Further, our data suggest that the quantity (and quality) of food supply, hydrodynamic conditions and sediment properties are the main drivers controlling the bathymetric distribution of benthic foraminiferal assemblages in both canyons. The middle – lower CANNW revealed unstable trophic conditions, related to terrigenous sediment input due to turbidity currents. In CAND, the foraminiferal assemblages ecology indicated sufficient organic matter supply that favors species establishment and diversity, indicating a more productive and less unstable environment than in CANWN.

58 GEOSCIENCES↗

Computer-aided Abnormality Detection in Chest Radiographs in a Clinical Setting via Domain-adaptation

Deep learning (DL) models are being deployed at medical centers to aid radiologists for diagnosis of lung conditions from chest radiographs. Such models are often trained on a large volume of publicly available labeled radiographs. These pre-trained DL models’ ability to generalize in clinical settings is poor because of the changes in data distributions between publicly available and privately held radiographs. In chest radiographs, the heterogeneity in distributions arises from the diverse conditions in X-ray equipment and their configurations used for generating the images. In the machine learning community, the challenges posed by the heterogeneity in the data generation source is known as domain shift, which is a mode shift in the generative model. In this work, we introduce a domain-shift detection and removal method to overcome this problem. Our experimental results show the proposed method’s effectiveness in deploying a pre-trained DL model for abnormality detection in chest radiographs in a clinical setting.

Dubey, Abhishek↗

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

Cross polarization in beam waveguide-fed Cassegrain reflector antennas

The sensitivity of the cross-polarization level to a deviation from the geometrical condition derived by Mizusawa and Kitsuregawa (1973) for different geometrical configurations is studied. This condition restricts the number of possible beam waveguide configurations for beam waveguide-fed Cassegrain reflector (BFCR) antennas. For a symmetrical feed, this condition results in a symmetrical aperture distribution with no cross-polarized component. By examining a number of beam waveguide configurations satisfying the condition, it was observed that for linearly polarized feed, the cross-polarization level is very sensitive to a deviation from this condition. For circularly polarized feed, deviation from this condition does not increase the cross-polarization level; however, it results in the squinting of the beam for BFCRs.

Houshmand, Bijan↗

Bimodal SLD Ice Accretion on a CRM Midspan Wing Section Model

An ice shape database has been created to document ice accretions on a 65 per cent scale version of the mid-span section of the Common Research Model resulting from an exposure to a Supercooled Large Drop (SLD) icing cloud with a bimodal drop size distribution. The ice shapes created were documented with photographs, laser scanned surface measurements over a section of the model span, and measurement of the ice mass over the same section of each accretion. Drop distribution effects were evaluated by using the same IRT icing conditions except with either monomodal or bimodal drop distributions. Ice shapes resulting from the bimodal distribution as well as from equivalent monomodal drop size distributions were obtained and compared. Results indicate that the ice shapes resulting from the monomodal and bimodal drop size distributions had equivalent mass and volume, as expected with this approach. Note: No audio portion for this record, additionally slide 9 has a playable video but must be downed in order to view.

Icing↗

Simulated Beam Extraction Performance Characterization of a 50-cm Ion Thruster Discharge

A 50 cm ion thruster is being developed to operate at >65 percent total efficiency at 11 kW, 2700 s Isp and over 25 kW, 4500 s Isp at a total efficiency of >75 percent. The engine is being developed to address the need for a multimode system that can provide a range of thrust-to- power to service national and commercial near-earth onboard propulsion needs such as station-keeping and orbit transfer. Operating characteristics of the 50 cm ion thruster were measured under simulated beam extraction. The discharge current distribution at the various magnet rings was measured over a range of operating conditions. The relationship between the anode current distribution and the resulting plasma uniformity and ion flux measured at the thruster exit plane is discussed. The thermal envelope will also be investigated through the monitoring of magnet temperatures over the range of discharge powers investigated. Discharge losses as a function of propellant utilization was also characterized at multiple simulated beam currents. Bulk plasma conditions such as electron temperature and electron density near engine centerline was measured over a range of operating conditions using an internal Langmuir probe. Sensitivity of discharge performance to chamber length is also discussed. This data acquired from this discharge study will be used in the refinement of a throttle table in anticipation for eventual beam extraction testing.

Foster, John E.↗

Successful Cleanroom Installation of PIP-II SSR2 Coupler Using Robotic Arm

The Fermilab Side-Coupled Linac accelerates H- beam from 116 MeV to 400 MeV through seven 805 MHz modules. Twelve wire scanners are present in the Side Coupled Linac and four are present in the transfer line between the Linac and the Booster synchrotron ring. These wire scanners act as important diagnostic instruments to directly collect information on the beam s transverse distribution. The manipulation of the conditions of wire scanner data collection enables further characterization of the beamline, such as calculating emittance and the Twiss parameters of the beam at select regions. Here we present the results of these studies and characterization of the non-Gaussian transverse beam distribution observed.

Narug, C.↗

Observing Changing Ecological Diversity in the Anthropocene

As the world enters the Anthropocene, the planet's environment is changing rapidly, putting critical ecosystem services at risk. Understanding and forecasting how ecosystems will change over the coming decades requires understanding the sensitivity of species to environmental change. The extant distribution of species and functional groups contains valuable information about the performance of different species in different environments. However, with high rates of environmental change, information inherent in ranges of many species will disappear, since that information exists only under quasi-equilibrium conditions. The information content of distributional data obtained now is greater than data obtained in the future. New remote sensing technologies can map chemical and structural traits of plant canopies and allow inference of trait and in many cases, species ranges. Current satellite remote sensing data can only produce relatively simple classifications, but new techniques have dramatically higher biological information content.

environmental changes↗

Early Failure of Lithium–Sulfur Batteries at Practical Conditions: Crosstalk between Sulfur Cathode and Lithium Anode

Lithium–sulfur (Li–S) batteries are one of the most promising next-generation energy storage technologies due to their high theoretical energy and low cost. However, Li–S cells with practically high energy still suffer from a very limited cycle life with reasons which remain unclear. Here, through cell study under practical conditions, it is proved that an internal short circuit (ISC) is a root cause of early cell failure and is ascribed to the crosstalk between the S cathode and Li anode. The cathode topography affects S reactions through influencing the local resistance and electrolyte distribution, particularly under lean electrolyte conditions. The inhomogeneous reactions of S cathodes are easily mirrored by the Li anodes, resulting in exaggerated localized Li plating/stripping, Li filament formation, and eventually cell ISC. Manipulating cathode topography is proven effective to extend the cell cycle life under practical conditions. The findings of this work shed new light on the electrode design for extending cycle life of high-energy Li–S cells, which are also applicable for other rechargeable Li or metal batteries.

25 ENERGY STORAGE↗

A Particle Method for the Multispecies Landau Equation

Abstract The multispecies Landau collision operator describes the two-particle, small scattering angle or grazing collisions in a plasma made up of different species of particles such as electrons and ions. Recently, a structure preserving deterministic particle method (Carrillo et al. in J. Comput. Phys. 7:100066, 2020) has been developed for the single species spatially homogeneous Landau equation. This method relies on a regularization of the Landau collision operator so that an approximate solution, which is a linear combination of Dirac delta distributions, is well-defined. Based on a weak form of the regularized Landau equation, the time dependent locations of the Dirac delta functions satisfy a system of ordinary differential equations. In this work, we extend this particle method to the multispecies case, and examine its conservation of mass, momentum, and energy, and decay of entropy properties. We show that the equilibrium distribution of the regularized multispecies Landau equation is a Maxwellian distribution, and state a critical condition on the regularization parameters that guarantees a species independent equilibrium temperature. A convergence study comparing an exact multispecies Bobylev-Krook-Wu (BKW) solution to the particle solution shows approximately 2nd order accuracy. Important physical properties such as conservation, decay of entropy, and equilibrium distribution of the particle method are demonstrated with several numerical examples.

Mathematics↗

Heavy ion irradiation effects on CrFeMnNi and AlCrFeMnNi high entropy alloys

Co-free but Al-included medium/high entropy alloys (M/HEAs) have gained increasing interests due to their lower cost and the potential to tune the multi-phase microstructure. The irradiation response of two Co-free HEAs, face-centered cubic (FCC) CrFeMnNi with limited Cr enriched α' phase and body-centered cubic (BCC) AlCrFeMnNi with B2s phase and nanoprecipitates were explored. Ion irradiations using 5 MeV Fe 2+ ions were performed at 500°C to a peak fluence of 50 and/or 100 displacements per atom (dpa). In dual-phase AlCrFeMnNi, there was no significant radiation induced segregation or chemical intermixing at the coherent matrix (FeCrMn-rich)/second phase (AlNi-rich) boundaries. In CrFeMnNi, limited voids were only detected at the peak damage location of ~ 50 dpa. On the other hand, voids were widely distributed in AlCrFeMnNi: under 50 and 100 dpa irradiation conditions, voids were found with larger dimension and denser distribution in the FeCrMn-rich matrix, smaller and slightly lower density in an AlNi-rich second phase. In addition, the diameter of the FeCMn-rich nanoprecipitates didn't reveal any tendency of dissolution or growth. This is correlated with their superior structural stability against irradiation. Significant radiation-induced hardening (increases from 3.8 ± 0.2 GPa to 4.7 ± 0.6 GPa) was measured in CrFeMnNi, but only ~ 4% hardness increase (from 7.4 ± 0.8 GPa to 7.7 ± 0.4 GPa) was noted in AlCrFeMnNi. Finally, in addition to the radiation-induced defects, such as voids, dislocation loops and point defects, other factors, such as chemical short-range ordering may play an important role.

36 MATERIALS SCIENCE↗

Predicting Emission Source Terms in a Reduced-Order Fire Spread Model—Part 1: Particulate Emissions

A simple, easy-to-evaluate, surrogate model was developed for predicting the particle emission source term in wildfire simulations. In creating this model, we conceptualized wildfire as a series of flamelets, and using this concept of flamelets, we developed a one-dimensional model to represent the structure of these flamelets which then could be used to simulate the evolution of a single flamelet. A previously developed soot model was executed within this flamelet simulation which could produce a particle size distribution. Executing this flamelet simulation 1200 times with varying conditions created a data set of emitted particle size distributions to which simple rational equations could be tuned to predict a particle emission factor, mean particle size, and standard deviation of particle sizes. These surrogate models (the rational equation) were implemented into a reduced-order fire spread model, QUIC-Fire. Using QUIC-Fire, an ensemble of simulations were executed for grassland fires, southeast U.S. conifer forests, and western mountain conifer forests. Resulting emission factors from this ensemble were compared against field data for these fire classes with promising results. Also shown is a predicted averaged resulting particle size distribution with the bulk of particles produced to be on the order of 1 μm in size.

54 ENVIRONMENTAL SCIENCES↗

Mechanically fastened composite laminates subjected to combined bearing-bypass and shear loading

Bolts and rivets provide a means of load transfer in the construction of aircraft. However, they give rise to stress concentrations and are often the source and location of static and fatigue failures. Furthermore, fastener holes are prone to cracks during take-off and landing. These cracks present the most common origin of structural failures in aircraft. Therefore, accurate determination of the contact stresses associated with such loaded holes in mechanically fastened joints is essential to reliable strength evaluation and failure prediction. As the laminate is subjected to loading, the contact region, whose extent is not known, develops between the fastener and the hole boundary through this contact region, which consists of slip and no-slip zones due to friction. The presence of the unknown contact stress distribution over the contact region between the pin and the composite laminate, material anisotropy, friction between the pin and the laminate, pin-hole clearance, combined bearing-bypass and shear loading, and finite geometry of the laminate result in a complex non-linear problem. In the case of bearing-bypass loading in compression, this non-linear problem is further complicated by the presence of dual contact regions. Previous research concerning the analysis of mechanical joints subjected to combined bearing-bypass and shear loading is non-existent. In the case of bearing-bypass loading only, except for the study conducted by Naik and Crews (1991), others employed the concept of superposition which is not valid for this non-linear problem. Naik and Crews applied a linear finite element analysis with conditions along the pin-hole contact region specified as displacement constraint equations. The major shortcoming of this method is that the variation of the contract region as a function of the applied load should be known a priori. Also, their analysis is limited to symmetric geometry and material systems, and frictionless boundary conditions. Since the contact stress distribution and the contact region are not known a priori, they did not directly impose the boundary conditions appropriate for modelling the contact and on-contact regions between the fastener and the hole. Furthermore, finite element analysis is not suitable for iterative design calculations for optimizing laminate construction in the presence of fasteners under complex loading conditions. In this study, the solution method developed by Madenci and Ileri (1992a,b) has been extended to determine the contact stresses in mechanical joints under combined bearing-bypass and shear loading, and bearing-bypass loading in compression resulting in dual contact regions.

Madenci, Erdogan↗

Multiarea Distribution System State Estimation via Distributed Tensor Completion

Here, this paper proposes a model-free distribution system state estimation method based on tensor completion using canonical polyadic decomposition. In particular, we consider a setting where the network is divided into multiple areas. The measured physical quantities at buses located in the same area are processed by an area controller. A three-way tensor is constructed to collect these measured quantities. The measurements are analyzed locally to recover the full state information of the network. A distributed closed-form iterative algorithm based on the alternating direction method of multipliers is developed to obtain the low-rank factors of the whole network state tensor where information exchange happens only between neighboring areas. The convergence properties of the distributed algorithm and the sufficient conditions on the number of samples for each smaller network that guarantee the identifiability of the factors of the state tensor are presented. To demonstrate the efficacy of the proposed algorithm and to check the identifiability conditions, numerical simulations are carried out using the IEEE 123-bus system and a large-scale real utility feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗