Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pulse shaping”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Understanding JET-C quiescent phases with edge harmonic magnetohydrodynamic activity and comparison with behaviour under ITER-like wall conditioning

An analysis of edge localised mode-free (quiescent) H-mode discharges exhibiting edge harmonic magnetoydrodynamic activity in the JET-carbon wall machine is presented. It is observed that the otherwise quiescent pulses with multiple-n harmonic oscillations are sustained until a threshold in pedestal electron density and collisionality is crossed. The macroscopic pedestal parameters associated with the quiescent phase are compared with those of a database of JET-ELMy discharges with both carbon and ITER-like wall (ILW). This comparison provides the identification of the existence regions in the relevant pedestal and global plasma parameters for edge harmonic oscillations (EHOs) in JET plasmas. Although the ELMy database scans pedestal collisionality and β values typical of ET-carbon quiescent operation, shaping and current are not simultaneously compatible with EHO existence. Nevertheless, ILW operation with JET-carbon quiescent-like parameters could in principle be achieved, and improved pedestal performance could be observed in more recent JET-ILW pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Analysis of Pulsed Airborne Lidar Measurements of Atmospheric CO2 Column Absorption from 3-13 km Altitudes

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS space mission [1]. It uses two pulsed laser transmitters allowing simultaneous measurement of a CO2 absorption line in the 1575 nm band, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are precisely stepped in wavelength across the CO2 line and an O2 line region during the measurement. The direct detection receiver measures the energies of the laser echoes from the surface along with the range profile of scattering in the path. The column densities for the CO2 and O2 gases are estimated from the ratio of the on- and off-line signals via the integrated path differential absorption (IPDA) technique. The time of flight of the laser pulses is used to estimate the height of the scattering surface and to reject laser photons scattered in the atmosphere. We developed an airborne lidar to demonstrate an early version of the CO2 measurement from the NASA Glenn Lear-25 aircraft. The airborne lidar stepped the pulsed laser's wavelength across the selected CO2 line with 20 wavelength steps per scan. The line scan rate is 450 Hz, the laser pulse widths are 1 usec, and laser pulse energy is 24 uJ. The time resolved laser backscatter is collected by a 20 cm telescope, detected by a NIR photomultiplier and is recorded on every other reading by a photon counting system [2]. During August 2009 we made a series of 2.5 hour long flights and measured the atmospheric CO2 absorption and line shapes using the 1572.33 nm CO2 line. Measurements were made at stepped altitudes from 3-13 km over locations in the US, including the SGP ARM site in Oklahoma, central Illinois, north-eastern North Carolina, and over the Chesapeake Bay and the eastern shore of Virginia. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes, and some measurements were made through thin clouds. The Oklahoma and east coast flights were coordinated with a LaRC/ITT CO2 lidar on the LaRC UC-12 aircraft, and in-situ measurements were made using its CO2 sensor and radiosondes. We have conducted an analysis of the ranging and IPDA lidar measurements from these four flights. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We used a cross-correlation approach to process the laser echo records. This was used to estimate the range to the scattering surface, to define the edges of the laser pulses and to determine echo pulse energy at each wavelength. We used a minimum mean square approach to fit an instrument response function and to solve for the best-fit CO2 absorption line shape. We then calculated the differential optical depth (DOD) of the fitted CO2 line. We computed its statistics at the various altitude steps, and compare them to the DODs calculated from spectroscopy based on HITRAN 2008 and the column conditions calculated from the airborne in-situ readings. The results show the lidar and in-situ measurements have very similar DOD change with altitude and greater than 10 segments per flight where the scatter in the lidar measurements are less than or equal to 1ppm. We also present the results from subsequent CO2 column absorption measurements, which were made with stronger detected signals during three flights on the NASA DC-8 over the southwestern US in during July 2010.

Abshire, James B.↗
DOE OSTI2023

Numerical and Experimental Study of an Aircraft Igniter Plasma Jet Discharge

The spark discharge of an aircraft plasma jet igniter is studied using high-fidelity numerical simulations and X-ray radiography measurements. The target problem here features the thermal expansion of hot gas introduced by the electric spark within a confined igniter cavity, which eventually evolves into a pulsed jet of a high-temperature kernel. A comprehensive set of models adapted from existing strategies for internal combustion engine spark plug discharge is extended to the target problem, including the modeling of energy deposition, plasma reactions, thermodynamic properties, and heat losses. A series of validation and parameter studies are performed and presented. The kernel size is found to be sensitive to heat losses arising from radiation and hot gas remained within the discharge cavity, rather than heat conduction to the wall in the discharge cavity. Depending on the enforced shape of the post-breakdown electric arc, the spark kernel can be off-centered, tilted, and considerably asymmetric. These features have been previously not considered when studying such igniter configurations and may have a first-order impact on the ignition process. Provided a proper setup of the heat loss models and electric arc shape, the numerical results are quantitatively comparable to the experimental results in terms of the kernel size, shape, and velocity throughout different stages after the spark discharge.

Tang, Yihao↗
DOE OSTI2021

Performance of Compact Pulsed Thermal Imaging System for In-Service Applications. Pulsed thermal tomography nondestructive examination of additively manufactured reactor materials and components

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of complex topology nuclear reactor parts from high-strength corrosion resistance alloys, such as stainless steel and Inconel. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which has the capability of melting metallic powder and net shaping the structures with relatively high precision. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to creep in high temperature nuclear reactor environment. Currently, there exist limited capabilities to evaluate actual AM structures nondestructively. Pulsed Thermography (PT) imaging provides a capability for non-destructive evaluation (NDE) of sub-surface defects in arbitrary size structures. The PT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. The data cube of PT measurements consists of surface temperature taken at sequential time intervals T(x,y,t). Material defects can be detected either by analyzing the thermograms T(x,y,t) data cube, or by using thermal tomography (TT) algorithm to obtain 3D spatial reconstruction of thermal effusivity e(x,y,z). To reduce the cost and enable in-service NDE in spatially constrained environment, it is highly desirable to develop PT with compact and inexpensive IR camera. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact PT system can also be used for in-service nondestructive evaluation (NDE) applications. However, data cube obtained with PT based on compact IR camera suffers from strong thermal noises and loss of features due to relatively low sampling rate. In this report we describe two unsupervised machine learning (ML) algorithms for enhancement of PT images obtained with compact IR camera. In one approach, we introduce Sparse Coding Discrete Cosine Transform (SC/DCT) algorithm to remove additive white Gaussian noise (AWGN) from spatial thermal effusivity reconstructions. In another approach we introduce a Spatial Temporal Denoised Thermal Source Separation (STDTSS) ML algorithm to process thermograms. The STDTSS algorithm consists of spatial and temporal denoising using Gaussian and Savitzky–Golay filtering, followed by the matrix decomposition using Principal Component Analysis (PCA), and Independent Component Analysis (ICA) to automatically detect flaws. In the work described in this report, we constructed a compact PT system using a relatively small and low-cost FLIR A65 camera, consisting on uncooled microbolometer detector. Performance of SC/DCT algorithm was demonstrated on enhancing TT images of Inconel 718 AM plate. Performance of the STDTSS methods was investigated using thermography data obtained from imaging stainless steel 316L specimens produced with LPBF method with imprinted calibrated porosity defects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Morphology Control of Self-Assembled Three-Phase Au-BaTiO 3 –ZnO Hybrid Metamaterial for Tunable Optical Properties

Microstructural control in metal-dielectric hybrid metamaterials presents enormous opportunities in tailoring the physical properties including the magnetic and optical properties. In this paper, we demonstrate a strong tunability achieved in the microstructure of self-assembled ordered three-phase Au-BaTiO 3 –ZnO hybrid metamaterial along with its optical properties, grown by a pulsed laser deposition method. Varying the growth temperature, deposition frequency, and template thickness evolves the microstructure by tuning the Au and ZnO pillar geometry as well as the shape and size of the Au nanoparticles capping the ZnO nanowires. The three-phase hybrid metamaterials exhibit unique optical properties, including enhanced nonlinear optical properties, hyperbolic dispersion in the visible and near-infrared wavelength region, and tuned epsilon-near-zero (ENZ) wavelength upon varying the deposition parameters. This study suggests that the three-phase hybrid metamaterials present great potential in the microstructure and optical property tuning that can also be applied to other two-phase and three-phase nanocomposite systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Spectral Properties of Globally Distributed ENA Fluxes across Diverse Regions of the Heliosphere

This study analyzes energetic neutral atom (ENA) spectral properties across distinct regions of globally distributed flux (GDF) sky maps, using Interstellar Boundary Explorer data from a full solar cycle, corrected for time dispersion. By time-shifting the data to the heliosheath using GDF source distances from D. B. Reisenfeld et al., we achieve a more accurate representation of heliosheath GDF energy spectra. We quantify ENA spectral characteristics, heliosheath line-of-sight-integrated proton pressure, and heliosheath proton temperature, comparing these to solar wind properties at 1 au and interplanetary scintillation-derived solar wind data. Our findings show that the spectral index is generally anticorrelated with heliosheath proton temperature and pressure, except in the central tail, where a partial positive correlation is observed. The lowest spectral index values occur when high-latitude heliosheath regions are dominated by fast solar wind from polar coronal holes. The south pole exhibits the flattest energy spectra due to plasma heating from both fast solar wind and a late-2014 pressure pulse. The central tail shows shorter variability (5–6 yr) for spectral index and heliosheath proton temperature, while proton pressure follows the 11 yr solar cycle. Most spectral shapes exhibit a “knee” distribution, peaking during solar maximum, with an “ankle” shape observed only at the south pole during solar cycle transitions. Asymmetry in proton pressure in the lobes is driven by the draping effect of the local interstellar magnetic field. This study provides insights into the energetic properties of GDF across the heliosphere, enhancing our understanding of the heliospheric environment.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

Short pulse radar used to measure sea surface wind speed and SWH

A joint airborne measurement program is being pursued by NRL and NASA Wallops Flight Center to determine the extent to which wind speed and sea surface significant wave height (SWH) can be measured quantitatively and remotely with a short pulse (2 ns), wide-beam (60 deg), nadir-looking 3-cm radar. The concept involves relative power measurements only and does not need a scanning antenna, Doppler filters, or absolute power calibration. The slopes of the leading and trailing edges of the averaged received power for the pulse limited altimeter are used to infer SWH and surface wind speed. The interpretation is based on theoretical models of the effects of SWH on the leading edge shape and rms sea-surface slope on the trailing-edge shape. The models include the radar system parameters of antenna beam width and pulsewidth.

Hammond, D. L.↗
NASA NTRS1983

Refinement and application of acoustic impulse technique to study nozzle transmission characteristics

An improved acoustic impulse technique was developed and was used to study the transmission characteristics of duct/nozzle systems. To accomplish the above objective, various problems associated with the existing spark-discharge impulse technique were first studied. These included (1) the nonlinear behavior of high intensity pulses, (2) the contamination of the signal with flow noise, (3) low signal-to-noise ratio at high exhaust velocities, and (4) the inability to control or shape the signal generated by the source, specially when multiple spark points were used as the source. The first step to resolve these problems was the replacement of the spark-discharge source with electroacoustic driver(s). These included (1) synthesizing on acoustic impulse with acoustic driver(s) to control and shape the output signal, (2) time domain signal averaging to remove flow noise from the contaminated signal, (3) signal editing to remove unwanted portions of the time history, (4) spectral averaging, and (5) numerical smoothing. The acoustic power measurement technique was improved by taking multiple induct measurements and by a modal decomposition process to account for the contribution of higher order modes in the power computation. The improved acoustic impulse technique was then validated by comparing the results derived by an impedance tube method. The mechanism of acoustic power loss, that occurs when sound is transmitted through nozzle terminations, was investigated. Finally, the refined impulse technique was applied to obtain more accurate results for the acoustic transmission characteristics of a conical nozzle and a multi-lobe multi-tube supressor nozzle.

Salikuddin, M.↗
NASA NTRS1999

Internal Heterogeneous Processes in Aluminum Combustion

This paper discusses the aluminum particle combustion mechanism which has been expanded by inclusion of gas dissolution processes and ensuing internal phase transformations. This mechanism is proposed based on recent normal and microgravity experiments with particles formed and ignited in a pulsed micro-arc. Recent experimental findings on the three stages observed in Al particle combustion in air and shows the burning particle radiation, trajectory (streak), smoke cloud shapes, and quenched particle interiors are summarized. During stage I, the radiation trace is smooth and the particle flame is spherically symmetric. The temperature measured using a three-color pyrometer is close to 3000 K. Because it exceeds the aluminum boiling point (2730 K), this temperature most likely characterizes the vapor phase flame zone rather than the aluminum surface. The dissolved oxygen content within particles quenched during stage I was below the detection sensitivity (about 1 atomic %) for Wavelength Dispersive Spectroscopy (WDS). After an increase in the radiation intensity (and simultaneous decrease in the measured color temperature from about 3000 to 2800 K) indicative of the transition to stage II combustion, the internal compositions of the quenched particles change. Both oxygen-rich (approx. 10 atomic %) and oxygen-lean (< 1 %) regions are identified within the particles using back-scattered electron imaging and WDS. During stage II, oscillations are observed in particle radiation and the flame and smoke cloud are distorted from their original spherically-symmetric shape. In stage III, particle radiation continues to exhibit oscillations, but its radiation intensity drops and remains at a nearly constant level. The measured temperature decreases to about 2300 K. Also, larger changes in particle velocities are observed, and oxide caps are found on quenched particle surfaces. While these results showed the correlation between the aluminum particle combustion behavior and the evolution of its internal composition, the change from the spherically symmetric to asymmetric flame shape occurring upon the transition from stage I to stage II combustion could not be understood based only on the fact that dissolved oxygen is detected in the particles. The connection between the two phenomena appeared even less significant because in earlier aluminum combustion studies carried in O2/Ar mixtures, flame asymmetry was not observed as opposed to experiments in air or O2/CO mixtures. It has been proposed that the presence of other gases, i.e., hydrogen, or nitrogen causes the change in the combustion regime.

Dreizin, E. L.↗
DOE OSTI2021

ATHENA: A unique radiation environment platform at the National Ignition Facility

This paper describes the ATHENA platform, an energy tuning assembly, which was developed to spectrally shape the National Ignition Facility (NIF) deuterium–tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with a capability to act as a short-pulse neutron source. This unique, otherwise inaccessible radiation environment complements existing experimental facilities and capabilities. Here, the flexible ATHENA irradiation positions were modeled using an ensemble of Monte Carlo simulations with stochastic sampling of the nuclear cross-sections to characterize the radiation environments and uncertainty for the platform. Validation of the internal neutron spectrum produced from fielding ATHENA at NIF occurred through neutron flux unfolding with 20 measured activation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Bus Design for the Poloidal Field Coils of the NSTX-Upgrade Fusion Device

The National Spherical Torus eXperiment (NSTX) has undergone a major upgrade to NSTX-U at Princeton Plasma Physics Laboratory (PPPL). NSTX upgrade (NSTX-U) will double the toroidal field, plasma current, and neutral beam injection heating power, as well as significantly increase the pulse duration. NSTX-U uses three poloidal field (PF) coils at the vessel top and three at the bottom near the divertor areas to control the local plasma shape there. These shaping coils operate at a maximum terminal voltage of 2 kV, corresponding to a maximum current about 20 kA. To supply the power to the PF coils, hard copper bus bars are typically used near the vacuum vessel, while water cooled flexible cables are used away from the vessel. The detailed design of the hard copper bus bars will be covered in this article. During operation, these hard bus bars are subject to high electromagnetic (EM) forces, thermal displacement loads, as well as plasma halo loads due to disruption. The EM, thermal, and structural analysis were performed, and the results revealed that, with the worst loads combined, the new design meets the NSTX-U thermal, structural, and fatigue cycle requirements. Finally, the manufacturing and installation process of the bus work will be discussed at the end of this article.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Athena-I CUBIT Journal Files

The Monte Carlo N-Particle (MCNP)1 transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into con structive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is a time-consuming and error prone process as the complexities of geometries increase. The UM capability was originally designed to work with UM models created with the Abaqus software and ASCII input files that it generates. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with version 6.3, MCNP can also process HDF5 mesh input files. External codes must be used to generate Abaqus input files for MCNP UM calculations. CUBIT, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files exported from CUBIT cannot be used for MCNP simulations because it lacks the proper syntax. A Python script was developed to convert an Abaqus file created by CUBIT to an Abaqus file format that MCNP can process. Creating UM models for complex geometries is not an easy task. The process of creating UM models in CUBIT for MCNP simulations is detailed in. CUBIT provides several user interface options including a graphical user interface (GUI) and a command line interface. A GUI provides an easy way to use CUBIT without learning the CUBIT command syntax. When using CUBIT with either interface option, command lines are written into an ASCII file known as a journal file; this journal file can be edited and archived so that it can be played back in CUBIT to automatically generate a UM model. This report describes the CUBIT journal files of the UM models developed for Athena-I. The Athena platform, an energy-tuning assembly, was developed to spectrally shape the National Ignition Facility (NIF) deuterium-tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with capability to act as a short-pulse neutron source. MCNP6 was used for the Athena experiment design analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Simulation of Acoustics for Ares I Scale Model Acoustic Tests

The Ares I Scale Model Acoustics Test (ASMAT) is a series of live-fire tests of scaled rocket motors meant to simulate the conditions of the Ares I launch configuration. These tests have provided a well documented set of high fidelity acoustic measurements useful for validation including data taken over a range of test conditions and containing phenomena like Ignition Over-Pressure and water suppression of acoustics. To take advantage of this data, a digital representation of the ASMAT test setup has been constructed and test firings of the motor have been simulated using the Loci/CHEM computational fluid dynamics software. Results from ASMAT simulations with the rocket in both held down and elevated configurations, as well as with and without water suppression have been compared to acoustic data collected from similar live-fire tests. Results of acoustic comparisons have shown good correlation with the amplitude and temporal shape of pressure features and reasonable spectral accuracy up to approximately 1000 Hz. Major plume and acoustic features have been well captured including the plume shock structure, the igniter pulse transient, and the ignition overpressure.

Putnam, Gabriel↗
NASA NTRS2011

Validation and Simulation of ARES I Scale Model Acoustic Test -1- Pathfinder Development

The Ares I Scale Model Acoustics Test (ASMAT) is a series of live-fire tests of scaled rocket motors meant to simulate the conditions of the Ares I launch configuration. These tests have provided a well documented set of high fidelity measurements useful for validation including data taken over a range of test conditions and containing phenomena like Ignition Over-Pressure and water suppression of acoustics. To take advantage of this data, a digital representation of the ASMAT test setup has been constructed and test firings of the motor have been simulated using the Loci/CHEM computational fluid dynamics software. Within this first of a series of papers, results from ASMAT simulations with the rocket in a held down configuration and without water suppression have then been compared to acoustic data collected from similar live-fire tests to assess the accuracy of the simulations. Detailed evaluations of the mesh features, mesh length scales relative to acoustic signals, Courant-Friedrichs-Lewy numbers, and spatial residual sources have been performed to support this assessment. Results of acoustic comparisons have shown good correlation with the amplitude and temporal shape of pressure features and reasonable spectral accuracy up to approximately 1000 Hz. Major plume and acoustic features have been well captured including the plume shock structure, the igniter pulse transient, and the ignition overpressure. Finally, acoustic propagation patterns illustrated a previously unconsidered issue of tower placement inline with the high intensity overpressure propagation path.

Putnam, G. C.↗
NASA NTRS2024

The Giant Out Burst of EXO 2030+375 II: Broadband Spectroscopy and Evolution

In 2021, the High-mass X-ray Binary EXO 2030+375 underwent a giant X-ray outburst, the first since 2006, which reached a peak flux of ∼600 mCrab (3–50 keV). The goal of this work is to study the spectral evolution over the course of the outburst, search for possible cyclotron resonance scattering features (CRSFs), and to associate spectral components with the emission pattern of the accretion column. We use broadband spectra taken with NuSTAR, NICER, and Chandra near the peak and during the decline phase of the outburst. We describe the data with established empirical continuum models and perform pulse-phase-resolved spectroscopy. We compare the spectral evolution with pulse phase with a proposed geometrical emission model. We find a significant spectral hardening toward lower luminosity, a behavior that is expected for super-critical sources. The continuum shape and evolution cannot be described by a simple power law model with exponential cutoff but requires additional absorption or emission components. We can confirm the presence of a narrow absorption feature at ∼10 keV in both NuSTAR observations. The absence of harmonics question the interpretation of this feature as a CRSF. The empirical spectral components cannot be directly associated with identified emission components from the accretion column.

R. Ballhausen↗
NASA NTRSDate not supplied

Surface Instability of Liquid Propellants in Microgravity During Pulsed Settling Operations

Pulsing reaction control system (RCS) thrusters, vent valves, or other propulsion devices can preserve propellant resources in space, but this operation also effectively introduces a vibration to the vehicle. When the vibration is perpendicular to the liquid propellant surface, Faraday waves may be generated at the liquid-vapor interface. These Faraday instabilities can perturb or break up the liquid surface of cryogenic tanks, leading to inefficiencies in thermal management or even ullage collapse. Drawing from theory and experiments, an engineering model defining the allowable design regions for pulsed settling in microgravity was assembled and verified with computational fluid dynamics (CFD) simulations. A traditional settling metric, the Bond number, was also overlaid in the model to indicate which duty cycles were insufficient to overcome surface tension and aggregate propellant. Mission planners and engineers can consult the tool to rapidly evaluate the stability of a liquid-vapor interface given the pulse frequency and the excitation acceleration. Expressions developed for Faraday waves induced by a sinusoidal forcing input at standard gravity were found to provide excellent predictive capabilities for pulsed, or rectangular, waveforms in the absence of a consistent gravitational acceleration. This study extends the usage of these equations to an alternative forcing function and microgravity environments for the purpose of estimating natural frequencies, surface mode shapes, surface wave amplitudes, and the onset of droplet ejection. CFD simulations with the Loci/STREAM-VoF (Volume of Fluid) solver were initially validated against experimental results in standard gravity. Discrete points on the design map were then investigated with CFD and confirmed that the engineering model reliably indicates surface stability and most Faraday wave characteristics without requiring higher-fidelity tools. The engineering model is highly extensible and can be adapted for various propellant fill fractions, fluid properties, and tank sizes.

Faraday Waves↗
NASA NTRSDate not supplied

Surface Instability of Liquid Propellants in Microgravity During Pulsed Settling Operations

Pulsing reaction control system (RCS) thrusters, vent valves, or other propulsion devices can preserve propellant resources in space, but this operation also effectively introduces a vibration to the vehicle. When the vibration is perpendicular to the liquid propellant surface, Faraday waves may be generated at the liquid-vapor interface. These Faraday instabilities can perturb or break up the liquid surface of cryogenic tanks, leading to inefficiencies in thermal management or even ullage collapse. Drawing from theory and experiments, an engineering model defining the allowable design regions for pulsed settling in microgravity was assembled and verified with computational fluid dynamics (CFD) simulations. A traditional settling metric, the Bond number, was also overlaid in the model to indicate which duty cycles were insufficient to overcome surface tension and aggregate propellant. Mission planners and engineers can consult the tool to rapidly evaluate the stability of a liquid-vapor interface given the pulse frequency and the excitation acceleration. Expressions developed for Faraday waves induced by a sinusoidal forcing input at standard gravity were found to provide excellent predictive capabilities for pulsed, or rectangular, waveforms in the absence of a consistent gravitational acceleration. This study extends the usage of these equations to an alternative forcing function and microgravity environments for the purpose of estimating natural frequencies, surface mode shapes, surface wave amplitudes, and the onset of droplet ejection. CFD simulations with the Loci/STREAM-VoF (Volume of Fluid) solver were initially validated against experimental results in standard gravity. Discrete points on the design map were then investigated with CFD and confirmed that the engineering model reliably indicates surface stability and most Faraday wave characteristics without requiring higher-fidelity tools. The engineering model is highly extensible and can be adapted for various propellant fill fractions, fluid properties, and tank sizes.

Faraday Waves↗