SEARCH · Search NASA
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.
The laws of fluid motion are examined systematically for the case where friction is assumed to be very slight. Calculations are carried out with the appropriate differential equation and practical investigations are illustrated.
The present FEM technique addresses both linear and nonlinear boundary value problems encountered in computational physics by handling general three-dimensional regions, boundary conditions, and material properties. The box finite elements used are defined by a Cartesian grid independent of the boundary definition, and local refinements proceed by dividing a given box element into eight subelements. Discretization employs trilinear approximations on the box elements; special element stiffness matrices are included for boxes cut by any boundary surface. Illustrative results are presented for representative aerodynamics problems involving up to 400,000 elements.
Recently, NASA, FAA, and other organizations have focused their attention upon the possible effects of rain on airfoil performance. Rhode carried out early experiments and concluded that the rain impacting the aircraft increased the drag. Bergrum made numerical calculation for the rain effects on airfoils. Luers and Haines did an analytic investigation and found that heavy rain induces severe aerodynamic penalties including both a momentum penalty due to the impact of the rain and a drag and lift penalty due to rain roughening of the airfoil and fuselage. More recently, Hansman and Barsotti performed experiments and declared that performance degradation of an airfoil in heavy rain is due to the effective roughening of the surface by the water layer. Hansman and Craig did further experimental research at low Reynolds number. E. Dunham made a critical review for the potential influence of rain on airfoil performance. Dunham et al. carried out experiments for the transport type airfoil and concluded that there is a reduction of maximum lift capability with increase in drag. There is a scarcity of published literature in analytic research of two-phase boundary layer around an airfoil. Analytic research is being improved. The following assumptions are made: the fluid flow is non-steady, viscous, and incompressible; the airfoil is represented by a two-dimensional flat plate; and there is only a laminar boundary layer throughout the flow region. The boundary layer approximation is solved and discussed.
Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.
Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.
Cavitation bubble collapse in venturi, discussing flow characteristics under large pressure gradient
The analytical prediction and description of transonic flow in turbomachinery is complicated by three fundamental effects: (1) the fluid equations describing the transonic regime are inherently nonlinear, (2) shock waves may be present in the flow, and (3) turbomachine blading is geometrically complex, possessing large amounts of curvature, stagger, and twist. A three-dimensional computation procedure for the study of transonic turbomachine fluid mechanics is described. The fluid differential equations and corresponding difference operators are presented, the boundary conditions for complex blade shapes are described, and the computational implementation and mapping procedures are developed. Illustrative results of a typical unthrottled transonic rotor are also presented.
An analytical solution is presented for the baroclinic stability problem of a Boussinesq fluid in a beta-plane channel with Ekman suction boundary conditions. All of the modes, stable and unstable, belonging to this problem are identified. It is found that an unstable mode exists for only a certain range of values of the Burger number. The value of the Burger number at the upper limit of this range increases as the Ekman number decreases. Beyond this upper limit only a damped mode exists. It is also found that this transition in parameter space from the unstable to the stable mode occurs in a discontinuous manner.
The major criterion for the Atmospheric General Circulation Experiment (AGCE) design is that it be possible to realize strong baroclinic instability in the spherical configuration chosen. A configuration was selected in which a hemispherical shell of fluid is subjected to latitudinal temperature gradients on its spherical boundaries and the latitudinal boundaries are insulators. Work in the laboratory with a cylindrical version of this configuration revealed more instabilities than baroclinic instability. Since researchers fully expect these additional instabilities to appear in the spherical configuration also, they decided to continue the laboratory cylindrical annulus studies. Four flow regimes were identified: an axisymmetric Hadley circulation, boundary layer convection, baroclinic waves and deep thermal convection. Regime diagrams were prepared.
An accurate finite difference scheme is used to study the performance of a newly developed multiple-scale turbulence model in the similarity region of plane/radial/round jet flows. Under the assumption that molecular viscosity can be neglected, the governing equations contain a hidden eigenvalue which corresponds to the sharp boundary of the flowfield. It is found that this eigenvalue is the locus of the free boundary separating the turbulent shear flow and the irrotational ambient fluid.
While 'no-slip' boundary conditions and the Navier-Stokes equations of continuum fluid mechanics have served the vapor transport community well until now, it is pointed out that transport conditions within highly nonisothermal ampoules are such that the nonisothermal side walls 'drive' the dominant convective flow, and the familiar Stokes-Fourier-Fick laws governing the molecular fluxes of momentum, energy, and (species) mass in the 'continuum' field equations will often prove to be inadequate, even at Knudsen numbers as small as 0.001. The implications of these interesting gas kinetic phenomena under microgravity conditions, and even under 'earth-bound' experimental conditions, are outlined here, along with a tractable approach to their systematic treatment.
The behavior of thermally conductive plasma flows in helmet-streamer coronal structures is investigated within the framework of the axisymmetric nonrotating one-fluid MHD model. Continuous subsonic-supersonic solutions satisfying observed boundary conditions at the sun as well as the vanishing of the temperature at infinity are obtained and presented. Special attention is paid to the combined effects of conductive flow (and corresponding thermal force) and rapidly diverging magnetic field on the critical points. In this, the heliocentric distance of the neutral point determining the separation between closed and open field lines (cusp) is treated as a free parameter. These thermally conductive solutions are contrasted with those provided by corresponding isothermal models.
An unsteady, 3D, implicit upwind Euler/Navier-Stokes algorithm is here used to compute the flutter characteristics of Wing 445.6, the AGARD standard aeroelastic configuration for dynamic response, with a view to the discrepancy between Euler characteristics and experimental data. Attention is given to effects of fluid viscosity, structural damping, and number of structural model nodes. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. The V-g analysis indicates that fluid viscosity has a significant effect on the supersonic flutter boundary for this wing.
The 'radial' transport of energy by internal ULF waves, stimulated by dayside magnetospheric boundary oscillations, is analyzed in the framework of one-fluid magnetohydrodynamics. (the term radial is used here to denote the direction orthogonal to geomagnetic flux surfaces.) The model for the inhomogeneous magnetospheric plasma and background magnetic field is axisymmetric and includes radial and parallel variations in the magnetic field, magnetic curvature, plasma density, and low but finite plasma pressure. The radial mode structure of the coupled fast and intermediate MHD waves is determined by numerical solution of the inhomogeneous wave equation; the parallel mode structure is characterized by a Wentzel-Kramer-Brillouin (WKB) approximation. Ionospheric dissipation is modeled by allowing the parallel wave number to be complex. For boudnary oscillations with frequencies in the range from 10 to 48 mHz, and using a dipole model for the background magnetic field, the combined effects of magnetic curvature and finite plasma pressure are shown to (1) enhance the amplitude of field line resonances by as much as a factor of 2 relative to values obtained in a cold plasma or box-model approximation for the dayside magnetosphere; (2) increase the energy flux delivered to a given resonance by a factor of 2-4; and (3) broaden the spectral width of the resonance by a factor of 2-3. The effects are attributed to the existence of an 'Alfven buoyancy oscillation,' which approaches the usual shear mode Alfven wave at resonance, but unlike the shear Alfven mode, it is dispersive at short perpendicular wavelengths. The form of dispersion is analogous to that of an internal atmospheric gravity wave, with the magnetic tension of the curved background field providing the restoring force and allowing radial propagation of the mode. For nominal dayside parameters, the propagation band of the Alfven buoyancy wave occurs between the location of its (field line) resonance and that of the fast mode cutoff that exists at larger radial distances.
Laser Doppler velocimetry (LDV) measurements and hydrogen bubble flow visualization techniques were used to examine the near-wall flow structure of 2D and 3D turbulent boundary layers (TBLs) over a range of low Reynolds numbers. The goals of this research were (1) an increased understanding of the flow physics in the near wall region of turbulent boundary layers,(2) to observe and quantify differences between 2D and 3D TBL flow structures, and (3) to document Reynolds number effects for 3D TBLs. The LDV data have provided results detailing the turbulence structure of the 2D and 3D TBLs. These results include mean Reynolds stress distributions, flow skewing results, and U and V spectra. Effects of Reynolds number for the 3D flow were also examined. Comparison to results with the same 3D flow geometry but at a significantly higher Reynolds number provided unique insight into the structure of 3D TBLs. While the 3D mean and fluctuating velocities were found to be highly dependent on Reynolds number, a previously defined shear stress parameter was discovered to be invariant with Reynolds number. The hydrogen bubble technique was used as a flow visualization tool to examine the near-wall flow structure of 2D and 3D TBLs. Both the quantitative and qualitative results displayed larger turbulent fluctuations with more highly concentrated vorticity regions for the 2D flow.
Spray cooling has high potential in thermal management and life support systems by overcoming the deleterious effect of microgravity upon two-phase heat transfer. In particular spray cooling offers several advantages in heat flux removal that include the following: 1. By maintaining a wetted surface, spray droplets impinge upon a thin fluid film rather than a dry solid surface 2. Most heat transfer surfaces will not be smooth but rough. Roughness can enhance conductive cooling, aid liquid removal by flow channeling. 3. Spray momentum can be used to a) substitute for gravity delivering fluid to the surface, b) prevent local dryout and potential thermal runaway and c) facilitate liquid and vapor removal. Yet high momentum results in high We and Re numbers characterizing the individual spray droplets. Beyond an impingement threshold, droplets splash rather than spread. Heat flux declines and spray cooling efficiency can markedly decrease. Accordingly we are investigating droplet impingement upon a) dry solid surfaces, b) fluid films, c) rough surfaces and determining splashing thresholds and relationships for both dry surfaces and those covered by fluid films. We are presently developing engineering correlations delineating the boundary between splashing and non-splashing regions.
Results of an experimental investigation on isolated synthetic jets in crossflow from rectangular orifices of different aspect ratio and orientation are presented. Three aspect ratios, AR = 4, 8, and 16, with pitch a = 90deg., were investigated, Additionally, the AR = 8 case was pitched at 20 deg.. The yaw angle, , was varied through 0 deg., 10 deg., 45 deg. and 90 deg.. All orifices had same exit area and the data were compared with synthetic as well as steady jet from a circular orifice of same area. Hotwire measurements were performed to obtain all three components of mean velocity and turbulent stresses. Data were acquired for momentum-flux ratio up to J = 50. Distributions of time- and phase-averaged data were obtained on the cross sectional plane at x/D = 0.5, 5 and 10, as well as on the axial plane of the symmetry. Qualitative flowfield similarity between synthetic and steady jets is observed. However, high-momentum 'cap' above the low-momentum 'dome', characteristic of steady jets, does not necessarily appear in the synthetic jet. The position and shape of the high-momentum region depend on the distance from the orifice, pitch, yaw as well as momentum-flux ratio. Consequently, the location of the minimum velocity in the 'dome' measured at the plane of symmetry, y(submin), is adopted as a reference for penetration estimate and trajectory comparison. For AR = 16, the dome is the largest in area with maximum velocity deficit. However, the penetration is somewhat higher for AR = 4. Increase in yaw reduces the spatial extent of the dome and the penetration height but augments the velocity deficit. At low J the dome is connected to the boundary layer and traces of the cap of high momentum fluid are visible above it. Increase in J lifts the dome and reorganizes the high-momentum fluid around its perimeter, eventually bringing it underneath. Phase-averaged data document dynamic topological changes within the cycle. Phase-averaged streamwise velocity contours on the cross-sectional plane exhibit behavior commensurate with that seen in time-averaged data at various J.
The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.