Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “generator parameter estimation”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Use of DAGMan in CRAB3 to improve the splitting of CMS user jobs

CRAB3 is a workload management tool used by CMS physicists to analyze data acquired by the Compact Muon Solenoid (CMS) detector at the CERN Large Hadron Collider (LHC). Research in high energy physics often requires the analysis of large collections of files, referred to as datasets. The task is divided into jobs that are distributed among a large collection of worker nodes throughout the Worldwide LHC Computing Grid (WLCG). Splitting a large analysis task into optimally sized jobs is critical to efficient use of distributed computing resources. Jobs that are too big will have excessive runtimes and will not distribute the work across all of the available nodes. However, splitting the project into a large number of very small jobs is also inefficient, as each job creates additional overhead which increases load on infrastructure resources. Currently this splitting is done manually, using parameters provided by the user. However the resources needed for each job are difficult to predict because of frequent variations in the performance of the user code and the content of the input dataset. As a result, dividing a task into jobs by hand is difficult and often suboptimal. In this work we present a new feature called “automatic splitting” which removes the need for users to manually specify job splitting parameters. We discuss how HTCondor DAGMan can be used to build dynamic Directed Acyclic Graphs (DAGs) to optimize the performance of large CMS analysis jobs on the Grid. We use DAGMan to dynamically generate interconnected DAGs that estimate the processing time the user code will require to analyze each event. This is used to calculate an estimate of the total processing time per job, and a set of analysis jobs are run using this estimate as a specified time limit. Some jobs may not finish within the alloted time; they are terminated at the time limit, and the unfinished data is regrouped into smaller jobs and resubmitted.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Determining Series Resistance for Equivalent Circuit Models of a PV Module

Literature describes various methods for determining a series resistance for a photovoltaic device from measured IV curves. We investigate use of these techniques to estimate the series resistance parameter for a single diode equivalent circuit model. With simulated IV curves we demonstrate that the series resistance values obtained by these techniques differ systematically from the known series resistance parameter values used to generate the curves, indicating that these methods are not suitable for determining the series resistance parameter for the single diode model equation. Furthermore, we present an alternative method to determine the series resistance parameter jointly with the other parameters for the single diode model equation, and demonstrate the accuracy and reliability of this technique in the presence of measurement errors.

42 ENGINEERING↗
DOE OSTI2025

Optimal Power Management for Large-Scale Battery Energy Storage Systems via Bayesian Inference

Large-scale battery energy storage systems (BESS) have found ever-increasing use across industry and society to accelerate clean energy transition and improve energy supply reliability and resilience. However, their optimal power management poses significant challenges: the underlying high-dimensional nonlinear nonconvex optimization lacks computational tractability in real-world implementation, and the uncertainty of the exogenous power demand makes exact optimization difficult. This paper presents a new solution framework to address these bottlenecks. The solution pivots on introducing power-sharing ratios to specify each cell’s power quota from the output power demand. To find the optimal power-sharing ratios, we formulate a nonlinear model predictive control (NMPC) problem to achieve power-loss-minimizing BESS operation while complying with safety, cell balancing, and power supply-demand constraints. We then propose a parameterized control policy for the power-sharing ratios, which utilizes only three parameters, to reduce the computational demand in solving the NMPC problem. This policy parameterization allows us to translate the NMPC problem into a Bayesian inference problem for the sake of 1) computational tractability, and 2) overcoming the nonconvexity of the optimization problem. We leverage the ensemble Kalman inversion technique to solve the parameter estimation problem. Concurrently, a low-level control loop is developed to seamlessly integrate our proposed approach with the BESS to ensure practical implementation. This low-level controller receives the optimal power-sharing ratios, generates output power references for the cells, and maintains a balance between power supply and demand despite uncertainty in output power. We conduct extensive simulations and experiments on a 20-cell prototype to validate the proposed approach.

Battery energy storage systems (BESSs)↗
DOE OSTI2016

A Dirichlet process mixture model for automatic {sup 18}F-FDG PET image segmentation: Validation study on phantoms and on lung and esophageal lesions

Purpose: The aim of this study was to implement a Dirichlet process mixture (DPM) model for automatic tumor edge identification on {sup 18}F-fluorodeoxyglucose positron emission tomography ({sup 18}F-FDG PET) images by optimizing the parameters on which the algorithm depends, to validate it experimentally, and to test its robustness. Methods: The DPM model belongs to the class of the Bayesian nonparametric models and uses the Dirichlet process prior for flexible nonparametric mixture modeling, without any preliminary choice of the number of mixture components. The DPM algorithm implemented in the statistical software package R was used in this work. The contouring accuracy was evaluated on several image data sets: on an IEC phantom (spherical inserts with diameter in the range 10–37 mm) acquired by a Philips Gemini Big Bore PET-CT scanner, using 9 different target-to-background ratios (TBRs) from 2.5 to 70; on a digital phantom simulating spherical/uniform lesions and tumors, irregular in shape and activity; and on 20 clinical cases (10 lung and 10 esophageal cancer patients). The influence of the DPM parameters on contour generation was studied in two steps. In the first one, only the IEC spheres having diameters of 22 and 37 mm and a sphere of the digital phantom (41.6 mm diameter) were studied by varying the main parameters until the diameter of the spheres was obtained within 0.2% of the true value. In the second step, the results obtained for this training set were applied to the entire data set to determine DPM based volumes of all available lesions. These volumes were compared to those obtained by applying already known algorithms (Gaussian mixture model and gradient-based) and to true values, when available. Results: Only one parameter was found able to significantly influence segmentation accuracy (ANOVA test). This parameter was linearly connected to the uptake variance of the tested region of interest (ROI). In the first step of the study, a calibration curve was determined to automatically generate the optimal parameter from the variance of the ROI. This “calibration curve” was then applied to contour the whole data set. The accuracy (mean discrepancy between DPM model-based contours and reference contours) of volume estimation was below (1 ± 7)% on the whole data set (1 SD). The overlap between true and automatically segmented contours, measured by the Dice similarity coefficient, was 0.93 with a SD of 0.03. Conclusions: The proposed DPM model was able to accurately reproduce known volumes of FDG concentration, with high overlap between segmented and true volumes. For all the analyzed inserts of the IEC phantom, the algorithm proved to be robust to variations in radius and in TBR. The main advantage of this algorithm was that no setting of DPM parameters was required in advance, since the proper setting of the only parameter that could significantly influence the segmentation results was automatically related to the uptake variance of the chosen ROI. Furthermore, the algorithm did not need any preliminary choice of the optimum number of classes to describe the ROIs within PET images and no assumption about the shape of the lesion and the uptake heterogeneity of the tracer was required.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Electrons for Neutrinos

The ability of current and next generation accelerator based neutrino oscillation measurementsto reach their desired sensitivity requires a high level of understanding of the neutrino-nucleus interactions. This includes precise estimation of the relevant cross sections for the incident neutrino energy from the measured final state particles. Incomplete understanding of these interactions can skew the reconstructed neutrino spectrum and thereby bias the extraction of fundamental oscillation parameters and searches for new physics. In this talk, I will present new results of wide phase-space electron scattering data, collected with the CLAS spectrometer at the Thomas Jefferson National Accelerator Facility (JLab), where the reconstruction of the incoming lepton energy from the measured final state is being tested. Disagreements with current event generators, used in the analysis of neutrino oscillation measurements, are observed which indicate underestimation of nuclear effects.

43 PARTICLE ACCELERATORS↗
NASA NTRS1990

Estimation of synoptic fields of middle atmosphere parameters from Nimbus-7 LIMS profile data

The most widely used version of the Nimbus-7 middle atmosphere dataset is the set of high quality, daily, and zonal Fourier coefficients that resolve information out to six wavenumbers at 12 UTC. A Kalman filter algorithm was applied to the original profile data in order to generate those fields for the data archive or LAMAT product. The characteristics and implementation of the algorithm are described in some detail, along with examples of the output for each of the LIMS parameters.

Remsberg, Ellis E.↗
NASA NTRS1992

Light element nucleosynthesis - A false clue?

It is proposed that the dynamically estimated value for the cosmological density parameter, Omega(dy) = 0.15 x 10 exp +/- 0.20, reflects the baryon density at decoupling, resulting in lower initial, primordial values of D and the He-3 than are observed. An early generation of massive stars, forming somewhat after decoupling, collapses to black holes with masses of about 10 exp 6.5 solar masses. If they later accrete gas and emit a quasarlike (X, gamma)-ray spectrum, then (gamma, (He-4)) photodisintegration reactions will increase D and He-3 to the observed range, leaving a high-energy background radiation field similar to that observed. The massive black holes become the dynamically observed dark matter galactic halos. This scenario obviates the need for nonbaryonic dark matter and provides a specific form for the requisite baryonic dark matter; it thus reduces the number of density parameters: Omega(0) - Omega(dy) = Omega(BBN) is approximately equal to 0.15.

Gnedin, N. IU.↗
NASA NTRS1990

Statistical circuit design for yield improvement in CMOS circuits

This paper addresses the statistical design of CMOS integrated circuits for improved parametric yield. The work uses the Monte Carlo technique of circuit simulation to obtain an unbiased estimation of the yield. A simple graphical analysis tool, the yield factor histogram, is presented. The yield factor histograms are generated by a new computer program called SPICENTER. Using the yield factor histograms, the most sensitive circuit parameters are noted, and their nominal values are changed to improve the yield. Two basic CMOS example circuits, one analog and one digital, are chosen and their designs are 'centered' to illustrate the use of the yield factor histograms for statistical circuit design.

Kamath, H. J.↗
NASA NTRS2010

Improved Determination of Surface and Atmospheric Temperatures Using Only Shortwave AIRS Channels: The AIRS Version 6 Retrieval Algorithm

AIRS was launched on EOS Aqua on May 4, 2002 together with ASMU-A and HSB to form a next generation polar orbiting infrared and microwave atmosphere sounding system (Pagano et al 2003). The theoretical approach used to analyze AIRS/AMSU/HSB data in the presence of clouds in the AIRS Science Team Version 3 at-launch algorithm, and that used in the Version 4 post-launch algorithm, have been published previously. Significant theoretical and practical improvements have been made in the analysis of AIRS/AMSU data since the Version 4 algorithm. Most of these have already been incorporated in the AIRS Science Team Version 5 algorithm (Susskind et al 2010), now being used operationally at the Goddard DISC. The AIRS Version 5 retrieval algorithm contains three significant improvements over Version 4. Improved physics in Version 5 allowed for use of AIRS clear column radiances (R(sub i)) in the entire 4.3 micron CO2 absorption band in the retrieval of temperature profiles T(p) during both day and night. Tropospheric sounding 15 micron CO2 observations were used primarily in the generation of clear column radiances (R(sub i)) for all channels. This new approach allowed for the generation of accurate Quality Controlled values of R(sub i) and T(p) under more stressing cloud conditions. Secondly, Version 5 contained a new methodology to provide accurate case-by-case error estimates for retrieved geophysical parameters and for channel-by-channel clear column radiances. Thresholds of these error estimates are used in a new approach for Quality Control. Finally, Version 5 contained for the first time an approach to provide AIRS soundings in partially cloudy conditions that does not require use of any microwave data. This new AIRS Only sounding methodology was developed as a backup to AIRS Version 5 should the AMSU-A instrument fail. Susskind et al 2010 shows that Version 5 AIRS Only sounding are only slightly degraded from the AIRS/AMSU soundings, even at large fractional cloud cover.

Susskind, Joel↗
DOE OSTI2022

Why are cell populations maintained via multiple compartments?

We consider the maintenance of ‘product’ cell populations from ‘progenitor’ cells via a sequence of one or more cell types, or compartments, where each cell’s fate is chosen stochastically. If there is only one compartment then large amplification, that is, a large ratio of product cells to progenitors comes with disadvantages. The product cell population is dominated by large families (cells descended from the same progenitor) and many generations separate, on average, product cells from progenitors. These disadvantages are avoided using suitably constructed sequences of compartments: the amplification factor of a sequence is the product of the amplification factors of each compartment, while the average number of generations is a sum over contributions from each compartment. Passing through multiple compartments is, in fact, an efficient way to maintain a product cell population from a small flux of progenitors, avoiding excessive clonality and minimizing the number of rounds of division en route. We use division, exit and death rates, estimated from measurements of single-positive thymocytes, to choose illustrative parameter values in the single-compartment case. We also consider a five-compartment model of thymocyte differentiation, from double-negative precursors to single-positive product cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Discrete Fracture Network Modeling to Estimate Upscaled Parameters for the Topopah Spring, Lava Flow, and Tiva Canyon Aquifers at Pahute Mesa, Nevada National Security Site

This report describes the results of Discrete Fracture Network (DFN) simulations for the Topopah Spring Aquifer (TSA), Lava Flow Aquifer, and Tiva Canyon Aquifer (TCA), at Pahute Mesa on the Nevada National Security Site (NNSS), formerly the Nevada Test Site. The research focuses on calculating upscaled groundwater flow and contaminant transport parameters using DFNs generated according to fracture characteristics observed in the TSA, LFA and TCA at Pahute Mesa. The highly fractured and heterogeneous nature of these aquifers makes them candidates for stochastic DFN modeling of radionuclide transport on a small scale with subsequent upscaling. One hundred independent DFN realizations are generated for each aquifer, and the upscaled parameters for continuum simulations of subsurface flow and transport in fractured media at Pahute Mesa are calculated. Our goal is to implement a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, such as channelization of flow and transport along a few well connected, large fractures. Additionally, to simulate advective and advective-diffusive transport through the fracture networks, the Time Domain Random Walk (TDRW) approach is applied to account for matrix diffusion into a finite half-space. Moreover, a novel approach to calculate dynamic (active) fracture surface area to reflect flow channeling is implemented. This work will improve the representation of radionuclide transport processes in largescale, regulatory-focused models by providing estimates of hard-to-measure flow and contaminant transport parameters at large scales. In this report, we (1) show recent results of flow and transport simulations on multiple DFN realizations of the TSA, LFA, TCA; (2) discuss the resulting distributions of estimated upscaled parameters; (3) describe the estimation of upscaled parameters for an equivalent parallel-plate continuum model and (4) present a comparison between simulated transport from the equivalent continuum model and an actual DFN.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Impacts of Instrument Limitations on Estimated Raindrop Size Distribution, Radar Parameters, and Model Microphysics during Mei-Yu Season in East China

Instrumentation limitations on measured raindrop size distributions (DSDs) and their derived relations and physical parameters are studied through a comparison of the DSD measurements during mei-yu season in east China by four collocated instruments, that is, a two-dimensional video disdrometer (2DVD), a vertically pointing Micro Rain Radar (MRR), and two laser-optical OTT Particle Size Velocity (PARSIVEL) disdrometers (first generation: OTT-1; second generation: OTT-2). Among the four instruments, the 2DVD provides the most accurate DSD and drop velocity measurements, so its measured rainfall amount has the best agreement with the reference rain gauge. Other instruments tend to miss more small drops ( D < 1 mm), leading to inaccurate DSDs and a lower rainfall amount. The low rainfall estimation becomes significant during heavy rainfall. The impacts of instrument limitations on the microphysical processes (e.g., evaporation and accretion rates) and convective storm morphology are evaluated. This is important especially for mei-yu precipitation, which is dominated by a high concentration of small drops. Furthermore, the instrument limitations need to be taken into account in both QPE and microphysics parameterization.

47 OTHER INSTRUMENTATION↗
NASA NTRS2006

Improved Soundings and Error Estimates using AIRS/AMSU Data

AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU A and HSB, to form a next generation polar orbiting infrared and microwave atmospheric sounding system. The primary products of AIRS/AMSU are twice daily global fields of atmospheric temperature-humidity profiles, ozone profiles, sea/land surface skin temperature, and cloud related parameters including OLR. The sounding goals of AIRS are to produce 1 km tropospheric layer mean temperatures with an rms error of 1 K, and layer precipitable water with an rms error of 20 percent, in cases with up to 80 percent effective cloud cover. The basic theory used to analyze AIRS/AMSU/HSB data in the presence of clouds, called the at-launch algorithm, and a post-launch algorithm which differed only in the minor details from the at-launch algorithm, have been described previously. The post-launch algorithm, referred to as AIRS Version 4.0, has been used by the Goddard DAAC to analyze and distribute AIRS retrieval products. In this paper we show progress made toward the AIRS Version 5.0 algorithm which will be used by the Goddard DAAC starting late in 2006. A new methodology has been developed to provide accurate case by case error estimates for retrieved geophysical parameters and for the channel by channel cloud cleared radiances used to derive the geophysical parameters from the AIRS/AMSU observations. These error estimates are in turn used for quality control of the derived geophysical parameters and clear column radiances. Improvements made to the retrieval algorithm since Version 4.0 are described as well as results comparing Version 5.0 retrieval accuracy and spatial coverage with those obtained using Version 4.0.

Susskind, Joel↗
NASA NTRS1979

The relationship of unsteadiness in downwash to the quality of parameter estimates

This paper investigates the relative importance of including unsteady effects in the lift and downwash in the longitudinal dynamics and parameter extraction algorithm. A simple vortex system has been used to model unsteady aerodynamic effects into the longitudinal equations of motion of an aircraft. Computer-generated data and flight data were used to demonstrate that inclusion of unsteady aerodynamics in the parameter-extraction algorithm produced aerodynamic parameters that were different from those extracted when unsteady aerodynamics were left out of the algorithm. The differences between derivatives associated with the two extraction algorithms (with and without unsteady aerodynamics) were related to acceleration derivatives which usually cannot be extracted individually.

Wells, W. R.↗
DOE OSTI2017

Understanding the Impact of Oligomeric Polystyrene Side Chain Arrangement on the All‐Polymer Solar Cell Performance

Abstract The introduction of oligomeric polystyrene (PS) side chains into the conjugated backbone is proven to enhance the processability and electronic properties of semiconducting polymers. Here, two series of donor and acceptor polymers are prepared with different molar percentages of PS side chains to elucidate the effect of their substitution arrangement on the all‐polymer solar cell performance. The observed device performance is lower when the PS side chains are substituted on the donor polymer and higher when on the acceptor polymer, indicating a clear arrangement effect of the PS side chain. The incorporation of PS side chains to the acceptor polymer contributes to the decrease in phase separation domain size in the blend films. However, the reduced domain size was still an order of magnitude larger than the typical exciton diffusion length. A detailed morphological study together with the estimation of solubility parameter of the pristine PS, donor, and acceptor polymers reveals that the relative value of solubility parameter of each component dominantly contributes to the purity of the phase separated domain, which strongly impacts the amount of generated photocurrent and overall solar cell performance. This study provides an understanding of the design strategies to improve the all‐polymer solar cells.

Kurosawa, Tadanori↗
DOE OSTI2023

Smart thermostat data-driven U.S. residential occupancy schedules and development of a U.S. residential occupancy schedule simulator

Occupancy schedule is one of the key inputs in Building Energy Modeling (BEM) to reflect the interaction between buildings and occupants. Over the past decades, standardized occupancy schedules, developed mainly by engineering rule-of-thumb, have been widely used in BEM due to its simplicity and lack of real measured occupancy data. However, the BEM community has recognized their association with uncertainty and reliability in simulation results from BEM. This study introduces representative occupancy schedules in the U.S. residential buildings, derived from a large smart thermostat dataset and time-series K-means clustering, and an open-source tool to generate a stochastic residential occupancy schedule. Over 90,000 residential occupancy schedules were estimated from the ecobee Donate Your Data dataset. Then, the representative occupancy schedules were identified through clustering. This study further investigated the impacts of three parameters (day, house type, and state) on residential occupancy schedules. Then, a tool, the Residential Occupancy Schedule Simulator (ROSS), is developed using the representative occupancy schedules derived in this study. Details of this tool are presented in this paper. In conclusion, the derived representative occupancy schedules and the ROSS tool can help improve the energy modeling of residential buildings.

42 ENGINEERING↗
DOE OSTI2017

Narrowing the window of inflationary magnetogenesis

We consider inflationary magnetogenesis where the conformal symmetry is broken by the term f {sup 2}(φ) F {sub αβ} F {sup αβ}. We assume that the magnetic field power spectrum today between 0.1 and 10{sup 4} Mpc is a power law, with upper and lower limits from observation. This fixes f to be close to a power law in conformal time in the window during inflation when the modes observed today are generated. In contrast to previous work, we do not make any assumptions about the form of f outside these scales. We cover all possible reheating histories, described by an average equation of state −1/3 < w-bar < 1. Requiring that strong coupling and large backreaction are avoided both at the background and perturbative level, we find the bound δ {sub B} {sub 0} < 5 ×10{sup −15} ( r /0.07){sup 1/2} κG for the magnetic field generated by inflation, where r is the tensor-to-scalar ratio and κ is a constant related to the form of f . This estimate has an uncertainty of one order of magnitude related to our approximations. The parameter κ is < 100, and values ∼> 1 require a highly fine-tuned form of f ; typical values are orders of magnitude smaller.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Frequency-domain full-waveform inversion with non-linear descent directions

Full-waveform inversion (FWI) is a highly non-linear inverse problem, normally solved iteratively, with each iteration involving an update constructed through linear operations on the residuals. Incorporating a flexible degree of non-linearity within each update may have important consequences for convergence rates, determination of low model wavenumbers and discrimination of parameters. In this study, we examine one approach for doing so, wherein higher order scattering terms are included within the sensitivity kernel during the construction of the descent direction, adjusting it away from that of the standard Gauss–Newton approach. These scattering terms are naturally admitted when we construct the sensitivity kernel by varying not the current but the to-be-updated model at each iteration. Linear and/or non-linear inverse scattering methodologies allow these additional sensitivity contributions to be computed from the current data residuals within any given update. We show that in the presence of pre-critical reflection data, the error in a second-order non-linear update to a background of s 0 is, in our scheme, proportional to at most (Δs/s 0 ) 3 in the actual parameter jump Δs causing the reflection. In contrast, the error in a standard Gauss–Newton FWI update is proportional to (Δs/s 0 ) 2 . For numerical implementation of more complex cases, we introduce a non-linear frequency-domain scheme, with an inner and an outer loop. A perturbation is determined from the data residuals within the inner loop, and a descent direction based on the resulting non-linear sensitivity kernel is computed in the outer loop. We examine the response of this non-linear FWI using acoustic single-parameter synthetics derived from the Marmousi model. Lastly, the inverted results vary depending on data frequency ranges and initial models, but we conclude that the non-linear FWI has the capability to generate high-resolution model estimates in both shallow and deep regions, and to converge rapidly, relative to a benchmark FWI approach involving the standard gradient.

58 GEOSCIENCES↗