Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Results”

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 631 records · Page 35

Estimation of Participation Factors Using the Synchrosqueezed Wavelet Transform

This paper proposes a data-driven approach for estimating participation factors for a power system using only simulation results on selected disturbances. The approach is purely response-based and does not need a linearized system model for eigen-analysis, which makes it applicable to systems whose detailed, complete mathematical models are not available. Considering the unavoidable nonlinearity as exhibited in the transient period of a system response, the Synchrosqueezed Wavelet Transform is applied to simulated responses for modal analysis to obtain participation factors. Based on simulations of Kundur's two-area system using both the electromagnetic transient model and phasor model, the participation factors estimated by the proposed approach are compared with two other signal processing tools, the Prony analysis and continuous wavelet transform, and are also benchmarked with conventional model-based participation factors.

participation factor↗

Photonic ring resonator notch filters for astronomical OH line suppression

Photonic ring resonator arrays used as notch filters are a promising novel solution to improve the signal-to-noise ratio of ground-based astronomical observations by suppressing OH emission lines in the near-infrared (NIR) wavelength range (0.9-2.5 μm). We aim to fabricate a series of ring resonators connected by a waveguide, each with its resonance wavelength and full-width-half-maximum (FWHM) matched with one of the OH emission lines. The fundamental structure of ring resonator-based notch filters is a ring resonator evanescently coupled to a waveguide. Due to the optical loss of the ring resonator, the transmission of light through the waveguide will be suppressed at the resonance wavelength of the ring resonator. The degree of suppression can be controlled by varying the coupling strength between the waveguide and the ring resonator, while the resonance wavelength can be controlled by the geometry and size of the ring. However, since a single ring resonator has multiple modes each with a different wavelength separated by its free spectral range (FSR), the additional suppression at other wavelengths may undermine the signal. Therefore, the ideal FSR off the ring resonators needs to be larger than the measurement window. For this application, we aim the FSR to be 1 larger than 30 nm . Our initial measurement results revealed two major challenges for this application. First, the FSR of the ring resonators is limited by the bending loss of ring resonators. In a standard 220 nm silicon-on-insulator material platform, the FSR of ring resonators cannot exceed 30 nm; Second, due to the wafer thickness variation and fabrication inaccuracy, the achievable accuracy in the resonance wavelength of a fabricated single ring resonator is only ~1.5 nm. In this work, we report our novel solutions to the above two challenges. For the first one, we have developed a novel structure based on the Vernier effect to greatly raise the FSR of ring resonators. Serially coupled Vernier ring structures is a typical approach to increase the FSR without encountering bending loss, but it requires restrictive matching of the coupling and optical loss of the rings and any imperfect matching can result in even noisier signals. On the contrary, our structure, as shown in shown in Fig. 1, does not require the restrictive matching. The signals of two output ports in our structure marked in Fig. 1 will be added together to form the final output signal, and Fig. 2 shows the finite-difference time-domain (FDTD) simulation results of the final output signal. As can be seen, the FSR of the entire structure is well-beyond 60 nm with little noise at other wavelengths. For the second one, we have developed an easy-to-implement method to permanently tune the resonance wavelengths of ring resonators within 0.1 nm accuracy. The tuning mechanism is based on the property of hydrogen silsesquioxane (HSQ) resist. As a negative tone electron beam resist, HSQ upon exposure will cross-link and eventually form SiO2. The refractive 2 index of HSQ resist and formed SiO2 are slightly different , enabling the tuning. By having a thin layer of HSQ resist on top of ring resonator structures followed by SiO2 cladding, the resonance wavelengths of ring resonators can be post-tuned Advances in Optical Astronomical Instrumentation 2019, edited by Simon Ellis, Céline d'Orgeville, Downloaded From: https://www.spiedigitallibrary.org/conference-proceedings-of-spie on 23 Sep 2020 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use through electron beam exposure. We have found that, within 5 nm tuning range, the resonance wavelength shift and dosage applied obeys a linear relationship. So a simple calibration beforehand can result in resonance wavelength tuning within 0.1 nm accuracy. The transmission spectra of a 3-ring device before and after tuning is shown in Fig. 3.

Astronomical optics↗

Analysis of predicted mean vote-based model predictive control for residential HVAC systems

Model Predictive Control (MPC) is an advanced process control method that has attracted much attention in building heating, ventilation, and air conditioning (HVAC) systems. Here, this paper analyzes the optimal precooling performance in residential buildings using MPC with two different comfort indices, namely, temperature and predicted mean vote (PMV). It first formulates, for each comfort index, an optimization problem that accounts for different factors, such as weather, home thermal condition, prediction horizon, time-of use (TOU) utility rate, and rated cooling capacity. The problem is then solved, resulting in an MPC strategy that determines the HVAC on/off control signal and minimizes energy cost over a receding time horizon while maintaining thermal comfort. The energy performance difference between temperature-based and PMV-based MPC strategies is subsequently investigated, especially in light of the interior wall surface temperature and under different combinations of the factors. Extensive simulation results demonstrated that the proposed MPC strategies are adaptive and their performances depend primarily on weather, home thermal condition, and prediction horizon, while the impact of TOU utility rate and rated cooling capacity is relatively small. Because the PMV-based MPC strategy can take advantage of the lower interior wall surface temperature due to precooling, it resulted in 8–45% cost savings for the scenarios investigated and an average increase of 0.042–0.113 in the absolute value of the PMV index compared to the temperature-based MPC strategy.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Influence of acid leaching surface treatment on indentation cracking of soda lime silicate glass

Past work has shown that water or acid soaking treatments can increase the mechanical strength of soda lime silicate (SLS) glasses. In this work, we show that acid leaching treatments result in an increase in the apparent crack resistance of the acid-leached surface of SLS glass during indentation. Vickers indentation tests in controlled environments show a humidity dependence of radial cracking, suggesting that the transport of water through the leached layer plays a critical role in the propagation of cracks to the glass surface. Here, molecular dynamics simulations with reactive force fields indicate that the leached surface layer can undergo pressure-induced mechanochemical reactions during indentation, which increases the bridging oxygen connectivity in the silica network of the leached layer. Such structural changes can hinder transport of water molecules from the environment to the subsurface crack tip. Based on experimental observations and simulation results, a new hypothesis is proposed that mechanochemical restructuring in the leached layer in response to the applied load may lower the transport kinetics of molecular water to critical flaws, resulting in an apparent enhancement in the crack resistance of the acid-leached surface of SLS glass.

36 MATERIALS SCIENCE↗

Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

The aim of the study is to analyze the stability of the energetic particle modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. Additionally, the analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n = 1 EPM and n = 2 global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis ($\rlap{-} \iota_{0}$). In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given $\rlap{-} \iota_{0}$, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n = 1 EPM, n = 2 toroidal AE (TAE) and n = 3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP β threshold with respect to the co- and no-ECCD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Predictions of Cell-to-Cell Propagation and Vent Gas Production in the Thermal Runaway of Lithium-Ion Battery Stacks

This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.

25 ENERGY STORAGE↗

Performance Evaluation of Intelligent Solar Control Software Through Hardware-in-the-Loop (CRADA Final Report)

Recent research has highlighted the potential for solar to act as a zero-marginal-cost and zero-emission flexibility resource on the bulk power system when operated with advanced control systems. To increase the performance of these systems, leading technologies, including machine learning (ML) and hierarchical inverter set point allocation, have been developed by Latimer Controls, Inc. to estimate the headroom of large PV plants for grid operation and control; however, these technologies lack comprehensive validation under real-world application scenarios. Latimer Controls, Inc. received two voucher awards for research at a national laboratory from the Department of Energy American Made Solar Prize Round 6. The National Renewable Energy Laboratory (NREL) was selected to collaborate with Latimer staff to conduct a performance evaluation of Latimer PV control software. The NREL team will develop a hardware-in-the-loop (HIL) testbed to perform testing and validation of the Latimer PV control technology in a de-risked yet realistic testbed environment. Latimer and NREL worked together to analyze the test data, draw conclusions from the results, and disseminate the resulting scientific findings. In this CRADA work, we propose to test and validate the real-world application of the Latimer Control solution in an HIL environment. We evaluate the performance of different flexible solar technologies in responding to automatic generation control signals in a closed-loop fashion. In particular, a data-driven potential high limit (PHL) estimation is developed for large solar plants to accurately estimate their headroom so that they have fast and short-time regulation and control capability to participate in grid services and respond to grid signals in real time (e.g., AGC). This PHL estimation algorithm is embedded in a hardware power plant controller (PPC) and tested with an IEEE-39 bus system model developed in RTDS. To account for the varying cloud conditions and diverse inverter dispatches, we developed a 135-MW PV plant with detailed modeling of 27 individual PV modules and inverters using RTDS. The real-world communications used in such big plants, such as ModBus TCP/IP for inverter level and DNP3 for plant level, were developed to emulate the real-world applications in big PV plants. The ML-based PHL estimation method is tested under nine separate weather scenarios against the ‘reference-control’ solution, hereafter referred to as the baseline solution. The baseline method reserves a subset of inverters (reference group) to operate at their PHL at all times and dispatches only the remaining inverters (control group) at curtailed levels to fulfill the flexibility need. Despite being successfully piloted by NREL in California in 2017 and Chile in 2020, there exist two gaps in the state of the art to fully unlock the flexibility of PV plants: a. There is a trade-off between the PHL estimation accuracy and the flexibility range. b. There lacks granularity in the PHL estimation to capture the variation across inverters. The Latimer solution seeks to address these gaps by applying machine learning methods to improve PHL estimation accuracy while accounting for variability at every inverter. Performance metrics were taken from the 2023 Georgia Power CARES utility-scale RFP. The results demonstrate that the ML-based approach outperforms the traditional baseline method in PHL estimation accuracy for 7 of 9 scenarios. The average PHL error across the nine scenarios was 7.40% for the ML-based method, 2.06% less than the 9.46% PHL error average across scenarios that was exhibited by the baseline method. Additionally, the PHL error was below 5% for at least 95% of the testing interval for 3 of 9 tested intervals with the ML approach, whereas it did not achieve this metric for any of the baseline tests. Overall, simulation results indicate the superior performance of an ML-based approach compared to the conventional baseline reference-control approach, showcasing its potential to support grid stability and operational efficiency. This laboratory HIL testing using real PPC, representative power system simulation models in real-time with detailed PV plant and inverter models, and real-world communication protocols gives us confidence that this machine learning based PHL estimation algorithm works well in the hardware PPC and therefore de-risks future field commissioning. The end goal of this project is to advance grid technology to address the grid operation challenges brought by solar plant’s variability and uncertainties in power generation.

14 SOLAR ENERGY↗

A New Efficient Quantitative Multi-component Phase Field: Lattice Boltzmann Model for Simulating Ti6Al4V Solidified Dendrite Under Forced Flow

Ti6Al4V is a widely used, multi-component alloy in additive manufacturing, during which the fluid flow in the molten pool significantly affects the solidified dendrites. To predict and further control the microstructure, modeling and simulating the microstructure evolution play a critical role. In this study, a newly developed, efficient, quantitative multi-component phase-field (PF) model is coupled with a lattice Boltzmann (LB) model to simulate Ti6Al4V solidified dendrite evolution under fluid flow. The accuracy and convergence behavior of the model is validated by the Gibbs–Thomson relation at the dendrite tip. Single and multiple two-dimensional (2D) equiaxed dendrite evolution cases under forced flow were simulated. Results show that the dendrite pattern is influenced remarkably by the fluid flow. Underlying mechanisms of the asymmetrical evolution are revealed by discussing the interaction among the flow, composition distribution and dendrite morphology, quantitatively. The dendrite kinetics are also derived, which ascertains the relationship between tip velocity and undercooling and inlet velocity and is the foundation for larger-scale simulation. We believe that the coupled quantitative multi-component PF–LB framework employed in this study helps in investigating the solidified dendrite morphology evolution in a deep and quantitate manner.

36 MATERIALS SCIENCE↗

Techno-economic analysis of cathode material production using flame-assisted spray pyrolysis

The cost of cathode materials contributes approximately 32.7% of the total cell construction cost of lithium-ion batteries, significantly affecting the price of battery packs. To reduce the cathode material manufacturing cost, a flame-assisted spray pyrolysis (FSP) method has been developed to utilize a sustainable solvent of glycerol to manufacture the LiNi 1/3 Mn 1/3 Co 1/3 O 2 (NMC333) cathode materials. The purpose of this study is to evaluate the minimum cathode material selling price (MCSP) of the FSP processes compared with a traditional carbonate co-precipitation pathway. Results show that the MCSP of the FSP is $\$$19.1/kg that is 17% lower than the traditional carbonate co-precipitation pathway as a result of lower fixed operating cost and variable overhead. Sensitivity analysis shows that when the new process is integrated with in-situ sintering and processing, the MCSP can be as low as $\$$15.6/kg. When all the material prices are decreased by 20%, the FSP process can synthesize NMC333 at a price of $\$$2.3/kg lower. According to the simulation result, LiNi 0·8 Mn 0·1 Co 0·1 O 2 (NMC811) has the best potential to meet the U.S. Department of Energy battery price target of $\$$125/kWh, demonstrating that the FSP process is an attractive manufacturing technology for NMC cathode powder material production.

25 ENERGY STORAGE↗

Data assimilation for combustion ignition delay time simulation using schlieren image velocimetry

This study sought to improve the accuracy of simulating spray penetration and combustion ignition delay by means of data assimilation (DA). The simulations were conducted using the Reynolds-averaged Navier–Stokes (RANS) equations and assimilating the schlieren image data. In DA, an ensemble square root filter (EnSRF) was used to build the statistical model, making the simulation results more accurate without any change in the governing equations. Recognizing that the spray-cone injection angle has a large effect on penetration, we created ensemble members with different injection angles. And we applied the two-component velocity distribution calculated via SIV and updated both velocity and temperature by using a DA statistical model derived from RANS ensemble simulations. The ignition delay time is generally known to vary even under the same experimental conditions because it is influenced by many factors. In this study, we attempted the transient DA-assisted RANS simulation to predict the ignition delay time even when the temporal resolution and accuracy of the observation data ware insufficient. Our trials offer an example of how a combination of techniques can be effectively used to assimilate experimental data obtained under restricted conditions.

Combustion simulation↗

Renewable Integration in Hybrid ac-dc Systems using Multi-port Autonomous Reconfigurable Solar power plant (MARS)

To increase inertia and improve primary frequency response of future grids, grid-forming inverters connecting PV to grid and energy storage systems (ESSs) are likely to play an important role. Moreover, high-voltage direct current (HVdc) links will also be an enabler to transfer remote PV generation and for improved grid stability. That is, with increased penetration of PV, discrete development of PV and ESS connecting to transmission ac grid and HVdc links is one of the solutions for stable operation of grid. Here, an integrated concept for integration of PV and ESS to transmission ac grid and HVdc links is proposed that is named as multi-port autonomous reconfigurable solar power plant (MARS). The proposed integrated development reduces the cost associated with respect to discrete development of PV-ESS systems connecting to transmission ac grid and HVdc links. The integrated development incorporates advanced control methods to provide inertial and primary frequency response, reactive power support, and transient stability to manage PV and ESS resources. In this paper, high-fidelity switched system model of the integrated system and grids are developed and detailed simulation results are provided to showcase the stable operation of the integrated system and provision of grid support functions.

14 SOLAR ENERGY↗

Attack-Resilient Weighted $\ell_{1}$ Observer with Prior Pruning

Security related questions for Cyber Physical Systems (CPS) have attracted much research attention in searching for novel methods for attack-resilient control and/or estimation. Specifically, false data injection attacks (FDIAs) have been shown to be capable of bypassing bad data detection (BDD), while arbitrarily compromising the integrity of state estimators and robust controller even with very sparse measurements corruption. Moreover, based on the inherent sparsity of pragmatic attack signals, ℓ1 -minimization scheme has been used extensively to improve the design of attack-resilient estimators. For this, the theoretical maximum for the percentage of compromised nodes that can be accommodated has been shown to be 50%. In order to guarantee correct state recoveries for larger percentage of attacked nodes, researchers have begun to incorporate prior information into the underlying resilient observer design framework. For the most pragmatic cases, this prior information is often obtained through some data-driven machine learning process. Existing results have shown strong positive correlation between the tolerated attack percentages and the precision of the prior information. In this paper, we present a pruning method to improve the precision of the prior information, given corresponding stochastic uncertainty characteristics of the underlying machine learning model. Then a weighted ℓ1 -minimization is proposed based on the pruned prior. The theoretical and simulation results show that the pruning method significantly improves the observer performance for much larger attack percentages, even when moderately accurate machine learning model used.

Resilient observer, Cyber-physical systems, prunin↗

Fixed Bit Rotary Drilling Failure Criteria Effect on Drilling Vibration

Polycrystalline diamond compact (PDC) bit use has been expanded to geothermal drilling applications due to their improved performance in drilling compared to roller cone bits. Multiple vibration models have been developed yet and there are limited use of field data to study the ability of numerical vibration model to predict recorded vibration in the field. The objective of this paper is to evaluate the capabilities of two different rock-bit interaction models with two different ROP formulations for predicting the actual behavior seen in the field. The first model (Model-I) considers the rock-bit interaction, incorporating the effect of subsurface formation characteristics. The second model (Model-II) is a PDC bit model that uses calculated threshold weight-on-cutter (WOCt) to study the interaction of the bit with the rock surface. The models, which couple the drill stem axial and torsional motions, are solved numerically using the fourth-order Runge-Kutta method and compared with recorded downhole drilling data from the Chocolate Mountain Aerial Gunnery Range (CMAGR) geothermal well located in California. An 8–½" PDC bit with 8 blades was modeled and vibration simulation was performed to investigate stick-slip vibrations at 1400 ft. depth. Simulation results using field data indicate that Model-II encounters a low range of stick-slip severity (SSV) as compared to Model-I, which is in the moderate range. Comparison with the field value for SSV shows Model-II is the best representative of the field case for this depth. The sensitivity analysis shows that results are sensitive to model constants making prediction of SSV from models uncertain.

15 GEOTHERMAL ENERGY↗

Mechanical behavior analysis using small-size graphite disc compression testing with digital image correlation methods

The ASTM International Standard Test Method for Tensile Stress-Strain of Carbon and Graphite (ASTM C749) requires dog-bone shaped specimens. Unfortunately, specimen sizes required in ASTM C749 are much too large for irradiation within a test reactor and are incompatible with irradiation capsule volumes, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed as an alternative that provides a means for testing tensile properties using smaller specimens when geometric constraints are present. However, the important stress–strain relation throughout the test loading history that could be recorded using ASTM C749 is lost by using ASTM D8289 because of the difficulties involved with using traditional strain gauges or extensometers on such a small specimen. In this work, a digital image correlation (DIC) method, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 splitting tensile test on eight small samples of graphite grade Mersen 2114 to track their deformation during compressive loading. Initial results showed that the DIC technique can measure the surface displacement/strain on these small (Ø6 × 3 mm) graphite specimens with high accuracy and good repeatability, indicating that the DIC technique, combined with the small-disc splitting tensile test in ASTM D8289, can capture the strain evolution in the test. In conclusion, further comparison of measured strains with analytical and numerical simulation results provides additional understanding of the mechanical behaviour of this nuclear graphite material in complex stress states.

36 MATERIALS SCIENCE↗

Isotope mass scaling and transport comparison between JET Deuterium and Tritium L-mode plasmas

The dimensionless isotope mass scaling experiment between pure Deuterium and pure Tritium plasmas with matched $\rho^*$, $\nu^*$ , β n , q and $T_e/T_i$ has been achieved in JET L-mode with dominant electron heating (NBI+ohmic) conditions. 28% higher scaled energy confinement time $B_t\tau_{E,\mathrm{th}}/A$ is found in favour of the Tritium plasma. This can be cast in the form of the dimensionless energy confinement scaling law as $\Omega_i \tau_{E,\mathrm{th}} \sim A^{0.48 \pm 0.16}$. This significant isotope mass scaling is consequently seen in the scaled one-fluid heat diffusion coefficient $A \chi_{\mathrm{eff}}/B_t$ which is around 50% lower in the Tritium plasma throughout the whole plasma radius. The isotope mass dependence in the particle transport channel is negligible, supported also by the perturbative particle transport analysis with gas puff modulation. The comparison of the edge particle fuelling or ionisation profiles from the EDGE2D-EIRENE simulations show that the absolute density differences that are necessary for the dimensionless match in the confined plasma dominate over any isotope mass dependencies of particle fuelling and ionization profiles at the plasma edge. Local GENE simulation results indicate a mild anti-gyroBohm effect at $\rho_\mathrm{tor}$ = 0.6 and thereby a small isotope mass dependence in favour of Tritium on heat transport and a negligible effect on particle transport. A significant fraction of the isotope scaling and reduced heat transport observed in the Tritium plasma is not captured in the GENE and ASTRA-TGLF-SAT2 simulations by simply changing the isotope mass for the same input profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Correlation-Aided Robust Decentralized Dynamic State Estimation of Power Systems with Unknown Control Inputs

This paper proposes a correlation-aided robust adaptive unscented Kalman filter for power system decentralized dynamic state estimation with unknown inputs, termed as robust AUKF-UI. The temporal and spatial correlations among the unknown inputs are used to derive a vector auto-regressive (VAR) model in an adaptive manner. This VAR model is further integrated together with state transition and measurement models for joint state and unknown inputs estimation. This allows taking into account the implicit cross-correlations between the states and the unknown inputs. As a result, the rank requirement for unknown input vector estimation is relaxed and the local generator frequency measurement is not required. The temporal correlations of time series innovation vectors, predicted state and input vectors are also leveraged by the robust AUKFUI to detect, identify and process bad data. Without these correlations, it is very challenging to address bad data with unknown inputs. Simulation results carried out on the IEEE 39-bus system demonstrate that the proposed robust AUKF-UI achieves much better results than other methods in the presence of low measurement redundancy, strong nonlinearity, and bad data.

Zhao, Junbo↗

Sensitivity and Effective Parameterization of a Multi-Scale Model of Proton-Exchange-Membrane Water Electrolysis

The development of proton-exchange-membrane water electrolysis (PEMWE) is crucial to establishing a green-hydrogen infrastructure. Active research in hydrogen technologies is focused on minimizing costs, improving efficiency, and mitigating safety risks. In this study, sensitivity analyses were performed using a physics-based model and used to identify critical parameters capable of reducing overpotentials and safety risks. These parameters include operational (i.e. temperature), structural (i.e. thickness of the components), and material properties (i.e. exchange current density). Simulated results suggest that better engineering of the membrane electrode assemblies (MEA) can offset intrinsic limitations. Membrane thickness and conductivity yielded the largest percent change on performance. Gas crossover was observed to be sensitive to diffusion, solubility, and bubble formation. Current density regimes were identified, where each crossover parameter was most influential. MEA parameterization serves as a framework for theoretical PEMWE optimization simulations and input for design criteria for future research targets.

Dizon, A↗

MEGASIM: Distribution and Detection of Earth Trojan Asteroids

Abstract Using N -body simulation results from the MEGASIM data set, we present spatial distributions of Earth Trojan Asteroids and assess the detectability of the population in current and next-generation ground-based astronomical surveys. Our high-fidelity Earth Trojan Asteroid (ETA) distribution maps show never-before-seen high-resolution spatial features that evolve over timescales up to 1 Gyr. The simulation was synchronized to start times and timelines of two observational astronomy surveys: (1) the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) and (2) the Zwicky Transient Facility (ZTF). We calculate upper limits for the number of ETAs potentially observable with both the ZTF and LSST surveys. Due to the Yarkovsky Effect, we find no stable ETAs on billion-year timescales likely to be detected by any ETA survey, as no C-type or S-type ETAs (with H < 22 and H < 24, respectively) are likely to be stable on billion-year timescales, and ETAs large enough to remain stable on billion-year timescales are very rare relative to the rest of the ETA population. We find that a twilight ETA survey will not drastically increase the likelihood of individual ETA detection, but it would provide orders of magnitude more observations of select ETA populations. The null detection to date from ZTF restricts the potential ETA population to hundreds of objects larger than 100 m (at H ≈ 22), while a null detection by LSST will further restrict the ETA population to tens of objects larger than 100 m.

79 ASTRONOMY AND ASTROPHYSICS↗