Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient dynamic 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 937 records · Page 52

Using the System Complexity Metric (SCM) to Compare CO2 Removal Systems

A fundamental cause of difficulty in large engineering projects is their inherent complexity. An impression of complexity occurs if a system is simply difficult to understand, where there is no obvious mental model that correctly predicts its behavior. Higher system complexity is usually associated with higher cost and higher failure rate. Complexity is perceived if a system has many diverse components, multiple interactions and feedback loops, transients and dynamic behavior, and unanticipated failure modes. Identifying and removing these signs of complexity should improve performance and reduce the cost and failure rate. Complexity can be directly measured by the number of components and their interactions. The System Complexity Metric (SCM) is defined as the sum of the number of parts in a system, N, plus the number of the one-way interconnections between them, I. SCM = N + I. The SCM is easily determined by direct inspection of the system block diagram. SCM can be used to compare systems and to guide their redesign to reduce cost and failure rate. Carbon dioxide removal systems are analyzed using SCM, cost, and failure rate. As in previous work, cost is directly proportional to SCM and that failure rate increases as a power of SCM for large differences in SCM. The SCM ranking of carbon dioxide removal systems is the same as their ranking in detailed analysis and practice.

Harry W. Jones↗

Using the System Complexity Metric (SCM) to Compare CO2 Reduction Systems

A fundamental cause of difficulty in large engineering projects is their inherent complexity. An impression of complexity occurs if a system is simply difficult to understand, where there is no obvious mental model that correctly predicts its behavior. Higher system complexity is usually associated with higher cost and higher failure rate. Complexity is perceived if a system has many diverse components, multiple interactions and feedback loops, transients and dynamic behavior, and unanticipated failure modes. Identifying and removing these signs of complexity should improve performance and reduce the cost and failure rate. Complexity can be directly measured by the number of components and their interactions. The System Complexity Metric (SCM) is defined as the sum of the number of parts in a system, N, plus the number of the one-way interconnections between them, I. SCM = N + I. The SCM is easily determined by direct inspection of the system block diagram. SCM can be used to compare systems and to guide their redesign to reduce cost and failure rate. Carbon dioxide reduction systems are analyzed using SCM, cost, and failure rate. As in previous work, cost is directly proportional to SCM and that failure rate increases as a power of SCM for large differences in SCM. The SCM ranking of carbon dioxide reduction systems is the same as their ranking in detailed analysis and practice.

Harry W. Jones↗

Earth-Atmospheric Coupling During Strong Earthquakes by Analyzing MODIS Data

Interactions between the Earth and the atmosphere during major earthquakes (M greater than 5) are the subject of this investigation. Recently a mechanism has been proposed predicting the build-up of positive ground potentials prior to strong earthquake activity. Connected phenomena include: transient conductivity of rocks, injection of currents, possibly also electromagnetic emission and light emission from high points at the surface of the Earth. To understand this process we analyze vertical atmospheric profiles, land surface and brightness (temperature) data, using the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard NASA's Terra satellite launched in December 1999. MODIS covers the entire Earth every 1-2 days in 36 wavelength bands (20 visible and 16 infrared) at different spatial resolutions (250 m, 500 m, and 1 km). Using MODIS data we look for correlations between the atmospheric dynamics and solid Earth processes for the January 2001 strong earthquakes in San Salvador and India. As part of the build-up of positive grounds potential, an IR luminescence is predicted to occur in the 8-12 micrometer band. We use the MODIS data to differentiate between true "thermal" signals and IR luminescence. Indeed, on the basis of a temporal and spatial distribution analysis, a thermal anomaly pattern is found that appears to be related to the seismic activity. Aerosol content and atmospheric instability parameters also change when ground charges build up causing ion emission and leading to a thin aerosol layer over land. We analyze the aerosol content, atmospheric pressure, moisture profile and lifted index. Anomalous trends have been identified in few days prior to the main shocks. The significance of this observation should be explored further using other data sets.

Ouzounov, Dimitar↗

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗

Neutrino and Electromagnetic Signatures from Superluminous Supernovae: A Case Study for SN 2017egm

Superluminous supernovae (SLSNe) are rare transients that are ∼10–100 times more luminous than ordinary stellar explosions, reaching peak optical luminosities ∼10 44 –10 45 erg s −1 . The energy source powering SLSNe remains uncertain. In this work, we explore the multiwavelength and multimessenger signatures of the scenario in which SLSNe are powered by a newly born millisecond magnetar. We model the dynamical evolution and emission from the coupled system comprised of the magnetar, wind, nebula, and supernova ejecta, consistently evaluating the pair multiplicity of the wind and nebula regions, and the bulk wind Lorentz factor governing the e + − e − injection spectra in the nebula. We compute the thermal and nonthermal electromagnetic signatures, neutrino signatures, and investigate their detection prospects. For SN 2017egm, the nearest observed SLSNe, our prediction for high-energy gamma rays matches the recent detection by the Fermi Large Area Telescope. For neutrinos, using SN 2017egm a canonical SLSNe, we find that in the era of the Vera C. Rubin Observatory, a stacking analysis with upcoming neutrino observatories can lead to 3σ detection significance of neutrino events from a population of SLSNe within a decade of operation.

Mukhopadhyay, Mainak [Fermi National Accelerator L↗

Autonomous Multi-Sensor Coordination: The Science Goal Monitor

Many dramatic earth phenomena are dynamic and coupled. In order to fully understand them, we need to obtain timely coordinated multi-sensor observations from widely dispersed instruments. Such a dynamic observing system must include the ability to Schedule flexibly and react autonomously to sciencehser driven events; Understand higher-level goals of a sciencehser defined campaign; Coordinate various space-based and ground-based resources/sensors effectively and efficiently to achieve goals. In order to capture transient events, such a 'sensor web' system must have an automated reactive capability built into its scientific operations. To do this, we must overcome a number of challenges inherent in infusing autonomy. The Science Goal Monitor (SGM) is a prototype software tool being developed to explore the nature of automation necessary to enable dynamic observing. The tools being developed in SGM improve our ability to autonomously monitor multiple independent sensors and coordinate reactions to better observe dynamic phenomena. The SGM system enables users to specify what to look for and how to react in descriptive rather than technical terms. The system monitors streams of data to identify occurrences of the key events previously specified by the scientisther. When an event occurs, the system autonomously coordinates the execution of the users' desired reactions between different sensors. The information can be used to rapidly respond to a variety of fast temporal events. Investigators will no longer have to rely on after-the-fact data analysis to determine what happened. Our paper describes a series of prototype demonstrations that we have developed using SGM and NASA's Earth Observing-1 (EO-1) satellite and Earth Observing Systems' Aqua/Terra spacecrafts' MODIS instrument. Our demonstrations show the promise of coordinating data from different sources, analyzing the data for a relevant event, autonomously updating and rapidly obtaining a follow-on relevant image. SGM was used to investigate forest fires, floods and volcanic eruptions. We are now identifying new Earth science scenarios that will have more complex SGM reasoning. By developing and testing a prototype in an operational environment, we are also establishing and gathering metrics to gauge the success of automating science campaigns.

Koratkar, Anuradha↗

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗

Adaptive computational methods for SSME internal flow analysis

Adaptive finite element methods for the analysis of classes of problems in compressible and incompressible flow of interest in SSME (space shuttle main engine) analysis and design are described. The general objective of the adaptive methods is to improve and to quantify the quality of numerical solutions to the governing partial differential equations of fluid dynamics in two-dimensional cases. There are several different families of adaptive schemes that can be used to improve the quality of solutions in complex flow simulations. Among these are: (1) r-methods (node-redistribution or moving mesh methods) in which a fixed number of nodal points is allowed to migrate to points in the mesh where high error is detected; (2) h-methods, in which the mesh size h is automatically refined to reduce local error; and (3) p-methods, in which the local degree p of the finite element approximation is increased to reduce local error. Two of the three basic techniques have been studied in this project: an r-method for steady Euler equations in two dimensions and a p-method for transient, laminar, viscous incompressible flow. Numerical results are presented. A brief introduction to residual methods of a-posterior error estimation is also given and some pertinent conclusions of the study are listed.

Oden, J. T.↗

Impeller leakage flow modeling for mechanical vibration control

HPOTP and HPFTP vibration test results have exhibited transient and steady characteristics which may be due to impeller leakage path (ILP) related forces. For example, an axial shift in the rotor could suddenly change the ILP clearances and lengths yielding dynamic coefficient and subsequent vibration changes. ILP models are more complicated than conventional-single component-annular seal models due to their radial flow component (coriolis and centrifugal acceleration), complex geometry (axial/radial clearance coupling), internal boundary (transition) flow conditions between mechanical components along the ILP and longer length, requiring moment as well as force coefficients. Flow coupling between mechanical components results from mass and energy conservation applied at their interfaces. Typical components along the ILP include an inlet seal, curved shroud, and an exit seal, which may be a stepped labyrinth type. Von Pragenau (MSFC) has modeled labyrinth seals as a series of plain annular seals for leakage and dynamic coefficient prediction. These multi-tooth components increase the total number of 'flow coupled' components in the ILP. Childs developed an analysis for an ILP consisting of a single, constant clearance shroud with an exit seal represented by a lumped flow-loss coefficient. This same geometry was later extended to include compressible flow. The objective of the current work is to: supply ILP leakage-force impedance-dynamic coefficient modeling software to MSFC engineers, base on incompressible/compressible bulk flow theory; design the software to model a generic geometry ILP described by a series of components lying along an arbitrarily directed path; validate the software by comparison to available test data, CFD and bulk models; and develop a hybrid CFD-bulk flow model of an ILP to improve modeling accuracy within practical run time constraints.

Palazzolo, Alan B.↗

User's guide to computer programs JET 5A and CIVM-JET 5B to calculate the large elastic-plastic dynamically-induced deformations of multilayer partial and/or complete structural rings

These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included.

Wu, R. W. H.↗

Excited-State Bond Contraction and Charge Migration in a Platinum Dimer Complex Characterized by X-ray and Optical Transient Absorption Spectroscopy

Interactions between metal centers in dimeric transition metal complexes (TMCs) play important roles in their excited-state energetics and pathways, and, thus, affect their photophysical properties relevant to their applications, e.g., photoluminescent materials and photocatalysis. In this study, we report electronic and nuclear structural dynamics studies of two photoexcited pyrazolate-bridged [Pt(ppy)(µ-R 2 pz)] 2 -type Pt(II) dimers (ppy = 2-phenylpyridine, µ-R 2 pz = 3,5-substituted pyrazolate): [Pt(ppy)(µ-H 2 pz)] 2 (1) and [Pt(NDI-ppy)(µ-Ph 2 pz)] 2 (2, NDI = 1,4,5,8-naphthalenediimide), both of which have distinct ground state Pt-Pt distances. X-ray transient absorption spectroscopy at the Pt L III -edge revealed a new d-orbital vacancy due to the one-electron oxidation of the Pt centers in 1 and 2. However, while a transient Pt-Pt contraction was observed in 2, such an effect was completely absent in 1, demonstrating how the excited states of these complexes are determined by the overlap of the Pt (d z 2 ) orbitals, which is tuned by steric bulk of the pyrazolate R-groups in the 3- and 5-positions. In tandem with analysis of the Pt-Pt distance structural parameter, we observed photoinduced charge separation in 2 featuring a covalently linked NDI acceptor on the ppy ligand. The formation and subsequent decay of the NDI radical anion absorption signals were detected upon photoexcitation using optical transient absorption spectroscopy. The NDI radical anion decayed on the same timescale, hundreds of picoseconds, as that of the d-orbital vacancy signal of the oxidized Pt-Pt core observed in the XTA measurements. The data indicated ultrafast formation of the charge-separated state and subsequent charge recombination to the original Pt(II-II) species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Balance of Plant Modeling and Real-Time Hardware-in-the-Loop Integration with the Microreactor Automated Control System

The advent of novel microreactor technology has driven a focused effort to explore safety and efficiency improvements that can be achieved through the use of automated system control. Development of control strategies, especially for initial demonstration, requires an adequate surrogate environment to safely research failure modes and control integration with realistic hardware delay. However, efficiency gains from control strategies are improved when the scope of controller action is expanded to include system-level dynamics such as downstream heat extraction and mass flow. For this reason, a balance-of-plant (BOP) model of a representative microreactor system has been developed using the TRANsient Simulation Framework of Reconfigurable Models library in Modelica. This model captures a reactor and primary NaK coolant loop that represent corresponding system components of the Microreactor Applications Research Validation and EvaLuation (MARVEL) design as well as a secondary coolant loop and heat extraction representative of the Microreactor Agile Non-Nuclear Experimental Test Bed (MAGNET). This model configuration allows for hardware-in-the-loop (HIL) integration with microreactor automated control system (MACS) hardware in real time through a Python-based gRPC client. Real-time simulation of model performance with emulated hardware and communication delay suggests that under independent proportional-integral-derivative control of BOP model drum dynamics and downstream heat extraction, stable power load following is achievable. A slight delay in load following, filtering of high-frequency dynamics, and localized temperature fluctation suggest room for improvement through the development of higher-level control strategies. The simulated coupling of the MAGNET facility lays the groundwork for future digital twin analysis with a coupled MACS-MAGNET HIL demonstration.

McConnell, Jono [ORNL] (ORCID:0000000238984741)↗

Gamma-ray burster recurrence timescales

Three optical transients have been found which are associated with gamma-ray bursters (GRBs). The deduced recurrence timescale for these optical transients (tau sub opt) will depend on the minimum brightness for which a flash would be detected. A detailed analysis using all available data of tau sub opt as a function of E(gamma)/E(opt) is given. For flashes similar to those found in the Harvard archives, the best estimate of tau sub opt is 0.74 years, with a 99% confidence interval from 0.23 years to 4.7 years. It is currently unclear whether the optical transients from GRBs also give rise to gamma-ray events. One way to test this association is to measure the recurrence timescale of gamma-ray events tau sub gamma. A total of 210 gamma-ray error boxes were examined and it was found that the number of observed overlaps is not significantly different from the number expected from chance coincidence. This observation can be used to place limits on tau sub gamma for an assumed luminosity function. It was found that tau sub gamma is approx. 10 yr if bursts are monoenergetic. However, if GRBs have a power law luminosity function with a wide dynamic range, then the limit is tau sub gamma 0.5 yr. Hence, the gamma-ray data do not require tau sub gamma and tau sub opt to be different.

Schaefer, B. E.↗

Detonation wave compression in gas turbines

A study was made of the concept of augmenting the performance of low pressure ratio gas turbines by detonation wave compression of part of the flow. The concept exploits the constant volume heat release of detonation waves to increase the efficiency of the Brayton cycle. In the models studied, a fraction of the compressor output was channeled into detonation ducts where it was processed by transient transverse detonation waves. Gas dynamic studies determined the maximum cycling frequency of detonation ducts, proved that upstream propagation of pressure pulses represented no problems and determined the variations of detonation duct output with time. Mixing and wave compression were used to recombine the combustor and detonation duct flows and a concept for a spiral collector to further smooth the pressure and temperature pulses was presented as an optional component. The best performance was obtained with a single firing of the ducts so that the flow could be re-established before the next detonation was initiated. At the optimum conditions of maximum frequency of the detonation ducts, the gas turbine efficiency was found to be 45 percent while that of a corresponding pressure ratio 5 conventional gas turbine was only 26%. Comparable improvements in specific fuel consumption data were found for gas turbines operating as jet engines, turbofans, and shaft output machines. Direct use of the detonation duct output for jet propulsion proved unsatisfactory. Careful analysis of the models of the fluid flow phenomena led to the conclusion that even more elaborate calculations would not diminish the uncertainties in the analysis of the system. Feasibility of the concept to work as an engine now requires validation in an engineering laboratory experiment.

Wortman, A.↗

The Impact of Carotenoid Energy Levels on the Exciton Dynamics and Singlet–Triplet Annihilation in Isolated Bacterial Light-Harvesting 2 Complexes

The light-harvesting 2 (LH2) complex of purple phototrophic bacteria plays a critical role in absorbing solar energy and distributing the excitation energy. Exciton dynamics within LH2 complexes are controlled by the structural arrangement and energy levels of the bacteriochlorophyll (BChl) and carotenoid (Car) pigments. However, there is still debate over the competing light-harvesting versus energy-dissipation pathways. In this work, we compared five variants of the LH2 complex from genetically modified strains of Rhodobacter sphaeroides, all containing the same BChls but different Cars with increasing conjugation: zeta-carotene (N = 7; LH2 Zeta ), neurosporene (N = 9; LH2 Neu ), spheroidene (N = 10; LH2 Spher ), lycopene (N = 11; LH2 Lyco ), and spirilloxanthin (N = 13; LH2 Spir ). Absorption measurements confirmed that the Car excited-state energy decreased with increasing conjugation. Similarly, fluorescence spectra showed that the B850 BChl emission peak had an increasing red shift from LH2 Zeta →(LH2 Neu /LH2 Spher )→LH2 Lyco →LH2 Spir . In contrast, time-resolved fluorescence and ultrafast transient absorption (fs-TA) revealed similar excited-state lifetimes (∼1 ns) for all complexes except LH2 Spir (∼0.7 ns). From fs-TA analysis, an additional ∼7 ps nonradiative dissipation step from B850 BChl was observed for LH2 Zeta . Further, singlet− singlet and singlet−triplet annihilation studies showed a ∼50% average fluorescence lifetime reduction in LH2Zeta at high laser power and high repetition rate, compared to ∼10−15% reductions in LH2 Neu /LH2 Spher /LH2 Lyco and minimal lifetime change in LH2 Spir . In LH2 Zeta , the fastest decay component (<50 ps) became prominent at high repetition rates, consistent with strong singlet−triplet annihilation. Nanosecond TA measurements revealed long-lived (>40 μs) BChl triplet states in LH2 Zeta and signs of damage caused by singlet oxygen, whereas other LH2s showed faster triplet quenching (∼18 ns) by Cars. These findings highlight a key design principle of LH2 complexes: the Car triplet energy must be significantly lower than the BChl triplet energy to efficiently quench BChl triplets that otherwise act as potent “trap states,” causing exciton annihilation in laser-based experiments or photodamage in native membranes.

Absorption↗

Impact Testing and Analysis of Composites for Aircraft Engine Fan Cases

The fan case in a jet engine is a heavy structure because of its size and because of the requirement that it contain a blade released during engine operation. Composite materials offer the potential for reducing the weight of the case. Efficient design, test, and analysis methods are needed to efficiently evaluate the large number of potential composite materials and design concepts. The type of damage expected in a composite case under blade-out conditions was evaluated using a subscale test in which a glass/epoxy composite half-ring target was impacted with a wedge-shaped titanium projectile. Fiber shearing occurred near points of contact between the projectile and target. Delamination and tearing occurred on a larger scale. These damage modes were reproduced in a simpler test in which flat glass/epoxy composites were impacted with a blunt cylindrical projectile. A surface layer of ceramic eliminated fiber shear fracture but did not reduce delamination. Tests on 3D woven carbon/epoxy composites indicated that transverse reinforcement is effective in reducing delamination. A 91 cm (36 in.) diameter full-ring sub-component was proposed for larger scale testing of these and other composite concepts. Explicit, transient, finite element analyses indicated that a full-ring test is needed to simulate complete impact dynamics, but simpler tests using smaller ring sections are adequate when evaluation of initial impact damage is the primary concern.

Roberts, Gary D.↗

Frequency Agile Solar Radiotelescope

The Frequency Agile Solar Radiotelescope (FASR) has been strongly endorsed as a top community priority by both Astronomy & Astrophysics Decadal Surveys and Solar & Space Physics Decadal Surveys in the past two decades. Although it was developed to a high state of readiness in previous years (it went through a CATE analysis and was declared “doable now”), the NSF has not had the funding mechanisms in place to fund this mid-scale program. Now it does, and the community must seize this opportunity to modernize the FASR design and build the instrument in this decade. The concept and its science potential have been abundantly proven by the pathfinding Expanded Owens Valley Solar Array (EOVSA), which has demonstrated a small subset of FASR’s key capabilities such as dynamically measuring the evolving magnetic field in eruptive flares, the temporal and spatial evolution of the electron energy distribution in flares, and the extensive coupling among dynamic components (flare, flux rope, current sheet). The FASR concept, which is orders of magnitude more powerful than EOVSA, is low-risk and extremely high reward, exploiting a fundamentally new research domain in solar and space weather physics. Utilizing dynamic broadband imaging spectropolarimetry at radio wavelengths, with its unique sensitivity to coronal magnetic fields and to both thermal plasma and nonthermal electrons from large flares to extremely weak transients, the ground-based FASR will make synoptic measurements of the coronal magnetic field and map emissions from the chromosphere to the middle corona in 3D. With its high spatial, spectral, and temporal resolution, as well as its superior imaging fidelity and dynamic range, FASR is poised to provide a system-wide perspective on myriad coupled phenomena. FASR will be a highly complementary and synergistic component of solar and heliospheric observing capabilities that is critically needed to support the next generation of solar science.

Dale E. Gary↗

Integrate FARM with PID controllers: IES Simulation Ecosystem Control System Development

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM was designed to support the HERON software module in the solution of the optimal dispatch problem for IES units. As the result of HERON-FARM dispatch simulation, the set-point trajectories are optimized to meet constraints on both the production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and the process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.) at a coarse time resolution (every 10 or 100 seconds) over long time horizons (several days or weeks). In case the operational constraints need to be met at finer time resolution, the computational burden of HERON-FARM would linearly increase with the sampling rate, and sub-optimal solutions might be obtained. System responses characterized by overshoots and damped oscillations temporarily violating the imposed constraints might occur during abrupt power transients. In this report, a hierarchical control system architecture for the operation of the Dynamic Energy Technology and Integration Laboratory (DETAIL) facility constructed at INL was proposed. First, the preliminary studies on the proposed control strategy for operating the facility and the designed PI controllers were reviewed. In particular, the current approach for generating the set-point trajectories was studied, and its limits were identified. To this aim, the inclusion of a Supervisory Control layer embedding a modified version of the FARM algorithm for preserving the system safe operation over both long and real-time horizons was proposed. In this way, FARM would be applied twice, i.e., the original version (“FARM Validator”) aiding the solution of the power dispatch problem, and the modified version (“FARM Supervisory” coordinating the PI controllers to address the real-time control tasks. Despite the kernel of the two modules is the same algorithm, their roles, tasks, and capabilities are quite different. A detailed description of the role of FARM at addressing low-level control tasks is provided, along with tentative operational procedures for training the models embedded into the algorithm by using the collected experimental data.

42 ENGINEERING↗