Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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 991 records · Page 55

Real-Time Wave Energy Converter Control Using Instantaneous Frequency

Wave Energy Converters (WECs) rely on effective Power Take-Off (PTO) control strategies to maximize energy absorption under dynamic sea conditions. Traditional hydrodynamic modeling techniques may require computationally intensive convolution calculations, making real-time control implementation challenging. This paper presents an alternative approach by leveraging instantaneous frequency estimation to dynamically adjust PTO damping in response to varying wave frequencies. Two real-time frequency estimation methods are explored: the Hilbert Transform (HT) and Phase-Locked Loop (PLL). The Hilbert Transform method provides accurate frequency tracking but introduces a delayed response due to its dependence on causal data. Conversely, the PLL approach demonstrates strong potential in frequency tracking but requires careful gain tuning, particularly in complex sea states. Comparative evaluations across multiple test cases—including sinusoidal variations, amplitude steps, frequency step changes, and real-world JONSWAP spectrum waves—highlight the strengths and limitations of each method. The two different PTO control techniques across the various frequency estimation methods were tested under real-sea states using a state-space model of a point-absorbing Wave Energy Converter. The Capture Width Ratio (CWR) is used as a performance metric, with results showing that the HT achieves a 10.6% improvement, while the PLL estimation yields a 0.9% improvement relative to the fixed parameter control baseline. These results highlight the effectiveness of real-time frequency estimation in improving energy absorption compared to static control parameters.

WEC control↗

An Investigation of Poorly Studied Open Cluster NGC 4337 Using Multicolor Photometric and Gaia DR2 Astrometric Data

We present a comprehensive analysis (photometric and kinematical) of the poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia DR2 databases. By determining the membership probabilities of stars, we identified the 624 most probable members with membership probability higher than 50% by using proper motion and parallax data taken from Gaia DR2. The mean proper motion of the cluster is obtained as μ{sub x}=−8.83±0.01 and μ {sub y} = 1.49 ± 0.006 mas yr{sup −1}. We find the normal interstellar extinction toward the cluster region. The radial distribution of members provides a cluster radius of 7.′75 (5.63 pc). The estimated age of 1600 ± 180 Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as 1.46±0.18 within the mass range 0.75−2.0 M{sub ⊙}, which is in fair agreement with Salpeter’s value (x = 1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates (A, D) by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are (A {sub conv}, D {sub conv}) = (96.°27 ± 0.°10, 13.°14 ± 0.°27) and (A {sub ◦}, D {sub ◦}) = (100.°282 ± 0.°10, 9.°577 ± 0.°323) respectively. We also computed the Velocity Ellipsoid Parameters, matrix elements (μ {sub ij}), direction cosines (l {sub j}, m {sub j}, n {sub j}), and the Galactic longitude of the vertex (l {sub 2}).

79 ASTRONOMY AND ASTROPHYSICS↗

Development of an open-source regional data assimilation system in PEcAn v. 1.7.2: application to carbon cycle reanalysis across the contiguous US using SIPNET

Abstract. The ability to monitor, understand, and predict the dynamics of the terrestrial carbon cycle requires the capacity to robustly and coherently synthesize multiple streams of information that each provide partial information about different pools and fluxes. In this study, we introduce a new terrestrial carbon cycle data assimilation system, built on the PEcAn model–data eco-informatics system, and its application for the development of a proof-of-concept carbon “reanalysis” product that harmonizes carbon pools (leaf, wood, soil) and fluxes (GPP, Ra, Rh, NEE) across the contiguous United States from 1986–2019. We first calibrated this system against plant trait and flux tower net ecosystem exchange (NEE) using a novel emulated hierarchical Bayesian approach. Next, we extended the Tobit–Wishart ensemble filter (TWEnF) state data assimilation (SDA) framework, a generalization of the common ensemble Kalman filter which accounts for censored data and provides a fully Bayesian estimate of model process error, to a regional-scale system with a calibrated localization. Combined with additional workflows for propagating parameter, initial condition, and driver uncertainty, this represents the most complete and robust uncertainty accounting available for terrestrial carbon models. Our initial reanalysis was run on an irregular grid of ∼ 500 points selected using a stratified sampling method to efficiently capture environmental heterogeneity. Remotely sensed observations of aboveground biomass (Landsat LandTrendr) and leaf area index (LAI) (MODIS MOD15) were sequentially assimilated into the SIPNET model. Reanalysis soil carbon, which was indirectly constrained based on modeled covariances, showed general agreement with SoilGrids, an independent soil carbon data product. Reanalysis NEE, which was constrained based on posterior ensemble weights, also showed good agreement with eddy flux tower NEE and reduced root mean square error (RMSE) compared to the calibrated forecast. Ultimately, PEcAn's new open-source regional data assimilation framework provides a scalable workflow for harmonizing multiple data constraints and providing a uniform synthetic platform for carbon monitoring, reporting, and verification (MRV) as well as accelerating terrestrial carbon cycle research.

54 ENVIRONMENTAL SCIENCES↗

Film Thickness and Composition Effects in Symmetric Ternary Block Copolymer/Homopolymer Blend Films: Domain Spacing and Orientation

We report thickness and homopolymer molar mass dependencies in thin film blends of a compositionally symmetric 75 kg/mol lamellar polystyrene-block-poly(methyl methacrylate) diblock copolymer (PS-b-PMMA) with near-equal volume fractions of PS and PMMA homopolymers on “neutral” surfaces that promote a vertical domain orientation. PS and PMMA homopolymers with equal molar masses of 1, 3, and 22 kg/mol are blended in a 1:1 mass ratio. To efficiently explore this parameter space, we prepared combinatorial samples with gradient homopolymer concentrations using electrospray deposition (ESD) and gradient film thicknesses using flow coating; these samples were subsequently characterized by grazing-incidence small-angle X-ray scattering and scanning electron microscopy, respectively. Our results demonstrate that the added homopolymers increase or reduce domain spacings with respect to the neat (unblended) block copolymer depending on their molar masses, in line with previous reports; a simple heuristic to estimate the domain scaling in ternary blends based on the homopolymer molar mass agrees well with data for the thinnest blend films studied here. More surprisingly, we observe that the measured domain spacing in blends also depends on film thickness, with thicker films exhibiting larger domain spacings than thinner films at the same homopolymer concentration. Furthermore, thickness gradients at fixed homopolymer concentrations reveal a change in lamellae orientation from vertical to horizontal as film thickness is increased beyond the nominal lamellar domain spacing. Coarse-grained molecular dynamics simulations indicate that this change in orientation is induced by homopolymer segregation to film surfaces. Finally, tying these effects together, we hypothesize that the observed dependence of domain spacing on film thickness is a consequence of changes in the vertical homopolymer composition profile in thicker films during lamellae reorientation and the conformational asymmetry between PMMA and PS segments.

36 MATERIALS SCIENCE↗

On interfacial properties of tetrahydrofuran: Atomistic and coarse-grained models from molecular dynamics simulation

We have determined the interfacial properties of tetrahydrofuran (THF) from direct simulation of the vapor-liquid interface. The molecules are modeled using six different molecular models, three of them based on the united-atom approach and the other three based on a coarse-grained (CG) approach. In the first case, THF is modeled using the transferable parameters potential functions approach proposed by Chandrasekhar and Jorgensen [J. Chem. Phys. 77, 5073 (1982)] and a new parametrization of the TraPPE force fields for cyclic alkanes and ethers [S. J. Keasler et al., J. Phys. Chem. B 115, 11234 (2012)]. In both cases, dispersive and coulombic intermolecular interactions are explicitly taken into account. In the second case, THF is modeled as a single sphere, a diatomic molecule, and a ring formed from three Mie monomers according to the SAFT-γ Mie top-down approach [V. Papaioannou et al., J. Chem. Phys. 140, 054107 (2014)]. Simulations were performed in the molecular dynamics canonical ensemble and the vapor-liquid surface tension is evaluated from the normal and tangential components of the pressure tensor along the simulation box. In addition to the surface tension, we have also obtained density profiles, coexistence densities, critical temperature, density, and pressure, and interfacial thickness as functions of temperature, paying special attention to the comparison between the estimations obtained from different models and literature experimental data. The simulation results obtained from the three CG models as described by the SAFT-γ Mie approach are able to predict accurately the vapor-liquid phase envelope of THF, in excellent agreement with estimations obtained from TraPPE model and experimental data in the whole range of coexistence. However, Chandrasekhar and Jorgensen model presents significant deviations from experimental results. We also compare the predictions for surface tension as obtained from simulation results for all the models with experimental data. The three CG models predict reasonably well (but only qualitatively) the surface tension of THF, as a function of temperature, from the triple point to the critical temperature. On the other hand, only the TraPPE united-atoms models are able to predict accurately the experimental surface tension of the system in the whole temperature range.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A method for the estimation of p-mode parameters from averaged solar oscillation power spectra

A new fitting methodology is presented that is equally well suited for the estimation of low-, medium-, and high-degree mode parameters from m-averaged solar oscillation power spectra of widely differing spectral resolution. This method, which we call the “Windowed, MuLTiple-Peak, averaged-spectrum” or WMLTP Method, constructs a theoretical profile by convolving the weighted sum of the profiles of the modes appearing in the fitting box with the power spectrum of the window function of the observing run, using weights from a leakage matrix that takes into account observational and physical effects, such as the distortion of modes by solar latitudinal differential rotation. We demonstrate that the WMLTP Method makes substantial improvements in the inferences of the properties of the solar oscillations in comparison with a previous method, which employed a single profile to represent each spectral peak. We also present an inversion for the internal solar structure, which is based upon 6366 modes that we computed using the WMLTP method on the 66 day 2010 Solar and Heliospheric Observatory/MDI Dynamics Run. To improve both the numerical stability and reliability of the inversion, we developed a new procedure for the identification and correction of outliers in a frequency dataset. We present evidence for a pronounced departure of the sound speed in the outer half of the solar convection zone and in the subsurface shear layer from the radial sound speed profile contained in Model S of Christensen-Dalsgaard and his collaborators that existed in the rising phase of Solar Cycle 24 during mid-2010.

79 ASTRONOMY AND ASTROPHYSICS↗

INITIAL SIZE DISTRIBUTION OF THE GALACTIC GLOBULAR CLUSTER SYSTEM

Despite the importance of their size evolution in understanding the dynamical evolution of globular clusters (GCs) of the Milky Way, studies that focus specifically on this issue are rare. Based on the advanced, realistic Fokker-Planck (FP) approach, we theoretically predict the initial size distribution (SD) of the Galactic GCs along with their initial mass function and radial distribution. Over one thousand FP calculations in a wide parameter space have pinpointed the best-fit initial conditions for the SD, mass function, and radial distribution. Our best-fit model shows that the initial SD of the Galactic GCs is of larger dispersion than today's SD, and that the typical projected half-light radius of the initial GCs is {approx}4.6 pc, which is 1.8 times larger than that of the present-day GCs ({approx}2.5 pc). Their large size signifies greater susceptibility to the Galactic tides: the total mass of destroyed GCs reaches 3-5 Multiplication-Sign 10{sup 8} M {sub Sun }, several times larger than previous estimates. Our result challenges a recent view that the Milky Way GCs were born compact on the sub-pc scale, and rather implies that (1) the initial GCs were generally larger than the typical size of the present-day GCs, (2) the initially large GCs mostly shrank and/or disrupted as a result of the galactic tides, and (3) the initially small GCs expanded by two-body relaxation, and later shrank by the galactic tides.

79 ASTRONOMY AND ASTROPHYSICS↗

Asymptotic Behavior of Memristive Circuits

The interest in memristors has risen due to their possible application both as memory units and as computational devices in combination with CMOS. This is in part due to their nonlinear dynamics, and a strong dependence on the circuit topology. We provide evidence that also purely memristive circuits can be employed for computational purposes. In the present paper we show that a polynomial Lyapunov function in the memory parameters exists for the case of DC controlled memristors. Such a Lyapunov function can be asymptotically approximated with binary variables, and mapped to quadratic combinatorial optimization problems. This also shows a direct parallel between memristive circuits and the Hopfield-Little model. In the case of Erdos-Renyi random circuits, we show numerically that the distribution of the matrix elements of the projectors can be roughly approximated with a Gaussian distribution, and that it scales with the inverse square root of the number of elements. This provides an approximated but direct connection with the physics of disordered system and, in particular, of mean field spin glasses. Using this and the fact that the interaction is controlled by a projector operator on the loop space of the circuit. We estimate the number of stationary points of the approximate Lyapunov function and provide a scaling formula as an upper bound in terms of the circuit topology only.

97 MATHEMATICS AND COMPUTING↗

Dynamic analysis of a circulating fluidized bed riser

A linear state model is proposed to analyze dynamic behavior of a circulating fluidized bed riser. Different operating regimes were attained with high density polyethylene beads at low and high system inventories. The riser was operated between the classical choking velocity and the upper transport velocity demarcating fast fluidized and transport regimes. At a given riser superficial gas velocity, the aerations fed at the standpipe were modulated resulting in a sinusoidal solids circulation rate that goes into the riser via L-valve. The state model was derived based on the mass balance equation in the riser. It treats the average solids fraction across the entire riser as a state variable. The total riser pressure drop was modeled using Newton’s second law of motion. The momentum balance equation involves contribution from the weight of solids and the wall friction caused by the solids to the riser pressure drop. The weight of solids utilizes the state variable and hence, the riser inventory could be easily calculated. The modeling problem boils down to estimating two parameters including solids friction coefficient and time constant of the riser. It has been shown that the wall friction force acts in the upward direction in fast fluidized regime which indicates that the solids were moving downwards on the average with respect to the riser wall. In transport regimes, the friction acts in the opposite direction. This behavior was quantified based on a sign of Fanning friction factor in the momentum balance equation. The time constant of the riser appears to be much higher in fast fluidized regime than in transport conditions.

42 ENGINEERING↗

Performance estimation of resonance-enhanced dual-buoy wave energy converter using coupled time-domain simulation

This article presents the modeling methodology and performance evaluation of the resonance-enhanced dual-buoy WEC (Wave Energy Converter) by HEM (hydrodynamic & electro-magnetic) fully-coupled-dynamics time-domain-simulation program. The numerical results are systematically compared with the authors’ 1/6-scale experiment. With a direct-drive linear generator, the WEC consists of dual floating cylinders and a moon-pool between the cylinders, which can utilize three resonance phenomena from moon-pool dynamics as well as heave motions of inner and outer buoys. The contact and friction between the two buoys observed in the experiment are also properly modeled in the time-domain simulation by the Coulomb-friction model. Moon-pool resonance peaks significantly exaggerated in linear potential theory are empirically adjusted through comparisons with measured values. A systematic comparative study between the simulations and experiments with and without PTO (power-take-off) is conducted, and the relative heave displacements/velocities and power outputs are well matched. Then, parametric studies are carried out with the simulation program to determine optimum generator parameters. The performance with various wave conditions is also assessed.

16 TIDAL AND WAVE POWER↗

Mass dependence of the activation enthalpy and entropy of unentangled linear alkane chains

The mass scaling of the self-diffusion coefficient D of polymers in the liquid state, D ∼ M{sup β}, is one of the most basic characteristics of these complex fluids. Although traditional theories such as the Rouse and reptation models of unentangled and entangled polymer melts, respectively, predict that β is constant, this exponent for alkanes has been estimated experimentally to vary from −1.8 to −2.7 upon cooling. Significantly, β changes with temperature T under conditions where the chains are not entangled and at temperatures far above the glass transition temperature T{sub g} where dynamic heterogeneity does not complicate the description of the liquid dynamics. Based on atomistic molecular dynamics simulations on unentangled linear alkanes in the melt, we find that the variation of β with T can be directly attributed to the dependence of the enthalpy ΔH{sub a} and entropy ΔS{sub a} of activation on the number of alkane backbone carbon atoms, n. In addition, we find a sharp change in the melt dynamics near a “critical” chain length, n ≈ 17. A close examination of this phenomenon indicates that a “buckling transition” from rod-like to coiled chain configurations occurs at this characteristic chain length and distinct entropy-enthalpy compensation relations, ΔS{sub a} ∝ ΔH{sub a}, hold on either side of this polymer conformational transition. We conclude that the activation free energy parameters exert a significant influence on the dynamics of polymer melts that is not anticipated by either the Rouse and reptation models. In addition to changes of ΔH{sub a} and ΔS{sub a} with M, we expect changes in these free energy parameters to be crucial for understanding the dynamics of polymer blends, nanocomposites, and confined polymers because of changes of the fluid free energy by interfacial interactions and geometrical confinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hyper-Velocity Nanoparticle Plasma Jet as Fast Probe for Runaway Electrons in Tokamak Disruptions

The large energy focalized deposition of the runaway electron (RE) beam on the first wall can potentially yield to destructive effects on the reactor vessel during the explosive instability of a major disruption in a large tokamak (e.g. RE beam estimation is multi-MA of ~10 MeV energy in ITER.) Impurity injection is seen as an effective way to suppress REs after thermal quench (TQ). Injected material in solid state is advantageous, but it must be in small size grain form, with large surface-to-mass ratio to facilitate its rapid ablation within plasma. However, the results using shattered pellet injection (SPI), a technique aiming to satisfy these condition, indicated that the relatively unexplored RE evolution phases and RE beam-plasma interaction are far from being properly diagnosed and fully understood. We proposed the idea to probe REs with a pulsed power hyper-velocity large-mass nanoparticle plasma jet (NPPJ). NPPJ provides a high specific surface area and swift injection to core, hence fast increase of three basic parameters: 1) the electron density ne (bound and free through ionization, if the cold post-TQ plasma has a still high enough T e ), 2) the effective charge (Z eff ) which rises both the dynamical friction force and plasma resistivity, and, finally, 3) the value of the (Dreicer and critical) threshold E-field above which electrons are accelerated and become REs. But the fast transient phenomena of RE ”seed” creation, RE ”avalanche” multiplication, and RE beam formation and interaction with the background plasma all take place dominantly in the core plasma, where the induction E-field has its maximum. Thus, the RE probe success requires a rapid (1 - 2 ms) injection of sufficient mass able to penetrate through the strong toroidal magnetic field (2 - 5 T) over a long distance (~1 - 2 m) with a large assimilation fraction in the core plasma. As a RE probe, NPPJ is a unique combination of fast response with trigger-to-delivery time of few milliseconds, large mass-velocity (>75 mg at several km/s) and deep injection into post-TQ central plasma. The Objectives of this Category #1 project (concept development) were as follows: 1: We investigated the essential basic physics aspects of creating a plasma with ~10 15 to 10 16 cm -3 boron nitride (BN) NP density suited for a pulsed power plasma gun acceleration to several km/s. We concluded that, for RE NPPJ probe, BN (diameter Ø10 to 100 nm NPs) are in all respects usable. 2: We advanced solutions to get BN NPPJ by adapting FAR-TECH’s technology for C 60 NPPJ. 3: We developed self-consistent models for RE current density j RE (R; t) evolution in the TQ aftermath, with fast and simultaneous increase of n e (R; t) and Z eff (R; t) > 1 upon the specific injection of C 60 (i.e. C ions) and BN NP (i.e. B and N ions), respectively. We showed that prompt post-TQ injection and high velocity are of essence in RE probe diagnostics and suppression. 4: We used the in-house HEM-2D code to simulate the C$_{60}^{q+}$ NPPJ transverse penetration through tokamak-level B-field (BT ~2 to 5 T) and distance (~1 to 2m), with fragmentation of C$_{60}^{q+}$ and developed the algorithm for BN Q+ NPPJ ablative sublimation of BN Q+ NP by the collisions with post-TQ plasma electrons. HEM-2D results indicated that fragmenting C$_{60}^{1+}$ NPPJ penetrates up to ~60 cm in 100 us when injecting at 1.6 ms after TQ and gradually released C 2 fragments and C ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Impact of Alfvénic Shear Flow on Magnetic Reconnection and Turbulence

Magnetic reconnection is a fundamental and omnipresent energy conversion process in plasma physics. Novel observations of fields and particles from Parker Solar Probe (PSP) have shown the absence of reconnection in a large number of current sheets in the near-Sun solar wind. Using near-Sun observations from PSP encounters 4–11 (2020 January–2022 March), we investigate whether reconnection onset might be suppressed by velocity shear. We compare estimates of the tearing mode growth rate in the presence of shear flow for time periods identified as containing reconnecting current sheets versus nonreconnecting times, finding systematically larger growth rates for reconnection periods. Upon examination of the parameters associated with reconnection onset, we find that 85% of the reconnection events are embedded in slow, non-Alfvénic wind streams. We compare with fast, slow non-Alfvénic, and slow Alfvénic streams, finding that the growth rate is suppressed in highly Alfvénic fast and slow wind, and reconnection is not seen in these wind types, as would be expected from our theoretical expressions. These wind streams have strong Alfvénic flow shear, consistent with the idea of reconnection suppression by such flows. This could help explain the frequent absence of reconnection events in the highly Alfvénic, near-Sun solar wind observed by PSP. Finally, we find a steepening of both the trace and magnitude magnetic field spectra within reconnection periods in comparison to ambient wind. We tie this to the dynamics of relatively balanced turbulence within these reconnection periods and the potential generation of compressible fluctuations.

slow solar wind↗

The triple decomposition of a fluctuating velocity field in a multiscale flow

A new method for the triple decomposition of a multiscale flow, which is based on the novel optimal mode decomposition (OMD) technique, is presented. OMD provides low order linear dynamics, which fits a given data set in an optimal way and is used to distinguish between a coherent (periodic) part of a flow and a stochastic fluctuation. The method needs no external phase indication since this information, separate for coherent structures associated with each length scale introduced into the flow, appears as the output. The proposed technique is compared against two traditional methods of the triple decomposition, i.e., bin averaging and proper orthogonal decomposition. This is done with particle image velocimetry data documenting the near wake of a multiscale bar array. It is shown that both traditional methods are unable to provide a reliable estimation for the coherent fluctuation while the proposed technique performs very well. The crucial result is that the coherence peaks are not observed within the spectral properties of the stochastic fluctuation derived with the proposed method; however, these properties remain unaltered at the residual frequencies. This proves the method’s capability of making a distinction between both types of fluctuations. The sensitivity to some prescribed parameters is checked revealing the technique’s robustness. Additionally, an example of the method application for analysis of a multiscale flow is given, i.e., the phase conditioned transverse integral length is investigated in the near wake region of the multiscale object array.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Accurate atmospheric parameters at moderate resolution using spectral indices: Preliminary application to the marvels survey

Studies of Galactic chemical, and dynamical evolution in the solar neighborhood depend on the availability of precise atmospheric parameters (effective temperature T {sub eff}, metallicity [Fe/H], and surface gravity log g) for solar-type stars. Many large-scale spectroscopic surveys operate at low to moderate spectral resolution for efficiency in observing large samples, which makes the stellar characterization difficult due to the high degree of blending of spectral features. Therefore, most surveys employ spectral synthesis, which is a powerful technique, but relies heavily on the completeness and accuracy of atomic line databases and can yield possibly correlated atmospheric parameters. In this work, we use an alternative method based on spectral indices to determine the atmospheric parameters of a sample of nearby FGK dwarfs and subgiants observed by the MARVELS survey at moderate resolving power (R ∼ 12,000). To avoid a time-consuming manual analysis, we have developed three codes to automatically normalize the observed spectra, measure the equivalent widths of the indices, and, through a comparison of those with values calculated with predetermined calibrations, estimate the atmospheric parameters of the stars. The calibrations were derived using a sample of 309 stars with precise stellar parameters obtained from the analysis of high-resolution FEROS spectra, permitting the low-resolution equivalent widths to be directly related to the stellar parameters. A validation test of the method was conducted with a sample of 30 MARVELS targets that also have reliable atmospheric parameters derived from the high-resolution spectra and spectroscopic analysis based on the excitation and ionization equilibria method. Our approach was able to recover the parameters within 80 K for T {sub eff}, 0.05 dex for [Fe/H], and 0.15 dex for log g, values that are lower than or equal to the typical external uncertainties found between different high-resolution analyses. An additional test was performed with a subsample of 138 stars from the ELODIE stellar library, and the literature atmospheric parameters were recovered within 125 K for T {sub eff}, 0.10 dex for [Fe/H], and 0.29 dex for log g. These precisions are consistent with or better than those provided by the pipelines of surveys operating with similar resolutions. These results show that the spectral indices are a competitive tool to characterize stars with intermediate resolution spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Fast estimation of dense gas dispersion from multiple continuous CO 2 surface leakage sources for risk assessment

Surface leakage of CO 2 , and associated potential impacts on health, safety, and the environment (HSE) are considered hazards of geologic carbon sequestration (GCS). There are two challenges associated with impact assessment of CO 2 surface dispersion. First, the fact that CO 2 is a dense gas makes its dispersion in air a complex process. Rigorous numerical solutions for modeling concentration distributions are relatively time-consuming. Second, impact assessment requires consideration of uncertainty, e.g., quantification of how much uncertainty is propagated through input parameters to model outputs by carrying out large numbers of model runs. In order to assess the potential consequences of surface leakage of CO 2 , it is useful to have a model that executes very quickly for repeated model calculations (e.g., in Monte Carlo mode) of the atmospheric dispersion of CO 2 (concentrations as a function of space and time). In addition, the model should be able to handle multiple surface leakage sources. In this work, we have extended the nomograph approach of Britter and McQuaid (1988) for estimating dense gas plume length from single leakage source to multiple leakage sources. The method is very fast and therefore amenable to general system-level GCS risk assessment including uncertainty quantification within the framework of the National Risk Assessment Partnership (NRAP) Integrated Assessment Model (IAM). The method is conservative in that it assumes the wind could be from any direction, and it handles multiple sources by a simple superposition approach. The method produces results in reasonable agreement with a sophisticated computational fluid dynamics (CFD) code, but runs in a small fraction of the time.

58 GEOSCIENCES↗

Enhancing e-waste estimates: Improving data quality by multivariate Input–Output Analysis

Highlights: • A multivariate Input–Output Analysis method for e-waste estimates is proposed. • Applying multivariate analysis to consolidate data can enhance e-waste estimates. • We examine the influence of model selection and data quality on e-waste estimates. • Datasets of all e-waste related variables in a Dutch case study have been provided. • Accurate modeling of time-variant lifespan distributions is critical for estimate. - Abstract: Waste electrical and electronic equipment (or e-waste) is one of the fastest growing waste streams, which encompasses a wide and increasing spectrum of products. Accurate estimation of e-waste generation is difficult, mainly due to lack of high quality data referred to market and socio-economic dynamics. This paper addresses how to enhance e-waste estimates by providing techniques to increase data quality. An advanced, flexible and multivariate Input–Output Analysis (IOA) method is proposed. It links all three pillars in IOA (product sales, stock and lifespan profiles) to construct mathematical relationships between various data points. By applying this method, the data consolidation steps can generate more accurate time-series datasets from available data pool. This can consequently increase the reliability of e-waste estimates compared to the approach without data processing. A case study in the Netherlands is used to apply the advanced IOA model. As a result, for the first time ever, complete datasets of all three variables for estimating all types of e-waste have been obtained. The result of this study also demonstrates significant disparity between various estimation models, arising from the use of data under different conditions. It shows the importance of applying multivariate approach and multiple sources to improve data quality for modelling, specifically using appropriate time-varying lifespan parameters. Following the case study, a roadmap with a procedural guideline is provided to enhance e-waste estimation studies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗