Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chaotic systems”

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 433 records · Page 24

Weather from the Stratosphere?

Is the stratosphere, the atmospheric layer between about 10 and 50 km, important for predicting changes in weather and climate? The traditional view is that the stratosphere is a passive recipient of energy and waves from weather systems in the underlying troposphere, but recent evidence suggests otherwise. At a workshop in Whistler, British Columbia (1), scientists met to discuss how the stratosphere responds to forcing from below, initiating feedback processes that in turn alter weather patterns in the troposphere. The lowest layer of the atmosphere, the troposphere, is highly dynamic and rich in water vapor, clouds, and weather. The stratosphere above it is less dense and less turbulent (see the figure). Variability in the stratosphere is dominated by hemispheric-scale changes in airflow on time scales of a week to several months. Occasionally, however, stratospheric air flow changes dramatically within just a day or two, with large-scale jumps in temperature of 20 K or more. The troposphere influences the stratosphere mainly through atmospheric waves that propagate upward. Recent evidence shows that the stratosphere organizes this chaotic wave forcing from below to create long-lived changes in the stratospheric circulation. These stratospheric changes can feed back to affect weather and climate in the troposphere.

Baldwin, Mark P.↗

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

The X-Ray Halo Scaling Relations of Supermassive Black Holes

We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, and masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being −, followed by −. The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the −) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondi-like models have inconsistent anticorrelation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as subdominant over most of cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion, the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev–Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

Supermassive black holes↗

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗

Secular Resonances in Nonrestricted Hierarchical Triple Systems

In this paper, the averaged Hamiltonian of a nonrestricted hierarchical triple system truncated at the third order is investigated. First, each secular resonant term is studied. For the well-studied secular quadrupole theory, it is analytically reformulated in a different manner in our work. The resonance width is numerically determined and displayed on the √(1−e{sub 1}{sup 2})−√(1−e{sub 2}{sup 2}) plane (also denoted as the e-tilde {sub 1}− e-tilde {sub 2} plane). In terms of the octupole terms, we show that for a near-planar configuration of the system, considerable variations of both the eccentricities of the inner and outer orbits can be generated by a single resonant term. The resonance width for every secular resonant angle from the octupole terms is also numerically determined and displayed on the e-tilde {sub 1}− e-tilde {sub 2} plane. The results show that an orbit flip with a near-perpendicular initial mutual inclination is possible for each secular resonance. By displaying the resonance widths of different resonant terms on the same e-tilde {sub 1}− e-tilde {sub 2} plane, we intuitively show the overlap of different secular resonances. Then, the full averaged Hamiltonian with both quadrupole and octupole terms is investigated using the Poincaré surface of section, with a special focus on the orbit flip. For the cases we exploited, we find that the near-planar flip of the inner orbit can be either regular or chaotic while the outer orbit flip is generally chaotic.

79 ASTRONOMY AND ASTROPHYSICS↗

Neurodynamic oscillators

Oscillation of electrical activity has been found in many nervous systems, from invertebrates to vertebrates including man. There exists experimental evidence of very simple circuits with the capability of oscillation. Neurons with intrinsic oscillation have been found and also neural circuits where oscillation is a property of the network. These two types of oscillations coexist in many instances. It is nowadays hypothesized that behind synchronization and oscillation there is a system of coupled oscillators responsible for activities that range from locomotion and feature binding in vision to control of sleep and circadian rhythms. The huge knowledge that has been acquired on oscillators from the times of Lord Rayleigh has made the simulation of neural oscillators a very active endeavor. This has been enhanced with more recent physiological findings about small neural circuits by means of intracellular and extracellular recordings as well as imaging methods. The future of this interdisciplinary field looks very promising; some researchers are going into quantum mechanics with the idea of trying to provide a quantum description of the brain. In this work we describe some simulations using neuron models by means of which we form simple neural networks that have the capability of oscillation. We analyze the oscillatory activity with root locus method, cross-correlation histograms, and phase planes. In the more complicated neural network models there is the possibility of chaotic oscillatory activity and we study that by means of Lyapunov exponents. The companion paper shows an example of that kind.

Espinosa, Ismael↗

Designing capture trajectories to unstable periodic orbits around Europa

The hostile environment of third body perturbations restricts a mission designer's ability to find well-behaved reproducible capture trajectories when dealing with limited control authority as is typical with low-thrust missions. The approach outlined in this paper confronts this shortcoming by utilizing dynamical systems theory and an extensive preexisting database of Restricted Three Body Problem (RTBP) periodic orbits. The stable manifolds of unstable periodic orbits are utilized to attract a spacecraft towards Europa. By selecting an appropriate periodic orbit, a mission designer can control important characteristics of the captured state including stability, minimum altitudes, characteristic inclinations, and characteristic radii among others. Several free parameters are optimized in the non-trivial mapping from the RTBP to a more realistic model. Although the ephemeris capture orbit is ballistic by design, low-thrust is used to target the state that leads to the capture orbit, control the spacecraft after arriving on the unstable quasi-periodic orbit, and begin the spiral down towards the science orbit. The approach allows a mission designer to directly target fuel efficient captures at Europa in an ephemeris model. Furthermore, it provides structure and controllability to the design of capture trajectories that reside in a chaotic environment.

low-thrust↗

Experimental Reacting Hydrogen Shear Layer Data at High Subsonic Mach Number

The flow in a planar shear layer of hydrogen reacting with hot air was measured with a two-component laser Doppler velocimeter (LDV) system, a schlieren system, and OH fluorescence imaging. It was compared with a similar air-to-air case without combustion. The high-speed stream's flow speed was about 390 m/s, or Mach 0.71, and the flow speed ratio was 0.34. The results showed that a shear layer with reaction grows faster than one without; both cases are within the range of data scatter presented by the established data base. The coupling between the streamwise and the cross-stream turbulence components inside the shear layers was low, and reaction only increased it slightly. However, the shear layer shifted laterally into the lower speed fuel stream, and a more organized pattern of Reynolds stress was present in the reaction shear layer, likely as a result of the formation of a larger scale structure associated with shear layer corrugation from heat release. Dynamic pressure measurements suggest that coherent flow perturbations existed inside the shear layer and that this flow became more chaotic as the flow advected downstream. Velocity and thermal variable values are listed in this report for a computational fluid dynamics (CFD) benchmark.

Chang, C. T.↗

Anomalous losses of energetic particles in the presence of an oscillating radial electric field in fusion plasmas

The confinement of energetic particles in nuclear fusion devices is studied in the presence of an oscillating radial electric field and an axisymmetric magnetic equilibrium. Here, it is shown that, despite the poloidal and toroidal symmetries, initially integrable orbits turn into chaotic regions that can potentially intercept the wall of the tokamak, leading to particle losses. It is observed that the losses exhibit algebraic time decay different from the expected exponential decay characteristic of radial diffusive transport. A dynamical explanation of this behaviour is presented, within the continuous time random walk theory. The central point of the analysis is based on the fact that, contrary to the radial displacement, the poloidal angle is not bounded and a proper statistical analysis can therefore be made, showing for the first time that energetic particle transport can be super-diffusive in the poloidal direction and characterised by asymmetric poloidal displacement. The connection between poloidal and radial positions ensured by the conservation of the toroidal canonical momentum, implies that energetic particles spend statistically more time in the inner region of the tokamak than in the outer one, which explains the observed algebraic decay. This indicates that energetic particles might be efficiently slowed down by the thermal population before leaving the system. Also, the asymmetric transport reveals a new possible mechanism of self-generation of momentum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Poor confinement in stellarators at high energy

High-energy particle resonances can modify particle distributions and even cause significant particle loss. Resonances can be present in any toroidal confinement device and can easily be found numerically. Many stellarators have weak magnetic shear so that large islands and large chaotic regions can be produced by resonant perturbations with small amplitudes. Additionally, while the choice of the field line helicity profile in the plasma can limit the presence of resonances at low particle energy, the resonance location is energy-dependent, and they can move into the plasma at higher energy. If resonances match the toroidal variation of the equilibrium, they can produce wide islands in the phase space of orbits even in the absence of perturbations due to instabilities. These islands increase in size with particle energy and can seriously affect the confinement of high-energy ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating Sulfur-Rich Mercury Analogs Exposed to Simulated Micrometeoroid Bombardment in the Laboratory

Space weathering (SW) continually alters the spectral, microstructural, and chemical characteristics of the surface of airless bodies across the solar system. The effects of SW vary depending on the heliocentric distance and the initial composition of the target surface. While SW on the Moon and S-type asteroids is well documented, our understanding of how this process affects Mercury is at an early stage. Mercury’s interplanetary environment is harsh, with the surface of the planet experiencing an intense solar wind flux as well as a higher flux and velocity of micrometeoroid impactors compared to the Moon and S-type asteroids. In addition, Mercury is also a geochemical endmember, with a surface composition low in Fe (<2 wt.%) and enriched in volatile components, such as sulfur (up to 4 wt.% in the low reflectance material (LRM)). These volatile components are thought to play a major role in the formation of hollows via their sublimation. Sulfur has been hypothesized to occur at the surface of Mercury as sulfide minerals (MgS, CaS) based on its correlation with Mg and Ca in remote sensing data. However, recent observations of chaotic terrains in the north polar area of Mercury and of glacier-like features at lower latitudes indicated that octasulfur (S 8 , elemental sulfur) is another likely constituent of a volatile-rich layer in Mercury’s crust. Its behavior on Mercury may result in a complex cycle of enrichment and depletion. While S is often depleted on small body surfaces, Mercury’s gravity could result in ejected S subsequently returning to the surface and coating regolith grains. Further, the reaction of reduced S-rich gas with glasses of a Mercury-like composition also produced S-rich coatings. However, the precise behavior and evolution of S-rich species exposed to the harsh SW on Mercury remains poorly understood and needs to be further investigated in the laboratory. Here, we present the results of our analyses of the spectral, microstructural, and chemical characteristics of S-rich Mercury analogs irradiated by pulsed laser to simulate the short duration, high temperature events associated with micrometeoroid impacts.

N Bott↗

2024 Nasa Lunabotics University Competition: Site Preparation With Bulk Regolith

Introduction: Lunabotics provides accredited institutions of higher learning students (vocational-technical, college, university) an opportunity to apply the NASA systems engineering process to design and build a prototype robot. This robot would be capable of performing the proposed operations on the Lunar surface in support of future Artemis mission goals. Lunabotics features a systems engineering design challenge to engage students in the next phase of hu-man space exploration supporting the Artemis missions. This two-semester event encourages students to design and build an autonomous or telerobotic robot designed to traverse the simulated Lunar surface and complete the assigned construction tasks. The number of teams accepted into this challenge is not predetermined but is based on the scores and overall quality of the Project Management Plans received and other factors. The culmination of the Lunabotics virtual challenge will be the design, build and operation of a functional prototype Lunar robot. Teams are required to submit the following: (1) Project Management Plan, (2) Systems Engineering Paper, a (3) STEM Engage-ment Report, and a (4) Proof of Life Video. This is an optional item, but to qualify for the grand prize a team must also submit a: (5) Presentation and Demonstration. Background: The NASA Lunabotics University Competition was first held in 2010 as a follow on to the NASA Regolith Excavation Challenge [1]. The high level of interest and participation from over 50 universities each year has led to a sustained annual competi-ion cadence: it has been held every May for the past 14 years [2,3]. Over 6,000 students have participated and been inspired to pursue Science, Technology, En-gineering and Mathematics careers (STEM). 2024 Competition: The necessary lunar surface tasks are evolving to meet the NASA Artemis Mission requirements. In the past Lunabotics challenges we gathered data to support Lunar mining for consuma-bles in the Lunar regolith. Now, in 2024, the task is to gather data on Lunar site preparation and construction by designing and building a robot that will traverse the chaotic Lunar terrain and construct a regolith-based berm. The goal is to build a berm structure which would be useful to the Artemis Mission for blast and ejecta protection during lunar landings and launches, shading cryogenic propellant tank farms, providing radiation protection around a nuclear power plant and other mission critical uses. Lunabotics will consist of three separate events this year. The first event is NASA’s Lunabotics Project Development Challenge, where teams submit various deliverables to be scored by judges. The second event will be the University of Central Florida (UCF) Lunabotics Qualification challenge, in the Exolith laboratory, where teams will put their de-signs to the test. The top ten scoring teams from the Qualification challenge are then invited to the third and final event, NASA’s Lunabotics On-Site Challenge at Kennedy Space Center in Florida. This presentation will summarize the results and lessons learned from the NASA Lunabotics University Competition held in May 2024.

Lunabotics↗

Replica symmetry breaking for the integrable two-site Sachdev–Ye–Kitaev model

We analyze a two-body non-Hermitian two-site Sachdev–Ye–Kitaev (SYK) model with the couplings of one site complex conjugated to the other site. This model, with no explicit coupling between the sites, shows an infinite number of second-order phase transitions, which is a consequence of the factorization of the partition function into a product over Matsubara frequencies. We calculate the quenched free energy in two different ways: first in terms of the single-particle energies and second by solving the Schwinger–Dyson equations of the two-site model. The first calculation can be done entirely in terms of a one-site model. The conjugate replica enters due to non-analyticities when Matsubara frequencies enter the spectral support of the coupling matrix. The second calculation is based on the replica trick of the two-site partition function. Both methods give the same result. The free-fermion partition function can also be rephrased as a matrix model for the coupling matrix. Up to minor details, this model is the random matrix model that describes the chiral phase transition of QCD, and the order parameter of the two-body model corresponds to the chiral condensate of QCD. Comparing to the corresponding four-body model, we are able to determine which features of the free energy are due to the chaotic nature of the four-body model. The high-temperature phase of both models is entropy dominated, and in both cases, the free energy is determined by the spectral density. The chaotic four-body SYK model has a low-temperature phase whose free energy is almost temperature-independent, signaling an effective gap of the theory even though the actual spectrum does not exhibit a gap. On the other hand, the low-temperature free energy of the two-body SYK model is not flat; in fact, it oscillates to arbitrarily low temperature. This indicates a less desirable feature that the entropy of the two-body model is not always positive in the low-temperature phase, which most likely is a consequence of the non-hermiticity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Basic studies of baroclinic flows

A fully nonlinear 3-dimensional numerical model (GEOSIM), previously developed and validated for several cases of geophysical fluid flow, has been used to investigate the dynamical behavior of laboratory experiments of fluid flows similar to those of the Earth's atmosphere. The phenomena investigated are amplitude vacillation, and the response of the fluid system to uneven heating and cooling. The previous year's work included hysteresis in the transition between axisymmetric and wave flow. Investigation is also continuing of the flows in the Geophysical Fluid Flow Cell (GFFC), a low-gravity Spacelab experiment. Much of the effort in the past year has been spent in validation of the model under a wide range of external parameters including nonlinear flow regimes. With the implementation of a 3-dimensional upwind differencing scheme, higher spectral resolution, and a shorter time step, the model has been found capable of predicting the majority of flow regimes observed in one complete series of baroclinic annulus experiments of Pfeffer and co-workers. Detailed analysis of amplitude vacillation has revealed that the phase splitting described in the laboratory experiments occurs in some but not all cases. Through the use of animation of the models output, a vivid 3-dimensional view of the phase splitting was shown to the audience of the Southeastern Geophysical Fluid Dynamics Conference in March of this year. A study on interannual variability was made using GEOSIM with periodic variations in the thermal forcing. Thus far, the model has not predicted a chaotic behavior as observed in the experiments, although there is a sensitivity in the wavenumber selection to the initial conditions. Work on this subject, and on annulus experiments with non-axisymmetric thermal heating, will continue. The comparison of GEOSIM's predictions will result from the Spacelab 3 GFFC experiments continued over the past year, on a 'back-burner' basis. At this point, the study (in the form of a draft of a journal article) is nearly completed. The results from GEOSIM compared very well with the experiments, and the use of the model allows the demonstration of flow mechanics that were not possible with the experimental data. For example, animation of the model output shows that the forking of the spiral bands is a transient phenomenon, due to the differential east-west propagation of convection bands from different latitudes.

Miller, Timothy L.↗

Lunar and Planetary Science XXXVI, Part 18

Topics discussed include: PoDS: A Powder Delivery System for Mars In-Situ Organic, Mineralogic and Isotopic Analysis Instruments Planetary Differentiation of Accreting Planetesimals with 26Al and 60Fe as the Heat Sources Ground-based Observation of Lunar Surface by Lunar VIS/NIR Spectral Imager Mt. Oikeyama Structure: First Impact Structure in Japan? Central Mounds in Martian Impact Craters: Assessment as Possible Perennial Permafrost Mounds (Pingos) A Further Analysis of Potential Photosynthetic Life on Mars New Insight into Valleys-Ocean Boundary on Mars Using 128 Pixels per Degree MOLA Data: Implication for Martian Ocean and Global Climate Change; Recursive Topography Based Surface Age Computations for Mars: New Insight into Surficial Processes That Influenced Craters Distribution as a Step Toward the Formal Proof of Martian Ocean Recession, Timing and Probability; Laser-induced Breakdown Spectroscopy: A New Method for Stand-Off Quantitative Analysis of Samples on Mars; Milk Spring Channels Provide Further Evidence of Oceanic, >1.7-km-Deep Late Devonian Alamo Crater, Southern Nevada; Exploration of Martian Polar Residual Caps from HEND/ODYSSEY Data; Outflow Channels Influencing Martian Climate: Global Circulation Model Simulations with Emplaced Water; Presence of Nonmethane Hydrocarbons on Pluto; Difference in Degree of Space Weathering on the Newborn Asteroid Karin; Circular Collapsed Features Related to the Chaotic Terrain Formation on Mars; A Search for Live (sup 244)Pu in Deep-Sea Sediments: Preliminary Results of Method Development; Some Peculiarities of Quartz, Biotite and Garnet Transformation in Conditions of Step-like Shock Compression of Crystal Slate; Error Analysis of Remotely-Acquired Mossbauer Spectra; Cloud Activity on Titan During the Cassini Mission; Solar Radiation Pressure and Transient Flows on Asteroid Surfaces; Landing Site Characteristics for Europa 1: Topography; and The Crop Circles of Europa.

Source record↗

Next-Cycle Optimal Dilute Combustion Control via Online Learning of Cycle-to-Cycle Variability Using Kernel Density Estimators

Dilute combustion using exhaust gas recirculation (EGR) presents a cost-effective method for increasing the efficiency of spark-ignition (SI) engines. However, the maximum amount of EGR that can be used at a given condition is limited by a rapid increment of cycle-to-cycle variability (CCV). This study describes a methodology to design a model-based stochastic optimal controller to adjust the cycle-to-cycle fuel injection quantity in order to reduce CCV and further extend the dilute limit. Given the complexity and chaotic nature of combustion events, the controller was enhanced with online learning in order to identify the statistical properties of combustion efficiency, which are needed to generate predictions for next-cycle events. This study showed that a kernel density estimator (KDE) can be used to learn the combustion properties in real time and can be incorporated into the feedback policy in order to calculate the optimal control command. Experimental results suggested that the dilute limit can be extended from 18.5% to 21% EGR fraction at an operating condition relevant for highway cruising. Additionally, the proposed controller can achieve a large CCV reduction with less fuel enrichment compared to previous methods, overall contributing to an increase in 0.2% indicated fuel conversion efficiency.

33 ADVANCED PROPULSION SYSTEMS↗

Visualization of Multi-Fidelity Approximations of Stochastic Economic Dispatch

As renewable energy generation deployment increases, the operation of electrical grids becomes more complex. Economic dispatch is part of a grid operator's regular decision process where the amount of energy to generate is determined based on the number of available generators and the actual level of energy demand. Renewable generators are inherently stochastic due to the chaotic nature of weather patterns, and thus, real-time decisions of economic dispatch become increasingly complex. Modeling efforts to assist in these decisions in the highest fidelity typically take hours to days to solve on leadership-class computers; too long for the 5-minute operational time-frame demanded of operators. Alternatively, multi-fidelity approximations can be used to predict generation levels quickly and with sufficient accuracy to be used for real-time operations. We have developed a visualization tool to demonstrate the utility of multi-fidelity approximations by displaying contextual results of economic dispatch approximations, comparisons across fidelity levels of generation levels and possible failures to meet demand, and meta-data on the modeling setup.

interactive visualization↗

Insights on the continuous representations of piecewise-smooth nonlinear systems: limits of applicability and effectiveness

Dynamical systems subject to intermittent contact are often modeled with piecewise-smooth contact forces. However, the discontinuous nature of the contact can cause inaccuracies in numerical results or failure in numerical solvers. Representing the piecewise contact force with a continuous and smooth function can mitigate these problems, but not all continuous representations may be appropriate for this use. Here, five representations used by previous researchers (polynomial, rational polynomial, hyperbolic tangent, arctangent, and logarithm-arctangent functions) are studied to determine which ones most accurately capture nonlinear behaviors including super- and subharmonic resonances, multiple solutions, and chaos. The test case is a single-DOF forced Duffing oscillator with freeplay nonlinearity, solved using direct time integration. This work intends to expand on past studies by determining the limits of applicability for each representation and what numerical problems may occur.

42 ENGINEERING↗