Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed method”

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 595 records · Page 33

Distributionally Robust Decision Making Leveraging Conditional Distributions

Distributionally robust optimization (DRO) is a powerful tool for decision making under uncertainty. It is particularly appealing because of its ability to leverage existing data. However, many practical problems call for decision-making with some auxiliary information, and DRO in the context of conditional distributions is not straightforward. We propose a conditional kernel distributionally robust optimization (CKDRO) method that enables robust decision making under conditional distributions through kernel DRO and the conditional mean operator in the reproducing kernel Hilbert space (RKHS). In particular, we consider problems where there is a correlation between the unknown variable y and an auxiliary observable variable x. Given past data of the two variables and a queried auxiliary variable, CKDRO represents the conditional distribution P(y|x) as the conditional mean operator in the RKHS space and quantifies the ambiguity set in the RKHS as well, which depends on the size of the dataset as well as the query point. To justify the use of RKHS, we demonstrate that the ambiguity set defined in RKHS can be viewed as a ball under a metric that is similar to the Wasserstein metric. The DRO is then dualized and solved via a finite dimensional convex program. The proposed CKDRO approach is applied to a generation scheduling problem and shows that the result of CKDRO is superior to common benchmarks in terms of quality and robustness.

Chen, Yuxiao↗

Methods for Incorporating Model Uncertainty into Exoplanet Atmospheric Analysis

A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of Calculated and Experimental Temperatures and Coolant Pressure Losses for a Cascade of Small Air-Cooled Turbine Rotor Blades

Average spanwise blade temperatures and cooling-air pressure losses through a small (1.4-in, span, 0.7-in, chord) air-cooled turbine blade were calculated and are compared with experimental nonrotating cascade data. Two methods of calculating the blade spanwise metal temperature distributions are presented. The method which considered the effect of the length-to-diameter ratio of the coolant passage on the blade-to-coolant heat-transfer coefficient and assumed constant coolant properties based on the coolant bulk temperature gave the best agreement with experimental data. The agreement obtained was within 3 percent at the midspan and tip regions of the blade. At the root region of the blade, the agreement was within 3 percent for coolant flows within the turbulent flow regime and within 10 percent for coolant flows in the laminar regime. The calculated and measured cooling-air pressure losses through the blade agreed within 5 percent. Calculated spanwise blade temperatures for assumed turboprop engine operating conditions of 2000 F turbine-inlet gas temperature and flight conditions of 300 knots at a 30,000-foot altitude agreed well with those obtained by the extrapolation of correlated experimental data of a static cascade investigation of these blades.

Stepka, Francis S↗

Satellite Remote Sensing Estimation of River Discharge: Application to the Yukon River Alaska

A methodology based on general hydraulic relations for rivers has been developed to estimate the discharge (flow rate) of rivers using satellite remote sensing observations. The estimates of discharge, flow depth, and flow velocity are derived from remotely observed water surface area, water surface slope, and water surface height, and demonstrated for two reaches of the Yukon River in Alaska, at Eagle (reach length 34.7 km) and near Stevens Village (reach length 38.3 km). The method is based on fundamental equations of hydraulic flow resistance in rivers, including the Manning equation and the Prandtl-von Karman universal velocity distribution equation. The method employs some new hydraulic relations to help define flow resistance and height of the zero flow boundary in the channel. Estimates are made both with and without calibration. The water surface area of the river reach is measured by using a provisional version of the U.S. Geological Survey (USGS) Landsat based product named Dynamic Surface Water Extent (DSWE). The water surface height and slope measurements require a self-consistent datum, and are derived from observations from the Jason-2 satellite altimeter mission. At both reach locations, the Jason-2 radar altimeter non-winter heights consistently tracked the stage recorded at USGS streamgages with a standard deviation of differences (error) during the non-winter periods of less than 7%. Part of the error may be due to differences in the gage and altimeter crossing locations with respect to the range of stage change and the response to changes in discharge at the upstream and downstream locations. For the non-winter periods, the radar derived slope estimates (mean = 0.0003) were constant over the mission lifetime, and in agreement with previously measured USGS water surface slopes and slopes determined from USGS topographic maps. The accuracy of the mean of the uncalibrated daily estimates of discharge varied between reaches, ranging from 13% near Stevens Village (N = 90) to −21% at Eagle (N = 246) based on the absolute error, and 5% to −6% based on the error of the log of the estimates. Calibrating to the mean of USGS daily discharge estimates from the streamflow rating for the same period of record at each streamgage resulted in mean absolute errors ranging from 1% to 2%, and log errors ranging from 1% or less. The error pattern of the estimates shows that without calibration, even though the mean is well simulated, the high and low end values over the range of estimates may have significant bias.

Bjerklie, David M.↗

Comparative Analysis of Type III Radio Bursts and Solar Flares: Spatial Localization and Correlation with Solar Flare Intensity

We present a comprehensive study of type III radio bursts and their association with solar flares of magnitude M1.0 and larger, as observed by four widely separated spacecraft (Parker Solar Probe, Solar Orbiter, STEREO-A, and Wind). Our main focus is the introduction and validation of two methods for localizing radio bursts using the available multispacecraft data. The first method utilizes intensity fitting with a circular Gaussian distribution, while the second method is based on the time arrival of radio bursts. We demonstrate the effectiveness of these methods through the analysis of a single type III burst event and compare their results with the traditional radio triangulation technique. Furthermore, we conduct a statistical study of 17 type III bursts associated with M- and X-class solar flares in years 2020–2022. Our findings suggest a possible correlation between solar flare intensities and longitudes, with east limb flares tending to be weaker than west limb flares. We also observe a systematic drift of radio burst longitudes toward the east, potentially explained by a poleward component of the local density gradient. Our results suggest a strong correlation between solar flare intensities and radio burst properties, enhancing our understanding of the relationship between solar flares and type III radio bursts.

Radio astronomy↗

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗

Assessment of NASA and RAE viscous-inviscid interaction methods for predicting transonic flow over nozzle afterbodies

The Langley Research Center of the National Aeronautics and Space Administration and the Royal Aircraft Establishment have undertaken a cooperative program to conduct an assessment of their patched viscous-inviscid interaction methods for predicting the transonic flow over nozzle afterbodies. The assessment was made by comparing the predictions of the two methods with experimental pressure distributions and boattail pressure drag for several convergent circular-arc nozzle configurations. Comparisons of the predictions of the two methods with the experimental data showed that both methods provided good predictions of the flow characteristics of nozzles with attached boundary layer flow. The RAE method also provided reasonable predictions of the pressure distributions and drag for the nozzles investigated that had separated boundary layers. The NASA method provided good predictions of the pressure distribution on separated flow nozzles that had relatively thin boundary layers. However, the NASA method was in poor agreement with experiment for separated nozzles with thick boundary layers due primarily to deficiencies in the method used to predict the separation location.

Putnam, L. E.↗

A Nonparametric Method for the Inference of Halo Occupation Distributions

The galaxy–halo connection traces processes by which galaxies form and evolve. The halo occupation distribution (HOD) describes the relationship between galaxies and their host dark matter haloes. Measurements of the galaxy two-point correlation function (2PCF) allow us to extract information about the HODs of observed galaxy samples. Several parametric HOD models have been proposed in the literature, but the choice of parameterization restricts the space of possible HODs. To resolve this issue, we introduce a nonparametric HOD fitting method in which we train an emulator to learn the mappings among the galaxy 2PCF, physical properties used to select galaxy samples, and the HOD, all obtained from simulated past light cones constructed with the Santa Cruz semianalytic model. Implementing this emulator within a likelihood analysis framework, we derive constraints on the HOD of a galaxy sample when provided with a measurement of its 2PCF. Using the emulator to accelerate likelihood evaluations, we test the nonparametric HOD approach on a set of 2PCFs for mock galaxy samples drawn from the TNG100-1 simulation and selected above threshold values of stellar mass and star formation rate. Our framework is able to recover TNG100-1 HODs within 0.2 dex. We use the TNG100-1 mocks to tune the reported uncertainties to estimate those expected in the analysis of observations. Comparing to parametric HOD modelling routines applied to the same mock galaxy samples, our approach consistently infers the HOD with comparable or greater precision and accuracy.

Kennedy, Jacob [Rutgers Univ., Piscataway, NJ (Uni↗

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗

Distribution of Snow and Maximum Snow Water Equivalent Obtained by LANDSAT Data and Degree Day Method

Maximum snow water equivalence and snowcover distribution are estimated using several LANDSAT data taken in snowmelting season over a four year period. The test site is Okutadami-gawa Basin located in the central position of Tohoku-Kanto-Chubu District. The year to year normalization for snowmelt volume computation on the snow line is conducted by year to year correction of degree days using the snowcover percentage within the test basin obtained from LANDSAT data. The maximum snow water equivalent map in the test basin is generated based on the normalized snowmelt volume on the snow line extracted from four LANDSAT data taken in a different year. The snowcover distribution on an arbitrary day in snowmelting of 1982 is estimated from the maximum snow water equivalent map. The estimated snowcover is compared with the snowcover area extracted from NOAA-AVHRR data taken on the same day. The applicability of the snow estimation using LANDSAT data is discussed.

Takeda, K.↗

Ceramic distribution members for solid state electrolyte cells and method of producing

A solid state electrolyte cells apparatus and method of producing is disclosed. The apparatus can be used for separating oxygen from an oxygen-containing feedstock or as a fuel cell for reacting fluids. Cells can be stacked so that fluids can be introduced and removed from the apparatus through ceramic distribution members having ports designed for distributing the fluids in parallel flow to and from each cell. The distribution members can also serve as electrodes to membranes or as membrane members between electrodes. The distribution member design does not contain any horizontal internal ports which allows the member to be thin. A method of tape casting in combination with an embossing method allows intricate radial ribs and bosses to be formed on each distribution member. The bosses serve as seals for the ports and allow the distribution members to be made without any horizontal internal ports.

Clark, Douglas J.↗

Cold dark matter. 1: The formation of dark halos

We use numerical simulations of critically closed cold dark matter (CDM) models to study the effects of numerical resolution on observable quantities. We study simulations with up to 256(exp 3) particles using the particle-mesh (PM) method and with up to 144(exp 3) particles using the adaptive particle-particle-mesh (P3M) method. Comparisons of galaxy halo distributions are made among the various simulations. We also compare distributions with observations, and we explore methods for identifying halos, including a new algorithm that finds all particles within closed contours of the smoothed density field surrounding a peak. The simulated halos show more substructure than predicted by the Press-Schechter theory. We are able to rule out all omega = 1 CDM models for linear amplitude sigma(sub 8) greater than or approximately = 0.5 because the simulations produce too many massive halos compared with the observations. The simulations also produce too many low-mass halos. The distribution of halos characterized by their circular velocities for the P3M simulations is in reasonable agreement with the observations for 150 km/s less than or = V(sub circ) less than or = 350 km/s.

Gelb, James M.↗

Non-destructive spatial analysis of phosphatase activity and total protein distribution in the rhizosphere using a root blotting method

Phosphorus (P) is an essential macronutrient for plant growth, but bioavailable P in soils is often limited due to immobilization resulting from pH and geochemical interactions. Understanding the dynamics of P in soils and elucidating the mechanisms by which plants access P from their environment are critical to evaluating productivity, particularly in nutrient poor environments. Phosphorus from organic matter can act as a major source of P for organisms in soil systems. Phosphatases, enzymes that liberate inorganic P from organic sources, are produced by both plants and microbes and are considered one of the most active classes of enzymes in soil. In this work we developed a root blotting method to spatially image phosphatase activity in the rhizosphere. Proteins from the rhizosphere are transferred to a nitrocellulose membrane while retaining their enzymatic activity and two-dimensional spatial distribution. Subsequent application of a fluorogenic phosphatase indicator, DDAO phosphate, enables visualization of the distribution of phosphatase activity in the sample. The proteins can then be fixed to the membrane and treated with a SYPRO® Ruby gel stain, a fluorescent total protein stain, allowing for visualization of total protein distribution. Taken together, the images of phosphatase activity and total protein localization can be mapped back to the root architecture and provide insight into factors affecting the spatial distribution of enzymatic activity and protein accumulation in the rhizosphere. Notably, this method can be applied to plants growing in rhizoboxes containing soil or soilless growth mixtures (e.g., sand or various potting mixes) and, because of the non-destructive nature of this approach, be performed over time to track changes. We anticipate that this fluorescent indicator imaging technique on root blots can be used in diverse plant-microbe-soil systems to better understand the role of phosphatases in P acquisition and soil P cycling.

59 BASIC BIOLOGICAL SCIENCES↗

Adsorption of argon on sintered tin dioxide analyzed by several methods

Argon adsorption measurements are presented over a wide range of temperature and coverage on a series of three progressively sintered SnO2 surfaces. These data are analyzed by mercury porosimetry, the BET method, the CAEDMON distribution analysis, and the Singleton-Halsey equation. Isosteric heats are computed, and the high-temperature virial expansion of the data presented. The advantages and disadvantages of each method are discussed with particular attention to the ability of physical adsorption to discriminate among surfaces beyond the measurement of surface area.

Hinman, D. C.↗

Multigrid methods and high order finite difference for flow in transition - Effects of isolated and distributed roughness elements

The high order finite difference and multigrid methods have been successfully applied to direct numerical simulation (DNS) for flow transition in 3D channels and 3D boundary layers with 2D and 3D isolated and distributed roughness in a curvilinear coordinate system. A fourth-order finite difference technique on stretched and staggered grids, a fully-implicit time marching scheme, a semicoarsening multigrid method associated with line distributive relaxation scheme, and a new treatment of the outflow boundary condition, which needs only a very short buffer domain to damp all wave reflection, are developed. These approaches make the multigrid DNS code very accurate and efficient. This makes us not only able to do spatial DNS for the 3D channel and flat plate at low computational costs, but also able to do spatial DNS for transition in the 3D boundary layer with 3D single and multiple roughness elements. Numerical results show good agreement with the linear stability theory, the secondary instability theory, and a number of laboratory experiments.

Liu, C.↗

Parameter identification in distributed systems

This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.

Baruh, H.↗

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗