Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power system oscillation”

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.

787 records · Page 44

Preliminary Planar Formation: Flight Dynamics Near Sun-Earth L2 Point

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of the Sun-Earth L2 libration point. Some of these projects will involve a distributed space system of telescope spacecraft acting together as a single telescope for high-resolution. The individual telescopes will be configured in a plane, surrounding a hub, where the telescope plane can be aimed toward various astronomical targets of interest. In preparation for these missions, it is necessary to develop an improved understanding of the dynamical behavior of objects in a planar configuration near L2. The classical circular restricted three body problem is taken as the basis for the analysis. At first order, the motion of such a telescope relative to the hub is described by a system of linear second order differential equations. These equations are identical to the circular restricted problem's linear equations describing the hub motion about L2. Therefore, the fundamental frequencies, both parallel to and normal to the ecliptic plane, are the same for the relative telescope motion as for the hub motion. To maintain the telescope plane for the duration necessary for the planned observations, a halo-type orbit of the telescopes about the hub is investigated. By using a halo orbit, the individual telescopes remain in approximately the same plane over the observation duration. For such an orbit, the fundamental periods parallel to and normal to the ecliptic plane are forced to be the same by careful selection of the initial conditions in order to adjust the higher order forces. The relative amplitudes of the resulting oscillations are associated with the orientation of the telescope plane relative to the ecliptic. As in the circular restricted problem, initial conditions for the linearized equations must be selected so as not to excite the convergent or divergent linear modes. In a higher order analysis, the telescope relative motion equations include the effects of the position of the hub relative to L2. In this paper, the differential equations are developed through second order in the distance of the hub from the libration point. A modified Lindstedt-Poincad perturbation method is employed to construct the solution of these differential equations through that same order of magnitude. In the course of the solution process, relationships are determined between the initial conditions of the telescopes, selected in order to avoid resonance excitation. As the differential equations include the hub position, it is necessary to simultaneously develop the solution for the hub. As has been done in past analyses of the circular restricted problem, the hub position is written in a power series formulation in terms of its distance from L2. Then, in order to be included in the telescope equations, the hub solution is cast in terms of the nonlinear frequency of the relative telescope motion. In the course of the analysis, it is determined that the hub should also maintain a halo orbit - about L2. Additionally, relationships are formed between the initial conditions of the telescopes and the hub. These relationships may be used to associate sets of initial conditions with particular orientations of the telescope plane. The accuracy of the analytical solution is verified through various simulations and comparison to numerical integration of the differential equations. The results of the simulations are presented, along with a graphical representation of the relationships between the initial conditions of the telescopes and hub.

Segerman, Alan M.↗

Characterization of a 100-kW Class Nested-Channel Hall Thruster

The next generation of electric propulsion consists of systems in excess of 300 kW of power. These systems enable a wide variety of missions, including crewed missions to near-Earth asteroids and Mars. Hall thrusters are a particularly attractive technology for these missions, but development and demonstration of 100-kW class devices has been limited to date. Here, a 100-kW class three-channel nested Hall thruster called the X3 was operated up to 102 kW total discharge power. The three channels of the X3 can be operated in any combination, providing seven unique configurations and a total throttling envelope of 2-200 kW. Previous testing of the X3 was limited to 30 kW and found that it was not providing state-of-the-art performance.Two low-power test campaigns were completed at the University of Michigan which identified potential mechanisms for this under-performance. Improvements to the thruster were made before a high-power performance characterization at NASA Glenn Research Center. There, the X3 was fired on xenon propellant from 5-102 kW total power. The thruster demonstrated stable operation in all seven channel combinations at discharge voltages from 300 V to 500 V and three different current densities. All seven channel combinations demonstrated similar performance at a given discharge voltage and current density. The largest thrust recorded was 5.4 N, and total efficiency and specific impulse ranged from 0.54 to 0.67 and 1800 seconds to 2650 seconds, respectively.In addition to the performance measurements, a suite of plasma diagnostics and a high-speed camera were used to study the operation of the thruster in greater detail. The probe results are compared against those in the literature and demonstrate that the X3, even in multi-channel operation, is producing similar charge, mass, current, and voltage utilization efficiencies as state of the art single-channel Hall thrusters. High-speed camera analysis identified that the X3 operated with similar discharge current oscillations at all conditions tested. Additionally, a preliminary calculation of cross-channel ingestion and its effect on thruster efficiency is made. This work represents the highest total power (102 kW), thrust (5.4 N), and discharge current (247 A) demonstrated by a Hall thruster to date, improvements of 6 percent, 64 percent, and 119 percent respectively over previous values. These results are discussed in the context of continued high-power Hall thruster development and future mission applications.

Hall, Scott James↗

Web-based wide-area monitoring platform for ringdown and clustering analytics in power systems

This paper introduces an open-source research platform for monitoring the Mexican interconnected power grid, allowing real-time processing and information extraction of the grid’s dynamic condition. Moreover, the platform is a Python-based development that embeds different ringdown and clustering analytics tools. In the case of ringdown analysis, the modal information can be extracted using some of the most known algorithms, i.e., Prony analysis, eigensystem realization algorithm (ERA), and matrix pencil (MP). For clustering analysis, the coherent behaviour of generator and non-generator buses is provided by applying recent state-of-the-art techniques such as affinity propagation, K-means, hierarchical agglomerative clustering, and typicality data analysis. The results of up to 93 PMUs show that this open-source platform suits researchers’ and engineers’ power system dynamic analysis requirements.

Clustering↗

Picture this SELFI: Submillimeter Enceladus Life Fundamentals Instrument

The Submillimeter Enceladus Life Fundamentals Instrument (SELFI) is a passive remote sensing submillimeter heterodyne spectrometer being developed at NASA GSFC under NASA's Maturation of Instruments for Solar System Exploration (MatISSE) program. SELFI will advance submillimeter receiver technology by 1) investigating the chemical and isotopic compositions and corresponding densities of Enceladus' plume material, their vertical thermal structures, and the transport mechanisms within the plumes, and 2) characterizing the source regions from which the plumes emerge. The Enceladus plumes are important in the context of life and habitability of its subsurface ocean environment. SELFI remote sensing measurements will 1) measure the spatial and temporal variabilities in the plume chemical compositions, 2) provide insight in to Enceladus' subsurface ocean environment by monitoring H2O, HDO, d18O, and d17O, 3) constrain the oxidation state of the subsurface ocean using H2O2 and O3, and 4) utilize the SO2 and H2S spectral signatures to constrain the impact arising from both the sea-floor volcanoes and pre-biotic molecules. Moreover, the detection of the remaining molecular species (14 in total) is vital to improving the current state of knowledge of Enceladus' subsurface ocean habitability this also permits us to explore the chemical alteration processes arising from primordial volatiles that have been observed in comets. Lastly, SELFI's continuum observations enable the correlation between observed variations in plume activity with surface temperatures.SELFI is currently being developed under a technology maturation program that will advance the RF-to-digital electronics of a submillimeter-wave heterodyne spectrometer to simultaneously observe fourteen molecular species with resonances between 530 GHz and 600 GHz. SELFI will have fine radiometric resolution, high spectral resolution (resolving power R > 106), multiple continuum channels and a high dynamical range, necessary to map Enceladus across its 30 K to 250 K temperature range. Under the MatISSE program, SELFI will advance from TRL 4 to 6 four key technologies of the RF-to-digital subsystem, which are: 1) the RF low noise amplifier design; 2) the single-sideband (SSB) mixer and local oscillator; 3) the IF assembly down-converter that maps the fourteen species to 2 x 500 MHz bandwidth; and 4) the digital spectrometer electronics.

Racette, Paul↗

Investigation of Theoretical Solutions to a Bottom-Raised Oscillating Surge Wave Energy Converter (OSWEC) Through Experimental and Parametric Studies

Experiments were conducted on a wave tank model of a bottom raised oscillating surge wave energy converter (OSWEC) model in regular waves. The OSWEC model shape was a thin rectangular flap, which was allowed to pitch in response to incident waves about a hinge located at the intersection of the flap and the top of the supporting foundation. Torsion springs were added to the hinge in order to position the pitch natural frequency at the center of the wave frequency range of the wave maker. The flap motion as well as the loads at the base of the foundation were measured. The OSWEC was modeled analytically using elliptic functions in order to obtain closed form expressions for added mass and radiation damping coefficients, along with the excitation force and torque. These formulations were derived and reported in a previous publication by the authors. While analytical predictions of the foundation loads agree very well with experiments, large discrepancies are seen in the pitch response close to resonance. These differences are analyzed by conducting a sensitivity study, in which system parameters, including damping and added mass values, are varied. The likely contributors to the differences between predictions and experiments are attributed to tank reflections, standing waves that can occur in long, narrow wave tanks, as well as the thin plate assumption employed in the analytical approach.

analytical↗

Satellite Altimetry, Ocean Circulation, and Data Assimilation

Ocean circulation is a critical factor in determining the Earth's climate. Satellite altimetry has been proven a powerful technique for measuring the height of the sea surface for the study of global ocean circulation dynamics. A major objective of my research is to investigate the utility of altimeter data for ocean circulation studies. The 6 years' data record of TOPEX/POSEIDON have been analyzed to study the spatial and temporal characteristics of large-scale ocean variability. A major result obtained in 1998 is the discovery of large-scale oscillations in sea level with a period of 25 days in the Argentine Basin of the South Atlantic Ocean (see diagram). They exhibit a dipole pattern with counterclockwise rotational propagation around the Zapiola Rise (centered at 45S and 317E), a small seamount in the abyssal plain of the basin. The peak-to-trough amplitude is about 10 cm over a distance of 500-1000 km. The amplitude of these oscillations has large seasonal-to-interannual variations. The period and rotational characteristics of these oscillations are remarkably similar to the observations made by two current meters deployed near the ocean bottom in the region. What TOPEX/POSEIDON has detected apparently are manifestations of the movement of the entire water column (barotropic motion). The resultant transport variation is estimated to be about 50 x 10(exp 6) cubic M/S, which is about 50% of the total water transport in the region. Preliminary calculations suggest that these oscillations are topographically trapped waves. A numerical model of the South Atlantic is used to investigate the nature of and causes for these waves. A very important property of sea surface height is that it is directly related to the surface geostrophic velocity, which is related to deep ocean circulation through the density field. Therefore altimetry observations are not only useful for determining the surface circulation but also for revealing information about the deep ocean. Another objective of my research is to investigate the extent of information one can infer about the deep ocean from altimetric observations. The approach is data assimilation by ocean models. The objective is to develop a modeling/data assimilation system that will produce estimates of the 3-dimensional state of the ocean for the entire duration of the TOPEX/POSEIDON Mission.

Fu, Lee-Lueng↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of Pulse Counting for the Mars Organic Mass Analyzer (MOMA) Ion Trap Detection Scheme

The Mars Organic Mass Analyzer is being developed at Goddard Space Flight Center to identify organics and possible biological compounds on Mars. In the process of characterizing mass spectrometer size, weight, and power consumption, the use of pulse counting was considered for ion detection. Pulse counting has advantages over analog-mode amplification of the electron multiplier signal. Some advantages are reduced size of electronic components, low power consumption, ability to remotely characterize detector performance, and avoidance of analog circuit noise. The use of pulse counting as a detection method with ion trap instruments is relatively rare. However, with the recent development of high performance electrical components, this detection method is quite suitable and can demonstrate significant advantages over analog methods. Methods A prototype quadrupole ion trap mass spectrometer with an internal electron ionization source was used as a test setup to develop and evaluate the pulse-counting method. The anode signal from the electron multiplier was preamplified. The an1plified signal was fed into a fast comparator for pulse-level discrimination. The output of the comparator was fed directly into a Xilinx FPGA development board. Verilog HDL software was written to bin the counts at user-selectable intervals. This system was able to count pulses at rates in the GHz range. The stored ion count nun1ber per bin was transferred to custom ion trap control software. Pulse-counting mass spectra were compared with mass spectra obtained using the standard analog-mode ion detection. Prelin1inary Data Preliminary mass spectra have been obtained for both analog mode and pulse-counting mode under several sets of instrument operating conditions. Comparison of the spectra revealed better peak shapes for pulse-counting mode. Noise levels are as good as, or better than, analog-mode detection noise levels. To artificially force ion pile-up conditions, the ion trap was overfilled and ions were ejected at very high scan rates. Pile-up of ions was not significant for the ion trap under investigation even though the ions are ejected in so-called 'ion-micro packets'. It was found that pulse counting mode had higher dynamic range than analog mode, and that the first amplification stage in analog mode can distort mass peaks. The inherent speed of the pulse counting method also proved to be beneficial to ion trap operation and ion ejection characterization. Very high scan rates were possible with pulse counting since the digital circuitry response time is so much smaller than with the analog method. Careful investigation of the pulse-counting data also allowed observation of the applied resonant ejection frequency during mass analysis. Ejection of ion micro packets could be clearly observed in the binned data. A second oscillation frequency, much lower than the secular frequency, was also observed. Such an effect was earlier attributed to the oscillation of the total plasma cloud in the ion trap. While the components used to implement pulse counting are quite advanced, due to their prevalence in consumer electronics, the cost of this detection system is no more than that of an analog mode system. Total pulse-counting detection system electronics cost is under $250

Van Amerom, Friso H.↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Raman-shifting an ArF excimer laser to generate new lines for obtaining optical diagnostic based information in flow fields

The application of tunable excimer lasers in combustion and flow diagnostics is almost routine nowadays. The properties of this laser system that enable density and temperature measurements in supersonic and hypersonic flow fields to be conducted are its high power, high repetition rate, and high spectral brightness. The limitation imposed by this system on these measurements is the paucity of lines in the wavelength region, the vacuum-ultraviolet, where species of interest, such as OH, N2, O2, H2, H2O, CO, NO, etc., are susceptible to electronic excitation to high-lying states. To circumvent this problem one normally resorts to nonlinear optical techniques such as frequency conversion via stimulated Raman scattering (SRS), more commonly known as Raman shifting or Raman mixing, to extend these nonintrusive and nonperturbing techniques to the shorter wavelengths in the VUV region and, for that matter, to longer wavelengths in the infrared region, if the need arises. The theoretical basis of SRS and its application are well documented in the literature. In essence, the Raman shift is a consequence of the inelastic scattering of the incident radiation by the sample. Most of the scattered radiation from the molecules of the sample is unchanged in frequency. However, a small fraction of the incident radiation is changed in frequency. This shift is a result of the fact that some of the incident photons on colliding with the molecules of the sample give up some of their energy and emerge with a lower energy resulting in the lower-frequency Stokes radiation. Other incident photons may increase their energy by colliding with the vibrationally excited molecules of the medium and emerge as higher-frequency antistokes radiation. The generation of the latter is the main objective of this project. The process, however, depends on several factors, including the beam quality of the pump laser, the cross-section of the gaseous medium, the gas pressure, and the ambient temperature of the gas near the focal region. Furthermore, since the Raman shifting process is polarization sensitive, it is necessary to have all of the laser energy in a single polarization. These factors were taken into consideration in the execution of the project. The implementation of the Raman shift was accomplished by focusing the 193 nm output of an ArF excimer laser (Lamda-Physik LPX 150) into a 1-meter long high pressure recirculating Raman cell filled with H2 gas. The laser system was modified in order to improve the mode quality of the pump beam to enhance the Raman shifting. To accomplish this feat, a prism beam expander and grating on the oscillator discharge provided wavelength tuning over the excimer gain profile. Furthermore, a triple-pass configuration, as opposed to unstable resonator optics, was employed in the operation of the amplifier cavity so that when the oscillator output radiation, focused by a 51-cm focal length fused silica lens through a 50 micron pinhole (serving as a spatial filter) and recollimated with a 25-cm focal length lens, was fed into the amplifier, it was injection locked, thereby providing tunable radiation with relatively low divergence. The forward scattered radiation emanating from the impingement of the modified pump beam on the Raman cell was detected using an energy meter after the latter had been separated from it using a dispersing prism.

Koker, Edmond B.↗

Mooring System Cost Estimates for Wave Energy Farms in Shared Mooring Arrays

As wave energy converters (WECs) become more advanced and cost-efficient, so too must their mooring systems. A key question in the development plans for WECs is the cost of the mooring system, particularly for large wave farms. WEC devices deployed in a WEC farm array, where each device can be connected by shared mooring lines in various array layouts, have potential to reduce mooring system costs significantly. This paper presents the modeling and designing of mooring systems for large WEC arrays and the calculation of the cost of each mooring system to determine the change in cost as the number of WECs in a farm increases. A baseline mooring system for a single floating oscillating water column (OWC) WEC was developed for this analysis. The mooring system utilizes four anchored mooring lines connected to a square assembly of wire rope mooring lines, supported by four floating buoys, and attached to the floating WEC by four polyester rope mooring lines. This assembly, referred to as a floating cell, can be tiled to form various rectangular WEC arrays. The objective of this analysis is to determine how the mooring system cost changes as more WECs are added to an array layout, each with their own interconnected floating cell. To do this, complete mooring systems need to be designed for each WEC array layout. To narrow down the design space of a WEC array mooring system, a couple assumptions were made. It was assumed that the floating cell parameters of the baseline design were to stay constant across all floating cells in the WEC array. It was also assumed that the anchored mooring lines would be of the baseline configuration, a predominantly chain mooring line with a short section of polyester rope near the fairlead, and a drag-embedment anchor. These anchored lines were assumed to extend from the outer edges of the WEC array, inline with the headings of the wire rope mooring lines of the floating cells, or diagonal if extending from a corner of the WEC array. Full mooring systems were designed for 2xN and NxN WEC array layouts and efficiently simulated in the mooring dynamics simulation tool, MoorDyn, to ensure all dynamic constraints were met. The system costs were calculated and then refined by shortening unnecessary chain line lengths and reducing the chain diameters of the downstream anchored mooring lines. It was found that, in general, mooring system costs per WEC decrease when WECs are installed in an array. Compared to the baseline mooring system for a single WEC, the 2x3 mooring system array had the lowest mooring system cost per WEC, reducing the cost per WEC by 59%. The 3x3 and 4x4 array mooring systems also saw significant reductions in cost per WEC but had negligible cost savings between the two designs, primarily because the larger 4x4 mooring system requires larger chain diameters, which increases cost. These results provide an interesting glimpse into modeling, designing, and calculating the cost of mooring systems for large WEC arrays.

cost↗

Real-time steerable frequency-stepped Doppler backscattering (DBS) system for local helicon wave electric field measurements on the DIII-D tokamak

A new frequency-stepped Doppler backscattering (DBS) system has been integrated into a real-time steerable electron cyclotron heating launcher system to simultaneously probe local background turbulence (f < 10 MHz) and high-frequency (20–550 MHz) density fluctuations in the DIII-D tokamak. The launcher allows for 2D steering (horizontally and vertically) over wide angular ranges to optimize probe location and wavenumber response. The vertical steering can be optimized during a discharge in real time. The new DBS system employs a programmable frequency synthesizer with adjustable dwell time as a source to launch either O or X-mode polarized millimeter waves. This system can step in real-time over the entire E-band frequency range (60–90 GHz). This combination of capabilities allows for the diagnosis of the complex internal spatial structure of high power (>200 kW) helicon waves (476 MHz) injected from an external antenna during helicon current drive experiments in DIII-D. Broadband density fluctuations around the helicon frequency are observed during real-time scans of measurement location and wavenumber during these experiments. Analysis indicates that these broadband high-frequency fluctuations are a result of backscattering of the DBS millimeter-wave probe beam from plasma turbulence modulated by the helicon wave. It is observed that background turbulence is effectively locally “tagged” with the helicon wave electric field, forming images of the turbulent spectrum in the overall density fluctuation spectrum that appear as high-frequency sidebands of the turbulence. These observations of background turbulence and high-frequency fluctuations open up the possibility of monitoring local helicon wave amplitude by comparing the high-frequency signal amplitude to the simultaneously measured background turbulence. In combination with the real-time measurement location and wavenumber scanning capabilities (offered by real-time frequency-stepping and steering), this allows rapid determination of the spatial distribution of the helicon wave power during steady-state plasma operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Physics-Based Methods of Failure Analysis and Diagnostics in Human Space Flight

The Integrated Health Management (IHM) for the future aerospace systems requires to interface models of multiple subsystems in an efficient and accurate information environment at the earlier stages of system design. The complexity of modern aeronautic and aircraft systems (including e.g. the power distribution, flight control, solid and liquid motors) dictates employment of hybrid models and high-level reasoners for analysing mixed continuous and discrete information flow involving multiple modes of operation in uncertain environments, unknown state variables, heterogeneous software and hardware components. To provide the information link between key design/performance parameters and high-level reasoners we rely on development of multi-physics performance models, distributed sensors networks, and fault diagnostic and prognostic (FD&P) technologies in close collaboration with system designers. The main challenges of our research are related to the in-flight assessment of the structural stability, engine performance, and trajectory control. The main goal is to develop an intelligent IHM that not only enhances components and system reliability, but also provides a post-flight feedback helping to optimize design of the next generation of aerospace systems. Our efforts are concentrated on several directions of the research. One of the key components of our strategy is an innovative approach to the diagnostics/prognostics based on the real time dynamical inference (DI) technologies extended to encompass hybrid systems with hidden state trajectories. The major investments are into the multiphysics performance modelling that provides an access of the FD&P technologies to the main performance parameters of e.g. solid and liquid rocket motors and composite materials of the nozzle and case. Some of the recent results of our research are discussed in this chapter. We begin by introducing the problem of dynamical inference of stochastic nonlinear models and reviewing earlier results. Next, we present our analytical approach to the solution of this problem based on the path integral formulation. The resulting algorithm does not require an extensive global search for the model parameters, provides optimal compensation for the effects of dynamical noise, and is robust for a broad range of dynamical models. In the following Section the strengths of the algorithm are illustrated illustrated by inferring the parameters of the stochastic Lorenz system and comparing the results with those of earlier research. Next, we discuss a number of recent results in application to the development of the IHM for aerospace system. Firstly, we apply dynamical inference approach to a solution of classical three tank problems with mixed unknown continuous and binary parameters. The problem is considered in the context of ground support system for filling fuel tanks of liquid rocket motors. It is shown that the DI algorithm is well suited for successful solution of a hybrid version of this benchmark problem even in the presence of additional periodic and stochastic perturbation of unknown strength. Secondly, we illustrate our approach by its application to an analysis of the nozzle fault in a solid rocket motor (SRM). The internal ballistics of the SRM is modelled as a set of one-dimensional partial differential equations coupled to the dynamics of the propellant regression. In this example we are specifically focussed on the inference of discrete and continuous parameters of the nozzle blocking fault and on the possibility of an application of the DI algorithm to reducing the probability of "misses" of an on-board FD&P for SRM. In the next section re-contact problem caused by first stage/upper stage separation failure is discussed. The reaction forces imposed on the nozzle of the upper stage during the re-contact and their connection to the nozzle damage and to the thrust vector control (TVC) signal are obtained. It is shown that transient impact induced torquean be modelled as a response of an effective damped oscillator. A possible application of the DI algorithm to the inference of damage parameters and predicting fault dynamics ahead of time using the actuator signal is discussed. Finally, we formulate Bayesian inferential framework for development of the IHM system for in-flight structural health monitoring (SHM) of composite materials. We consider the signal generated by piezoelectric actuator mounted on composite structure generating elastic waves in it. The signal received by the sensor is than compared with the baseline signal. The possibility of damage inference is discussed in the context of development of the SHM.

Smelyanskiy, Vadim N.↗