SEARCH · Engineering Papers

Results for “Conditional distributions”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

NASA NTRS1985

The cleft ion fountain - A two-dimensional kinetic model

The transport of ionospheric ions from a source in the polar cleft ionosphere through the polar magnetosphere is investigated using a two-dimensional, kinetic, trajectory-based code. The transport model includes the effects of gravitation, longitudinal magnetic gradient force, convection electric fields, and parallel electric fields. Individual ion trajectories as well as distribution functions and resulting bulk parameters of density, parallel average energy, and parallel flux for a presumed cleft ionosphere source distribution are presented for various conditions to illustrate parametrically the dependences on source energies, convection electric field strengths, ion masses, and parallel electric field strengths. The essential features of the model are consistent with the concept of a cleft-based ion fountain supplying ionospheric ions to the polar magnetosphere, and the resulting plasma distributions and parameters are in general agreement with recent low-energy ion measurements from the DE 1 satellite.

Horwitz, J. L.↗

DOE OSTI2020

Online Optimization for Networked Distributed Energy Resources With Time-Coupling Constraints

This paper proposes a Lyapunov optimization-based online distributed (LOOD) algorithmic framework for active distribution networks (ADNs) with numerous photovoltaic inverters and inverter air conditionings (IACs). In the proposed scheme, ADNs can track an active power setpoint reference at the substation in response to transmission-level requests while concurrently minimizing the social utility loss and ensuring the security of voltages. Conventional distributed optimization methods are rarely feasible to track the optimal solutions in fast variable environments using a fine-grained sampling interval where the underlying optimization problem evolves with the iterations of the algorithms. In contrast, based on the framework of online convex optimization (OCO), the developed approach uses a distributed algebraic update to compute the next round decisions relying on the current feedback of measurements. Notably, the time-coupling constraints of IACs are decoupled for online implementation with Lyapunov optimization technique. An incentive scheme is tailored to coordinate the customer-owned assets in lieu of the direct control from network operators. Optimality and convergency are characterized analytically. Finally, we corroborate the proposed method on a modified version of 33-node test feeder. Benchmark tests show that the proposed method is computationally and economically efficient, and outperforming existing algorithms.

active distribution networks↗

DOE OSTI2024

3002 Humidified Tandem Differential Mobility Analyzer (HTDMA) Instrument Handbook

The Brechtel Manufacturing Inc. (BMI) Humidified Tandem Differential Mobility Analyzer (HT-DMA Model 3002) (Brechtel and Kreidenweis 2000a,b, Henning et al. 2005, Xerxes et al. 2014) measures how aerosol particles of different initial dry sizes grow or shrink when exposed to changing relative humidity (RH) conditions. It uses two different mobility analyzers (DMA) and a humidification system to make the measurements. One DMA selects a narrow size range of dry aerosol particles, which are exposed to varying RH conditions in the humidification system. The second (humidified) DMA scans the particle size distribution output from the humidification system. Scanning a wide range of particle sizes enables the second DMA to measure changes in size or growth factor (growth factor = humidified size/dry size), due to water uptake by the particles. A Condensation Particle Counter (CPC) downstream of the second DMA counts particles as a function of selected size in order to obtain the number size distribution of particles exposed to different RH conditions.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

Analysis of hypersonic arcjet flow fields and surface heating of blunt bodies

A Gauss-Seidel implicit aerothermodynamic Navier-Stokes computational CFD code with thermochemical surface conditions (GIANTS) and with appropriate boundary conditions was developed for the analysis of hypersonic arcjet flows over large angle blunt bodies. Two-dimensional axisymmetric multicomponent full Navier-Stokes equations were solved, together with boundary conditions, to account for the surface slip, surface catalysis, and energy and mass conservation at the body surface. It is shown that the predictions for isothermal, noncatalytic, and no-slip surface conditions which are made using the GIANTS code with simple boundary conditions agree with those of the LAURA code, and the predicted normalized heating distributions over the forebody surface with boundary conditions for various cone angles are consistent with the trends of data from arcjet experiments.

Chen, Y.-K.↗

DOE OSTI2015

Molecular structure, morphology and growth mechanisms and rates of 5-hydroxymethyl furfural (HMF) derived humins

We apply ATR-FTIR spectroscopy, Scanning Electron Microscopy (SEM) and Dynamic Light Scattering (DLS) experiments to investigate the molecular structure, morphology and growth mechanism of 5-hydroxymethyl furfural (HMF) derived humins as a function of HMF conversion. Our FTIR data support a reaction pathway in which humins form either through a ring opening mechanism and/or through substitution at the α or β position via nucleophilic attack. The addition of DMSO as a co-solvent leads to significant changes in the FTIR spectra of humins. We find that the nucleophilic attack pathway is suppressed in the presence of DMSO co-solvent and rationalizes the very small humin particles (~100 nm) observed in SEM images contrary to the large particles (with multimodal size distribution and largest particles of up to 3–4 μm) observed in neat water. DLS experiments under several reaction conditions further confirm the particle size distribution observed via SEM. A plausible reaction network for humin formation, which rationalizes qualitatively our experimental results as well as those reported in the literature, is also postulated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

MatPhase: Material phase prediction for Li-ion Battery Reconstruction using Hierarchical Curriculum Learning

Li-ion Batteries (LIB), one of the most efficient energy storage devices, are used extensively in many industrial applications. These batteries consist of electrodes that are put together with heterogeneous material compositions. Imaging data of these battery electrodes obtained from X-ray tomography can explain the distribution of material constituents and allow reconstructions to study electron transport pathways. Such reconstructions of material constituents help quantify various associated properties of electrodes (e.g., volume-specific surface area, porosity) which determine the performance of batteries. These images often suffer from low image contrast between multiple material constituents, hence making it difficult for humans to distinguish and characterize these constituents through visual inspection. A minor error in detecting distributions of the material constituents can lead to magnified errors in the calculated parameters of material properties (e.g., porosity). We present MatPhase, a novel hierarchical curriculum learning technique to address the complex task of estimating material constituent distribution in battery electrodes. MatPhase comprises three modules: (i) an uncertainty-aware global model trained to yield inferences conditioned upon global knowledge of material distribution, (ii) a local model to capture relatively more fine-grained (local) distributional signals, (iii) an aggregator model to appropriately fuse the local and global effects towards obtaining the final distribution. On average, MatPhase improves prediction up to 8.5% relative to other sophisticated modeling pipelines and state-of-the-arts (SOTA) object detection models employed in the performance comparison.

Tabassum, Anika↗

NASA NTRS1996

The Global Perspective on the Evolution of Solids in a Protoplanetary Disk

It is currently thought that planets around solar-type stars form by the accumulation of solid matter entrained in a gaseous, turbulent protoplanetary disk. We have developed a model designed to simulate the part of this process that starts from small particles suspended in the gaseous disk at the end of the formation stage, and ends up with most of the solid material aggregated into 1-10-km planetesimals. The major novelty of our approach is its emphasis on the global, comprehensive treatment of the problem, as our model simultaneously keeps track of the evolution of gas and solid particles due to gas-solid coupling, coagulation, sedimentation, and evaporation/condensation. The result of our calculations is the radial distribution of solid material circumnavigating a star in the form of a planetesimal swarm. Such a distribution should well approximate the radial apportionment of condensed components of the planets spread over the radial extent of the mature planetary system. Therefore we view our calculations as an attempt to predict the large-scale architecture of planetary systems and to assess their potential diversity. In particular, we have found that some initial conditions lead to all solids being lost to the star, but we can also identify initial conditions leading to a radial distribution of solid material quite reminiscent of what is found in our solar system.

Stepinski, T. F.↗

DOE OSTI2014

Gas-kinetic unified algorithm for hypersonic flows covering various flow regimes solving Boltzmann model equation in nonequilibrium effect

Based on the Gas-Kinetic Unified Algorithm (GKUA) directly solving the Boltzmann model equation, the effect of rotational non-equilibrium is investigated recurring to the kinetic Rykov model with relaxation property of rotational degrees of freedom. The spin movement of diatomic molecule is described by moment of inertia, and the conservation of total angle momentum is taken as a new Boltzmann collision invariant. The molecular velocity distribution function is integrated by the weight factor on the internal energy, and the closed system of two kinetic controlling equations is obtained with inelastic and elastic collisions. The optimization selection technique of discrete velocity ordinate points and numerical quadrature rules for macroscopic flow variables with dynamic updating evolvement are developed to simulate hypersonic flows, and the gas-kinetic numerical scheme is constructed to capture the time evolution of the discretized velocity distribution functions. The gas-kinetic boundary conditions in thermodynamic non-equilibrium and numerical procedures are studied and implemented by directly acting on the velocity distribution function, and then the unified algorithm of Boltzmann model equation involving non-equilibrium effect is presented for the whole range of flow regimes. The hypersonic flows involving non-equilibrium effect are numerically simulated including the inner flows of shock wave structures in nitrogen with different Mach numbers of 1.5-Ma-25, the planar ramp flow with the whole range of Knudsen numbers of 0.0009-Kn-10 and the three-dimensional re-entering flows around tine double-cone body.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1990

Pitch angle distributions of energetic ions in the lobes of the distant geomagnetic tail

A model is presented which accounts for observed pitch angle distributions of energetic ions in the lobes of the deep geomagnetic tail. Energetic (greater than 35 keV) ion data from the ISEE-3 spacecraft obtained during 1982-1983, when the spacecraft made a series of transversals of the distant geomagnetic tail (X sub GSE greater than -238 R sub E) is considered. Analysis indicates that under usual conditions the lobe pitch angle distribution is highly peaked at pitch angles closely perpendicular to the magnetic field direction while also showing a small downtail anisotropy. It is noted, however, that during periods of solar particle enhancements the lobe distributions become much more clearly isotropic. The proposed model successfully reproduces these features of the data and also accounts for the more isotropic ion population observed in the lobe during solar particle events, when the terrestrial component of the magnetosheath source may be considered negligeable in comparison to the enhanced solar component.

Owen, C. J.↗

DOE OSTI2024

Machine learning assisted bayesian inference of mix and hot-spot conditions in NIF implosions

Experiments on the National Ignition Facility (NIF) have provided clear evidence of ablator material mixing into the Hot-Spot, leading to degraded performance. However, inferring the amount of mix and Hot-Spot conditions from typical experimental observations (e.g. x-ray spectra and images) is highly challenging. Here, we have developed an analysis method that utilizes machine learning assisted Bayesian inference to find the probability distributions of the Hot-Spot and mix conditions. This approach uses a neural network, trained on an idealized 2-dimensional representation of the Hot-Spot and mix distribution, and Bayesian inference to find the statistical distributions of Hot-Spot conditions that provide a match with observations. We have tested this method with synthetic data from simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Kinetics Measurements in Resistive Electrolytes Using Ring-Disk Electrode: Ring as Current “Shield” Enables Uniform Disk Current Distribution

Rotating disk electrodes are commonly used for electrochemical kinetics measurements. A major disadvantage of these types of electrodes is their nonuniform secondary current distribution, especially when performing electroanalytical measurements in resistive electrolytes. Such nonuniform current distribution can render the values of kinetics constants extracted from the disk electrode to be highly inaccurate. Furthermore, one emerging class of electrolytes that suffer from low ionic conductivities is deep eutectic solvents (DES). DES are a promising class of electrolytes for various emerging applications; however, due to their resistive nature, the secondary current distribution when using them is typically highly nonuniform. For example, the Wagner number when measuring Cu²⁺/Cu⁺ kinetics in choline chloride–ethylene glycol DES (1:4 molar ratio of ChCl:EG) is very low (<0.1), indicating highly nonuniform current distribution over the disk electrode. We show here that the Cu²⁺/Cu⁺ exchange current density measured using disk electrodes is very inaccurate due to the aforementioned nonuniform current distribution. To obtain uniform disk current distribution, we employ here a coplanar concentric rotating ring-disk electrode (RRDE), where the ring serves the function of a current “shield.” Specifically, we show using modeling that the ring minimizes the current distribution nonuniformity at the disk by effectively shielding the disk against current spikes near the disk edge. This enables improved precision in electrode kinetics measurements for the Cu²⁺/Cu⁺ couple in resistive DES. To enable broad applicability of this technique, an analytical expression based on the Wagner number is integrated into an iterative algorithm to help users identify ring conditions to achieve uniform current distribution and thus improved electroanalytics at the disk.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗

DOE OSTI2016

Modeling and control of fuel distribution in a dual-fuel internal combustion engine leveraging late intake valve closings

In internal combustion engines, cycle-to-cycle and cylinder-to-cylinder variations of the combustion process have been shown to negatively impact the fuel efficiency of the engine and lead to higher exhaust emissions. The combustion variations are generally tied to differences in the composition and condition of the trapped mass throughout each cycle and across individual cylinders. Thus, advanced engines featuring exhaust gas recirculation, flexible valve actuation systems, advanced fueling strategies, and turbocharging systems are prone to exhibit higher variations in the combustion process. In this study, the cylinder-to-cylinder variations of the combustion process in a dual-fuel internal combustion engine leveraging late intake valve closing are investigated and a model to predict and address one of the root causes for these variations across cylinders is developed. The study is conducted on an inline six-cylinder heavy-duty dual-fuel engine equipped with exhaust gas recirculation, a variable geometry turbocharger, and a fully flexible variable intake valve actuation system. The engine is operated with late intake valve closure timings in a dual-fuel combustion mode in which a high reactivity fuel is directly injected into the cylinders and a low reactivity fuel is port injected into the cylinders. The cylinder-to-cylinder variations observed in the study have been associated with the maldistribution of the port-injected fuel, which is exacerbated at late intake valve timings. The resulting difference in indicated mean effective pressure between the cylinders ranges from 9% at an intake valve closing of 570° after top dead center to 38% at an intake valve closing of 620° after top dead center and indicates an increasingly uneven fuel distribution. The study leverages both experimental and simulation studies to investigate the distribution of the port-injected fuel and its impact on cylinder-to-cylinder variation. The effects of intake valve closing as well as the impact of intake runner length on fuel distribution were quantitatively analyzed, and a model was developed that can be used to accurately predict the fuel distribution of the port-injected fuel at different operating conditions with an average estimation error of 1.5% in cylinder-specific fuel flow. A model-based control strategy is implemented to adjust the fueling at each port and shown to significantly reduce the cylinder-to-cylinder variations in fuel distribution.

Kassa, Mateos↗

DOE OSTI2016

Minimizing Uncertainties Impact in Decision Making with an Applicability Study for Economic Power Dispatch

This report summaries the work performed under the LDRD project on the preliminary study on knowledge automation, where specific focus has been made on the investigation of the impact of uncertainties of human decision making onto the optimization of the process operation. At first the statistics on signals from the Brain-Computing Interface (BCI) is analyzed so as to obtain the uncertainties characterization of human operators during the decision making phase using the electroencephalogram (EEG) signals. This is then followed by the discussions of an architecture that reveals the equivalence between optimization and closed loop feedback control design, where it has been shown that all the optimization problems can be transferred into the control design problem for closed loop systems. This has led to a “closed loop” framework, where the structure of the decision making is shown to be subjected to both process disturbances and controller’s uncertainties. The latter can well represent the uncertainties or randomness occurred during human decision making phase. As a result, a stochastic optimization problem has been formulated and a novel solution has been proposed using probability density function (PDF) shaping for both the cost function and the constraints using stochastic distribution control concept. A sufficient condition has been derived that guarantees the convergence of the optimal solution and discussions have been made for both the total probabilistic solution and chanced constrained optimization which have been well-studied in optimal power flows (OPF) area. A simple case study has been carried out for the economic dispatch of powers for a grid system when there are distributed energy resources (DERs) in the system, and encouraging results have been obtained showing that a significant savings on the generation cost can be expected.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

DOE OSTI2022

Measurement of Triple-Differential Inclusive Muon Neutrino Charged-Current Cross Sections on Argon with the MicroBooNE Detector

We report the first measurement of the differential cross section $d^2σ(E_ν)/dcos(θ_µ)dP_µ$for inclusive muon neutrino charged-current scattering on argon. This measurement utilizes data from6.5×10 20 protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located on-axis in the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8 GeV. In this note, the mapping from reconstructed kinematics to truth quantities, particularly from reconstructed to true neutrino energy, is validated by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after conditional constraint from the muon momentum distribution in data. The success of this validation gives confidence that the missing energy in the MicroBooNE detector is well-modeled within its uncertainties in simulation, enabling the unfolding to an energy-dependent triple-differential measurement over muon momentum, muon scattering angle, and neutrino energy. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future LArTPC experiments measuring neutrino oscillations. Comparisons with model predictions show the best agreement with Neut 5.4.0.1 at low energy, and with NuWro 19.02.01 at higher energies, particularly at forward muon scattering angles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1997

Achieving Zero Current for Polar Wind Outflow on Open Flux Tubes Subjected to Large Photoelectron Fluxes

In this study we investigate how the condition of zero current on open flux tubes with polar wind outflow, subjected to large photoelectron fluxes, can be achieved. We employ a steady state collisionless semikinetic model to determine the density profiles of O(+), H(+), thermal electrons and photoelectrons coming from the ionosphere along with H(+), ions and electrons coming from the magnetosphere. The model solution attains a potential distribution which both satisfies the condition of charge neutrality and zero current. For the range of parameters considered in this study we find that a 45-60 volt discontinuous potential drop may develop to reflect most of the photoelectrons back toward the ionosphere. This develops because the downward flux of electrons from the magnetosphere to the ionosphere on typical open flux tubes (e.g. the polar rain) appears to be insufficient to balance the photoelectron flux from the ionosphere.

Wilson, G. R.↗

DOE OSTI2021

Development and Evaluation of Chemistry‐Aerosol‐Climate Model CAM5‐Chem‐MAM7‐MOSAIC: Global Atmospheric Distribution and Radiative Effects of Nitrate Aerosol

Abstract An advanced aerosol treatment, with a focus on semivolatile nitrate formation, is introduced into the Community Atmosphere Model version 5 with interactive chemistry (CAM5‐chem) by coupling the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) with the 7‐mode Modal Aerosol Module (MAM7). An important feature of MOSAIC is dynamic partitioning of all condensable gases to the different fine and coarse mode aerosols, as governed by mode‐resolved thermodynamics and heterogeneous chemical reactions. Applied in the free‐running mode from 1995 to 2005 with prescribed historical climatological conditions, the model simulates global distributions of sulfate, nitrate, and ammonium in good agreement with observations and previous studies. Inclusion of nitrate resulted in ∼10% higher global average accumulation mode number concentrations, indicating enhanced growth of Aitken mode aerosols from nitrate formation. While the simulated accumulation mode nitrate burdens are high over the anthropogenic source regions, the sea‐salt and dust modes respectively constitute about 74% and 17% of the annual global average nitrate burden. Regional clear‐sky shortwave radiative cooling of up to −5 W m −2 due to nitrate is seen, with a much smaller global average cooling of −0.05 W m −2 . Significant enhancements in regional cloud condensation nuclei (at 0.1% supersaturation) and cloud droplet number concentrations are also attributed to nitrate, causing an additional global average shortwave cooling of −0.8 W m −2 . Taking into consideration of changes in both longwave and shortwave radiation under all‐sky conditions, the net change in the top of the atmosphere radiative fluxes induced by including nitrate aerosol is −0.7 W m −2 .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Universal method for the optimization of HDC coating uniformity on non-planar, non-stationary substrates for inertial confinement fusion targets

The thickness uniformity of chemical vapor deposited (CVD) diamond coatings on non-planar, non-stationary substrates depends on both the intrinsic instantaneous coating thickness distribution (ICTD) of the coating conditions used and, if applicable, on the frequency of substrate reorientation. While important for many CVD diamond applications, the relative impact of the ICTD and substrate reorientation on the coating thickness uniformity has not been studied. In this work, we systematically investigate the effect of these factors for microwave-plasma chemical vapor deposition (MPCVD) of diamond (referred to as high density carbon (HDC) in the inertial confinement fusion (ICF) community) coatings on spherical, rolling substrates. This coating technique is used to fabricate capsules for ICF experiments, which require extreme coating uniformity with <0.3 % thickness variation (so-called Mode 1 or M1) to ensure symmetric compression of imploding targets. To extract the otherwise unobservable reorientation timescale (Δt), Monte Carlo simulations were performed using experimental ICTD data as input. This combined approach confirms scaling relationships between the substrate reorientation timescale as well as coating thickness and coating uniformity, as expected from a 3D random walk. Simulations confirm that M1 is Rayleigh-distributed and scales as (Δt) 1/2 , consistent with the randomization of two angles that determine orientation of a sphere. We also demonstrate that, under the conditions studied, Δt is the dominant factor in determining thickness uniformity while the intrinsic ICTD has minimal impact. Finally, experiments show that Δt can be affected by total batch size under constant agitation conditions due to space constraints that limit the capsule reorientation kinetics. In conclusion, this study highlights the utility of a combined experiment-simulation approach as a general methodology for understanding and improving coating uniformity on non-planar, non-stationary substrates.

Capsule↗