Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dense flow”

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 415 records · Page 23

Effect of surface defects on damage and failure of dynamically driven cylindrical shells

Here, in this work, expanding cylinder experiments are performed using a gas gun on the aluminum 6061 alloy to understand the dynamic deformation and failure behavior of this material in the presence of surface defects in the form of scratches. The experiments are then modeled with the modified Tepla (for tensile plasticity) model, which is equipped with the damage yield function of Gurson [“Continuum theory of ductile rupture by void nucleation and growth: Part I—Yield criteria and flow rules for porous ductile media,” (1977)] and Tvergaard and Needleman [Acta Metall. 32(1), 157–169 (1984)], over-stress viscosity [F. Addessio and J. Johnson, J. Appl. Phys. 74(3), 1640–1648 (1993)], and the underlying viscoplastic model of Preston et al. [J. Appl. Phys. 93(1), 211–220 (2003)] (Preston–Tonks–Wallace model) for the dense material. The Tepla model parameters were calibrated using the Bayesian approach against experimentally measured free-surface velocity data obtained from plate impact experiments on the aluminum 6061 alloy. The evolution of radial expansion and shapes of the deformed cylinder from the simulations on a pristine cylinder wall and a cylinder wall with a longitudinal surface defect shows excellent agreement with the corresponding experimental results. Finally, the role of surface defect size and cylinder wall geometry on the dynamic strength of the cylinder wall was also investigated through a series of simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗

Volcanic recycling of carbonates on Mars

Thermal erosion of carbonate deposits by turbulently-flowing lava is investigated as a means of recycling carbon dioxide back into the atmosphere of Mars. Erosion rates of several meters/day are found, implying that up to hundreds of meters of carbonate could be removed over the lifetime of a flow. A large fraction of the northern plains and other parts of Mars were covered by lava during the Hesperian, and may have released the carbon dioxide trapped in carbonate deposits. This period of time, several times 10 exp 8 yrs, is comparable to that for the redeposition of such carbonate deposits. Therefore, there could have existed a relatively dense atmosphere, and enhanced weathering and erosion, after the Noachian era. This may help explain the apparent observational evidence for late fluvial and lacustrine activity on Mars.

Schaefer, M. W.↗

Physiological water model development

The water of the human body can be categorized as existing in two main compartments: intracellular water and extracellular water. The intracellular water consists of all the water within the cells and constitutes over half of the total body water. Since red blood cells are surrounded by plasma, and all other cells are surrounded by interstitial fluid, the intracellular compartment has been subdivided to represent these two cell types. The extracellular water, which includes all of the fluid outside of the cells, can be further subdivided into compartments which represent the interstitial fluid, circulating blood plasma, lymph, and transcellular water. The interstitial fluid surrounds cells outside of the vascular system whereas plasma is contained within the blood vessels. Avascular tissues such as dense connective tissue and cartilage contain interstitial water which slowly equilibrates with tracers used to determine extracellular fluid volume. For this reason, additional compartments are sometimes used to represent these avascular tissues. The average size of each compartment, in terms of percent body weight, has been determined for adult males and females. These compartments and the forces which cause flow between them are presented. The kidneys, a main compartment, receive about 25 percent of the cardiac output and filters out a fluid similar to plasma. The composition of this filtered fluid changes as it flows through the kidney tubules since compounds are continually being secreted and reabsorbed. Through this mechanism, the kidneys eliminate wastes while conserving body water, electrolytes, and metabolites. Since sodium accounts for over 90 percent of the cations in the extracellular fluid, and the number of cations is balanced by the number of anions, considering the renal handling sodium and water only should sufficiently describe the relationship between the plasma compartment and kidneys. A kidney function model is presented which has been adapted from a previous model of normal renal function in man. To test the validity of the proposed kidney model, results predicted by the model will be compared to actual data involving injected or ingested fluids and subsequent urine flow rates. Comparison of the model simulation to actual data following the ingestion of 1 liter of water is shown. The model simulation is also shown with actual data following the intravenous infusion of hypertonic saline.

Doty, Susan↗

Validation and Analysis of SRTM and VCL Data Over Tropical Volcanoes

The focus of our investigation was on the application of digital topographic data in conducting first-order volcanological and structural studies of tropical volcanoes, focusing on the Java, the Philippines and the Galapagos Islands. Kilauea volcano, Hawaii, served as our test site for SRTM data validation. Volcanoes in humid tropical environments are frequently cloud covered, typically densely vegetated and erode rapidly, so that it was expected that new insights into the styles of eruption of these volcanoes could be obtained from analysis of topographic data. For instance, in certain parts of the world, such as Indonesia, even the regional structural context of volcanic centers is poorly known, and the distribution of volcanic products (e.g., lava flows, pyroclastic flows, and lahars) are not well mapped. SRTM and Vegetation Canopy Lidar (VCL) data were expected to provide new information on these volcanoes. Due to the cancellation of the VCL mission, we did not conduct any lidar studies during the duration of this project. Digital elevation models (DEMs) such as those collected by SRTM provide quantitative information about the time-integrated typical activity on a volcano and allow an assessment of the spatial and temporal contributions of various constructional and destructional processes to each volcano's present morphology. For basaltic volcanoes, P_c?w!m-d and Garbed (2000) have shown that gradual slopes (less than 5 deg.) occur where lava and tephra pond within calderas or in the saddles between adjacent volcanoes, as well as where lava deltas coalesce to form coastal plains. Vent concentration zones (axes of rift zones) have slopes ranging from 10 deg. to 12 deg. Differential vertical growth rates between vent concentration zones and adjacent mostly-lava flanks produce steep constructional slopes up to 40". The steepest slopes (locally approaching 90 deg.) are produced by fluvial erosion, caldera collapse, faulting, and catastrophic avalanches, all of which are usually identifiable. Due to the delay in the release of the SRTM data following the February 2000 flight, a significant part of our effort was devoted to the analog studies of the SRTM topographic data using topographic data from airborne interferometric radars. As part of the original SRTM Science Team, we proposed four study sites (Kilauea, Hawaii; Mt. Pinatubo, Philippines; Cerro Am1 and Femandina volcanoes, Galapagos Islands; and Tengger caldera, Java) where we could conduct detailed geologic studies to evaluate the uses of SRTM data for the analysis of lava flows, lahars, erosion of ash deposits, and an evaluation of the structural setting of the volcanoes. Only near the end of this project was one of these SRTM Science Team products (Luzon Island, the Philippines) released to the community, and we only had limited time to work on these data.

Mouginis-Mark, Peter J.↗

Ionization fronts and shocked flows - The structure of the Orion Nebula at 0.1 arcsec

Initial results are presented from images of the Orion Nebula obtained with the HST Wide-Field Camera with a resolution of better than about 10 to the 15th cm. The morphology of the forbidden S II emission is markedly different from that seen in H-beta and forbidden O II. The forbidden S II is dominated by filamentary features with widths between 0.1 and 1 arcsesc which sharply highlight ionization fronts moving into dense neutral material. While the surfaces into which these fronts are moving are textured, subarcsec clumps with high density contrast are uncommon. The fields contains a number of HH objects whose structure is revealed by the HST, and which are described.

Hester, J. Jeff↗

Development of a Test-Bed for Testing and Refining EarthEn’s Supercritical CO 2 Based Energy Storage System

EarthEn’s energy storage concept leverages supercritical carbon dioxide (sCO 2 ) as a working fluid and relies on compact, high-performance components operating at elevated pressures and temperatures. To accelerate component development and reduce technical risk prior to larger-scale demonstrations, Oak Ridge National Laboratory (ORNL) developed a 100 kW-scale sCO 2 test-bed under a Cooperative Research and Development Agreement with EarthEn (CRADA NO. NFE-24-10050). The objective of the work was to design and construct a flexible experimental facility capable of reproducing key thermodynamic state points and heat-transfer conditions relevant to EarthEn’s thermal energy storage (TES) cycle, with particular emphasis on enabling development and evaluation of next-generation heat exchangers and TES concepts. The test-bed consists of a closed-loop sCO 2 circulation system housed within an open-topped enclosure. In its as-installed configuration, dense-phase sCO 2 is recirculated through a printed circuit recuperator, an electrically heated section, a throttling device used to simulate turbine expansion, and a water-cooled printed circuit heat exchanger that rejects heat to the building chilled-water system before returning to the pump. The pump is driven by a variable frequency drive, enabling controlled adjustment of flow and operating point. A comprehensive instrumentation suite was integrated to support both safe operation and high-quality data collection. Installed sensors include Coriolis flow meters for sCO 2 flow rate and density, resistance temperature detectors and thermocouples distributed throughout the loop (including the heated section and key heat exchanger ports), and pressure transducers for absolute and differential pressure measurements. The facility was designed to support high-pressure (19 MPa nominal) and high-temperature (575°C nominal) operation with credited overpressure protection provided by a rupture disk. Nominal operating conditions were selected to support 100 kW-class testing while maintaining flexibility for non-heated and heated shakedown, control development, and future integration of advanced TES test sections. In parallel with facility development, a system-level thermal-hydraulic model was created using Modelica-based tools to support component sizing, anticipate performance over targeted test conditions, and establish a framework for future model calibration against experimental data. At the conclusion of the project performance period, the facility was in final assembly, and the pressure boundary was nearly completed. However, several practical challenges associated with high-pressure/high-temperature systems and specialized component procurement impacted schedule and prevented initial pump-driven operation and full commissioning within the available resources. This report documents the as-built design, operating capabilities, and instrumentation, and it summarizes key lessons learned related to heater fabrication and testing, first-of-a-kind assembly factors, specialty flange supply constraints, and fill pump corrective actions. Finally, it outlines a phased plan for future commissioning and experimental campaigns, including control and instrumentation shakedown, heater characterization, model calibration, and testing at state points representative of EarthEn’s TES cycle.

25 ENERGY STORAGE↗

Single-Phase Dielectric Fluid Thermal Management for Power-Dense Automotive Power Electronics

This paper describes the design and performance of a dielectric fluid cooling concept for automotive power electronics. The concept combines a low-thermal-resistance package (which eliminates metalized ceramic substrates) with a high-performance convective cooling strategy (slot jets impinging on finned surfaces). Modeling was first used to design the cooling system to maximize thermal performance and minimize pumping power. Additionally, a prototype was then fabricated, and experiments were conducted to validate the model predictions using three fluids at various fluid flow rates (16.7 cm 3 /s [1 L/min] to 68.3 cm 3 /s [4.1 L/min]) and inlet temperatures (30C and 70C). The final design was compact (120-cm 3 total volume, including heat exchanger and conceptual power modules) and cooled 12 devices (e.g., silicon carbide). The validated model was then used to predict the junction-to-fluid thermal resistance and pumping power for various conditions, including 40C fluid temperature. The results predict thermal resistance values as low as 19 mm 2 K/W are possible using the dielectric fluid cooling approach. The dielectric fluid cooling system is predicted to provide thermal resistance and pumping power values that are approximately 56% and 90% lower, respectively, compared to an automotive power electronics cooling system.

30 DIRECT ENERGY CONVERSION↗

Optimizing aluminum oxide passivation layers—Laser-induced plasmas for layer formation, real-time diagnostics, and layer evaluation

Molten salts have beneficial thermophysical and electrolytic properties, but the aggressive nature of molten salts requires corrosion mitigation strategies, such as structural material surface treatment. In this study, a laser-induced breakdown spectroscopy (LIBS) system was used to simultaneously ablate and monitor de-excitation spectra to form and identify dense alumina passivation layers. Additionally, high-frequency (kHz) LIBS imaging allowed rapid elemental mapping and depth profiling for surface O/Al ratios via minimally destructive analysis. The paired t-test rejected the null hypothesis that air and Ar cover gases were equivalent (p = 0.00012), and the analysis of variance (ANOVA) showed that overlap between shots (p = 0.03112) and cover gas flow rate (p = 0.02729) were statistically significant. At 75% overlap, the O/Al ratio increased to 10% ± 2%, and 90% overlap rose further to 31% ± 2%. Regardless of overlap, tuning the gas flow from 0.5 to 2.5 L min −1 enhanced the O/Al ratio by 24% ± 4%. LIBS depth profiling showed a 3× increase in layer thickness from 75% to 90% overlap. Molecular band peaks of AlO, a precursor to Al 2 O 3 found during treatment, suggest that band intensity could be used for real-time Al 2 O 3 optimization. Ultimately, by adjusting laser spot overlap and cover gas flow rate, a 266 nm laser was shown to form Al 2 O 3 layers on an Al6061 sample, with cover gas flow rate primarily affecting the O/Al ratio and layer thickness correlating with laser spot overlap. This combination of molecular spectra collection, imaging, and depth profiling demonstrated the robust layer analysis capabilities of LIBS.

Corrosion↗

Structure and dynamics of magneto-inertial, differentially rotating laboratory plasmas

We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at ${\approx }100\,{\rm km}\,{\rm s}^{-1}$ . The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$ . As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value ${\approx }20\,{\rm km}\,{\rm s}^{-1}$ . We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements.

Valenzuela-Villaseca, V. (ORCID:0000000167721441)↗

Low-energy plasma observations at synchronous orbit

The University of California at San Diego Auroral Particles Experiment on the ATS 6 satellite in synchronous orbit has detected a low-energy plasma population which is separate and distinct from both the ring current and the plasma sheet populations. The density and temperature of this low-energy population are highly variable, with temperatures in the range kT = 1-30 eV and densities ranging from less than 1 per cu cm to more than 10 per cu cm. The occurrence of a dense low-energy plasma is most likely in the afternoon and dusk local time sectors, whereas n greater than 1 per cu cm is seen in the local night sector only during magnetically quiet periods. These observations suggest that this plasma is the outer zone of the plasmasphere. During magnetically active periods this low-energy plasma is often observed flowing sunward. In the dusk sector, strong sunward plasma flow is often observed for 1-2 hours prior to the onset of a substorm-associated particle injection.

Lennartsson, W.↗

Ambient Effects on Basalt and Rhyolite Lavas under Venusian, Subaerial, and Subaqueous Conditions

Both subaerial and subaqueous environments have been used as analog settings for Venus volcanism. To assess the merits of this, the effects of ambient conditions on the physical properties of lava on Venus, the seafloor, and land on Earth are evaluated. Rhyolites on Venus and on the surface of Earth solidify before basalts do because of their lower eruption temperatures. Rhyolite crust is thinner than basalt crust at times less than about an hour, especially on Venus. At later times, rhyolite crust is thicker because of its lower latent heat relative to basalt. The high pressure on the seafloor and Venus inhibits the exsolution of volatiles in lavas. Vesicularity and bulk density are proportional, so that lavas of the same composition should be more dense on the seafloor and less dense on land. Because viscosity depends partly upon the fraction of unvesiculated water in a melt, basalts with the same initial volatile abundance will be least viscous on the seafloor and most viscous on land. Assuming the same preeruptive H2O contents, molten rhyolites on Venus will have viscosities approx. 10% that of rhyolites on land. Despite lower expected viscosities, under-water flows are more buoyant and should have heights like subaerial and Venusian lavas of the same composition and extrusive history. In cases where the influence of crust is insignificant, a volume of rhyolite will have a higher aspect ratio than the same volume of basalt, no matter what the environment. If flow rheology is dominated by the presence of strong crust, aspect ratios differ little among environments or between compositions. These analyses support a rhyolitic interpretation for the composition of Venusian festooned flows and a basaltic interpretation for the composition of Venusian steep-sided domes. Although ambient effects are significant, extrusion rate and eruption history must also be considered to explain analogous volcanic landforms on Earth and Venus.

Bridges, Nathan T.↗

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗

Ocean Variability Beneath the Filchner-Ronne Ice Shelf Inferred From Basal Melt Rate Time Series

Fourteen phase-sensitive radars (ApRES) were deployed on the Filchner-Ronne Ice Shelf (FRIS) to measure variability in its basal melt rate. Melt rates from sites along the Ronne Depression vary seasonally, consistent with the dynamics of the propagation of seasonal dense water from the western ice front into the cavity. Several sites at the back of the FRIS cavity feature a signal with two seasonal maxima. Sub-ice shelf oceanographic data available from one of the sites indicate that this signal is caused by two different pathways followed by the same source water. Inter-annual variability is strongest along a direct flow pathway between western Ronne Ice Front and western Berkner Island. Highest melting occurred in 1999 and 2018, following anomalously low summer sea-ice concentrations in front of the ice shelf. Inter-annual melt rate variability at the back of the FRIS cavity is limited. If present, it is expressed as a suppression or delay in the arrival of the seasonal melt rate minimum, which can be understood in terms of inter-annual stratification changes and variable inflow pathways toward the sites. Long term mean ApRES melt rates agree with estimates from satellite data over eastern FRIS. However, the satellite estimates overstate the area of active basal freezing in the western part of the ice shelf. Furthermore, the temporal melt rate variability from the satellite estimates exaggerates the range of variability at both seasonal and inter-annual time scales with any correspondence between the in-situ and remotely derived inter-annual variability being limited to a single site.

54 ENVIRONMENTAL SCIENCES↗

Viscosity of Xenon Examined in Microgravity

Why does water flow faster than honey? The short answer, that honey has a greater viscosity, merely rephrases the question. The fundamental answer is that viscosity originates in the interactions between a fluid s molecules. These interactions are so complicated that, except for low-density gases, the viscosity of a fluid cannot be accurately predicted. Progress in understanding viscosity has been made by studying moderately dense gases and, more recently, fluids near the critical point. Modern theories predict a universal behavior for all pure fluids near the liquid-vapor critical point, and they relate the increase in viscosity to spontaneous fluctuations in density near this point. The Critical Viscosity of Xenon (CVX) experiment tested these theories with unprecedented precision when it flew aboard the Space Shuttle Discovery (STS-85) in August 1997. Near the critical point, xenon is a billion times more compressible than water, yet it has about the same density. Because the fluid is so "soft," it collapses under its own weight when exposed to the force of Earth s gravity - much like a very soft spring. Because the CVX experiment is conducted in microgravity, it achieves a very uniform fluid density even very close to the critical point. At the heart of the CVX experiment is a novel viscometer built around a small nickel screen. An oscillating electric field forces the screen to oscillate between pairs of electrodes. Viscosity, which dampens the oscillations, can be calculated by measuring the screen motion and the force applied to the screen. So that the fluid s delicate state near the critical point will not be disrupted, the screen oscillations are set to be both slow and small.

Zimmerli, Gregory A.↗

Multiwell Fiber Optic Sensing Reveals Effects of CO2 Flow on Triggered Seismicity

Induced seismicity is one of the main risks for gigaton-scale geological storage of carbon dioxide (⁠CO 2 ⁠). Thus, passive seismic monitoring is often recommended as a necessary component of the monitoring systems for CO 2 storage projects, with a particular forcus on risk mitigation. We present the first field study, CO2CRC Otway Project Stage 3 (Victoria, Australia), where distributed acoustic sensing (DAS) enabled high-precision tracking of the induced seismicity triggered by a small CO 2 injection and also informed the reservoir models. In 610 days of passive seismic monitoring of the Stage 3 injection, we detected 17 microseismic events (maximum moment magnitude Mw 0.1) using five deep boreholes equipped with enhanced-sensitivity optical fiber. The DAS array has sensitivity sufficient for detection and location of induced events with M w ~ ₋2 in a monitoring borehole located up to 1500 m away. Thanks to the dense spatial sampling by the DAS, we were able to estimate the focal mechanisms for events with M w > ₋1.5; although the monitoring boreholes provided very limited angular coverage. The main cluster of the events has the same location and source mechanism as the one triggered by the previous CO 2 injection at the Otway Project site, Stage 2C. Surprizingly, the Stage 2C and Stage 3 events closely followed the actual movement of the CO 2 saturation plume front (not the pressure front), as observed using controlled-source reflection seismic images. The nature of the plume-fault interaction remains unclear, but some alteration of the fault gouge by CO 2 might be responsible for the faults’ reactivation by the pressure perturbation. Importantly, the seismogenic fault could not be identified in the seismic images and was only revealed by DAS observations, which also demonstrated the signature of fluid–rock interaction, that may control the CO 2 flow.

58 GEOSCIENCES↗

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Enceladus and Europa: How Does Hydrothermal Activity Begin at the Surface?

The question of how the surface hydrothermal activity (e.g., eruptive plumes and heat flow) is initiated can be addressed within the frame-work of our "Perrier Ocean" model. This model delivers the necessary heat and chemicals to support the heat flow and plumes observed by Cassini in Enceladus' South Polar Region. The model employs closed-loop circulation of water from a sub-surface ocean. The ocean is the main reservoir of heat and chemicals, including dissolved gases. As ocean water moves up toward the surface, pressure is re-duced and gases exsolve forming bubbles. This bub-bly mixture is less dense than the icy crust and the buoyant ocean-water mixture rises toward the surface. Near the surface, heat and chemicals, including some volatiles, are delivered to the chambers in which plumes form and also to shallow reservoirs that keep the surface ice "warm". (Plume operations, per se, are as described by Schmidt et al. and Postberg et al. and are adopted by us.) After transferring heat, the water cools, bubbles contract and dissolve, and the mixture is now relatively dense. It descends through cracks in the crust and returns to the ocean. Once the closed-loop circulation has started it is self-sustaining. Loss of water via the erupting plumes is relatively negligible compared to the amount needed to maintain the heat flow.We note that the activity described herein for the the "Perrier-Ocean" model could, a priori, apply to all small icy bodies that sheltered an interior ocean at some point in their history.

Enceladus↗