Search NASA⌕ Search

SEARCH · Search NASA

Results for “Modulational instability”

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.

319 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

A Cloud-Resolving Simulation of Hurricane Bob (1991): Storm Structure and Eyewall Buoyancy

A numerical simulation of Hurricane Bob (1991) is conducted using the Penn State University-National Center for Atmospheric Research mesoscale model MM5 with a horizontal grid spacing of 1.3 Km on the finest nested mesh The model produces a realistic hurricane that intensifies slowly during the period of fine-scale simulation. Time-averaged results reveal the effects of storm motion. vertical shear, beta gyres and deformation forcing on the structure of radial inflow, vertical motion, and precipitation. Instantaneous model fields show that radial inflow in the eyewall is very intense near the surface but transitions to strong low-level outflow near the top of the boundary layer. The low-level structure is modulated by a wavenumber 2 disturbance that rotates around the eyewall at half the speed of the maximum tangential winds and is consistent with a vortex Rossby edge wave. The statistical distribution of vertical velocity in the eyewall indicates that the eyewall is composed of a small number of intense updrafts that account for the majority of the upward mass flux rather than a more gradual and symmetric eyewall circulation, consistent with the concept of hot towers. Tongues of high equivalent potential temperature, Theta(sub e), are seen along the inner edge of the eyewall updraft and within the low-level outflow. This air originates from outside of the eyewall with the highest theta(sub e) air coming from the layer closest to the surface after penetrating closest to the center. Occasionally, high Theta(sub e), air within the eye is drawn into the eyewall updrafts. The high Theta(sub e), air rising within the eyewall is shown to be associated with positive eyewall buoyancy with sufficient convective available potential energy along its path to produce relatively strong convective updrafts. Although the requirements for conditional symmetric instability are met within the eyewall and the air parcel trajectories follow slanted paths, the radial displacement of air parcels in the low-level outflow moves the air parcel sufficiently far away from the upper- warm core that the air becomes unstable to vertical displacements. Hence, convective instability rather than symmetric instability accounts for the stronger updrafts in the eyewall.

Braun, Scott A.↗
NASA NTRS2011

AXIOM: Advanced X-Ray Imaging of the Magnetosphere

Planetary plasma and magnetic field environments can be studied in two complementary ways by in situ measurements, or by remote sensing. While the former provide precise information about plasma behaviour, instabilities and dynamics on local scales, the latter offers the global view necessary to understand the overall interaction of the magnetospheric plasma with the solar wind. Some parts of the Earth's magnetosphere have been remotely sensed, but the majority remains unexplored by this type of measurements. Here we propose a novel and more elegant approach employing remote X-ray imaging techniques, which are now possible thanks to the relatively recent discovery of solar wind charge exchange X-ray emissions in the vicinity of the Earth's magnetosphere. In this article we describe how an appropriately designed and located X-ray telescope, supported by simultaneous in situ measurements of the solar wind, can be used to image the dayside magnetosphere, magnetosheath and bow shock, with a temporal and spatial resolution sufficient to address several key outstanding questions concerning how the solar wind interacts with the Earth's magnetosphere on a global level. Global images of the dayside magnetospheric boundaries require vantage points well outside the magnetosphere. Our studies have led us to propose AXIOM: Advanced X-ray Imaging Of the Magnetosphere, a concept mission using a Vega launcher with a LISA Pathfinder-type Propulsion Module to place the spacecraft in a Lissajous orbit around the Earth Moon L1 point. The model payload consists of an X-ray Wide Field Imager, capable of both imaging and spectroscopy, and an in situ plasma and magnetic field measurement package. This package comprises a Proton-Alpha Sensor, designed to measure the bulk properties of the solar wind, an Ion Composition Analyser, to characterize the minor ion populations in the solar wind that cause charge exchange emission, and a Magnetometer, designed to measure the strength and direction of the solar wind magnetic field. We also show simulations that demonstrate how the proposed X-ray telescope design is capable of imaging the predicted emission from the dayside magnetosphere with the sensitivity and cadence required to achieve the science goals of the mission.

Branduardi-Raymont, G.↗
NASA NTRS2012

AXIOM: Advanced X-ray Imaging of the Magnetosphere

Planetary plasma and magnetic field environments can be studied in two complementary ways - by in situ measurements, or by remote sensing. While the former provide precise information about plasma behaviour, instabilities and dynamics on local scales, the latter offers the global view necessary to understand the overall interaction of the magnetospheric plasma with the solar wind. Some parts of the Earth's magnetosphere have been remotely sensed, but the majority remains unexplored by this type of measurements. Here we propose a novel and more elegant approach employing remote X-ray imaging techniques. which are now possible thanks to the relatively recent discovery of solar wind charge exchange X-ray emissions in the vicinity of the Earth's magnetosphere. In this article we describe how an appropriately designed and located. X-ray telescope, supported by simultaneous in situ measurements of the solar wind, can be used to image the dayside magnetosphere, magnetosheath and bow shock. with a temporal and spatial resolution sufficient to address several key outstanding questions concerning how the solar wind interacts with the Earth's magnetosphere on a global level. Global images of the dayside magnetospheric boundaries require vantage points well outside the magnetosphere. Our studies have led us to propose 'AXIOM: Advanced X-ray Imaging Of the Magnetosphere', a concept mission using a Vega launcher with a LISA Pathfinder-type Propulsion Module to place the spacecraft in a Lissajous orbit around the Earth - Moon Ll point. The model payload consists of an X-ray Wide Field Imager, capable of both imaging and spectroscopy, and an in situ plasma and magnetic field measurement package. This package comprises a Proton-Alpha Sensor, designed to measure the bulk properties of the solar wind, an Ion Composition Analyser, to characterize the minor ion populations in the solar wind that cause charge exchange emission, and a Magnetometer, designed to measure the strength and direction of the solar wind magnetic field. We also show simulations that demonstrate how the proposed X-ray telescope design is capable of imaging the predicted emission from the dayside magnetosphere with the sensitivity and cadence required to achieve the science goals of the mission.

Branduardi-Raymont, G.↗
DOE OSTI2024

Polymeric ionic liquids containing copper(I) and copper(II) ions as gas chromatographic stationary phases for olefin separations

Copper(I) ions (Cu + ) are used in olefin separations due to their olefin complexing ability and low cost, but their instability in the presence of water and gases limits their widespread use. Ionic liquids (ILs) have emerged as stabilizers of Cu + ions and prevent their degradation, providing high olefin separation efficiency. There is limited understanding into the role that polymeric ionic liquids (PILs), which possess similar structural characteristics to ILs, have on Cu + ion-olefin interactions. Moreover, copper ions with diverse oxidation states, including Cu + and Cu 2+ ions, have been rarely employed for olefin separations. In this study, gas chromatography (GC) is used to investigate the interaction strength of olefins to stationary phases composed of the 1-hexyl-3-methylimidazolium bis[(trifluoromethyl)sulfonyl]imide ([C 6 MIM + ][NTf 2 – ]) IL and the poly(1-hexyl-3-vinylimidazolium [NTf 2 – ]) (poly([C 6 VIM + ][NTf 2 – ])) PIL containing monovalent and divalent copper salts (i.e., [Cu + ][NTf 2 – ] and [Cu 2+ ]2[NTf 2 – ]). The chromatographic retention of alkenes, alkynes, dienes, and aromatic compounds was examined. Incorporation of the [Cu 2+ ]2[NTf 2 – ] salt into a stationary phase comprised of poly(dimethylsiloxane) resulted in strong retention of olefins, while its addition to the [C 6 MIM + ][NTf 2 – ] IL and poly([C 6 VIM + ][NTf 2 – ]) PIL allowed for the interaction strength to be modulated. Olefins exhibited greater affinities toward IL and PIL stationary phases containing the [Cu 2+ ]2[NTf 2 – ] salt compared to those with the [Cu + ][NTf 2 – ] salt. Elimination of water from both copper salts was observed to be an important factor in promoting olefin interactions, as evidenced by increased olefin retention upon exposure of the stationary phases to high temperatures. Furthermore, to evaluate the long-term thermal stability of the stationary phase, chromatographic retention of probes was measured on the [Cu 2+ ]2[NTf 2 – ]/[C 6 MIM + ][NTf 2 – ] IL stationary phase after its exposure to helium at a temperature of 110°C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Particle-in-cell simulations of the critical ionization velocity effect in finite size clouds

The critical ionization velocity (CIV) mechanism in a finite size cloud is studied with a series of electrostatic particle-in-cell simulations. It is observed that an initial seed ionization, produced by non-CIV mechanisms, generates a cross-field ion beam which excites a modified beam-plasma instability (MBPI) with frequency in the range of the lower hybrid frequency. The excited waves accelerate electrons along the magnetic field up to the ion drift energy that exceeds the ionization energy of the neutral atoms. The heated electrons in turn enhance the ion beam by electron-neutral impact ionization, which establishes a positive feedback loop in maintaining the CIV process. It is also found that the efficiency of the CIV mechanism depends on the finite size of the gas cloud in the following ways: (1) Along the ambient magnetic field the finite size of the cloud, L (sub parallel), restricts the growth of the fastest growing mode, with a wavelength lambda (sub m parallel), of the MBPI. The parallel electron heating at wave saturation scales approximately as (L (sub parallel)/lambda (sub m parallel)) (exp 1/2); (2) Momentum coupling between the cloud and the ambient plasma via the Alfven waves occurs as a result of the finite size of the cloud in the direction perpendicular to both the ambient magnetic field and the neutral drift. This reduces exponentially with time the relative drift between the ambient plasma and the neutrals. The timescale is inversely proportional to the Alfven velocity. (3) The transvers e charge separation field across the cloud was found to result in the modulation of the beam velocity which reduces the parallel heating of electrons and increases the transverse acceleration of electrons. (4) Some energetic electrons are lost from the cloud along the magnetic field at a rate characterized by the acoustic velocity, instead of the electron thermal velocity. The loss of energetic electrons from the cloud seems to be larger in the direction of plasma drift relative to the neutrals, where the loss rate is characterized by the neutral drift velocity. It is also shown that a factor of 4 increase in the ambient plasma density, increases the CIV ionization yield by almost 2 orders of magnitude at the end of a typical run. It is concluded that a larger ambient plasma density can result in a larger CIV yield because of (1) larger seed ion production by non-CIV mechanisms, (2) smaller Alfven velocity and hence weak momentum coupling, and (3) smaller ratio of the ion beam density to the ambient ion density, and therefore a weaker modulation of the beam velocity. The simulation results are used to interpret various chemical release experiments in space.

Moghaddam-Taaheri, E.↗
NASA NTRS2000

Mesosphere Dynamics with Gravity Wave Forcing: Planetary Waves - 2

We present results from a non-linear, 3D, time dependent numerical spectral model (NSM) which extends from the ground up into the thermosphere and incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Our focal point is the mesosphere where wave interactions are playing a dominant role. We discuss planetary waves in the present paper and diurnal and semi-diurnal tides in the companion paper. Without external time dependent energy or momentum sources, planetary waves (PWs) are generated in the model for zonal wavenumbers 1 to 4, which have amplitudes in the mesosphere above 50 km as large as 30 m/s and periods between 2 and 50 days. The waves are generated primarily during solstice conditions, which indicates that the baroclinic instability (associated with the GW driven reversal in the latitudinal temperature gradient) is playing an important role. Results from a numerical experiment show that GWs are also involved directly in generating the PWs. For the zonal wavenumber m = 1, the predominant wave periods in summer are around 4 days and in winter between 6 and 10 days. For m = 2, the periods are in summer and close to 2.5 and 3.5 days respectively For m = 3, 4 the predominant wave periods are in both seasons close to two days. The latter waves have the characteristics of Rossby gravity waves with meridional winds at equatorial latitudes. A common feature of the PWs (m = 1 to 4) generated in summer and winter is that their vertical wavelengths throughout the mesosphere are large which indicates that the waves are not propagating freely but are generated throughout the region. Another common feature is that the PWs propagate preferentially westward in summer and eastward in winter, being launched from the westward and eastward zonal winds that prevail respectively in summer and winter altitudes below 80 km. During spring and fall, for m = 1 and 2 eastward propagating long period PWs are generated that are launched from the smaller eastward zonal winds that prevail in these seasons. PWs generated in the model produce large amplitude modulations of the diurnal tides at altitudes above 80 km and contribute to their seasonal variations.

Mayr, H. G.↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗
NASA NTRS2023

Heat Transfer and Interfacial Flow Physics of Microgravity Flow Boiling in Single-Side-Heated Rectangular Channel with Subcooled Inlet Conditions – Experiments Onboard the International Space Station

This study is the culmination of a long-term collaborative effort between researchers from the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center to investigate gravitational effects on flow boiling and flow condensation. The science and design concepts for this large-scale effort were initiated in 2011 and included several studies detailing various aspects of two-phase fluid physics in both Earth gravity and microgravity, culminating in construction of the large-scale experimental facility named “Flow Boiling and Condensation Experiment (FBCE)”. The experiment was launched to the International Space Station (ISS) in August 2021. Following the successful installation of FBCE, equipped with the Flow Boiling Module (FBM), onboard the ISS and completion of several safety checks, flow boiling experiments were performed for five months from February 2022 until July 2022. This resulted in a large flow boiling database covering broad ranges of operating parameters and heating configurations spanning several research objectives. This study investigates microgravity flow boiling of n-perfluorohexane with subcooled inlet in a single-side-heated rectangular channel of dimensions 114.6-mm heated length, 2.5-mm heated width, and 5.0-mm height. Key operating parameters investigated are mass velocity (199.90 – 3200.13 kg/m 2 s), inlet subcooling (0.10 – 45.76°C), and inlet pressure (113.30 – 164.29 kPa). Images and image sequences acquired via high-speed-video are presented to elucidate the interfacial flow physics. To analyze and explain the effects of various parameters in microgravity, heat transfer results are presented as flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric trends of local and averaged heat transfer coefficient. Mass velocity and inlet subcooling significantly influenced most of the aforementioned aspects of flow boiling, whereas effects of inlet pressure were comparatively insignificant. Although the data and observed flow physics might be different, the parametric effects and trends in microgravity are similar to vertical upflow in Earth gravity. Some cases, especially low mass velocities, high heat fluxes, and large degrees of inlet subcooling, experienced temporally anomalous flow behaviors caused by two-phase flow instabilities manifesting as flow reversals and resulted in deviations in overall trends. Severe thermodynamic non-equilibrium is observed throughout the channel. Overall, FBCE’s ISS experiments were successful for subcooled inlet with single-sided heating of rectangular channel, and the collected data well established the various effects on flow boiling physics in highly controlled long-term microgravity conditions.

flow boiling↗
NASA NTRS2023

Effects of Heating Configuration and Operating Parameters on Heat Transfer and Interfacial Physics of Microgravity Flow Boiling With Subcooled Inlet Conditions –Experiments Onboard the International Space Station

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m 2 s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m 2 s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

Microgravity↗
NASA NTRS2007

Inter-annual and Long-term Temperature Variations in the Mesopause Region at High Latitudes Generated by the Stratospheric QBO

The Numerical Spectral Model (NSM) simulates the Quasi-biennial Oscillation (QBO) that dominates the zonal circulation of the lower stratosphere at low latitudes. In the model, the QBO is generated with parameterized small-scale gravity waves (GW), which are partially augmented in 3D with planetary waves owing to baroclinic instability. Due to GW filtering, the QBO extends into the upper mesosphere, evident in UARS zonal wind and TIMED temperature measurements. While the QBO zonal winds are confined to equatorial latitudes, even in simulations with latitude-independent wave source, the associated temperature variations extend to high latitudes. The meridional circulation redistributes some of the QBO energy to focus it partially onto the Polar Regions. The resulting QBO temperature variations away from the equator tend to increase at higher altitudes to produce inter-annual variations that can exceed 5 K in the polar mesopause region -- and our 3D model simulations show that the effect is variable from year to year and can produce large differences between the two hemispheres, presumably due to interactions involving the seasonal variations. Modeling studies with the NSM have shown that long-term variations can also be generated by the QBO interacting with the seasonal cycles through OW node-filtering. A 30-month QBO, optimally synchronized by the 6-month Semi-Annual Oscillation (SAO), thus produces a 5-year or semi-decadal (SD) oscillation -- and observational evidence for that has been provided by a recent analysis of stratospheric NCEP data. In a simulation with the 2D version of the NSM, this SD oscillation extends into the upper mesosphere, and we present results to show that the related temperature variations could contribute significantly to the long-term variations of the polar mesopause region. Quasi-decadal variations could furthermore arise from the modeled solar cycle modulations of the QBO and 12-month annual oscillation. Our numerical results are discussed in the context of the observed low summer temperatures reproduced by the model, to demonstrate that the above interannual and long-term variations could contribute significantly to the climatology of Polar Mesospheric Clouds (PMC) investigated by the Aeronomy of Ice in the Mesosphere (AIM) mission.

Mayr, Hans G.↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
DOE OSTI2024

Updraft Width Modulates Ambient Atmospheric Controls on Convective Cloud Depth

Abstract The depth of convective clouds affects vertical transport of atmospheric constituents, influencing downstream weather and climate. Atmospheric controls on the maximum depth reached by moist convection are investigated with radar‐tracked convective cells tagged with sounding‐derived atmospheric parameters from a field campaign in central Argentina. Regression analyses show that narrow (<12‐km diameter) and wide (>16‐km diameter) cell depths respond to disparate factors, where cell areas are defined using composite reflectivity signatures. Undiluted lifted parcel indices including convective available potential energy (CAPE) and level of neutral buoyancy (LNB) are top predictors of wide cell maximum depth while mid‐tropospheric relative humidity is the top predictor of narrow cell maximum depth. Because narrow cells are more numerous than wide cells, the overall outcome of the full cell population does not strongly correlate with CAPE and LNB conditions. Tracked cells and atmospheric conditions in a simulation with 3‐km grid spacing covering the field campaign produce similar results to those observed. Narrow cells that are relatively deep have a cooler and moister mid‐troposphere with weaker free tropospheric subsidence, while relatively deep wide cells have much warmer and moister lower tropospheric conditions. These atmospheric differences are present 1 hr before cell initiation at both a fixed observing site and variable cell initiation locations. Simulated narrow cell maximum equivalent potential temperature decreases with height at a rate similar to the ambient vertical gradient, causing these cells to fall short of their LNB and supporting the view that entrainment‐driven dilution is a dominant control on their depth.

54 ENVIRONMENTAL SCIENCES↗