Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear response”

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 595 records · Page 33

Additive friction stir deposition induced stress ripples in aluminum alloy

In additive manufacturing, modification of the viscosity of softened or molten metal plays a significant role in the development of mechanical properties distribution. Similar to soft polymers, molten or softened metal usually exhibits vastly different viscoelastic properties. During scratching tests, viscoelastic polymers undergo spatiotemporal stress development and relaxation, causing stress ripples to form on their surfaces. Here, it is observed that friction stir-based solid-state manufacturing can yield similar stress ripples on softened aluminum alloys due to the linear motion of the deposition tool. Softened aluminum alloys exhibit viscoelastic and viscoplastic responses similar to polymeric material during the solid-state deposition process. In the final product, residual stress ripples appear on the surface of the metal alloy due to periodic residual stress relaxation, leading to oscillatory behavior from the bulk mechanical properties of the material. The orientation of the in-built stress pattern formed in the bulk of the additively manufactured metal alloy underneath the visually evident onion ring-like surface pattern formation is diametrically in the opposite direction. The experimentally observed results are consistent with analytical calculation and numerical simulations in terms of behaviors.

36 MATERIALS SCIENCE↗

Data for "Soil carbon dynamics during drying vs. rewetting: Importance of antecedent moisture conditions"

This dataset contains data used for the paper "Soil carbon dynamics during drying vs. rewetting: importance of antecedent moisture conditions". The Related References field will be updated with a full citation when available.Soil moisture influences soil carbon dynamics, including microbial growth and respiration. The response of such ‘soil respiration’ to moisture changes is generally assumed to be linear and reversible, i.e. to depend only on the current moisture state. Current models thus do not account for antecedent soil moisture conditions when determining soil respiration or the available substrate pool. We conducted a laboratory incubation to determine how the antecedent conditions of drought and flood influenced soil organic matter (SOM) chemistry, bioavailability, and respiration. We sampled soils from an upland coastal forest, Beaver Creek, WA USA, and subjected them to drying and rewetting treatments. For the drying treatment, field moist soils were saturated and then dried to 75, 50, 35, and 5 % saturation. In the rewetting treatment, field moist soils were air-dried and then rewet to 35, 50, 75, and 100 % saturation. We measured respiration and water extractable organic carbon (WEOC) concentrations and used 1H-NMR and FT-ICR-MS to characterize the WEOC pool across the treatments. The drying vs. wetting treatment strongly influenced SOM bioavailability, as rewet soils (with antecedent drought) had greater WEOC concentrations and respiration fluxes compared to the drying soils (with antecedent flood). In addition, air-dry soils had the highest WEOC concentrations, and the NMR-resolved peaks showed a strong contribution of protein groups in these soils. Both NMR and FT-ICR-MS analyses indicated increased contribution of complex aromatic groups/molecules in the rewet soils, compared to the drying soils. We suggest that drying introduced organic matter into the WEOC pool via desorption of aromatic molecules and/or by microbial cell lysis, and this stimulated microbial mineralization rates. Our work indicates that even short-term shifts in antecedent moisture conditions can strongly influence soil C dynamics at the core scale. The predictive uncertainties in current soil models may be reduced by a more accurate representation of soil water and C persistence that includes a mechanistic and quantitative understanding of the impact of antecedent moisture conditions.This dataset contains a compressed (.zip) archive of the data and R scripts used for this manuscript. The dataset includes files in .csv and .txt format, which can be accessed and processed using MS Excel or R. NMR data are provided as raw output data (accessed in Bruker TopSpin or MestreNova) as well as the MestreNova-processed files. This archive can also be accessed on GitHub at https://github.com/kaizadp/hysteresis_and_soil_carbon (DOI: 10.5281/zenodo.4432885).

1H-NMR↗

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Nuclear Data Adjustment for Nonlinear Applications in the OECD/NEA WPNCS SG14 Benchmark -- A Bayesian Inverse UQ-based Approach for Data Assimilation

The Organization for Economic Cooperation and Development (OECD) Working Party on Nuclear Criticality Safety (WPNCS) proposed a benchmark exercise to assess the performance of current nuclear data adjustment techniques applied to nonlinear applications and experiments with low correlation to applications. This work introduces Bayesian Inverse Uncertainty Quantification (IUQ) as a method for nuclear data adjustments in this benchmark, and compares IUQ to the more traditional methods of Generalized Linear Least Squares (GLLS) and Monte Carlo Bayes (MOCABA). Posterior predictions from IUQ showed agreement with GLLS and MOCABA for linear applications. When comparing GLLS, MOCABA, and IUQ posterior predictions to computed model responses using adjusted parameters, we observe that GLLS predictions fail to replicate computed response distributions for nonlinear applications, while MOCABA shows near agreement, and IUQ uses computed model responses directly. We also discuss observations on why experiments with low correlation to applications can be informative to nuclear data adjustments and identify some properties useful in selecting experiments for inclusion in nuclear data adjustment. Performance in this benchmark indicates potential for Bayesian IUQ in nuclear data adjustments.

FOS: Computer and information sciences↗

Impact of ITCZ width on global climate: ITCZ-MIP

Abstract. The width of the Inter-Tropical Convergence Zone (ITCZ) affects tropical rainfall, Earth's albedo, large-scale circulation, and climate sensitivity. To better understand the ITCZ width and its effects on global climate, we present a protocol to force changes in ITCZ width in climate models. Starting from an aquaplanet configuration with a slab ocean, adding surface heat fluxes in the deep tropics forces the ITCZ to narrow, and subtracting them causes it to widen. The protocol successfully generates changes in ITCZ width in all four climate models used in this study. ITCZ width in each model responds linearly to forcing magnitude and sign. Comparing across the four climate models, a response to varying ITCZ width that is remarkably consistent among models is the ITCZ strength, which is greater the narrower the ITCZ. On the other hand, the effect of varying ITCZ width on climate sensitivity is divergent among our four models, varying even in sign. Results from this pilot study highlight the connections between surface fluxes, ITCZ width, and the wider climate. A comprehensive model intercomparison project (MIP) has the potential to advance understanding of both the physical processes shaping ITCZ width and its influence on remote atmospheric circulations and global climate.

Geology↗

Time-dependent SOLPS-ITER simulations of the tokamak plasma boundary for model predictive control using SINDy *

Abstract Time-dependent SOLPS-ITER simulations have been used to identify reduced models with the sparse identification of nonlinear dynamics (SINDy) method and develop model-predictive control of the boundary plasma state using main ion gas puff actuation. A series of gas actuation sequences are input into SOLPS-ITER to produce a dynamic response in upstream and divertor plasma quantities. The SINDy method is applied to identify reduced linear and nonlinear models for the electron density at the outboard midplane n e , s e p O M P and the electron temperature at the outer divertor T e , s e p d i v . Note that T e , s e p d i v is not necessarily the peak value of T e along the divertor. The identified reduced models are interpretable by construction (i.e. not black box), and have the form of coupled ordinary differential equations. Despite significant noise in T e , s e p d i v , the reduced models can be used to predict the response over a range of actuation levels to a maximum deviation of 0.5% in n e , s e p O M P and 5%–10% in T e , s e p d i v for the cases considered. Model retraining using time history data triggered by a preset error threshold is also demonstrated. A model predictive control strategy for nonlinear models is developed and used to perform feedback control of a SOLPS-ITER simulation to produce a setpoint trajectory in n e , s e p O M P using the integrated plasma simulator framework. The developed techniques are general and can be applied to time-dependent data from other boundary simulations or experimental data. Ongoing work is extending the approach to model identification and control for divertor detachment, which will present transient nonlinear behavior from impurity seeding, including realistic latency and synthetic diagnostic signals derived from the full SOLPS-ITER output.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust Output Feedback Control Design for Inertia Emulation by Wind Turbine Generators

Wind generation has gained widespread use as a renewable energy source. Most wind turbines and other renewables connected to the grid through converters result in a reduction in the natural inertial response to grid frequency changes. The doubly-fed induction generator (DFIG) can be controlled to compensate for this reduction and, in fact, provide faster response than traditional synchronous machines. This paper proposes to design observer based output feedback linear quadratic regulator (LQR) and H control laws to realize the inertia emulation function and deliver fast frequency support. Furthermore, the aim is to track the reference speed by a diesel synchronous generator (DSG) in order to reach the desired inertia. The control signal is computed based on a reduced order model using the balanced truncation technique. A comparison with selective modal analysis (SMA) and balanced truncation model reduction techniques is presented. Comprehensive results show the effective emulation of synthetic inertia by implementing the control laws on a nonlinear three- phase diesel-wind system. The proposed technique is analyzed for different short circuit ratio (SCR) scenarios.

17 WIND ENERGY↗

Late day measurement of excised branches results in uncertainty in the estimation of two stomatal parameters derived from response curves in Populus deltoides Bartr. x Populus nigra L.

Many terrestrial biosphere models depend on an understanding of the relationship between stomatal conductance and photosynthesis. However, unlike the measurement of photosynthetic parameters, such as the maximum carboxylation capacity, where standard methods (e.g., CO 2 response or ACi curves) are widely accepted, a consensus method for empirically measuring parameters representing stomatal response has not yet emerged. Most models of stomatal response to environment represent stomatal conductance as being bounded by a lower intercept parameter (g 0 ), and linearly scaled based on a multivariate term described by the stomatal slope parameter (g 1 ). Here we employ the widely used Unified Stomatal Optimization model, to test whether g 1 and g 0 parameters are impacted by the choice of measurement method, either on an intact branch, or a cut branch segment stored in water. We measured paired stomatal response curves on intact and excised branches of a hybrid poplar clone (Populus deltoides Bartr. x Populus nigra L. OP367), measured twice over a diurnal period. We found that predawn branch excision did not significantly affect measured g 0 and g 1 when measured within 4 h of excision. Measurement in the afternoon resulted in significantly higher values of g 1 and lower values of g0, with values changing by 55% and 56% respectively. Excision combined with afternoon measurement resulted in a marked effect on parameter estimates, with g1 increasing 89% from morning to afternoon and a 25% lower g 1 for cut branches than those measured in situ. We also show that in hybrid poplar the differences in parameter estimates obtained from plants measured under different conditions can directly impact models of canopy function, reducing modeled transpiration by 18% over a simulated 12.5-hour period. While these findings are only for a single isohydric woody species, our findings suggest that stomatal optimality parameters may not remain constant throughout the day.

54 ENVIRONMENTAL SCIENCES↗

Mixed-strategy approach to band-edge analysis and modeling in semiconductors

In semiconductor physics, the Urbach rule describes an exponential dependence of the UV-Vis absorption coefficient on photon energy at the band edge. The rule is so ubiquitously applied that a single common mechanism, variability in the local potential, is generally considered to be the cause of exponential band edges. Variability in the local potential may be caused by alloy inhomogeneity, structural defects, impurities, thermal excitation, or any other mechanism that interrupts ideal crystallinity and creates disorder. A crucial point of the Urbach rule is its development for band edges which approach linearity when plotted on a logarithmic scale. The rule was not developed in the context of nonlinear band edges with complex fine structure. Even so, this model has been applied extensively to band edges with such features. Using differential analysis and deconvolution, we develop an analytical technique capable of separating the components responsible for the fine structure observed. This approach reduces challenging signal-artifact convolution problems to a simple linear combination. Furthermore, we develop a model that can be applied to the nonlinear band-edge signals extracted. In this work, we show that in the limiting case of linearity, our model is consistent with and produces the same results as the linear Urbach model. Our model is capable of bypassing “tailing” as an indirect measure of the spread in the local potential and gives a direct measure of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Two-Stage Reinforcement Learning Policy Search for Grid-Interactive Building Control

This paper develops an intelligent grid-interactive building controller, which optimizes building operation during both normal hours and demand response (DR) events. To avoid costly on-demand computation and to adapt to non-linear building models, the controller utilizes reinforcement learning (RL) and makes real-time decisions based on a near-optimal control policy. Learning such a policy typically amounts to solving a hard non-convex optimization problem. We propose to address this problem with a novel global-local policy search method. In the first stage, an RL algorithm based on zero-order gradient estimation is leveraged to search for the optimal policy globally, due to its scalability and the potential to escape some poor performing local optima. The obtained policy is then fine-tuned locally to bring the first-stage solution closer to that of the original unsmoothed problem. Experiments on a simulated five-zone commercial building demonstrate the advantages of the proposed method over existing learning approaches. They also show that the learned control policy outperforms a pragmatic linear model predictive controller (MPC) and approaches the performance of an oracle MPC in testing scenarios. Using a state-of-the-art advanced computing system, we demonstrate that the controller can be learned and deployed within hours of training.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗

Collective excitations and quantum incompressibility in electron-hole bilayers

In this study, we apply quantum continuum mechanics to the calculation of the excitation spectrum of a coupled electron-hole bilayer. The theory expresses excitation energies in terms of ground-state intra- and interlayer pair correlation functions, which are available from quantum Monte Carlo calculations. The final formulas for the collective modes deduced from this approach coincide with the formulas obtained in the “quasilocalized particle approximation” by Kalman et al. [G. Kalman, V. Valtchinov, and K. I. Golden, Phys. Rev. Lett. 82 , 3124 (1999) ], and likewise the theory predicts the existence of gapped excitations in the charged channel, with the gap arising from electron-hole correlation. An immediate consequence of the gap is that the static density-density response function of the charged channel vanishes as q 2 for wave vector q → 0 , rather than linearly in q , as commonly expected. In this sense, the system is incompressible . This feature, which has no analog in the classical electron-hole plasma, is consistent with the existence of an excitonic ground state and implies the existence of a discontinuity in the chemical potential of electrons and holes when the numbers of electrons and holes are equal. It should be experimentally observable by monitoring the densities of electrons and holes in response to potentials that attempt to change these densities in opposite directions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

(U) PRAD0697 & PRAD0698: Complex Loading of CeO₂ Powder

Cerium(IV) oxide (CeO₂) powder is shock compressed using the Precision High Energy-density Liner Implosion eXperiment (PHELIX) platform. Experimental results are compared against several modeling approaches. Compaction behavior is best captured with a P-∝ model, which calculates CeO₂ powder bulk densities within 80-99% of experimental values but overpredicts densi cation at the cylindrical target's outer radius and center by up to 20%. Preliminary calculations suggest that accuracy could be increased with the inclusion of a coupled strength model. Several common computational modeling approaches for the shock compression response of granular materials and the magnetohydrodynamic (MHD) force upon the impactor/liner in pulsed power compression experiments are investigated and analyzed for their validity. The Bi-linear Ramp, P-∝ PACXP, and P-∝ Menikoff-Kober continuum compaction models are calibrated to planar impact Hugoniot data for CeO₂ powder and used to predict the powder's shock compaction response under non-planar shock wave compression. MHD calculations of the PHELIX pulsed power driver are performed using an idealized resistor-inductor-capacitor (RLC) circuit calibrated to previous experiments. All simulations are performed using the LANL code FLAG. Two validation experiments are computationally designed using the calibrated compaction and circuit models, executed using the PHELIX platform on CeO₂ targets with initial porous densities of 3.95 and 4.03 g/cm³, measured with proton radiography, and analyzed against the model predictions. The two P-∝ models more accurately describe CeO₂ powder densi cation than the Bi-linear Ramp model. However, the two P-∝ models overpredict bulk density of the shock compressed CeO₂ powder by up to 20% when the appropriate impact velocities are applied. MHD calculations for both validation experiments underpredict liner impact velocities by 4-11% when using the idealized RLC circuit model calibrated to previous experiments. Compensating underpredictions of impact velocity and overpredictions of powder densication lead to a false accuracy in pre-shot calculations compared to experimental data. To improve correlation between simulations and experiments, the following improvements are suggested: 1. A coupled strength model for CeO₂ powder that updates strength as a function of porosity and applied stress. 2. An improved MHD circuit model that more accurately captures the PHELIX machine.

36 MATERIALS SCIENCE↗

Nuclear Data Adjustment for Nonlinear Applications in the OECD/NEA WPNCS SG14 Benchmark—A Bayesian Inverse UQ-Based Approach for Data Assimilation

The Organisation for Economic Co-operation and Development Working Party on Nuclear Criticality Safety has proposed a benchmark exercise to assess the performance of current nuclear data adjustment techniques applied to nonlinear applications and experiments with low correlation to applications. This work introduces Bayesian inverse uncertainty quantification (IUQ) employing scientific machine learning surrogate models as a method for nuclear data adjustments in this benchmark, and compares IUQ to the more traditional methods of generalized linear least squares (GLLS) and Monte Carlo Bayes (MOCABA). Posterior predictions from IUQ showed agreement with GLLS and MOCABA for linear applications. Here, when comparing GLLS, MOCABA, and IUQ posterior predictions to computed model responses using adjusted parameters, we observe that the GLLS predictions failed to replicate the computed response distributions for nonlinear applications, while MOCABA showed near agreement, and IUQ used the computed model responses directly. We also discuss observations on why experiments with low correlation to applications can be informative to nuclear data adjustments and identify some properties useful in selecting experiments for inclusion in nuclear data adjustment. Performance in this benchmark indicates potential for Bayesian IUQ in nuclear data adjustments.

Bayesian calibration↗

EPR Investigation of Temporal and Thermal Stability of Irradiated Sweeteners

Abstract Electron paramagnetic resonance was used to study the dosimetric properties of commercial sweeteners (sucrose, saccharin, and stevia). The results concluded that sugar has dosimetric properties of interest, which makes it suitable for dose reconstruction, emergency response, and nuclear forensics applications. The radiation dose response was investigated by irradiating samples up to 10 Gy and found to be linear. The post-irradiation stability of samples was observed to be useful for the intended applications. The results showed that the post-irradiation intensities vary from one sample to another, with the general trend demonstrating a small increase in intensities shortly after irradiation followed by an approximately constant leveling out. Finally, the effect of isochronal annealing was studied by heating samples to various temperatures. The results showed an overall stable behavior of the signal below melting temperature.

Environmental Sciences & Ecology↗

Factors Affecting Override Behavior during Demand Flexibility from High-Resolution Smart Thermostat Data

In this study, we delineate key weather and demographic predictors of override behavior in residential buildings during connected thermostats demand response events. Anticipating and reducing overrides is critical to demand flexibility (DF) program success. We use high- dimensional fixed effects linear regression techniques on a large ecobee dataset for about 5,000 enrolled households in the United States. We identify critical weather (indoor and outdoor temperature), building (housing type), and occupant (previous DF overrides and previous DF event exposure) factors influencing the override patterns of customers during thermostat demand response. We also differentiate DF events by season to understand the different indoor and outdoor conditions that might influence seasonal override rates. We found significant differences in override rates between building types, with single-family semi-detached homes generally having the highest overrides. Having a history of overrides was additionally a critical factor in predicting occupant response to future DF events. Understanding these override differences is necessary for rural electric cooperatives and emerging DF programs without access to large DR historical data. Overall, we provide critical information on technical and local demographic characteristics that may be correlated to building type to influence strategies to reduce overrides and improve the adoption of DF technologies.

connected thermostats (CTs)↗

Colossal Cryogenic Electro‐Optic Response Through Metastability in Strained BaTiO 3 Thin Films

The search for thin film electro-optic materials that can retain superior performance under cryogenic conditions has become critical for quantum computing. Barium titanate thin films show large linear electro-optic coefficients in the tetragonal phase at room temperature, which is severely degraded down to ≈200 pm V −1 in the rhombohedral phase at cryogenic temperatures. There is immense interest in manipulating these phase transformations and retaining superior electro-optic properties down to liquid helium temperature. Utilizing the thermodynamic theory of optical properties, a large low-temperature electro-optic response is designed by engineering the energetic competition between different ferroelectric phases, leading to a low-symmetry monoclinic phase with a massive electro-optic response. The existence of this phase is demonstrated in a strain-tuned BaTiO 3 thin film that exhibits a linear electro-optic coefficient of 2516 ± 100 pm V −1 at 5 K, which is an order of magnitude higher than the best reported performance thus far. Importantly, the electro-optic coefficient increases by 100 × during cooling, unlike the conventional films, where it degrades. Further, at the lowest temperature, significant higher order electro-optic responses also emerge. These results represent a new framework for designing materials with property enhancements by stabilizing highly tunable metastable phases with strain.

barium titanate↗