Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Analyses of Model-Based Predictive Control for a S-CO 2 Brayton Cycle Power Converter

The supercritical carbon dioxide (S-CO 2 ) Brayton cycle is a promising alternative to steam cycles as an energy converter for advanced nuclear reactors as well as fossil-fired and concentrated solar power plants. For advanced reactors, such as Sodium-cooled Fast Reactors (SFRs), the S-CO 2 cycle offers the benefits of high cycle efficiency, very compact turbomachinery, and elimination of sodium-water reactions, all of which may enable improved economics. The efficiency gains calculated for the S-CO 2 cycle are particularly significant if advantage can be taken of the sharp changes in the carbon dioxide thermo-physical properties near the critical point. On the other hand, the same property variations give rise to challenges in the cycle design, analysis, and control. The onset of two-phase flow (where the minimum cycle temperature and pressure drop drop below the critical value) can degrade compressor performance and potentially cause damage to the compressor blades. In previous work, the feasibility of the Model-based Predictive Control (MPC) approach was assessed for improving regulation of the precooler outlet temperature. A linear reduced-order model for the precooler was generated and for this simplified standalone system operating with an MPC controller the response of outlet temperature was evaluated through simulations. These evaluations included the presence of temperature measurement time delay representing a degraded sensor. It was demonstrated that the MPC controller maintained stable precooler temperature for those delay cases where a proportional-integral controller became unstable. In the present work the simulations have been extended to evaluate the MPC controller operating under more realistic conditions, namely in the split-flow recompression cycle non-linear plant. This compares with the previous standalone reduced-order linear system described above. In the present work the MPC-based control scheme was implemented in a GPASS 1-D systems code representation of the split-flow recompression cycle. The performance of the controller was again investigated for degraded sensor conditions and compared to that of a proportional-integral controller. The results were qualitatively similar to the previous results with the MPC controller and again providing superior stability characteristics. This work provides a sound starting point for the development of model-based control systems. In particular, the MPC control approach applied to Multiple Input Multiple Output systems can be used to synchronize the operation of the different system components to improve control system dynamic performance.

20 FOSSIL-FUELED POWER PLANTS↗
NASA NTRSDate not supplied

Recent Improvements and Verification of A Full Body Model in Opensim

BACKGROUND: The dynamic feasibility [1,2] criterion, that the subject’s Center of Pressure (COP) be located within the Base of Support (BOS) which outlines the feet, has aided in assessing the stability of human motion recorded on earth while performing the recorded tasks in lunar gravity or as countermeasures exercises on a vibration isolation and stabilization system in microgravity. The convex hull of the BOS on the platform under the subject’s feet was estimated using virtual markers on the feet of the scaled subject model. The COP was calculated using the ground reaction forces and moments determined from motion capture data with biomechanical modeling tools [3]. Occasionally, large-amplitude oscillatory spikes or “artifacts” were observed in the subject’s linear and angular momentum derivatives, affecting some COP data derived from motion capture. The purpose of this investigation was to assess and improve the accuracy of model scaling and BOS estimation as well as to determine the efficacy of model adjustments in mitigating artifacts influencing motion capture-derived ground reaction force and COP results. METHODS: To aid evaluation of proposed process and model updates, motion capture data were collected for two subjects during unit test and range of motion trials, lunar tasks, and countermeasures exercise motions. Markers were added to the full body Plug-in Gait marker set [4] during data collection. New markers were placed on the front, back, sides, and top of the head to improve scaling using distances between marker pairs. Medial elbow markers were added to stabilize the upper arm during OpenSim Inverse Kinematics (IK) [5]. Finally, markers were added on the outer edge of the heels and on the outside edges of the first and last toes on each foot. These additional foot markers were made available to test new automated foot scaling techniques and to calculate the error between the subject’s estimated and recorded BOS. The modified unscaled OpenSim Full Body Rajagopal Model [6,7] was adjusted using some previously investigated techniques [8] to mitigate rapid shifts in joint angles occurring during IK, as these were found to cause the spike artifacts observed in subsequent stages of analysis. Since OpenSim models use Euler angles and rotation sequences, the arm axes of rotation were adjusted, and the pelvis order of rotation was changed to minimize the likelihood of encountering “gimbal lock” during common human motion. The model clavicle, arm, elbow, wrist, pelvis, and ankle angle limits were adjusted to better accommodate the full human range of motion seen in exercise and lunar data. The shoulder joint center was calculated using a “pivoting” algorithm [9], and both shoulder joint center and upper arm markers were included during IK to provide additional shoulder stability on a case-by-case basis. The quality of IK results was assessed by three criteria: minimizing error between recorded motion capture markers and model markers, checking for reasonable rates of change in joint angles between fames (i.e., no IK artifacts), and ensuring the absence of spikes in the inertial forces and angular momentum derivatives calculated using a custom OpenSim plugin [10]. RESULTS: The additional markers placed on the subject during data collection allowed the head and feet to be scaled more accurately using distances between new marker pairs. Scaling with BOS markers placed on the subject and removing the limit on subtalar angle resulted in more accurate BOS determination. Unrealistically large changes in joint angles between frames could be reduced by including clavicle, sternum, and medial elbow markers during IK. In cases with large arm ranges of motion, results could be further improved by running IK using medial elbow and virtual shoulder joint center markers. Model adjustments significantly improved the IK results affecting COP calculation and increased the accuracy of BOS estimation.

C A Bell↗
NASA NTRS2001

Analysis and visualization of single-trial event-related potentials

In this study, a linear decomposition technique, independent component analysis (ICA), is applied to single-trial multichannel EEG data from event-related potential (ERP) experiments. Spatial filters derived by ICA blindly separate the input data into a sum of temporally independent and spatially fixed components arising from distinct or overlapping brain or extra-brain sources. Both the data and their decomposition are displayed using a new visualization tool, the "ERP image," that can clearly characterize single-trial variations in the amplitudes and latencies of evoked responses, particularly when sorted by a relevant behavioral or physiological variable. These tools were used to analyze data from a visual selective attention experiment on 28 control subjects plus 22 neurological patients whose EEG records were heavily contaminated with blink and other eye-movement artifacts. Results show that ICA can separate artifactual, stimulus-locked, response-locked, and non-event-related background EEG activities into separate components, a taxonomy not obtained from conventional signal averaging approaches. This method allows: (1) removal of pervasive artifacts of all types from single-trial EEG records, (2) identification and segregation of stimulus- and response-locked EEG components, (3) examination of differences in single-trial responses, and (4) separation of temporally distinct but spatially overlapping EEG oscillatory activities with distinct relationships to task events. The proposed methods also allow the interaction between ERPs and the ongoing EEG to be investigated directly. We studied the between-subject component stability of ICA decomposition of single-trial EEG epochs by clustering components with similar scalp maps and activation power spectra. Components accounting for blinks, eye movements, temporal muscle activity, event-related potentials, and event-modulated alpha activities were largely replicated across subjects. Applying ICA and ERP image visualization to the analysis of sets of single trials from event-related EEG (or MEG) experiments can increase the information available from ERP (or ERF) data. Copyright 2001 Wiley-Liss, Inc.

Non-NASA Center↗
DOE OSTI2023

Nonlinear analysis and vibro-impact characteristics of a shaft-bearing assembly

Here this study investigates the nonlinear frequency response of a shaft-bearing assembly with vibro-impacts occurring at the bearing clearances. The formation of nonlinear behavior as system parameters change is examined, along with the effects of asymmetries in the nominal, inherently symmetric system. The primary effect of increasing the forcing magnitude or decreasing the contact gap sizes is the formation of grazing-induced chaotic solution branches occurring over a wide frequency range near each system resonance. The system's nominal setup has very hard contact stiffness and shows no evidence of isolas or superharmonic resonances over the frequency ranges of interest. Moderate contact stiffnesses cause symmetry breaking and introduce superharmonic resonance branches of primary resonances. Even if some primary resonances are not present due to the system's inherent symmetry, their superharmonic resonances still manifest. Branches of quasiperiodic isolas (isolated resonance branches) are also discovered, along with a cloud of isolas near a high-frequency resonance. Parameter asymmetries are found to produce a few significant changes in behavior: asymmetric linear stiffness, contact stiffness, and gap size could affect the behavior of primary resonant frequencies and isolas.

42 ENGINEERING↗
DOE OSTI2015

SU-E-T-118: Analysis of Variability and Stability Between Two Water Tank Phantoms Utilizing Water Tank Commissioning Procedures

Purpose: The commissioning criteria of water tank phantoms are essential for proper accuracy and reproducibility in a clinical setting. This study outlines the results of mechanical and dosimetric testing between PTW MP3-M water tank system and the Standard Imaging Doseview 3D water tank system. Methods: Measurements were taken of each axis of movement on the tank using 30 cm calipers at 1, 5, 10, 50, 100, and 200 mm for accuracy and reproducibility of tank movement. Dosimetric quantities such as percent depth dose and dose profiles were compared between tanks using a 6 MV beam from a Varian 23EX LINAC. Properties such as scanning speed effects, central axis depth dose agreement with static measurements, reproducibility of measurements, symmetry and flatness, and scan time between tanks were also investigated. Results: Results showed high geometric accuracy within 0.2 mm. Central axis PDD and in-field profiles agreed within 0.75% between the tanks. These outcomes test many possible discrepancies in dose measurements across the two tanks and form a basis for comparison on a broader range of tanks in the future. Conclusion: Both 3D water scanning phantoms possess a high degree of spatial accuracy, allowing for equivalence in measurements regardless of the phantom used. A commissioning procedure when changing water tanks or upon receipt of a new tank is nevertheless critical to ensure consistent operation before and after the arrival of new hardware.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2021

Hypersonic Boundary-Layer Transition on Blunted Cones at Angle of Attack

Experimental studies of cones at several high-speed facilities have demonstrated that, for small nosetip bluntness, transition onset over a circular cone moves upstream along the leeward side and downstream along the windward side, but this trend may be reversed at large bluntness values, where transition onset moves downstream along the leeward side and upstream along the windward side. A theoretical and numerical investigation is performed to characterize the effects of nose bluntness on disturbance amplification over the circular cone for several angles of attack, with the goal of understanding the potential physical mechanisms behind the experimental observations. The three-dimensional laminar basic states over a 1.5 m long, 7-degree half-angle cone with 9.525 mm nosetip radius are computed for selected angles of attack values and freestream conditions that are selected to match the Mach 10 experiments conducted within the Hypervelocity Wind Tunnel 9 at the Arnold Engineering Development Complex (AEDC). The solutions at a freestream unit Reynolds number of 17.1 million per meter are used to perform detailed instability analyses for angles of attack equal to 0, 1, 3, and 5 degrees. Results indicate that the linear amplification of stationary crossflow waves along inflection lines may begin to influence transition along the acreage of the cone for angles of attack equal to or larger than 5 degrees. The measured trend in transition front with respect to increasing angle of attack is found to be consistent with the predicted increase in the amplification factors for Mack mode disturbances along the streamline trajectories. The increase in Mack mode amplification along the windward ray for higher angles of attack is shown to be the result of a progressively earlier entropy-layer swallowing. Computations also indicate that the transition amplification factor along the windward ray is not constant and increases with the angle of attack and that the transition #-factors along the leeward ray are rather small. The nonmodal analysis for zero degrees angle of attack shows that entropy-layer disturbances with appreciably strong energy growth can coexist with Mack mode instabilities at the measured transition location.

Boundary Layer Transition↗
DOE OSTI2017

Discrete Adjoint Sensitivity Analysis of Hybrid Dynamical Systems With Switching [Discrete Adjoint Sensitivity Analysis of Hybrid Dynamical Systems]

Sensitivity analysis is an important tool for describing power system dynamic behavior in response to parameter variations. It is a central component in preventive and corrective control applications. The existing approaches for sensitivity calculations, namely, finite-difference and forward sensitivity analysis, require a computational effort that increases linearly with the number of sensitivity parameters. In this paper, we investigate, implement, and test a discrete adjoint sensitivity approach whose computational effort is effectively independent of the number of sensitivity parameters. The proposed approach is highly efficient for calculating sensitivities of larger systems and is consistent, within machine precision, with the function whose sensitivity we are seeking. This is an essential feature for use in optimization applications. Moreover, our approach includes a consistent treatment of systems with switching, such as dc exciters, by deriving and implementing the adjoint jump conditions that arise from state-dependent and time-dependent switchings. The accuracy and the computational efficiency of the proposed approach are demonstrated in comparison with the forward sensitivity analysis approach. In conclusion, this paper focuses primarily on the power system dynamics, but the approach is general and can be applied to hybrid dynamical systems in a broader range of fields.

42 ENGINEERING↗
NASA NTRS2017

Aerodynamic Coefficients from Aeroballistic Range Testing of Deployed- and Stowed-SIAD SFDT Models

This report documents a ballistic-range test campaign conducted in 2012 in order to estimate the aerodynamic stability characteristics of two configurations of the Supersonic Flight Dynamics Test (SFDT) vehicle prior to its initial flight in 2014. The SFDT vehicle was a test bed for demonstrating several new aerodynamic decelerator technologies then being developed under the Low-Density Supersonic Decelerator (LDSD) Project. Of particular interest here is the Supersonic Inflatable Aerodynamic Decelerator (SIAD), an inflatable attached torus used to increase the drag surface area of an entry vehicle during the supersonic portion of the entry trajectory. Two model configurations were tested in the ballistic range: one representing the SFDT vehicle prior to deployment of the SIAD, and the other representing the nominal shape with the SIAD inflated. Both models were fabricated from solid metal, and therefore, the effects of the flexibility of the inflatable decelerator were not considered. The test conditions were chosen to match, as close as possible, the Mach number, Reynolds number, and motion dynamics expected for the SFDT vehicle in flight, both with the SIAD stowed and deployed. For SFDT models with the SIAD stowed, 12 shots were performed covering a Mach number range of 3.2 to 3.7. For models representing the deployed SIAD, 37 shots were performed over a Mach number range of 2.0 to 3.8. Pitch oscillation amplitudes covered a range from 0.7 to 20.6 degrees RMS. Portions of this report (data analysis approach, aerodynamic modeling, and resulting aerodynamic coefficients) were originally published as an internal LDSD Project report [1] in 2012. In addition, this report provides a description of the test design approach, the test facility, and experimental procedures. Estimated non-linear aerodynamic coefficients, including pitch damping, for both model configurations are reported, and the shot-by-shot trajectory measurements, plotted in comparison with calculated trajectories based on the derived non-linear aerodynamic coefficients, are provided as appendices. Since the completion of these tests, two full-scale SFDT flights have been successfully conducted: one in June 2014 [2, 3], and one in June 2015 [3].

Ground Test↗
DOE OSTI2026

Performance and transport in the ARC tokamak

The ARC TM tokamak, a high-field (𝐵 𝑇 = 11.4 T) fusion power plant, under development by Commonwealth Fusion Systems, is studied using a suite of integrated modelling tools to predict its fusion power generation (𝑃𝑓⁡𝑢⁢𝑠), transport and confinement properties. Analysis is based off an ARC operational point scoped first with zero-dimensional (0-D) plasma operational contour (POPCON) modelling to produce 1.13 GW of fusion power. A suite of integrated modelling tools (TRANSP, ASTRA and TORAX) were applied to predict the performance and kinetic profiles of the ARC design point, yielding a range of predicted performance spanning from ∼900 to 1300 MW in rough quantitative agreement with POPCON predictions. The sensitivity of these results to uncertain modelling inputs was probed using scans of pedestal boundary conditions around EPED-predicted values (total pressure and temperature ratios), tungsten concentration and seperatrix density around their nominal assumptions. Pedestal pressure and pedestal top (𝑇 𝑖 /𝑇 𝑒 ) play a large role in 1.5-dimensional performance predictions, able to modify the predicted 𝑃 𝑓⁡𝑢⁢𝑠 by a factor of 2 within reasonable assumptions. High-fidelity core nonlinear gyrokinetic profile predictions, performed using CGYRO (Candy et al. 2016 J. Comput. Phys., vol. 324, pp. 73–93) coupled with the PORTALS (Rodriguez-Fernandez et al. 2024 Nucl. Fusion, vol. 64, 076034; Phys. Plasmas, vol. 31, 2024, 062501) framework, yield substantially lower performance (𝑃 𝑓⁡𝑢⁢𝑠 =677 MW) compared with 0-D and medium-fidelity modelling for nominal assumptions, showing that there is non-negligible uncertainty between models and that future work on SPARC may help resolve discrepancies. Lower overall performance results from significantly reduced volume-averaged densities and temperatures, along with reduced levels of density and temperature peaking. Turbulence and transport are largely dominated by ion temperature gradient across the profile, confirmed by both linear stability and the response of the nonlinear fluxes to changes in gradients, with some impact of kinetic ballooning modes in the deep core. This work represents one of the most complete scoping of potential fusion power plant conditions performed to date. The extensive integrated modelling provides confidence in ARC performance approaching 1 GW, while nonlinear gyrokinetic modelling results in open questions into the physics of density and temperature peaking in fusion-power-plant-relevant operational space. A discussion of results and the role that the SPARC tokamak (Creely et al. 2020 J. Plasma Phys., vol. 86, 865860502) will play in informing ARC design, performance and operation is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2005

Forecasting Lightning at Kennedy Space Center/Cape Canaveral Air Force Station, Florida

The Applied Meteorology Unit (AMU) developed a set of statistical forecast equations that provide a probability of lightning occurrence on Kennedy Space Center (KSC) I Cape Canaveral Air Force Station (CCAFS) for the day during the warm season (May September). The 45th Weather Squadron (45 WS) forecasters at CCAFS in Florida include a probability of lightning occurrence in their daily 24-hour and weekly planning forecasts, which are briefed at 1100 UTC (0700 EDT). This information is used for general scheduling of operations at CCAFS and KSC. Forecasters at the Spaceflight Meteorology Group also make thunderstorm forecasts for the KSC/CCAFS area during Shuttle flight operations. Much of the current lightning probability forecast at both groups is based on a subjective analysis of model and observational data. The objective tool currently available is the Neumann-Pfeffer Thunderstorm Index (NPTI, Neumann 1971), developed specifically for the KSCICCAFS area over 30 years ago. However, recent studies have shown that 1-day persistence provides a better forecast than the NPTI, indicating that the NPTI needed to be upgraded or replaced. Because they require a tool that provides a reliable estimate of the daily thunderstorm probability forecast, the 45 WS forecasters requested that the AMU develop a new lightning probability forecast tool using recent data and more sophisticated techniques now possible through more computing power than that available over 30 years ago. The equation development incorporated results from two research projects that investigated causes of lightning occurrence near KSCICCAFS and over the Florida peninsula. One proved that logistic regression outperformed the linear regression method used in NPTI, even when the same predictors were used. The other study found relationships between large scale flow regimes and spatial lightning distributions over Florida. Lightning, probabilities based on these flow regimes were used as candidate predictors in the equation development. Fifteen years (1 989-2003) of warm season data were used to develop the forecast equations. The data sources included a local network of cloud-to-ground lightning sensors called the Cloud-to-Ground Lightning Surveillance System (CGLSS), 1200 UTC Florida synoptic soundings, and the 1000 UTC CCAFS sounding. Data from CGLSS were used to determine lightning occurrence for each day. The 1200 UTC soundings were used to calculate the synoptic-scale flow regimes and the 1000 UTC soundings were used to calculate local stability parameters, which were used as candidate predictors of lightning occurrence. Five logistic regression forecast equations were created through careful selection and elimination of the candidate predictors. The resulting equations contain five to six predictors each. Results from four performance tests indicated that the equations showed an increase in skill over several standard forecasting methods, good reliability, an ability to distinguish between non-lightning and lightning days, and good accuracy measures and skill scores. Given the overall good performance the 45 WS requested that the equations be transitioned to operations and added to the current set of tools used to determine the daily lightning probability of occurrence.

Lambert, Winfred↗
NASA NTRS2023

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
NASA NTRSDate not supplied

A Path to Improving Simulated Properties of Low Clouds over the Beaufort Sea using Airborne In Situ Observations of Subgrid-Scale Variability

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. One of the main reasons for this problem are the possible unrealistic assumptions that models/reanalyses make about the subgrid-scale (SGS) variability of meteorological properties, as well as the relationship between SGS variability and grid-scale (GS) cloud properties. We utilize cloud and thermodynamic data of low level (primarily) liquid clouds collected from two aircraft campaigns conducted over the Beaufort Sea to better understand and characterize this problem. Examining data from the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign and the 1998 First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)–Arctic Cloud Experiment (ACE) reveals that GS cloud water variability is closely related with SGS distribution of total water (i.e. water vapor + cloud water). We examine the influence of three assumed SGS parameters on the estimation of GS cloud water: the width of the SGS PDF of total water (ΔQT), the shape of the SGS PDF, and the critical saturation ratio (Scrit), which partitions the SGS PDF into water that is available or not for conversion to cloud water. Both Scrit and ΔQT influence predicted GS cloud water strongly. ARISE and FIRE-ACE disagree somewhat in the details, but both campaigns support the possibility that more realistic representations of these two parameters may lead to more realistic GS cloud water. In particular, both parameters are sensitive to GS relative humidity, and so we examine the effects of allowing these quantities to vary as an empirically-derived linear function of GS relative humidity. In contrast, the shape of the PDF has little effect on the predicted cloud water – a surprising result that merits additional investigation.

J Brant Dodson↗
DOE OSTI2023

An Analysis of the Effects of Renewable Energy Intermittency on the 2030 Korean Electricity Market

Republic of Korea has unique geographical characteristics similar to those of an island, resulting in an isolated power system. For this reason, securing sufficient operating reserves for the system’s stability and reliability in the face of the intermittency of increasing variable renewable energy (VRE) is paramount, and this will pave the way to achieving the nation’s decarbonization target and carbon neutrality. However, the current reserve-operation method in Republic of Korea does not take into account energy-system conditions, such as the intermittency of the VRE. Therefore, this paper presents an analysis of the impact of changes in reserve-operation methods on the electricity market in the future Republic of Korean power system, with the increased levels of VRE that are currently envisioned. Specifically, three reserve-operation methods, including Korea’s current reserve-power-operation standards, were applied to the two power-system plans announced by the Korean government to analyze the annual generator operation and costs. The analysis results show that securing reserves proportional to the VRE would exert negative effects, such as increased power-generation costs and the curtailment of nuclear and VRE generation. These results can contribute to the estimation of operational reserves needed for high levels of VRE and to the design of new the Korean reserve market, to be introduced in 2025.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Understanding Self-Assembly and the Stabilization of Liquid/Liquid Interfaces: The Importance of Ligand Tail Branching and Oil-Phase Solvation

Organophosphorus-based ligands represent a versatile set of solvent extraction reagents whose chemical makeup plays an important role in extraction mechanism. We hypothesize that the branching of the extractant hydrophobic tail and its oil-phase solvation affect the liquid/liquid interfacial structure. Understanding the structure mediated adsorption and interfacial ordering becomes key in designing ligands with enhanced selectivity and efficiency for targeted extractions. We employed vibrational sum frequency generation spectroscopy and interfacial tension measurements to extract thermodynamic adsorption energies, map interfacial ordering, and rationalize disparate behaviors of model di-(2-ethylhexyl) phosphoric acid and dioctyl phosphoric acid ligands at the hexadecane water interface. With increased surface loading, ligands with branched hydrophobic tails formed stable interfaces at much lower concentrations than those observed for ligands with linear alkyl tails. The lack of an oil phase and associated solvation results in markedly different interfacial properties, and thus measurements made at air/liquid surfaces cannot be assumed to correlate with the processes occurring at buried liquid/liquid interfaces. Here, we attribute these differences in the surface mediated self-assembly to key variations in hydrophobic interactions and tail solvation taking place in the oil phase demonstrating that interactions in both the polar and nonpolar phases are essential to understand self-assembly and function.

36 MATERIALS SCIENCE↗
NASA NTRS2023

Automated Fluidics Device for Extraction and Quantification of miRNA Biomarkers From Blood

Radiation Assessment DuRing Exposure And long-Duration Spaceflight (RADREADS) demonstrates space-compatible point-of-care technology for quantitative biological monitoring of blood miRNA biomarkers in response to long-term low dose radiation exposure. This individualized monitoring approach will inform targeted treatment strategies to maximize medical resource utilization by accounting for individual susceptibility to radiation-related illnesses. As human spaceflight progresses beyond Earth’s magnetic shielding, radiation exposure poses a significant risk to astronaut health and safety. Extended operation in this environment comes with an increased risk of radiation exposure, leading to higher risks of radiation sickness, cancer, central nervous system effects, and degenerative diseases. While conventional physical dosimetry techniques capture radiation dose, individualistic susceptibility to radiation damage is varied. Multiple characteristics, including age, body weight, sex, genetics, and immune status, have been found to influence radiosensitivity (Liu et al. 2011, and Bouffler 2016). This differential response necessitates individualized monitoring and targeted treatment strategies to maximize medical resource utilization; however, a practical diagnostic platform for quantifying long-term, low dose radiation-induced tissue damage does not currently exist. MicroRNAs (miRNAs) are a class of small, non-coding RNAs that regulate gene expression by mediating the degradation of messenger RNA. The levels of particular miRNAs are influenced by biological processes such as inflammation and serve as biomarkers for a variety of conditions including cancer (Singh et al. 2017). MicroRNAs are found in various bodily fluids and are amenable to collection via liquid biopsies, providing a minimally invasive and easily quantifiable readout for a variety of radiosensitive reporters. A preliminary signature of 15 spaceflight sensitive miRNA has been identified in rodent and human studies, including miR-21-5p, miR-24-3p, miR-92a-3p, miR-17-5p, miR-16a-3p, miR-34a-3p, and miR-223-3p. These targets generally increased expression with radiation dose and linear energy transfer, though variation between individuals is not yet described. Current gaps in the field include a lack of understanding of longitudinal biological responses to long-term, low dose radiation exposure and the absence of space-compatible point-of-care technology for quantitative biological monitoring. In this body of work, we aim to develop an automated bleed-to-read system to process whole blood for the detection of miRNA biomarkers in order to monitor individualistic responses to radiation exposure. This will be achieved via separating serum (or plasma) from whole blood, followed by extraction, amplification, and quantification of the miRNA using a RT-qPCR reaction. Previously, the WetLab-2 hardware enabled execution of a RT-qPCR reaction aboard ISS; however, it is a manual system that requires crew manipulation and bulky components (Parra et al. 2017). To address these issues, automated fluid handling hardware was developed for each stage of sample preparation. Extraction of total RNA is achieved by sequentially pumping reagents through an off-the-shelf nucleic acid binding column (miRNeasy Serum/Plasma Advanced Kit, Qiagen). This approach eliminates several manual pipetting and centrifuging steps and limits the use of toxic chemicals commonly found in other sample processing techniques. The resulting elution will then be automatically dispensed for RT-qPCR analysis using a compact rotary qPCR (Mic qPCR Cycler, Bio Molecular Systems) that will improve spaceflight compatibility by removing bubbles from the detection region, another challenge highlighted by WetLab-2 (Parra et al. 2017). Efforts are also being made to simplify the RT-qPCR reaction to a 1-step air-dryable mix to improve long-term reagent stability at room temperature and reduce system complexity. By automating the RT-qPCR processes via microfluidic manipulation, RADREADS will reduce crewmember hands-on time and enable the personalized detection of radiation-induced tissue damage during long duration missions. Minimally invasive, longitudinal monitoring of individual’s response to radiation exposure will inform how the physiological system responds to long-term low dose space radiation and enables development of targeted countermeasures by the medical team. Ultimately, this portable technology will require minimal technical expertise and can also be used to monitor miRNA biomarkers associated with other diseases.

Tristen Head↗
DOE OSTI2020

Coupling between interfacial charge and mechanical deformation at high temperatures in ceramics. Final report

Basic research, under this project, has led to new fundamental discoveries in the effect of electric field (and temperature) on diffusion kinetics, electroluminescence and electronic conductivity in ceramic materials. For example, yttria stabilized zirconia sinters in just seconds well below 1000°C, while, at the same time, exhibiting a highly non-linear rise in conductivity, and intense electroluminescence. The phenomenon has been shown to occur in nearly all ceramics whether they are ionic or electronic conductors, or insulators. The rise in conductivity leads to Joule heating but in-situ measurements with a platinum standard at APS and Brookhaven (NSLS-II) have shown that the temperatures are much below what would have been needed to sinter in a few seconds. We have explored the hypothesis that huge concentrations of Frenkel pairs of vacancies and interstitials are produced. While these point defects accelerate sintering (by the interstitials being transported to the pores and the vacancies to the grain boundaries), they may also produce electron-hole pairs which can impart high electrical conductivity or recombine to emit photons. Many measurements, including in-situ measurements of lattice expansion are consistent with this hypothesis. However, the energy barrier for generating Frenkels in zirconia for example is much larger than can be expected at temperature and electrical fields in these experiments. Most recent work suggests that short wavelength phonons play a critical role in the genesis of the “flash” phenomenon. For example, the Debye temperature has been shown to be the lower bound for initiating flash. Molecular dynamics simulations have shown that large concentrations of Frenkel pairs can form under the proliferation of “optical” phonons at the edge of the Brillouin zone, that is with a wavelength at or below the lattice parameter. This project has led to the development of new experimental techniques as well as a software hardware interface for process and control and data analysis in real time, with a time scale of less than one second. This approach spells a new paradigm in materials science research. It potential impact on high rate manufacturing is self evident.

36 MATERIALS SCIENCE↗