SEARCH · Engineering Papers

Results for “rotating fluid”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

NASA NTRS2004

Numerical, Analytical, Experimental Study of Fluid Dynamic Forces in Seals Volume 6: Description of Scientific CFD Code SCISEAL

The objectives of the program were to develop computational fluid dynamics (CFD) codes and simpler industrial codes for analyzing and designing advanced seals for air-breathing and space propulsion engines. The CFD code SCISEAL is capable of producing full three-dimensional flow field information for a variety of cylindrical configurations. An implicit multidomain capability allow the division of complex flow domains to allow optimum use of computational cells. SCISEAL also has the unique capability to produce cross-coupled stiffness and damping coefficients for rotordynamic computations. The industrial codes consist of a series of separate stand-alone modules designed for expeditious parametric analyses and optimization of a wide variety of cylindrical and face seals. Coupled through a Knowledge-Based System (KBS) that provides a user-friendly Graphical User Interface (GUI), the industrial codes are PC based using an OS/2 operating system. These codes were designed to treat film seals where a clearance exists between the rotating and stationary components. Leakage is inhibited by surface roughness, small but stiff clearance films, and viscous pumping devices. The codes have demonstrated to be a valuable resource for seal development of future air-breathing and space propulsion engines.

Athavale, Mahesh↗

NASA NTRS2004

Numerical, Analytical, Experimental Study of Fluid Dynamic Forces in Seals: Description of Seal Dynamics Code DYSEAL and Labyrinth Seals Code KTK - Volume 5

Shapiro, Wilbur↗

DOE OSTI2025

Mechanical Resonant Sensing of Spin Texture Dynamics in a 2D Antiferromagnet

Abstract The coupling between the spin degrees of freedom and macroscopic mechanical motions, including striction, shearing, and rotation, has attracted wide interest with applications in actuation, transduction, and information processing. Experiments so far have established the mechanical responses to the long‐range ordered or isolated single spin states. However, it remains elusive whether mechanical motions can couple to a different type of magnetic structure, the non‐collinear spin textures, which exhibit nanoscale spatial variations of spin (domain walls, skyrmions,etc.) and are promising candidates to realize high‐speed computing devices. Here, collective spin texture dynamics is detected with nanoelectromechanical resonators fabricated from 2D antiferromagnetic (AFM) MnPS 3 with 10 −9 strain sensitivity. By examining radio frequency mechanical oscillations under magnetic fields, new magnetic transitions are identified with sharp dips in resonant frequency. They are attributed to collective AFM domain wall motions as supported by the analytical modeling of magnetostriction and large‐scale spin‐dynamics simulations. Additionally, an abnormally large modulation in the mechanical nonlinearity at the transition field infers a fluid‐like response due to ultrafast domain motion. The work establishes a strong coupling between spin texture and mechanical dynamics, laying the foundation for electromechanical manipulation of spin texture and developing quantum hybrid devices.

Chemistry↗

DOE OSTI2020

Characterizing magnetized plasmas with dynamic mode decomposition

Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Design optimization of a brush turbine with a cleaner/water based solution

Recently, a fluid turbine which has a brush attached to it has been designed and tested with water as fluid. The purpose of the turbine-brush is to clean up fouling in a tube. The Montreal Protocol prohibits the use of CFC products from refrigeration industry or from industry in general as a cleanser in 1996. Alternatives for the cleansers, devices or a combination of alternative devices with a cleanser should be found. One of the methods is to develop a device which cleans fouling with a cleaning medium. In this paper, we describe a turbine connected with a brush. However, the turbine with the brush should be simple and easy to install. This device is a combined small liquid turbine with a brush. The turbine is activated by the liquid flowing through the tube. Then the turbine turns the brush cleaning fouling along the tube. Based on the energy conservation and the Bernoulli equation along with an empirical relationship of drag force obtained from an experimental apparatus, a relationship of the rotational speed, the number of blades, and geometric variables of the turbine-brush was obtained. The predicted rotational speeds were compared with the experimental observations. Further work was recommended for improvements.

Kim, Rhyn H.↗

DOE OSTI2026

Numerical simulation of vortex-induced vibration response of a single IEA 10-MW wind turbine blade

Three-dimensional simulation of vortex-induced vibration (VIV) of a single International Energy Agency (IEA) 10-MW reference wind turbine blade with a length of 97.325 m is performed using the ExaWind stack, an open-source suite of codes. This study aims to illustrate the spanwise VIV response characteristics and cross-validate the results with an existing commercial framework. Five near-body meshes and three time steps are selected for the convergence study. To improve computational efficiency, several VIV triggering methods are also compared to shorten the VIV development period. The ExaWind-based VIV simulation strategy for a single IEA 10-MW blade is determined. First, the modal shape is validated against published results. Then, spanwise VIV responses of four blade configurations under a fixed and varied incoming flow velocity are analyzed. Results show that the VIV response is dominated by the first edgewise (second overall) mode. Little first-mode contributions appear near the second-mode node, producing a pi phase jump, and a higher harmonics response occurs near the blade root. Rotational degrees of freedom are minor compared with translational motion. The response versus reduced velocity is analyzed, showing a two-branch behavior similar to that of VIV for a bluff cylinder. Across all tested cases, the dominant frequency remains locked to the natural frequency of the second mode with no observed desynchronization. A mild deviation is observed for the case of 90-degree pitch and 310-degree azimuth rotation near a reduced velocity of 6, which will be examined with additional cases in future work. These findings indicate that severe VIV responses can arise under specific configurations and flow conditions, thereby increasing the potential for VIV fatigue damage and requiring greater attention during operation.

17 WIND ENERGY↗

NASA NTRSDate not supplied

Differential Gene Expression in A Cross-Feeding Two-Species Model Microbial Community Under Simulated Microgravity and Deep-Space Radiation

A long-term goal of space biology is to understand interspecies microbial interactions in space. Presently, little is known about the combined effect of microgravity and ionizing radiation on bacterial community response when species are interdependent through exchange of metabolites in fluid medium (cross-feeding). Microgravity is expected to slow interspecies mass transfer and growth in cross-feeding communities in the low-shear, diffusion-limited environment, while ionizing radiation may influence stress response to direct (DNA damage) and indirect damage (ROS). Using a well-understood, two-species (Escherichia coli and Salmonella enterica) microbial community engineered to be a model for studying cross-feeding, we simulated galactic cosmic rays (GCRsim) and microgravity to test the hypothesis: exposure to ionizing radiation causes cell damage or stress, altering transcriptomic community responses in metabolically interdependent cells, which is exacerbated by microgravity. We expect to see differential gene expression between cross-feeding and non-cross-feeding communities. We measured GCRsim effects on growth and gene expression in well-mixed versus simulated-microgravity conditions and in cross-feeding and non-cross-feeding medium. Microbial cultures were inoculated into liquid medium in rotating wall vessels (RWV) with different rotation rates: 5 RPM (simulated microgravity) and 50 RPM (well-mixed). The E. coli-S. enterica consortium, under simulated microgravity, were exposed to 500 mGy of Simplified 5-ion Galactic Cosmic Ray Simulation for 2 hours at Brookhaven National Lab. We harvested samples 40 minutes after irradiation for extraction and sequencing (NASA GeneLab). Here we present the differential gene expression analysis results, which reveal altered transcriptomic community responses, even where growth rate differences are not observed. Gene expression of these actively metabolizing microbial communities in GCRsim may illuminate molecular mechanisms of microbial interactions in space. Understanding how microbial community gene expression, metabolism, and other cellular processes are influenced by spaceflight stressors can inform the use of microbes in human life support for low Earth orbit missions and beyond.

microgravity↗

NASA NTRS1993

Supersonic shock wave/vortex interaction

Although shock wave/vortex interaction is a basic and important fluid dynamics problem, very little research has been conducted on this topic. Therefore, a detailed experimental study of the interaction between a supersonic streamwise turbulent vortex and a shock wave was carried out at the Penn State Gas Dynamics Laboratory. A vortex is produced by replaceable swirl vanes located upstream of the throat of various converging-diverging nozzles. The supersonic vortex is then injected into either a coflowing supersonic stream or ambient air. The structure of the isolated vortex is investigated in a supersonic wind tunnel using miniature, fast-response, five-hole and total temperature probes and in a free jet using laser Doppler velocimetry. The cases tested have unit Reynolds numbers in excess of 25 million per meter, axial Mach numbers ranging from 2.5 to 4.0, and peak tangential Mach numbers from 0 (i.e., a pure jet) to about 0.7. The results show that the typical supersonic wake-like vortex consists of a non-isentropic, rotational core, where the reduced circulation distribution is self similar, and an outer isentropic, irrotational region. The vortex core is also a region of significant turbulent fluctuations. Radial profiles of turbulent kinetic energy and axial-tangential Reynolds stress are presented. The interactions between the vortex and both oblique and normal shock waves are investigated using nonintrusive optical diagnostics (i.e. schlieren, planar laser scattering, and laser Doppler velocimetry). Of the various types, two Mach 2.5 overexpanded-nozzle Mach disc interactions are examined in detail. Below a certain vortex strength, a 'weak' interaction exists in which the normal shock is perturbed locally into an unsteady 'bubble' shock near the vortex axis, but vortex breakdown (i.e., a stagnation point) does not occur. For stronger vortices, a random unsteady 'strong' interaction results that causes vortex breakdown. The vortex core reforms downstream of the rear stagnation point, and the reduced circulation distribution once again becomes self-similar in this region. A-new model of this interaction is proposed. Finally, a curve defining the approximate limits of supersonic vortex breakdown is presented.

Settles, G. S.↗

NASA NTRS1993

Foil bearings

The rolling element bearings (REB's) which support many turbomachinery rotors offer high load capacity, low power requirements, and durability. Two disadvantages of REB's are: (1) rolling or sliding contact within the bearing has life-limiting consequences; and (2) REB's provide essentially no damping. The REB's in the Space Shuttle Main Engine (SSME) turbopumps must sustain high static and dynamic loads, at high speeds, with a cryogenic fluid as lubricant and coolant. The pump end ball bearings limit the life of the SSME high pressure oxygen turbopump (HPOTP). Compliant foil bearing (CFB) manufacturers have proposed replacing turbopump REB's with CFB's CFB's work well in aircraft air cycle machines, auxiliary power units, and refrigeration compressors. In a CFB, the rotor only contracts the foil support structure during start up and shut down. CFB damping is higher than REB damping. However, the load capacity of the CFB is low, compared to a REB. Furthermore, little stiffness and damping data exists for the CFB. A rotordynamic analysis for turbomachinery critical speeds and stability requires the input of bearing stiffness and damping coefficients. The two basic types of CFB are the tension-dominated bearing and the bending-dominated bearing. Many investigators have analyzed and measured characteristics of tension-dominated foil bearings, which are applied principally in magnetic tape recording. The bending-dominated CFB is used more in rotating machinery. This report describes the first phase of a structural analysis of a bending-dominated, multileaf CFB. A brief discussion of CFB literature is followed by a description and results of the present analysis.

Elrod, David A.↗

NASA NTRS1996

Large deflection analysis of a tension-foil bearing

The rolling element bearings (REB's) which support many turbomachinery rotors offer high load capacity, low power requirements, and durability. Two disadvantages of REB's are as follows: rolling or sliding contact within the bearing has life-limiting consequences; and REB's provide essentially no damping. The REB's in the Space Shuttle Main Engine (SSME) turbopumps must sustain high static and dynamic loads, at high speeds, with a cryogenic fluid as lubricant and coolant. The pump end ball bearings limit the life of the SSME high pressure oxygen turbopump (HPOTP). Compliant foil bearing (CFB) manufacturers have proposed replacing turbopump REB's with CFB's. CFB's work well in aircraft air cycle machines, auxiliary power units, and refrigeration compressors. In a CFB, the rotor only contacts the foil support structure during start up and shut down. CFB damping is higher than REB damping. However, the load capacity of the CFB is low, compared to a REB. Furthermore, little stiffness and damping data exist for the CFB. A rotordynamic analysis for turbomachinery critical speeds and stability requires the input of bearing stiffness and damping coefficients. The two basic types of CFB are the tension-dominated bearing and the bending-dominated bearing. Many investigators have analyzed and measured characteristics of tension-dominated foil bearings, which are applied principally in magnetic tape recording. The bending-dominated CFB is used more in rotating machinery. Recently, a new tension-foil bearing configuration has been proposed for turbomachinery applications.

Elrod, David A.↗

NASA NTRS1993

Navier-Stokes analysis of an oxidizer turbine blade with tip clearance with and without a mini-shroud

The Gas Generator Oxidizer Turbine (GGOT) Blade is being analyzed by various investigators under the NASA MSFC-sponsored Turbine Stage Technology Team design effort. The present work concentrates on the tip clearance region flow and associated losses; however, flow details for the passage region are also obtained in the simulations. The present calculations simulate the rotor blade row in a rotating reference frame with the appropriate coriolis and centrifugal acceleration term included in the momentum equations. The upstream computational boundary is located about one axial chord from the blade leading edge. The boundary conditions at this location have been determined by Pratt & Whitney using an Euler analysis without the vanes to obtain approximately the same flow profiles at the rotor as were obtained with the Euler stage analysis including the vanes. Inflow boundary layer profiles are then constructed assuming the skin friction coefficient at both the hub and the casing. The downstream computational boundary is located about one axial chord from the blade trailing edge, and the circumferentially averaged static pressure at this location was also obtained from the P&W Euler analysis. Results obtained for the 3-D baseline GGOT geometry at the full scale design Reynolds number show a region of high loss in the region near the casing. Particle traces in the near tip region show vortical flow behavior of the fluid which passes through the clearance region and exits at the downstream edge of the gap. In an effort to reduce clearance flow losses, the mini-shroud concept was proposed by the Pratt & Whitney design team. Calculations were performed on the GGO geometry with the mini-shroud. Results of these calculations indicate that the mini-shroud does not significantly affect the flow in the passage region, and although the tip clearance flow is different, the mini-shroud does not seem to prevent the above-mentioned vortical flow behavior. Since both flow distortion and total pressure losses are similar for both geometries, the addition of the mini-shroud does not seem to reduce the tip clearance flow effects.

Chan, Tony↗

NASA NTRS1993

The Influence of Swirl Brakes and a Tip Discharge Orifice on the Rotordynamic Forces Generated by Discharge-to-Suction Leakage Flows in Shrouded Centrifugal Pumps

Recent experiments conducted in the Rotor Force Test Facility at the California Institute of Technology have examined the effects of a tip leakage restriction and swirl brakes on the rotordynamic forces due to leakage flows on an impeller undergoing a prescribed circular whirl. The experiments simulate the leakage flow conditions and geometry of the Alternate Turbopump Design (ATD) of the Space Shuttle High Pressure Oxygen Turbopump and are critical to evaluating the pump's rotordynamic instability problems. Previous experimental and analytical results have shown that discharge-to-suction leakage flows in the annulus of a shrouded centrifugal pump contribute substantially to the fluid induced rotordynamic forces. Also, previous experiments have shown that leakage inlet (pump discharge) swirl can increase the cross-coupled stiffness coefficient and hence increase the range of positive whirl for which the tangential force is destabilizing. In recent experimental work, the present authors demonstrated that when the swirl velocity within the leakage path is reduced by the introduction of ribs or swirl brakes, then a substantial decrease in both the destabilizing normal and tangential forces could be achieved. Motivation for the present research is that previous experiments have shown that restrictions such as wear rings or orifices at pump inlets affect the leakage forces. Recent pump designs such as the Space Shuttle Alternate Turbopump Design (ATD) utilize tip orifices at discharge for the purpose of establishing axial thrust balance. The ATD has experienced rotordynamic instability problems and one may surmise that these tip discharge orifices may also have an important effect on the normal and tangential forces in the plane of impeller rotation. The present study determines if such tip leakage restrictions contribute to undesirable rotordynamic forces. Additional motivation for the present study is that the widening of the leakage path annular clearance and the installation of swirl brakes in the ATD has been proposed to solve its instability problems. The present study assesses the effect of such a design modification on the rotordynamic forces. The experimental apparatus consists of a solid or dummy impeller, a housing instrumented for pressure measurements, a rotating dynamometer and an eccentric whirl mechanism. The solid impeller is used so that leakage flow contributions to the forces are measured, but the main throughflow contributions are not experienced. The inner surface of the housing has been modified to accommodate meridional ribs or swirl brakes within the leakage annulus. In addition, the housing has been modified to accommodate a discharge orifice that qualitatively simulates one side of the balance piston orifice of the Space Shuttle ATD. Results indicate the detrimental effects of a discharge orifice and the beneficial effects of brakes. Plots of the tangential and normal forces versus whirl ratio show a substantial increase in these forces along with destabilizing resonances at some positive whirl ratios when a discharge orifice is added. When brakes are added, some of the detrimental effects of the orifice are reduced. For the tangential force, a plot versus whirl ratio shows a significant reduction and a destabilizing resonance appears to be eliminated. For the normal force, although the overall force is not reduced, again a destabilizing resonance appears to be eliminated.

Sivo, Joseph M.↗

NASA NTRSDate not supplied

GT2024-139262: Development of a Computational Fluid Dynamics Model for a High-Speed Centrifugal Compressor

Computational fluid dynamics (CFD) has become widely used in the design and analysis of turbomachinery components such as centrifugal compressors. However, CFD is only a limited representation of experimental cases and struggles to model complex flows or lack model details to increase computation speed. The High Efficiency Centrifugal Compressor (HECC) was designed by United Technologies Research Center for NASA to investigate the difficulties behind improving centrifugal compressor technology and to provide an open case to the turbomachinery research community. The CFD technology seen in previous numerical models of the compressor have improved since then. In this study, a high-fidelity computational model for the vaneless configuration of HECC operating at its design point was developed and compared to an equivalent experimental case. This model incorporated GPU versions of mesh generation code WAND and solver code LEO from AeroDynamic Solutions, Inc. Blade sections extracted from the solid model formed the basis of the impeller mesh, while the remainder of the stage was modeled as ducts. Standard day aerodynamic conditions, design rotational speed, design mass flow rate, and inlet flow angles served as additional inputs. A mesh sensitivity study was performed prior to using the model for analysis. The accuracy of the model’s predictions of impeller and stage performance were assessed by comparing the total pressure ratio, total temperature rise ratio, and adiabatic efficiency to the corresponding parameters observed in the experimental case. Results showed a general overprediction of compressor performance by the model. Predictions were more accurate for the impeller performance compared to those for the full compressor stage. Spanwise analyses showed that the model tended to follow flow behavior patterns exhibited in the experiment. Beyond this study, this model was used to examine HECC aerodynamics and performance in research following thereafter.

centrifugal compressor↗

DOE OSTI2024

Corrosion-Resistant Non-Carbon Electrocatalyst Supports for PEFCS

Insufficient durability of commercial catalysts is a persistent issue for PEFC applications.1 Herein, Pt on Nb-doped-TiO 2 is reported as a highly durable electrocatalyst on which the direct 4-electron reduction of oxygen to water is more facile compared to Pt/C. All performance metrics are reported comparing Pt/Nb-TiO2 with a 15% Pt loading, against a commercial Pt/C catalyst (46.5% Pt loading, Tanaka, K. K.). Nb-doped-TiO 2 with high surface area and high electronic conductivity was synthesized using the supercritical fluid method. Initially, the durability of the catalyst was characterized using accelerated stability tests (ASTs) involving 10,000 high potential cycles (DOE/FCCJ protocol) and the Pt/Nb-TiO 2 was found to retain 78% of its initial electrochemically active surface area (ECSA) compared to the 57.6 % retained by Pt/C. These observations were in excellent agreement with previous reports that the Pt particle size of Pt/C grew from 2nm to 8nm during the AST protocol along with severe corrosion and amorphization of the carbon surface.2, 3 TEM and XPS studies of the Pt-Nb-doped-TiO 2 catalyst showed that the Pt particle size grew from 3nm to 6nm and the Nb(IV) and Ti(III) in the support were oxidized to Nb(V) and Ti(IV) after the durability test. Thus, the improvement in Pt/Nb-TiO 2 ECSA retention was attributed to the lower extent of particle growth and lack of oxidative support loss upon oxidation as compared to Pt/C. The oxygen reduction reaction (ORR) performance was characterized by linear polarization using a rotating disk electrode (RDE). The electrochemical surface areas of Pt/Nb-TiO 2 and Pt/C were found to be 48m2·g -1 and 83m2·g -1 respectively, and the mass activity for the ORR at 0.9V vs. RHE were found to be 150 mA·mg-1Pt and 124 mA·mg-1Pt respectively. The improved mass activity on Pt/Nb-TiO 2 was attributed to strong metal support interaction (SMSI) between the Nb-TiO 2 support and the Pt catalyst based on the 625 meV decrease in the binding energy of the Pt4f x-ray photo-electron spectroscopy (XPS) peaks of Pt/Nb-TiO 2 compared to Pt/C. To quantify the impact of the SMSI, a kinetic model was applied to calculate the elementary reaction rate constants for the various steps of the ORR on both catalysts. The reaction rate constant (k1) for the direct 4-electron transfer pathway to produce H 2 O was significantly larger in Pt/Nb-TiO 2 as compared to Pt/C. Thus, the reduction in the electron binding energy as observed in the XPS was found to aid the facile filling of the higher energy 2p orbitals of O 2 thereby leading to improved 4-electron transfer kinetics and improved overall activity.

08 HYDROGEN↗

NASA NTRS2011

Computational Fluid Dynamics Demonstration of Rigid Bodies in Motion

The Design Analysis Branch (NE-Ml) at the Kennedy Space Center has not had the ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom (6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations were based on a conceptual Space Launch System (SLS). The first simulation that was prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from its core stage. To reduce computational time during the development of the simulation, only half of the physical domain with respect to the symmetry plane was simulated. Then a full solution was prepared and computed. The second simulation was a model of the SLS as it departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The simulation was then constrained to only allow translations.

Camarena, Ernesto↗

NASA NTRS2002

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗

DOE OSTI2022

Confinement-Driven Heterogeneous Benzene Crystallization in Silica Nanopores

Nanoconfinement alters the thermodynamics, dynamics, and kinetics of fluids hosted in nanoscale solid nanopores to an extent that depends on the characteristics of the confining space and the chemistry of the confined fluids. Confinement-induced alterations in the phase behavior of confined energetic fluids under high pressure or low temperature are highly relevant to subsurface and subsea phenomena such as fluid flow in porous media, hydrate formation and dissociation, and gas storage capacity. Although extensive efforts have been directed toward understanding the phase behavior of confined fluids, the role of solid-liquid interfaces in the phase transitions of organic liquids has not been resolved yet. Here, we explore the onset and growth of benzene crystallization confined in 6 nm sized SBA-15 silica nanopores in the temperature range 300-200 K using in situ extended range small-angle and wide-angle X-ray scattering (SAXS/WAXS) measurements and atomistic classical molecular dynamics (MD) simulations. The crystallization onset of confined benzene depresses to 265 K compared to the freezing point of bulk benzene ( ~278 K), followed by the continuous growth of the emerged crystals in the pore space with complete crystallization at 200 K. The orientation of the emerged benzene crystals is dominated by parallel (π-π stacking) and perpendicular (T-shape stacking) orientations along the cylindrical pore radius and pore length, respectively. The onset of benzene crystals occurs heterogeneously on the pore surface and grows continuously toward the pore center. Further, confined benzene undergoes a dynamical crossover from fragile to strong dynamics behavior, inferred from the rotational and translational diffusion. The insights provided by this study have significant implications for the phase transitions of confined organic liquids that are relevant to a wide range of applications in the biological, geological, environmental, and chemical fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗