Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution parameters”

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 1,135 records · Page 63

Enhancing Even Gas Distribution in Porous Media with Radial Flow

Packed beds, such as those used in thermal energy storage (TES) systems, typically use flow from one end to the other. This axial flow configuration leads to inefficiencies due to dispersion effects and high pressure drop. Instead of axial flow, this work proposes radial flow in packed beds, where a central tube provides flow that transports fluid from the center to the bed’s wall. Radial flow could be a promising solution to increase the efficiency of charging/discharging processes in TES systems. For instance, studies have been conducted on the thermal behavior during the charging process using radial flow and axial flow, and it was found that radial flow is better than axial flow in terms of thermal performance, where more energy can be stored in the storage tank during the charging process. However, the design of the radial pipe should be optimized to improve even flow distribution. In this work, numerical investigations on different designs were analyzed to enhance even flow gas distribution in the porous media when using this radial technique. This work shows the impact of different parameters on the even flow distribution into the packed bed for two different designs: 1) one radial tube at the center to provide the radial flow along with an annular tube at the wall to receive the flow, and 2) one radial tube at the center and four radial tubes at the wall. Air was used as fluid and 6 mm alumina beads as packing materials. Computational fluid dynamics (CFD) models in COMSOL Multiphysics were used to simulate the behavior of air flow through the piping and packed bed. It was found that the flow into the bed from the radial tube can be affected by different parameters: space between the holes, size/diameter of the holes, number of segments/zones, number of the holes in each zone and the length of each zone.

Beck, David M.↗

Sea Ice Rheology Experiment (SIREx): 1. Scaling and Statistical Properties of Sea-Ice Deformation Fields

As the sea-ice modeling community is shifting to advanced numerical frameworks, developing new sea-ice rheologies, and increasing model spatial resolution, ubiquitous deformation features in the Arctic sea ice are now being resolved by sea-ice models. Initiated at the Forum for Arctic Modeling and Observational Synthesis, the Sea Ice Rheology Experiment (SIREx) aims at evaluating state-of-the-art sea-ice models using existing and new metrics to understand how the simulated deformation fields are affected by different representations of sea-ice physics (rheology) and by model configuration. Part 1 of the SIREx analysis is concerned with evaluation of the statistical distribution and scaling properties of sea-ice deformation fields from 35 different simulations against those from the RADARSAT Geophysical Processor System (RGPS). For the first time, the viscous-plastic (and the elastic-viscous-plastic variant), elastic-anisotropic-plastic, and Maxwell-elasto-brittle rheologies are compared in a single study. We find that both plastic and brittle sea-ice rheologies have the potential to reproduce the observed RGPS deformation statistics, including multi-fractality. Model configuration (e.g., numerical convergence, atmospheric representation, spatial resolution) and physical parameterizations (e.g., ice strength parameters and ice thickness distribution) both have effects as important as the choice of sea-ice rheology on the deformation statistics. It is therefore not straightforward to attribute model performance to a specific rheological framework using current deformation metrics. In light of these results, we further evaluate the statistical properties of simulated Linear Kinematic Features in a SIREx Part 2 companion paper.

54 ENVIRONMENTAL SCIENCES↗

Revisiting Experimental Signatures of the Ponderomotive Force

The classical theory of single-electron dynamics in focused laser pulses is the foundation of both the relativistic ponderomotive force (RPF), which underlies models of laser-collective-plasma dynamics, and the discovery of novel strong-field radiation dynamics. Despite this bedrock importance, consensus eludes the community as to whether acceleration of single electrons in vacuum has been observed in experimental conditions. We analyze an early experiment on the RPF with respect to several features that were neglected in modeling and that can restore consistency between theory predictions and experimental data. The right or wrong pulse profile function, laser parameters, or initial electron distribution can each make or break the agreement between predictions and data. The laser phase at which the electron’s interaction with the pulse begins has a large effect, explaining why much larger energies are achieved by electrons liberated in the focal region by photoionization from high-Z atoms and by electrons ejected from a plasma mirror. Finally, we compute the difference in a typical electron spectrum arising from fluctuating focal spot size in state-of-the-art ultra-relativistic laser facilities. Our results emphasize the importance of thoroughly characterizing laser parameters in order to achieve quantitatively accurate predictions and the precision required for discovery science.

Hegelich, Bjorn Manuel↗

Learning with Adaptive Conservativeness for Distributionally Robust Optimization: Incentive Design for Voltage Regulation: Preprint

Information asymmetry between the Distribution System Operator (DSO) and Distributed Energy Resource Aggregators (DERAs) obstructs designing effective incentives for voltage regulation. To capture this effect, we employ a Stackelberg game-theoretic framework, where the DSO seeks to overcome the information asymmetry and refine its incentive strategies by learning from DERA behavior over multiple iterations. We introduce a model-based online learning algorithm for the DSO, aimed at inferring the relationship between incentives and DERA responses. Given the uncertain nature of these responses, we also propose a distributionally robust incentive design model to control the probability of voltage regulation failure and then reformulate it into a convex problem. This model allows the DSO to periodically revise distribution assumptions on uncertain parameters in the decision model of the DERA. Finally, we present a gradient-based method that permits the DSO to adaptively modify its conservativeness level, measured by the size of a Wasserstein metric-based ambiguity set, according to historical voltage regulation performance. The effectiveness of our proposed method is demonstrated through numerical experiments.

distribution system operator↗

Assessment of Internal Fabric of Particulate Materials

Particle arrangement and distribution within a soil matrix has long been recognized as having significant influence on the mechanical behavior of cohesionless soils. It is well known that two soil specimens having the same grain type (e.g., quartz, feldspar, etc.), same grain size distribution and relative density (or void ratio) can display completely different mechanical behavior. Because of the different fabric configurations in the otherwise similar specimens, they are likely to have different mechanical properties such as stress-strain response, anisotropy, dilatancy, etc. Soil Fabric is defined as the arrangement of particles, particle groups and associated pore space. In the literature, fabric analysis techniques are mainly classified as destructive (e.g., specimen stabilization, thin-sectioning, and microscopy), and nondestructive techniques (e.g., magnetic resonance imaging, ultrasonic testing, x-ray radiography, and computed tomography). Quantifying the void ratio and its distribution is the main parameter used to describe the fabric of particulate materials. This paper presents a comprehensive literature review of fabric analysis techniques applied to particulate materials. In addition, the results of a comprehensive investigation to quantify void ratio of sand specimens will be presented and discussed. The sand used in the experiments in a natural, uniform rounded to sub-rounded silica sand known as F-75 banding sand with mean particle size of 0.22 mm. Uniform specimens and specimens subjected to different axial-strain levels tested under triaxial and biaxial conditions are examined to evaluate void ratio evolution and its distribution using destructive thin-sectioning and nondestructive Computed Tomography (CT) techniques. Details of a new innovative polygon generation technique called Voronoi tessellation used to quantify void ratio of microscopic images of sand grains will be presented and compared to classical Oda's method. Finally, frequency distribution of void ratio calculated using CT analysis will be presented and discussed.

Alshibi, Khalid A.↗

Learning with Adaptive Conservativeness for Distributionally Robust Optimization: Incentive Design for Voltage Regulation

Information asymmetry between the Distribution System Operator (DSO) and Distributed Energy Resource Aggregators (DERAs) obstructs designing effective incentives for voltage regulation. To capture this effect, we employ a Stackelberg game-theoretic framework, where the DSO seeks to overcome the information asymmetry and refine its incentive strategies by learning from DERA behavior over multiple iterations. We introduce a model-based online learning algorithm for the DSO, aimed at inferring the relationship between incentives and DERA responses. Given the uncertain nature of these responses, we also propose a distributionally robust incentive design model to control the probability of voltage regulation failure and then reformulate it into a convex problem. This model allows the DSO to periodically revise distribution assumptions on uncertain parameters in the decision model of the DERA. Finally, we present a gradient-based method that permits the DSO to adaptively modify its conservativeness level, measured by the size of a Wasserstein metric-based ambiguity set, according to historical voltage regulation performance. The effectiveness of our proposed method is demonstrated through numerical experiments.

adaptation models↗

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗

Probabilistic Modeling of the Renal Stone Formation Module

The Integrated Medical Model (IMM) is a probabilistic tool, used in mission planning decision making and medical systems risk assessments. The IMM project maintains a database of over 80 medical conditions that could occur during a spaceflight, documenting an incidence rate and end case scenarios for each. In some cases, where observational data are insufficient to adequately define the inflight medical risk, the IMM utilizes external probabilistic modules to model and estimate the event likelihoods. One such medical event of interest is an unpassed renal stone. Due to a high salt diet and high concentrations of calcium in the blood (due to bone depletion caused by unloading in the microgravity environment) astronauts are at a considerable elevated risk for developing renal calculi (nephrolithiasis) while in space. Lack of observed incidences of nephrolithiasis has led HRP to initiate the development of the Renal Stone Formation Module (RSFM) to create a probabilistic simulator capable of estimating the likelihood of symptomatic renal stone presentation in astronauts on exploration missions. The model consists of two major parts. The first is the probabilistic component, which utilizes probability distributions to assess the range of urine electrolyte parameters and a multivariate regression to transform estimated crystal density and size distributions to the likelihood of the presentation of nephrolithiasis symptoms. The second is a deterministic physical and chemical model of renal stone growth in the kidney developed by Kassemi et al. The probabilistic component of the renal stone model couples the input probability distributions describing the urine chemistry, astronaut physiology, and system parameters with the physical and chemical outputs and inputs to the deterministic stone growth model. These two parts of the model are necessary to capture the uncertainty in the likelihood estimate. The model will be driven by Monte Carlo simulations, continuously randomly sampling the probability distributions of the electrolyte concentrations and system parameters that are inputs into the deterministic model. The total urine chemistry concentrations are used to determine the urine chemistry activity using the Joint Expert Speciation System (JESS), a biochemistry model. Information used from JESS is then fed into the deterministic growth model. Outputs from JESS and the deterministic model are passed back to the probabilistic model where a multivariate regression is used to assess the likelihood of a stone forming and the likelihood of a stone requiring clinical intervention. The parameters used to determine to quantify these risks include: relative supersaturation (RS) of calcium oxalate, citrate/calcium ratio, crystal number density, total urine volume, pH, magnesium excretion, maximum stone width, and ureteral location. Methods and Validation: The RSFM is designed to perform a Monte Carlo simulation to generate probability distributions of clinically significant renal stones, as well as provide an associated uncertainty in the estimate. Initially, early versions will be used to test integration of the components and assess component validation and verification (V&V), with later versions used to address questions regarding design reference mission scenarios. Once integrated with the deterministic component, the credibility assessment of the integrated model will follow NASA STD 7009 requirements.

kidney stones↗

Model-Free Voltage Control of Active Distribution System with PVs Using Surrogate Model-Based Deep Reinforcement Learning

Accurate knowledge of the distribution system topology and parameters is required to achieve good voltage control performance, but this is difficult to obtain in practice. This paper proposes a physical-model-free voltage control method based on a surrogate-model-enabled deep reinforcement learning approach. Specifically, a surrogate model is trained in a supervised manner using the recorded limited number of historical data to learn the relationship between the power injections and voltage fluctuations of each node. Then, the deep reinforcement learning algorithm is applied to learn an optimal control strategy from the experiences obtained by continuous interactions with the surrogate model. The proposed method can achieve physical-model-free control of unbalanced distribution network and inform real-time decisions to deal with fast voltage fluctuations caused by the rapid variation of PV generation. Simulation results on an unbalance IEEE 123-bus system show that the proposed method can achieve similar performance as that of perfect physical-model-based approaches while being advantageous over other traditional methods.

active distribution network↗

Effect of structure on current and potential distributions in porous electrode

Porous electrodes generally contain constricted macropores and localized micropores. The effects of the macropore constrictions on the resistance of a capillary were studied and an analytical model was developed for predicting the current distribution in a constricted macropore which directly includes constriction effects and does not require an empirical tortuosity parameter. The current and concentration distributions in localized micropores were also investigated and it was shown that the microporous area is fully accessible to charge and mass transfer processes. From these analyses it was concluded that the micropores primarily affect the kinetics of the interfacial processes by contributing to the interfacial area, while the macropores impose ohmic and mass transport limitations through the volume of the porous electrode.

Lanzi, Oscar↗

An overview of controls research on the NASA Langley Research Center grid

The NASA Langley Research Center has assembled a flexible grid on which control systems research can be accomplished on a two-dimensional structure that has many physically distributed sensors and actuators. The grid is a rectangular planar structure that is suspended by two cables attached to one edge so that out of plane vibrations are normal to gravity. There are six torque wheel actuators mounted to it so that torque is produced in the grid plane. Also, there are six rate gyros mounted to sense angular motion in the grid plane and eight accelerometers that measure linear acceleration normal to the grid plane. All components can be relocated to meet specific control system test requirements. Digital, analog, and hybrid control systems capability is provided in the apparatus. To date, research on this grid has been conducted in the areas of system and parameter identification, model estimation, distributed modal control, hierarchical adaptive control, and advanced redundancy management algorithms. The presentation overviews each technique and presents the most significant results generated for each area.

Montgomery, Raymond C.↗

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Investigating the large-scale environment of wide-angle tailed radio galaxies in the local Universe

We present a statistical analysis of the large-scale (up to two Mpc) environment of a homogeneous and complete sample, both in radio and optical selection, of wide-angle tailed radio galaxies (WATs) in the local Universe (i.e., with redshifts z ≲ 0.15). The analysis is carried out using the parameters obtained from cosmological neighbors within two Mpc of the target source. Results regarding the large-scale environments of the WATs are then compared with those of Fanaroff–Riley type I (FR Is) and type II (FR IIs) radio galaxies listed in two other homogeneous and complete catalogs and selected with the same criterion adopted for the WATs catalog. We obtain indication that at low redshift, WATs inhabit environments with a larger number of galaxies than those of FR Is and FR IIs. In the explored redshift range, the physical size of the galaxy group or cluster in which the WATs reside appears to be almost constant with respect to FR Is and FR IIs, being around one Mpc. From the distribution of the concentration parameter, defined as the ratio between the number of cosmological neighbors lying within 500 kpc and within one Mpc, we conclude that WATs tend to inhabit the central region of the group or cluster in which they reside, in agreement with the general paradigm that WATs are the brightest cluster or group galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Variability and associated uncertainty in image analysis for soiling characterization in solar energy systems

The accumulation of soiling on photovoltaic modules and on the mirrors of concentrating solar power systems causes non-negligible energy losses with economic consequences. These challenges can be mitigated, or even prevented, through appropriate actions if the magnitude of soiling is known. Particle counting analysis is a common procedure to characterize soiling, as it can be easily performed on micrographs of glass coupons or solar devices that have been exposed to the environment. Particle counting does not, however, yield invariant results across institutions. The particle size distribution analysis is affected by the operator of the image analysis software and the methodology utilized. The results of a round-robin study are presented in this work to explore and elucidate the uncertainty related to particle counting and its effect on the characterization of the soiling of glass surfaces used in solar energy conversion systems. An international group of soiling experts analyzed the same 8 micrographs using the same open-source ImageJ software package. The variation in the particle analyses results were investigated to identify specimen characteristics with the lowest coefficient of variation (CV) and the least uncertainty among the various operators. The mean particle diameter showed the lowest CV among the investigated characteristics, whereas the number of particles exhibited the largest CV. Additional parameters, such as the fractional area coverage by particles and parameters related to the distribution's shape yielded intermediate CV values. These results can provide insights on the magnitude inter-lab variability and uncertainty for optical and microscope-based soiling monitoring and characterization.

14 SOLAR ENERGY↗

Analysis of plasma measurements for the Geotail mission

The CPI plasma measurements from the Geotail spacecraft are currently used by a number of scientists in support of varied research projects. Measurements from the CPI hot plasma analyzer have been processed to compute one-minute averages of plasma densities, temperatures, and velocities for the period 1 Oct. 1992 - 31 July 1993. These parameters are used in a preliminary survey of the magnetotail for distances from earth of 10 to 210 earth radii. The cold plasmas within the magnetotail drift towards the midplane and are thought to be a principal source for the hot plasma sheet. A remarkable result from Geotail is the observation of cold ion beams coexisting as distinct components in the presence of hot plasma-sheet plasmas. The development and evolution of plasmoids is a topic of considerable interest in studies of the magnetotail and magnetospheric substorms. A number of possible plasmoids have been identified in the Geotail data set based upon reversals of the Z component of the magnetic field as observed with the MGF instrumentation. An intensive study of the electron and ion velocity distributions and the plasma parameters for several of these events reveals unexpected features. Collaborative work with several Geotail instrument teams and other researchers is also ongoing.

Frank, Louis A.↗

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electron-beam energy reconstruction for neutrino oscillation measurements

Neutrinos exist in one of three types or ‘flavours’—electron, muon and tau neutrinos—and oscillate from one flavour to another when propagating through space. This phenomena is one of the few that cannot be described using the standard model of particle physics (reviewed in ref. 1), and so its experimental study can provide new insight into the nature of our Universe (reviewed in ref. 2). Neutrinos oscillate as a function of their propagation distance (L) divided by their energy (E). Therefore, experiments extract oscillation parameters by measuring their energy distribution at different locations. As accelerator-based oscillation experiments cannot directly measure E, the interpretation of these experiments relies heavily on phenomenological models of neutrino–nucleus interactions to infer E. Here we exploit the similarity of electron–nucleus and neutrino–nucleus interactions, and use electron scattering data with known beam energies to test energy reconstruction methods and interaction models. We find that even in simple interactions where no pions are detected, only a small fraction of events reconstruct to the correct incident energy. More importantly, widely used interaction models reproduce the reconstructed energy distribution only qualitatively and the quality of the reproduction varies strongly with beam energy. This shows both the need and the pathway to improve current models to meet the requirements of next-generation, high-precision experiments such as Hyper-Kamiokande (Japan) and DUNE (USA).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The effects of correlation on goodness of fit

Autocorrelated normal random variates were generated via computer and the effects of various levels of correlation on goodness of fit problems were studied. The results are useful in determining the distribution or estimating the parameters of populations that correlate observations such as wind speeds and temperature. The model used to generate the autocorrelated data is an autoregressive process of order 1. The Kolmogorov-Smirnov and chi-square statistics are used in the analysis. It was observed in the simulation that high positive correlations tend to shift the sample mean away from the population mean and negative correlations tend to shift the sample mean towards the population mean. In many cases, it was observed that positive and negative correlations tend to decrease the standard deviation. However, since this did not occur in all cases, no definite conclusion can be made regarding the standard deviation. Since the autoregressive process is a linear transformation, it is not surprising that normality was preserved. However, a possible extension of this problem could be to generate non-normal data and observe how the distribution is affected by correlation. Another extension might utilize another model such as autoregressive of order k or a moving average process of order k.

Kitchens, L. J.↗