Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Critical Fluids”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Validation of Transient Spacecraft Refueling Model with Gateway Breadboard Test Data

The ability to refuel spacecraft on-orbit is a key enabling technology for deep space human exploration. On-orbit refueling requires priming a pressurized liquid propellant from a source tank into fluid transfer lines at vacuum conditions, which could result in excessive surge pressures beyond the system’s qualified operating limits and potentially damage hardware. Another source of potentially damaging surge pressures is the water hammer transient that occurs when the isolation valves are closed to pause or end the refueling operation. Therefore, it is imperative to employ numerical models during the development of flight hardware to parametrically characterize system performance over a wide range of operating conditions to minimize the risk of exceeding pressure limits during flight. Since rapid pressure transients can be difficult to model accurately, it is important to validate those models with test data. As part of the Gateway bipropellant refueling system development, a simplified fluidic breadboard system was created using water as a simulant for fuel to support early program risk reduction of refueling operations. The objectives of breadboard testing were to gather performance data to characterize and demonstrate critical refueling operations and to validate numerical models that can be extended to predicting flight system performance. A description of the breadboard test article and the results of numerical model validation is presented. The model is shown to have very good agreement with test data.

Refueling↗
DOE OSTI2025

Learning and discovering multiple solutions using physics-informed neural networks with random initialization and deep ensemble

In this work we explore the capability of physics-informed neural networks (PINNs) to discover multiple solutions. Many real-world phenomena governed by nonlinear differential equations (DEs), such as fluid flow, exhibit multiple solutions under the same conditions, yet capturing this solution multiplicity remains a significant challenge. A key difficulty lies in providing appropriate initial conditions or guesses, as widely used time-marching schemes and Newton’s method are highly sensitive to these choices when solving complex computational problems. While machine learning models, particularly PINNs, have shown promise in solving DEs, their ability to capture multiple solutions remains underexplored. In this work, we propose a simple and practical approach using PINNs to learn and discover multiple solutions. We first demonstrate that PINNs, when combined with random initialization and deep ensemble method—originally developed for uncertainty quantification—can effectively uncover multiple solutions to nonlinear ordinary and partial DEs. Although training large ensembles of PINNs may appear computationally demanding, this can be done efficiently using vectorization techniques supported by modern deep learning frameworks, allowing many networks to be trained simultaneously. Our approach highlights the critical role of initialization in shaping solution diversity, addressing an often-overlooked aspect of machine learning for scientific computing. Furthermore, we propose utilizing PINN-generated solutions as initial conditions or initial guesses for conventional numerical solvers to enhance accuracy and efficiency in capturing multiple solutions. Extensive numerical experiments, including the Allen–Cahn equation and cavity flow, where our approach successfully identifies both stable and unstable solutions, validate the effectiveness of our method. These findings establish a general and efficient framework for addressing solution multiplicity in nonlinear DEs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Perturbations of the Richardson number field by gravity waves

An analytic solution is presented for a stratified fluid of arbitrary constant Richardson number. By computer aided analysis the perturbation fields, including that of the Richardson number can be calculated. The results of the linear analytic model were compared with nonlinear simulations, leading to the following conclusions: (1) the perturbations in the Richardson number field, when small, are produced primarily by the perturbations of the shear; (2) perturbations of in the Richardson number field, even when small, are not symmetric, the increase being significantly larger than the decrease (the linear analytic solution and the nonlinear simulations both confirm this result); (3) as the perturbations grow, this asymmetry increases, but more so in the nonlinear simulations than in the linear analysis; (4) for large perturbations of the shear flow, the static stability, as represented by N2, is the dominating mechanism, becoming zero or negative, and producing convective overturning; and (5) the convectional measure of linearity in lee wave theory, NH/U, is no longer the critical parameter (it is suggested that (H/u sub 0) (du sub 0/dz) takes on this role in a shearing flow).

Wurtele, M. G.↗
NASA NTRS1993

Microwave properties and characterization of co-evaporated BSCCO thin films

An extensive characterization of Bi-Sr-Ca-Cu-O (BSCCO) thin films deposited by co-evaporation on LaAlO3 and SrTiO3 substrates was performed. The films had a T(sub c) (R = O) of approximately 78 K, and were predominantly c-axis oriented, with critical current densities (J(sub c)) at 4.5 K of 1.6 x 10(exp 6) and 1.1 x 10(exp 6) A cm(sup -2), for the samples on SrTiO3 and LaAlO3, respectively. The microwave properties of the films were examined by three techniques. The complex conductivity sigma(sub *) = sigma(sub 1) - j(sigma(sub 2)) and the magnetic penetration depth (A) were measured by power transmission at 30.6 GHz; the surface resistance (R(sub s)) was measured using a cavity resonator at 58.9 GHz, and the transmission line losses were determined by measuring the quality factor (Q) of a linear microstrip resonator at 10.4 and 20.2 GHz. The complex conductivity for the film on LaAlO3 was determined to be (2.0-j10) x 10(exp 5) S/m at 77 K. It was observed that in the superconducting state sigma(sub 1) deviates from both the Bardeen-Cooper-Schrieffer (BCS) theory and the two-fluid model. Values of lambda were found to be approximately 2.0 and 1.1 microns at 77 K and 20 K respectively, and were obtained for the film on LaAlO3. The value of lambda at 20 K was approximately three times larger than that of BSCCO single crystals. R(sub s) values of 865 and 1391 mOmega were obtained for the films on SrTiO3 and LaAlO3, respectively, at 77 K and 58.9 GHz. Unloaded Q factors at 20 K of approximately 1100 and 800 at 10.4 and 20.2 GHz respectively, were measured for the BSCCO resonator. Unloaded Q values of 290 and 405 at 20 K were obtained at 10.4 GHz and 20.2 GHz respectively, for an all gold (Au) resonator.

Miranda, F. A.↗
NASA NTRS2002

Aspects of Supercritical Turbulence: Direct Numerical Simulation of O2/H2 and C7H16/N2 Temporal Mixing Layers

Results from Direct Numerical Simulations of temporal, supercritical mixing layers for two species systems are analyzed to elucidate species-specific turbulence aspects. The two species systems, O2/H2 and C7HG16/N2, have different thermodynamic characteristics; thus, although the simulations are performed at similar reduced pressure (ratio of the pressure to the critical pressure), the former system is dose to mixture ideality and has a relatively high solubility with respect to the latter, which exhibits strong departures from mixture ideality Due to the specified, smaller initial density stratification, the C7H16/N2 layers display higher growth and increased global molecular mixing as well as larger turbulence levels. However, smaller density gradients at the transitional state for the O2/H2 system indicate that on a local basis, the layer exhibits an enhanced mixing, this being attributed to the increased solubility and to mixture ideality. These thermodynamic features are shown to affect the irreversible entropy production (i.e. the dissipation), which is larger for the O2/H2 layer and is primarily concentrated in high density-gradient magnitude regions that are distortions of the initial density stratification boundary. In contrast, the regions of largest dissipation in the C7H16/N2 layer are located in high density-gradient magnitude regions resulting from the mixing of the two fluids.

Bellan, J.↗
DOE OSTI2024

A versatile pressure-cell design for studying ultrafast molecular-dynamics in supercritical fluids using coherent multi-pulse x-ray scattering

Supercritical fluids (SCFs) can be found in a variety of environmental and industrial processes. They exhibit an anomalous thermodynamic behavior, which originates from their fluctuating heterogeneous micro-structure. Characterizing the dynamics of these fluids at high temperature and high pressure with nanometer spatial and picosecond temporal resolution has been very challenging. The advent of hard x-ray free electron lasers has enabled the development of novel multi-pulse ultrafast x-ray scattering techniques, such as x-ray photon correlation spectroscopy (XPCS) and x-ray pump x-ray probe (XPXP). These techniques offer new opportunities for resolving the ultrafast microscopic behavior in SCFs at unprecedented spatiotemporal resolution, unraveling the dynamics of their micro-structure. However, harnessing these capabilities requires a bespoke high-pressure and high-temperature sample system that is optimized to maximize signal intensity and address instrument-specific challenges, such as drift in beamline components, x-ray scattering background, and multi-x-ray-beam overlap. We present a pressure cell compatible with a wide range of SCFs with built-in optical access for XPCS and XPXP and discuss critical aspects of the pressure cell design, with a particular focus on the design optimization for XPCS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A constitutive model for sheared dense suspensions of rough particles

In a concentrated suspension, particles come into contact due to the presence of asperities on their surfaces. As a result, the contact forces and interparticle friction become one of the important factors governing the rheology of rigid particle suspensions at high concentrations. We show that a load-dependent friction model can be used to reproduce the experimentally observed shear thickening [ST—continuous and discontinuous (DST)] behavior with an increase in the asperity size. Increasing the particle surface roughness size leads to a decrease in the critical shear rate above which shear-thickening takes place, a reduction in the critical volume fraction for DST and an increase in the viscosity jump across non-shear-thickened to shear-thickened regime. In this paper, we propose a constitutive model to quantify the effect of increasing the roughness size on the rheology of dense ST suspensions as well as on the critical shear rate for ST and the critical volume fraction for DST. We fit this model to our simulation data for stress controlled shear flow of dense rough particle suspensions. Once the fitting is complete, these equations are used to predict exact volume fractions and shear stress values for transitions between three regimes on the shear stress-shear rate flow state diagram for different roughness values. Furthermore, the results of this study can be used to tune the particle surface roughness for manipulating the dense suspension rheology according to different applications.

42 ENGINEERING↗
NASA NTRS1998

Numerical Modeling of Physical Vapor Transport in Contactless Crystal Growth Geometry

Growth from the vapor under conditions of limited contact with the walls of the growth ampoule is beneficial for the quality of the growing crystal due to reduced stress and contamination which may be caused by interactions with the growth container. The technique may be of a particular interest for studies on crystal growth under microgravity conditions: elimination of some factors affecting the crystal quality may make interpretation of space-conducted processes more conclusive and meaningful. For that reason, and as a part of our continuing studies on 'contactless' growth technique, we have developed a computational model of crystal growth process in such system. The theoretical model was built, and simulations were performed using the commercial computational fluid dynamics code, (CFD) ACE. The code uses an implicit finite volume formulation with a gray discrete ordinate method radiation model which accounts for the diffuse absorption and reflection of radiation throughout the furnace. The three-dimensional model computes the heat transfer through the crystal, quartz, and gas both inside and outside the ampoule, and mass transport from the source to the crystal and the sink. The heat transport mechanisms by conduction, natural convection, and radiation, and mass transport by diffusion and convection are modeled simultaneously and include the heat of the phase transition at the solid-vapor interfaces. As the thermal boundary condition, temperature profile along the walls of the furnace is used. For different thermal profiles and furnace and ampoule dimensions, the crystal growth rate and development of the crystal-vapor and source-vapor interfaces (change of the interface shape and location with time) are obtained. Super/under-saturation in the ampoule is determined and critical factors determining the 'contactless' growth conditions are identified and discussed. The relative importance of the ampoule dimensions and geometry, the furnace dimensions and its temperature, and the properties of the grown material are analyzed. The results of the simulations are compared with related experimental results on growth of CdTe, CdZnTe, ZnTe, PbTe, and PbSnTe crystals by this technique.

Palosz, W.↗
DOE OSTI2021

Techno-Economic Analysis of Lithium Extraction from Geothermal Brines

The United States has a large, domestic source of lithium in geothermal fluids, especially at the Salton Sea region of southern California, where estimates of lithium pass-through at geothermal plants exceed 24,000 metric tons per year, based on 2019 geothermal plant operations. Lithium extraction from geothermal brines offers the potential to provide the United States with a secure, domestic supply of lithium to meet the increasing demands of electric vehicles, grid energy storage, portable electronics, and other end-use applications. Additionally, the use of direct extraction technologies allows for a more sustainable lithium supply relative to current evaporative brine and hardrock mining operations in terms of land use, water use, time to market with lithium products, and carbon intensity of operations. This report is part of an effort to assess geothermal brines as a source of commercial lithium supply for the United States. In this study, the National Renewable Energy Laboratory (NREL) reviews and summarizes public technoeconomic analyses of lithium extraction technologies. The work was coordinated with the Critical Minerals Institute at the Colorado School of Mines who focused on supply chain analysis of lithium.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Nanosecond Repetitively Pulsed (NRP) Plasmas: Relationship Between Induced Flow and Plasma Characteristics at Atmospheric Pressure (Final Technical Report)

The local flow fields induced by plasma discharges at atmospheric pressure can play a critical role in aerodynamic flow and combustion control and have significant implications for plasma-based technologies in medicine and environmental engineering. However, the coupling between the plasma characteristics, the geometry and the induced flow field is not well understood. The primary goal of this work is to determine the relationship between the plasma characteristics and the induced flow field for nanosecond pulsed discharges at atmospheric pressure. A comprehensive experimental investigation of plasmas produced by nanosecond high-voltage pulses in a pin-to-pin electrode configuration was conducted, employing several measurement techniques including microwave and laser Rayleigh scattering, optical emission spectroscopy, and high-speed imaging. The flow field induced by the nanosecond plasma was measured and characterized using measurements of the velocity (particle image velocimetry) and density (background oriented schlieren) with unprecedented spatiotemporal resolution. The experimental measurements were complemented by development of a one-dimensional model of the plasma discharge and a high-fidelity computational fluid dynamics (CFD) simulation of the subsequent plasma-induced flow field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Computational Fluid Dynamics (CFD) Simulation of Hypersonic Turbine-Based Combined-Cycle (TBCC) Inlet Mode Transition

Methods of computational fluid dynamics were applied to simulate the aerodynamics within the turbine flowpath of a turbine-based combined-cycle propulsion system during inlet mode transition at Mach 4. Inlet mode transition involved the rotation of a splitter cowl to close the turbine flowpath to allow the full operation of a parallel dual-mode ramjet/scramjet flowpath. Steady-state simulations were performed at splitter cowl positions of 0deg, -2deg, -4deg, and -5.7deg, at which the turbine flowpath was closed half way. The simulations satisfied one objective of providing a greater understanding of the flow during inlet mode transition. Comparisons of the simulation results with wind-tunnel test data addressed another objective of assessing the applicability of the simulation methods for simulating inlet mode transition. The simulations showed that inlet mode transition could occur in a stable manner and that accurate modeling of the interactions among the shock waves, boundary layers, and porous bleed regions was critical for evaluating the inlet static and total pressures, bleed flow rates, and bleed plenum pressures. The simulations compared well with some of the wind-tunnel data, but uncertainties in both the windtunnel data and simulations prevented a formal evaluation of the accuracy of the simulation methods.

Slater, John W.↗
NASA NTRS2011

The Effects of Microgravity on Seated Height (Spinal Elongation)

ABSTRACT Many physiological factors, such as spinal elongation, fluid shifts, bone atrophy, and muscle loss, occur during an exposure to a microgravity environment. Spinal elongation is just one of the factors that can also affect the safety and performance of a crewmember while in space. Spinal elongation occurs due to the lack of gravity/compression on the spinal column. This allows for the straightening of the natural spinal curve. There is a possible fluid shift in the inter-vertebral disks that may also result in changes in height. This study aims at collecting the overall change in seated height for crewmembers exposed to a microgravity environment. During previous Programs, Apollo-Soyuz Test Project (ASTP) and Skylab, spinal elongation data was collected from a small number of subjects in a standing posture but were limited in scope. Data from these studies indicated a quick increase in stature during the first few days of weightlessness, after which stature growth reached a plateau resulting in up to a 3% increase of the original measurement [1-5]. However, this data was collected only for crewmembers in standing posture and not in a seated posture. Seated height may have a different effect than standing height due to a change in posture as well as due to a compounded effect of wearing restraints and a potential compression of the gluteal area. Seated height was deemed as a critical measurement in the design of the Constellation Program s (CxP) Crew Exploration Vehicle (CEV), called Orion which is now the point-of-departure vehicle for the Multi-Purpose Crew Vehicle (MPCV) Program; therefore a better understanding of the effects of microgravity on seated height is necessary. Potential changes in seated height that may not have impacted crew accommodation in previous Programs will have significant effects on crew accommodation due to the layout of seats in the Orion.. The current and existing configuration is such that the four crewmembers are stacked two by two with the commander and pilot seats on the top and the two remaining seats underneath, thereby limiting the amount of clearance for the crewmembers seated in the bottom seat. The inner mold line of these types of vehicles are fixed due to other design constraints; therefore, it is essential that all seats incorporate additional clearance to account for adequate spinal growth thereby ensuring that the crew can safely ingress the seat and be strapped in prior to its return to earth. If there is not enough clearance to account for spinal growth deltas between seats then there is the potential that crewmembers will not be able to comfortably and safely fit into their seats. The crewmember in the bottom stacked seat may even have negative clearance with the seat above him or her which could lead to potential ingress/egress issues or potentially injury of the crewmember during landing. These impacts are specific to these types of vehicles with stacked seat configuration. Without proper knowledge of the amount of spinal elongation, or growth, which occurs due to microgravity and space flight, the design of future vehicle(s) or suits may cause injury, discomfort, and limit crew accommodation and crew complements. The experiment primarily aimed to collect seated height data for subjects exposed to microgravity environments, and feed new information regarding the effect of elongation of the spine forward into the design of the Orion. The data collected during the experiment included, two seated height measurement and two digital pictures of seated height pre-, in-, and post-flight. In addition to seated height, crewmembers had an optional task of collecting stature , standing height. Seated height data was obtained from 29 crewmembers that included 8 ISS increment crew (2 females and 6 males) and 21 Shuttle crew (1 female, 20 males), and whose mean age was 48 years ( 4 years). This study utilized the last six Shuttle flights, STS-128 to STS-134. The results show that partipating crewmembers experienced growth up to 6% in seated height and up to 3% in stature. Based on the worst case statistical analysis of the subject data, the recommended seated height growth of 6% will be provided to the designers as the necessary seated height adjustment.

Young, K. S.↗
DOE OSTI2023

A Simulation Framework for Analyzing the Impact of Stochastic Occupant Behaviors on Demand Flexibility in Typical Commercial Buildings

As one of the primary users of the electric grid, buildings and building equipment, including heating, ventilation, and air conditioning (HVAC) systems, can be leveraged to provide the flexible demand needed to balance the grid. Typical strategies to achieve demand flexibility are to reduce electricity use during peak or critical periods by shutting down equipment or relaxing system setpoints, which will inevitably impact the occupants’ comfort. When occupants feel uncomfortable, they may take actions to regain their comfort, and some of those actions (such as turning on a personal fan) may have a negative impact on meeting the demand response goal. Therefore, it is important to incorporate occupant behaviors into the assessment ofthe building demand flexibility potential. In this study, a simulation framework that includes simulation of zone thermal loads, an HVAC system, and occupant behaviors, was developed to investigate the impact of occupant behaviors on demand flexibility. A case study was conducted using a small office model from the U.S. Department of Energy (DOE) Commercial Prototype Building Models to simulate the building envelope and zone loads. An agent-based occupant thermal behavior model was adapted to forecast occupants’ thermal comfort and their resulting thermal behaviors. An artificial neural network (ANN) based airflow model trained from a computational fluid dynamics (CFD) model of the zone was adopted to better predict the ambient environment of each occupant. An air-source heat pump simulation model that was calibrated from a real two-stage air-source heat pump system was used as the HVAC system. A typical load shedding event during peak hours was studied. Repeated simulations were conducted to capture the stochastic effects of occupant behaviors. The interplay between the demand flexibility, occupant comfort and behavior were analyzed by evaluating key performance indicators, including the energy use, occupant discomfort duration, and occupant behavior duration during the peak period. The results suggest that this framework can be used to analyze typical commercial buildings and their HVAC systems in terms ofdemand flexibility potential under the impact of occupant behaviors.

Chen, Zhelun↗
NASA NTRS1980

Thermal energy storage systems using fluidized bed heat exchangers

A rotary cement kiln and an electric arc furnace were chosen for evaluation to determine the applicability of a fluid bed heat exchanger (FBHX) for thermal energy storage (TES). Multistage shallow bed FBHX's operating with high temperature differences were identified as the most suitable for TES applications. Analysis of the two selected conceptual systems included establishing a plant process flow configuration, an operational scenario, a preliminary FBHX/TES design, and parametric analysis. A computer model was developed to determine the effects of the number of stages, gas temperatures, gas flows, bed materials, charge and discharge time, and parasitic power required for operation. The maximum national energy conservation potential of the cement plant application with TES is 15.4 million barrels of oil or 3.9 million tons of coal per year. For the electric arc furnance application the maximum national conservation potential with TES is 4.5 million barrels of oil or 1.1 million tons of coal per year. Present time of day utility rates are near the breakeven point required for the TES system. Escalation of on-peak energy due to critical fuel shortages could make the FBHX/TES applications economically attractive in the future.

Weast, T.↗
NASA NTRS2010

High-Altitude Hydration System

Three methods are being developed for keeping water from freezing during high-altitude climbs so that mountaineers can remain hydrated. Three strategies have been developed. At the time of this reporting two needed to be tested in the field and one was conceptual. The first method is Passive Thermal Control Using Aerogels. This involves mounting the fluid reservoir of the climber s canteen to an inner layer of clothing for better heat retention. For the field test, bottles were mounted to the inner fleece layer of clothing, and then aerogel insulation was placed on the outside of the bottle, and circumferentially around the drink straw. When climbers need to drink, they can pull up the insulated straw from underneath the down suit, take a sip, and then put it back into the relative warmth of the suit. For the field test, a data logger assessed the temperatures of the water reservoir, as well as near the tip of the drink straw. The second method is Passive Thermal Control with Copper-Shielded Drink Straw and Aerogels, also mounted to inner layers of clothing for better heat retention. Braided wire emanates from the inside of the fleece jacket layer, and continues up and around the drink straw in order to use body heat to keep the system-critical drink straw warm enough to keep water in the liquid state. For the field test, a data logger will be used to compare this with the above concept. The third, and still conceptual, method is Active Thermal Control with Microcontroller. If the above methods do not work, microcontrollers and tape heaters have been identified that could keep the drink straw warm even under extremely cold conditions. Power requirements are not yet determined because the thermal environment inside the down suit relative to the external environment has not been established. A data logger will be used to track both the external and internal temperatures of the suit on a summit day.

Parazynski, Scott E.↗
DOE OSTI2024

Identification of Dynamic Force Coefficients for an Additively Manufactured Hermetic Squeeze Film Bearing Support Damper Utilizing a Pass-Through Channel

Abstract The following paper presents breakthrough experimental results for a new hermetic squeeze film damper (HSFD) concept that is integrally designed within an externally pressurized tilting-pad radial gas bearing support. The flexibly damped gas bearing module was designed for a 7.2″ (183 mm) diameter shaft and fabricated using direct metal laser melting (DMLM); also known as additive manufacturing. The bearing and HSFD were sized based on ongoing studies for oil-free supercritical carbon dioxide (sCO2) power turbines in the 8.5 MW–10 MW power range. The development of the new damper concept was motivated by past dynamic testing on HSFD, which generated frequency-dependent stiffness and damping force coefficients. In efforts to eliminate the frequency dependency, a new HSFD architecture was conceived that adds accumulator volumes and a pass-through channel to previously conceived HSFD flow network designs. The other motivation for the work is the need to develop a cost-effective and reliable oil-free bearing technology that is scalable to large power turbomachinery applications. There were several objectives for the following work. The first objective was to successfully design and fabricate a single piece bearing-damper using additive manufacturing, while dimensionally controlling critical design features. The paper discusses the manufacturing steps and shows cut-ups that reveal adequate clearance control capability with internal damper clearances. The second objective was to perform experimental testing with the new HSFD design in efforts to extract stiffness and damping coefficients for excitation frequencies within 20–160 Hz and peak vibration amplitudes between 0.25 mils (6.35 microns) to 1 mil (25.4 microns). The test results for a single HSFD bearing module indicated that the design modifications to the HSFD architecture were successful in eliminating nearly all the frequency dependencies for the stiffness force coefficient. The dynamic tests yielded a stiffness coefficient that varied between 112 klb/in. (19.6 MN/m) and 96 klb/in. (16.8 MN/m). The damping force coefficient however, exhibited relatively more variation with frequency with values residing between 175 lb-s/in. (31 kN-s/m) to 214 lb-s/in. (37 kN-s/m). Finally, the paper advances a three-dimensional fluid-structure interaction (FSI) model using transient finite element analysis (FEA) coupled to a computational fluid dynamics (CFD) model. The FSI analysis performed between 20 Hz and 80 Hz was used to predict the stiffness and damping of the HSFD using a quarter-section model of the damper. The FSI analysis was able to support test results by showing only a 6–7.4% change in the magnitude of force coefficients. Stiffness predictions agree reasonably well with experiments whereas damping is underpredicted.

Engineering↗
NASA NTRS1993

Computational Fluid Dynamics (CFD) applications in rocket propulsion analysis and design

Computational Fluid Dynamics (CFD) has been used in recent applications to affect subcomponent designs in liquid propulsion rocket engines. This paper elucidates three such applications for turbine stage, pump stage, and combustor chamber geometries. Details of these applications include the development of a high turning airfoil for a gas generator (GG) powered, liquid oxygen (LOX) turbopump, single-stage turbine using CFD as an integral part of the design process. CFD application to pump stage design has emphasized analysis of inducers, impellers, and diffuser/volute sections. Improvements in pump stage impeller discharge flow uniformity have been seen through CFD optimization on coarse grid models. In the area of combustor design, recent CFD analysis of a film cooled ablating combustion chamber has been used to quantify the interaction between film cooling rate, chamber wall contraction angle, and geometry and their effects of these quantities on local wall temperature. The results are currently guiding combustion chamber design and coolant flow rate for an upcoming subcomponent test. Critical aspects of successful integration of CFD into the design cycle includes a close-coupling of CFD and design organizations, quick turnaround of parametric analyses once a baseline CFD benchmark has been established, and the use of CFD methodology and approaches that address pertinent design issues. In this latter area, some problem details can be simplified while retaining key physical aspects to maintain analytical integrity.

Mcconnaughey, P. K.↗