Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic State 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 703 records · Page 39

Recent Rates of Forest Harvest and Conversion in North America

Incorporating ecological disturbance into biogeochemical models is critical for estimating current and future carbon stocks and fluxes. In particular, anthropogenic disturbances, such as forest conversion and wood harvest, strongly affect forest carbon dynamics within North America. This paper summarizes recent (2000.2008) rates of extraction, including both conversion and harvest, derived from national forest inventories for North America (the United States, Canada, and Mexico). During the 2000s, 6.1 million ha/yr were affected by harvest, another 1.0 million ha/yr were converted to other land uses through gross deforestation, and 0.4 million ha/yr were degraded. Thus about 1.0% of North America fs forests experienced some form of anthropogenic disturbance each year. However, due to harvest recovery, afforestation, and reforestation, the total forest area on the continent has been roughly stable during the decade. On average, about 110 m3 of roundwood volume was extracted per hectare harvested across the continent. Patterns of extraction vary among the three countries, with U.S. and Canadian activity dominated by partial and clear ]cut harvest, respectively, and activity in Mexico dominated by conversion (deforestation) for agriculture. Temporal trends in harvest and clearing may be affected by economic variables, technology, and forest policy decisions. While overall rates of extraction appear fairly stable in all three countries since the 1980s, harvest within the United States has shifted toward the southern United States and away from the Pacific Northwest.

Masek, Jeffrey G.↗

Data-driven modeling and control of dynamical systems using Koopman and Perron-Frobenius operators

This dissertation studies the data-driven modeling and control problem of nonlinear systems by exploiting the linear operator theoretic framework involving Koopman and Perro-Frobenius operator. A systematic linear-operator based controller design procedure has been established, which can be used to solve a variety of nonlinear control problems, including feedback stabilization using control Lyapunov functions, optimal quadratic regulation using Koopman eigenfunctions and convex optimization formulation of optimal control problem using P-F and Koopman operator approximation. As the core of data-driven modeling, we first propose a new algorithm for the finite-dimensional approximation of the linear transfer Koopman and Perron-Frobenius operator from time-series data. We argue that the existing approach for the finite-dimensional approximation of these transfer operators such as Dynamic Mode Decomposition (DMD) and Extended Dynamic Mode Decomposition (EDMD) do not capture two important properties of these operators, namely positivity and Markov property. The algorithm we propose preserves these two properties. We call the proposed algorithm as naturally structured DMD (NSDMD) since it retains the inherent properties of these operators. Naturally structured DMD algorithm leads to a better approximation of the steady-state dynamics of the system regarding computing Koopman and Perron- Frobenius operator eigenfunctions and eigenvalues. However, preserving positivity property is critical for capturing the real transient dynamics of the system. This positivity property of the transfer operators and it's finite-dimensional approximation play an important role for controller and estimator design of nonlinear systems. To solve the feedback stabilization problem for nonlinear control systems, we tried to take advantage of the Koopman operator framework. The Koopman operator approach provides a linear representation for a nonlinear dynamical system and a bilinear representation for a nonlinear control system. The problem of feedback stabilization of a nonlinear control system is then transformed to the stabilization of a bilinear control system. We propose a control Lyapunov function (CLF)-based approach for the design of stabilizing feedback controllers for the bilinear system. The search for finding a CLF for the bilinear control system is formulated as a convex optimization problem. This leads to a schematic procedure for designing CLF-based stabilizing feedback controllers for the bilinear system and hence the original nonlinear system. Another advantage of the proposed controller design approach outlined in this dissertation is that it does not require explicit knowledge of system dynamics. In particular, the bilinear representation of a nonlinear control system in the Koopman eigenfunction space can be obtained from time-series data. Next, we study the optimal quadratic regulation problem for nonlinear systems. The linear operator theoretic framework involving the Koopman operator is used to lift the dynamics of nonlinear control system to an infinite-dimensional bilinear system. The optimal quadratic regulation problem for nonlinear system is formulated in terms of the finite-dimensional approximation of the bilinear system. A convex optimization-based approach is proposed for solving the quadratic regulator problem for bilinear system. We applied a variety of examples and compared the simulation results between our framework and conventional LQR control using linearized model. For more general optimal control problems, we provide a density-function based convex formulation for the optimal control problem of the nonlinear system. The convex formulation relies on the duality result in the stability theory of a dynamical system involving density function and Perron-Frobenius operator. The optimal control problem is formulated as an infinite-dimensional convex optimization program. The finite-dimensional approximation of the optimization problem relies on the recent advances made in the data-driven computation of the Koopman operator, which is dual to the Perron-Frobenius operator. Simulation results are presented to demonstrate the application of the developed framework.

Huang, Bowen↗

Leveraging Cell Expansion Sensing in State of Charge Estimation: Practical Considerations

Measurements such as current and terminal voltage that are typically used to determine the battery’s state of charge (SOC) are augmented with measured force associated with electrode expansion as the lithium intercalates in its structure. The combination of the sensed behavior is shown to improve SOC estimation even for the lithium ion iron phosphate (LFP) chemistry, where the voltage–SOC relation is flat (low slope) making SOC estimation using measured voltage difficult. For the LFP cells, the measured force has a non-monotonic F–SOC relationship. This presents a challenge for estimation as multiple force values can correspond to the same SOC. The traditional linear quadratic estimator can be driven to an incorrect SOC value. To address these difficulties, a novel switching estimation gain is used based on determining the operating region that corresponds to the actual SOC. Moreover, a drift in the measured force associated with a shift of the cell SOC–expansion behavior over time is addressed with a bias estimator for the force signal. The performance of Voltage-based (V) and Voltage and Force-based (V&F) SOC estimation algorithms are then compared and evaluated against a desired ±5% absolute error bound of the SOC using a dynamic stress test current protocol that tests the proposed estimation scheme across wide range of SOC and current rates.

25 ENERGY STORAGE↗

Accelerating uncertainty quantification in incremental dynamic analysis using dimension reduction-based surrogate modeling

We propose a surrogate modeling framework based on dimension reduction to facilitate the quantification of seismic risk of structural systems in performance-based earthquake engineering. The framework adopts incremental dynamic analysis (IDA) for addressing hazard variability, and promotes significant computational efficiency improvement for propagating epistemic uncertainties associated with the structural models. It utilizes both linear and nonlinear dimension reduction approaches, equipped with inverse mappings, to learn a functional between the input parameter space (e.g., the epistemic uncertainties of the structure) to the high-dimensional output space created through the IDA implementation across different ground motions and seismic intensity levels. Polynomial chaos expansion is adopted as the surrogate model to learn this functional in the reduced space. A nine-story steel moment-resisting frame with uncertain structural properties is used as a testbed. Furthermore, we select the seismic fragility curves as a measure of the structure’s seismic performance, since it provides an estimate of the probability of entering specified damage states for given levels of ground shaking.

42 ENGINEERING↗

Open-Source Ambient-Signal-Based Inertia Monitoring Tool for Power systems [SWR-25-59]

This script, developed by NREL, implements an inertia estimation algorithm using a sliding window approach. The inertia (H) is estimated based on frequency deviations and active power deviations, utilizing the swing equation from power system dynamics. An example application is provided using the WECC 240-bus system.

Peng, Jiangkai [National Renewable Energy Laborato↗

Results from a wake-steering experiment at a commercial wind plant: investigating the wind speed dependence of wake-steering performance

Wake steering is a wind farm control strategy in which upstream wind turbines are misaligned with the wind to redirect their wakes away from downstream turbines, thereby increasing the net wind plant power production and reducing fatigue loads generated by wake turbulence. In this paper, we present results from a wake-steering experiment at a commercial wind plant involving two wind turbines spaced 3.7 rotor diameters apart. During the 3-month experiment period, we estimate that wake steering reduced wake losses by 5.6% for the wind direction sector investigated. After applying a long-term correction based on the site wind rose, the reduction in wake losses increases to 9.3%. As a function of wind speed, we find large energy improvements near cut-in wind speed, where wake steering can prevent the downstream wind turbine from shutting down. Yet for wind speeds between 6–8 m/s, we observe little change in performance with wake steering. However, wake steering was found to improve energy production significantly for below-rated wind speeds from 8–12 m/s. By measuring the relationship between yaw misalignment and power production using a nacelle lidar, we attribute much of the improvement in wake-steering performance at higher wind speeds to a significant reduction in the power loss of the upstream turbine as wind speed increases. Additionally, we find higher wind direction variability at lower wind speeds, which contributes to poor performance in the 6–8 m/s wind speed bin because of slow yaw controller dynamics. Further, we compare the measured performance of wake steering to predictions using the FLORIS (FLOw Redirection and Induction in Steady State) wind farm control tool coupled with a wind direction variability model. Although the achieved yaw offsets at the upstream wind turbine fall short of the intended yaw offsets, we find that they are predicted well by the wind direction variability model. When incorporating the expected yaw offsets, estimates of the energy improvement from wake steering using FLORIS closely match the experimental results.

17 WIND ENERGY↗

Allometric relationships between primary size measures and sapwood area for six common tree species in snow-dependent ecosystems in the Southwest United States

High-elevation, snow-dependent, semiarid ecosystems across southwestern United States are expected to be vulnerable to climate change, including drought and fire, with implications for various aspects of the water cycle. To that end, much less is known about the dynamics of transpiration, an important component of the water cycle across this region. At the individual-tree scale, transpiration is estimated by scaling mean sap flux density by the hydroactive sapwood area (SA). SA also remains a key factor in effectively scaling individual tree water-use to stand level. SA across large spatial scales is normally established by relating SA of a few trees to primary size measures, e.g., diameter at breast height (DBH), tree height ( H ), or canopy diameter (CD). Considering the importance of SA in scaling transpiration, the primary objective of this study was therefore to establish six species-specific (aspen, maple, white fir, ponderosa pine, Douglas fir, Englemann spruce) allometric relationships between SA and three primary size measures (DBH, CD, or H ) across two high-elevation, snow-dependent, semiarid ecosystems in New Mexico and Arizona. Based on multiple statistical criteria (coefficient of determination, index of agreement, Nash–Sutcliffe efficiency) and ease of measurement in the forest, we identified DBH as the primary independent variable for estimating SA across all sites. Based on group regression analysis, we found allometric relationships to be significantly ( p < 0.05) different for the same species (ponderosa pine, Douglas-fir) across different sites. Overall, our allometric relationships provide a valuable database for estimating transpiration at different spatial scales from sap flow data in some of our most vulnerable ecosystems.

Forestry↗

Online and Offline Identification of False Data Injection Attacks in Battery Sensors Using a Single Particle Model

The cells in battery energy storage systems are monitored, protected, and controlled by battery management systems whose sensors are susceptible to cyberattacks. False data injection attacks (FDIAs) targeting batteries’ voltage sensors affect cell protection functions and the estimation of critical battery states like the state of charge (SoC). Inaccurate SoC estimation could result in battery overcharging and over discharging, which can have disastrous consequences on grid operations. This paper proposes a three-pronged online and offline method to detect, identify, and classify FDIAs corrupting the voltage sensors of a battery stack. To accurately model the dynamics of the series-connected cells a single particle model is used and to estimate the SoC, the unscented Kalman filter is employed. FDIA detection, identification, and classification was accomplished using a tuned cumulative sum (CUSUM) algorithm, which was compared with a baseline method, the chi-squared error detector. Online simulations and offline batch simulations were performed to determine the effectiveness of the proposed approach. Throughout the batch simulations, the CUSUM algorithm detected attacks, with no false positives, in 99.83% of cases, identified the corrupted sensor in 97% of cases, and determined if the attack was positively or negatively biased in 97% of cases.

25 ENERGY STORAGE↗

Vertical Axis Wind Turbine (VAWT) with Thermoplastic Composite Blades

Rising concerns about global climate change, planet pollution and increasing use of energy have accelerated wider utilization of renewable energy sources. A recent trend in energy system development is decentralization of power generation systems. The opportunity for individuals or organizations to generate electricity locally decreases losses and part of the electricity price caused by the long-distance power transmission. Vertical axis wind turbines (VAWTs) cannot compete with traditional horizontal axis wind turbines (HAWTs) on the basis of aerodynamic performance or maximum power output generated for a given footprint. However, their ability to perform at a wide range of wind speeds, indifference to wind direction, simplicity of fabrication, transportability and maintainability make them ideal candidates for decentralized installations in urban and residential areas. The project team including Steelhead Composites (SHC), Colorado State University (CSU), National Renewable Energy Laboratory (NREL) and Arkema Inc. designed and fabricated a VAWT rotor assembly with thermoplastic composite blades using novel infusion and fabrication techniques. This full-scale VAWT rotor assembly was designed for a rated power output of 0.5-1 kW, as estimated through analytical computational fluid dynamics (CFD) methods. Those CFD simulations were combined with structural finite element analysis (FEA) to optimize the shape and composite layup of the wind blades. Much of the emphasis of this program was to characterize the recyclable thermoplastic material and composite structural design, in order to take full advantage of the material properties of such composites. Additionally, this project examined the potential of thermoplastic resin systems to transform the way VAWT rotor assemblies are joined both in the factory in in the field, by utilizing thermally welded joints to minimize field failure. The infused composite blades manufactured during this program are an excellent example of these benefits made possible by the Arkema thermoplastic resin system. This prototype blade manufacturing also provided the groundwork to inform decisions for larger scale manufacturing techniques applicable to a commercialization strategy. There is specific commercial interest for small scale VAWT’s that are competitively priced, perform to expectations, and have environmentally conscious end-of-life characteristics. Economic analysis was also conducted to compare designs and material usage, providing an estimated levelized cost of energy (LCOE) for the VAWT. The next steps for commercialization would include testing of the rotor assembly to verify predicted performance, specification of the balance of plant components (generator, inverter, etc.), and identification of a manufacturing partner with capability to scale-up production. This program has showed the potential for a commercially viable product through detailed analysis, material characterization and prototype fabrication.

17 WIND ENERGY↗

Reconstruction and robust reduced-order observation of flexible variables

This paper discusses reconstruction and estimation of flexible variables from multiple strain measurements for use in state feedback control of flexible manipulators. Reconstruction is proposed for obtaining flexible mode amplitudes from the measurements, and estimation for the modal velocities. Design of the observer for estimation of the velocities is discussed with regard to robust implementation. The performance of the observer is examined experimentally for several specifications of the error dynamics.

Hastings, Gordon G.↗

In-Time UAV Flight-Trajectory Estimation and Tracking Using Bayesian Filters

Rapid increase of UAV operation in the next decade in areas of on-demand delivery, medical transportation services, law enforcement, traffic surveillance and several others pose potential risks to the low altitude airspace above densely populated areas. Safety assessment of airspace demands the need for a novel UAV traffic management (UTM) framework for regulation and tracking of the vehicles. Particularly for low-altitude UAV operations, quality of GPS measurements feeding into the UAV is often compromised by loss of communication link caused by presence of trees or tall buildings in proximity to the UAV flight path. Inaccurate GPS locations may yield to unreliable monitoring and inaccurate prognosis of remaining battery life and other safety metrics which rely on future expected trajectory of the UAV. This work therefore proposes a generalized monitoring and prediction methodology for autonomous UAVs using in-time GPS measurements. Firstly, a typical 4D smooth trajectory generation technique from a series of waypoint locations with associated expected times-of-arrival based on B-spline curves is presented. Initial uncertainty in the vehicle's expected cruise velocity is quantified to compute confidence intervals along the entire flight trajectory using error interval propagation approach. Further, the generated planned trajectory is considered as the prior knowledge which is updated during its flight with incoming GPS measurements in order to estimate its current location and corresponding kinematic profiles. Estimation of position is denoted in dicrete state-space representation such that position at a future time step is derived from position and velocity at current time step and expected velocity at the future time step. A linear Bayesian filtering algorithm is employed to efficiently refine position estimation from noisy GPS measurements and update the confidence intervals. Further, a dynamic re-planning strategy is implemented to incorporate unexpected detour or delay scenarios. Finally, critical challenges related to uncertainty quantification in trajectory prognosis for autonomous vehicles are identified, and potential solutions are discussed at the end of the paper. The entire monitoring framework is demonstrated on real UAV flight experiments conducted at the NASA Langley Research Center.

Banerjee, Portia↗

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit &#x039B; on each local quantum number, and if the initial state has low local quantum numbers, then an error at most &#x03F5; can be achieved by choosing &#x039B; to scale polylogarithmically with &#x03F5; &#x2212; 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on &#x039B; that achieves accuracy &#x03F5; , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy &#x03F5; by truncating local quantum numbers at &#x039B; = polylog ( &#x03F5; &#x2212; 1 ) .

Tong, Yu↗

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Along-Trajectory Acoustic Signal Variations Observed During the Hypersonic Re-Entry of the OSIRIS-REx Sample Return Capsule

The re-entry of the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) sample return capsule (SRC) on 24 September 2023 presented a rare opportunity to study atmospheric entry dynamics through a dense network of ground-based infrasound sensors. As the first interplanetary capsule to re-enter over the United States since Stardust in 2006, this event allowed for unprecedented observations of infrasound signals generated during hypersonic descent. We deployed 39 single-sensor stations across Nevada and Utah, strategically distributed to capture signals from distinct trajectory points. Infrasound data were analyzed to examine how signal amplitude and period vary with altitude and propagation path for a nonablating hypersonic object with well-defined physical and aerodynamic properties. Raytracing simulations incorporated atmospheric specifications from the ground-2-space model to estimate source altitudes for observed signals. Results confirmed ballistic arrivals at all stations, with source altitudes ranging from 44 to 62 km along the trajectory. Signal period and amplitude exhibited strong dependence on source altitude, with higher altitudes corresponding to lower amplitudes, longer periods, and reduced high-frequency content. Regression analysis demonstrated strong correlations between signal characteristics and both altitude and propagation geometry. Our results suggest, when attenuation is considered, the amplitude is primarily determined by the source, with the propagation path playing a secondary role over the distances examined. These findings emphasize the utility of controlled SRC re-entries for advancing our understanding of natural meteoroid dynamics, refining atmospheric entry models, and improving methodologies for planetary defense. The OSIRIS-REx SRC campaign represents the most comprehensive infrasound study of a hypersonic re-entry to date, showcasing the potential of coordinated geophysical observational networks for high-energy atmospheric phenomena, including space debris re-entries.

58 GEOSCIENCES↗

Overview of Ice-Sheet Mass Balance and Dynamics from ICESat Measurements

The primary purpose of the ICESat mission was to determine the present-day mass balance of the Greenland and Antarctic ice sheets, identify changes that may be occurring in the surface-mass flux and ice dynamics, and estimate their contributions to global sea-level rise. Although ICESat's three lasers were planned to make continuous measurements for 3 to 5 years, the mission was re-planned to operate in 33-day campaigns 2 to 3 times each year following failure of the first laser after 36 days. Seventeen campaigns were conducted with the last one in the Fall of 2009. Mass balance maps derived from measured ice-sheet elevation changes show that the mass loss from Greenland has increased significantly to about 170 Gt/yr for 2003 to 2007 from a state of near balance in the 1990's. Increased losses (189 Gt/yr) from melting and dynamic thinning are over seven times larger'than increased gains (25 gt/yr) from precipitation. Parts of the West Antarctic ice sheet and the Antarctic Peninsula are losing mass at an increasing rate, but other parts of West Antarctica and the East Antarctic ice sheet are gaining mass at an increasing rate. Increased losses of 35 Gt/yr in Pine Island, Thwaites-Smith, and Marie-Bryd.Coast are more than balanced by gains in base of Peninsula and ice stream C, D, & E systems. From the 1992-2002 to 2003-2007 period, the overall mass balance for Antarctica changed from a loss of about 60 Gt/yr to near balance or slightly positive.

Zwally, H. Jay↗

Dynamic Light-Scattering Probe Used for the Very Early Detection of Cataracts and to Measure Response to Therapy (or Treatment)

Half of all blindness worldwide is due to cataracts, and 34 million Americans over the age of 65 have cataracts. This figure is expected to rise to 70 million by the year 2030. Currently, no medical treatment is available to prevent or halt the progression of a cataract; nor is there any way to reverse a cataract once it has been detected by conventional methods. The only known treatment is surgical removal of the lens. It is estimated that over $5 billion will be spent this year for the treatment of cataract patients in the United States alone. According to Carl Kupfer, M.D., former Director of the National Eye Institute at the National Institutes of Health, "A delay in cataract formation of about 10 years would reduce the prevalence of visually disabling cataracts by about 45 percent." However, this vision remains a dream because of the lack of objective methods that can be used to detect cataracts very early and, hence, can be used for screening potential anticataract drugs. This is about to change. The National Eye Institute and NASA entered into a formal interagency agreement in late 1996 to develop and test a dynamic light-scattering (DLS) device to detect the earliest changes due to cataract formation. The new DLS probe developed by Dr. Ansari at the NASA Glenn Research Center is several orders of magnitude more sensitive than conventional clinical systems. This ultrahigh sensitivity enables the probe to detect the onset of a cataract before it has any effect on vision and allows researchers to test the effectiveness of anticataract drugs in reversing cataracts. This device (see the following illustration) was originally designed to study transport phenomena in microgravity fluid physics experiments onboard the space shuttle orbiters and space station. It will be very useful in testing anticataract drugs to halt or reverse the progression of cataracts during longitudinal clinical trials.

Ansari,Rafat R.↗

Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation

Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quantifying the Opportunity Limits of Automatic Residential Electric Load Shaping

Electric utility residential demand response programs typically reduce load a few times a year during periods of peak energy use. In the future, utilities and consumers may monetarily and environmentally benefit from continuously shaping load by alternatively encouraging or discouraging the use of electricity. One way to shape load and introduce elasticity is to broadcast forecasts of dynamic electricity prices that orchestrate electricity supply and demand in order to maximize the efficiency of conventional generation and the use of renewable resources including wind and solar energy. A binary control algorithm that influences the on and off states of end uses was developed and applied to empirical time series data to estimate price-based instantaneous opportunities for shedding and adding electric load. To overcome the limitations of traditional stochastic methods in quantifying diverse, non-Gaussian, non-stationary distributions of observed appliance behaviour, recent developments in wavelet-based analysis were applied to capture and simulate time-frequency domain behaviour. The performance of autoregressive and spectral reconstruction methods was compared, with phase reconstruction providing the best simulation ensembles. Results show spatiotemporal differences in the amount of load that can be shed and added, which suggest further investigation is warranted in estimating the benefits anticipated from the wide-scale deployment of continuous automatic residential load shaping. Empirical data and documented software code are included to assist in reproducing and extending this work.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗