Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Calculation of hypersonic shock structure using flux-split algorithms

There exists an altitude regime in the atmosphere that is within the continuum domain, but wherein the conventional Navier-Stokes equations cease to be accurate. The altitude limits for this so called continuum transition regime depend on vehicle size and speed. Within this regime the thickness of the bow shock wave is no longer negligible when compared to the shock stand-off distance and the peak radiation intensity occurs within the shock wave structure itself. For this reason it is no longer valid to treat the shock wave as a discontinuous jump and it becomes necessary to compute through the shock wave itself. To accurately calculate hypersonic flowfields, the governing equations must be capable of yielding realistic profiles of flow variables throughout the structure of a hypersonic shock wave. The conventional form of the Navier-Stokes equations is restricted to flows with only small departures from translational equilibrium; it is for this reason they do not provide the capability to accurately predict hypersonic shock structure. Calculations in the continuum transition regime, therefore, require the use of governing equations other than Navier-Stokes. Several alternatives to Navier-Stokes are discussed; first for the case of a monatomic gas and then for the case of a diatomic gas where rotational energy must be included. Results are presented for normal shock calculations with argon and nitrogen.

Eppard, W. M.↗
NASA NTRS1993

A theoretical and numerical study of polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum

In this paper, theoretical and numerical results of the polarimetric scattering and emission from random rough surfaces with anisotropic directional spectrum are presented for the remote sensing of ocean and soil surfaces. The polarimetric scattered field for rough dielectric surfaces is derived to the second order by the small perturbation method (SPM). It is found that the second-order scattered field is coherent in nature, and its coefficients for different polarizations present the lowest-order corrections to the Fresnel reflection coefficients of the surfaces. In addition, the cross-polarized (HV and VH) components of the coherent fields are reciprocal and not zero for surfaces with anisotropic directional spectrum when the azimuth angle of the incident direction is not aligned with the symmetry directions of surfaces. In order to verify the energy conservation condition of the theoretical results, which is important if the theory is to be applied to the passive polarimetry of rough surfaces, a Monte Carlo simulation is performed to numerically calculate the polarimetric reflectivities of one-dimensional random rough surfaces which are generated with a prescribed power-law spectrum in the spectral domain and transformed to the spatial domain by the FFT. The surfaces simulated by this approach are periodic with the period corresponding to the low-wavenumber cutoff. To calculate the scattering from periodic dielectric surfaces, the authors present a new numerical technique which applies the Floquet theorem to reduce the problem to one period and does not require the evaluation of one-dimensional periodic Green's function used in the conventional method of moment formulation. Once the scattering coefficients are obtained, the polarimetric Stokes vectors for the emission from the random surfaces are then calculated according to the Kirchhoff's law and are illustrated as functions of relative azimuth observation and row directions. The second-order SPM is also reduced to this one-dimensional case, and theoretical results for all four Stokes parameters are in excellent agreement with the numerical results obtained from the Monte Carlo simulation. Moreover, the second order coherent fields are found to be necessary for the theoretical evaluation of the third and fourth Stokes parameters. Without that, the reflectivities of random rough surfaces are significantly over-estimated, and the signs of the third and the fourth Stokes parameters are found to be incorrect, if calculated by using the first-order incoherent fields alone.

Yueh, S. H.↗
NASA NTRS1993

Prediction of incidence and surface roughness effects on turbine performance

The results of a Navier-Stokes analysis for predicting the change in turbine efficiency due to a change in either incidence or surface roughness is discussed. It was experimentally determined by Boynton, Tabibzadeh, and Hudson that polishing the SSME high pressure fuel turbine blades improved turbine efficiency by about 2 points over a wide range of operating conditions. These conditions encompassed the range of incidence seen by the turbine blading during flight. It is also necessary to be able to predict turbine performance at various operating points for future rocket turbopump applications. The code RVCQ3D, developed by Rod Chima, was used to determine the effects of changes in incidence angle on turbine blade row efficiency. The midspan Navier-Stokes results were used in conjunction with an inviscid flow analysis code to predict the efficiency of the two stage SSME over a wide range of operating conditions for smooth and rough turbine blades. The use of the Navier-Stokes analysis to predict changes in turbine efficiency due to variation in incidence angles was found to be superior to other incidence loss correlations available in the literature. The sensitivity of the Navier-Stokes results to grid parameters is discussed. The effects of the surface roughness were accounted for using the Cebeci-Chang rough wall turbulence model. This model was implemented in the code RVCQ3D. The implementation of this model for predicting the change in efficiency is also discussed.

Boyle, R. J.↗
NASA NTRS1994

Unique determination of solutions to the Burnett equations

Recent success in applying the Burnett equations to the one-dimensional shock-structure problem has raised the issue of whether the full Burnett equations can be used replace the Navier-Stokes equations for solving boundary-value problems in rarefied gasdynamics. As is familiar from the classical rarefied gasdynamics literature, the Burnett equations, if not solved as a successive approximation to the Navier-Stokes equations for a small Knudsen number, would require more boundary conditions than those in the Navier-Stokes system, owing to the presence of the higher-order derivatives. In this paper, this issue is examined with concrete solution examples for the steady Couette flows, addressing specifically whether solutions to the full Burnett equations can be uniquely determined without adding more boundary conditions than those in the Navier-Stokes system. The analysis, supported by detailed numerical solutions, confirms that additional boundary conditions are needed as long as the Knudsen number is not identically zero, lest the solution to the Burnett equations is not unique.

Lee, C. J.↗
NASA NTRS1997

Heat Transfer Computations of Internal Duct Flows With Combined Hydraulic and Thermal Developing Length

This study investigated the Navier-Stokes computations of the surface heat transfer coefficients of a transition duct flow. A transition duct from an axisymmetric cross section to a non-axisymmetric cross section, is usually used to connect the turbine exit to the nozzle. As the gas turbine inlet temperature increases, the transition duct is subjected to the high temperature at the gas turbine exit. The transition duct flow has combined development of hydraulic and thermal entry length. The design of the transition duct required accurate surface heat transfer coefficients. The Navier-Stokes computational method could be used to predict the surface heat transfer coefficients of a transition duct flow. The Proteus three-dimensional Navier-Stokes numerical computational code was used in this study. The code was first studied for the computations of the turbulent developing flow properties within a circular duct and a square duct. The code was then used to compute the turbulent flow properties of a transition duct flow. The computational results of the surface pressure, the skin friction factor, and the surface heat transfer coefficient were described and compared with their values obtained from theoretical analyses or experiments. The comparison showed that the Navier-Stokes computation could predict approximately the surface heat transfer coefficients of a transition duct flow.

Wang, C. R.↗
NASA NTRS1959

Structure of Weak Shock Waves in a Monatomic Gas

The profiles and thicknesses of normal shock waves in argon at Mach numbers of 1.335, 1.454, 1.576, and 1-713 were determined experimentally by means of a free-molecule probe whose equilibrium temperature is related by kinetic theory to the local flow properties and their gradients. Comparisons were made between the experimental shock profiles and the theoretical profiles calculated from the Navier-Stokes equations, the Grad 13-moment equations, and the Burnett equations. New, very accurate numerical integrations of the Burnett equations were obtained for this purpose with results quite different from those found by Zoller, to whom the solution of this problem is frequently attributed. The experimental shock profiles were predicted with approximately equal success by the Navier-Stokes and Burnett theories, while the 13-moment method was definitely less satisfactory. A surprising feature of the theoretical results is the relatively small difference in predictions between the Navier-Stokes and Burnett theories in the present range of shock strengths and the contrastingly large difference between predictions of Burnett and the 13-moment theories. It is concluded that the Navier-Stokes equations are correct for weak shocks and that within the present shock strength range the Burnett equations make no improvement which merits the trouble of solving them. For shocks of noticeably greater strength, say with a shock Mach number of more than 2.5, it remains fundamentally doubtful that any of these theories can be correct.

Sherman, F. S.↗
NASA NTRS1999

Parametric Studies of the Ejector Process within a Turbine-Based Combined-Cycle Propulsion System

Performance characteristics of the ejector process within a turbine-based combined-cycle (TBCC) propulsion system are investigated using the NPARC Navier-Stokes code. The TBCC concept integrates a turbine engine with a ramjet into a single propulsion system that may efficiently operate from takeoff to high Mach number cruise. At the operating point considered, corresponding to a flight Mach number of 2.0, an ejector serves to mix flow from the ramjet duct with flow from the turbine engine. The combined flow then passes through a diffuser where it is mixed with hydrogen fuel and burned. Three sets of fully turbulent Navier-Stokes calculations are compared with predictions from a cycle code developed specifically for the TBCC propulsion system. A baseline ejector system is investigated first. The Navier-Stokes calculations indicate that the flow leaving the ejector is not completely mixed, which may adversely affect the overall system performance. Two additional sets of calculations are presented; one set that investigated a longer ejector region (to enhance mixing) and a second set which also utilized the longer ejector but replaced the no-slip surfaces of the ejector with slip (inviscid) walls in order to resolve discrepancies with the cycle code. The three sets of Navier-Stokes calculations and the TBCC cycle code predictions are compared to determine the validity of each of the modeling approaches.

Georgiadis, Nicholas J.↗
NASA NTRS2002

Data Analysis for the SOLIS Vector Spectromagnetograph

The National Solar Observatory's SOLIS Vector Spectromagnetograph (VSM), which will produce three or more full-disk maps of the Sun's photospheric vector magnetic field every day for at least one solar magnetic cycle, is in the final stages of assembly. Initial observations, including cross-calibration with the current NASA/NSO spectromagnetograph (SPM) will soon be carried out at a test site in Tucson. This paper discusses data analysis techniques for reducing the raw data, calculation of line-of-sight magnetograms and both quick-look and high-precision inference of vector fields from Stokes spectral profiles. Existing SPM algorithms, suitably modified to accomodate the cameras, scanning pattern, and polarization calibration optics for the VSM, will be used to "clean" the raw data and to process line-of-sight, magnetograms. A recent. version of the High Altitude Observatory Milne-Eddington (HAO-ME) inversion code (Skumanich and Lites; 1987, 11)J 322, p. 473) will he used for high-precision vector fields since the algorithm has been extensively tested, is well understood, and is fast enough to complete data analysis within 24 hours of data acquisition. The simplified inversion algorithm of Auer, Heasley. arid House (1977, Sol. Phys. 55, p. 47) forms the initial guess for this version of the HAO-ME code and will be used for quick-look vector analysis of VSM data since its performance on simulated Stokes profiles is better than other candidate methods. Improvements (e.g., principal components analysis or neural networks) are under consideration and will be straightforward to implement. However, current resources are sufficient to store the original Stokes profiles only long enough for high-precision analysis. Retrospective reduction of Stokes data with improved methods will not be possible, and modifications will only be introduced when the advantages of doing so are compelling enough to justify discontinuity in the long-term data stream.

Jones, Harrison P.↗
NASA NTRS2004

CFD Approaches for Simulation of Wing-Body Stage Separation

A collection of computational fluid dynamics tools and techniques are being developed and tested for application to stage separation and abort simulation for next-generation launch vehicles. In this work, an overset grid Navier-Stokes flow solver has been enhanced and demonstrated on a matrix of proximity cases and on a dynamic separation simulation of a belly-to-belly wing-body configuration. Steady cases show excellent agreement between Navier-Stokes results, Cartesian grid Euler solutions, and wind tunnel data at Mach 3. Good agreement has been obtained between Navier-Stokes, Euler, and wind tunnel results at Mach 6. An analysis of a dynamic separation at Mach 3 demonstrates that unsteady aerodynamic effects are not important for this scenario. Results provide an illustration of the relative applicability of Euler and Navier-Stokes methods to these types of problems.

Buning, Pieter G.↗
NASA NTRS2004

Transitional Bubble in Periodic Flow Phase Shift

One particular characteristic observed in unsteady shear layers is the phase shift relative to the main flow. In attached boundary layers this will have an effect both on the instantaneous skin friction and heat transfer. In separation bubbles the contribution to the drag is dominated by the pressure distribution. However, the most significant effect appears to be the phase shift on the transition process. Unsteady transition behaviour may determine the bursting of the bubble resulting in an un-recoverable full separation. An early analysis of the phase shift was performed by Stokes for the incompressible boundary layer of an oscillating wall and an oscillating main flow. An amplitude overshoot within the shear layer as well as a phase shift were observed that can be attributed to the relatively slow diffusion of viscous stresses compared to the fast change of pressure. Experiments in a low speed facility with the boundary layer of a flat plate were evaluated in respect to phase shift. A pressure distribution similar to that on the suction surface of a turbomachinery aerofoil was superimposed generating a typical transitional separation bubble. A periodically unsteady main flow in the suction type wind tunnel was introduced via a rotating flap downstream of the test section. The experiments covered a range of the three similarity parameters of momentum-loss-thickness Reynolds-number of 92 to 226 and Strouhal-number (reduced frequency) of 0.0001 to 0.0004 at the separation point, and an amplitude range up to 19 %. The free stream turbulence level was less than 1% .Upstream of the separation point the phase shift in the laminar boundary layer does not appear to be affected significantly bay either of the three parameters. The trend perpendicular to the wall is similar to the Stokes analysis. The problem scales well with the wave velocity introduced by Stokes, however, the lag of the main flow near the wall is less than indicated analytically. The separation point comes closest to the Stokes analysis but the phase is still 20 degrees lower at the wall.

Talan, M.↗
NASA NTRS2010

The Primordial Inflation Polarization Explorer (PIPER)

The Primordial Inflation Polarization Explorer (PIPER) is it balloon-borne instrument designed to search for the faint signature of inflation in the polarized component of the cosmic microwave background (C-N-113). Each flight will be configured for a single frequency, but in order to aid in the removal of the polarized foreground signal due to Galactic dust, the filters will be changed between flights. In this way, the CMB polarization at a total of four different frequencies (200, 270, 350, and 600 GHz) will be, measured on large angular scales. PIPER consists of a pair of cryogenic telescopes, one for measuring each of Stokes Q and U in the instrument frame. Each telescope receives both linear orthogonal polarizations in two 32 x 40 element planar arrays that utilize Transition-Edge Sensors (TES). The first element in each telescope is a variable-delay polarization modulator (VPM) that fully modulates the linear Stokes parameter to which the telescope is sensitive. There are several advantages to this architecture. First, by modulating at the front of the optics, instrumental polarization is unmodulated and is therefore cleanly separated from source polarization. Second, by implementing this system with the appropriate symmetry, systematic effects can be further mitigated. In the PIPER design, many of the. systematics are manifest in the unmeasured linear Stokes parameter for each telescope and this can be separated from the desired signal. Finally, the modulation cycle never mixes the Q and U linear Stokes parameters, and thus residuals in the modulation do not twist the observed polarization vector. This is advantageous because measuring the angle of linear polarization is critical for separating the inflationary signal from other polarized components.

Chuss, David T.↗
NASA NTRS2011

Establishing Approaches to Modeling the Ares I-X and Ares I Roll Control System with Free-stream Interaction

Approaches were established for modeling the roll control system and analyzing the jet interactions of the activated roll control system on Ares-type configurations using the USM3D Navier-Stokes solver. Components of the modeling approach for the roll control system include a choice of turbulence models, basis for computing a dynamic equivalence of the real gas rocket exhaust flow in terms of an ideal gas, and techniques to evaluate roll control system performance for wind tunnel and flight conditions. A simplified Ares I-X configuration was used during the development phase of the roll control system modeling approach. A limited set of Navier-Stokes solutions was obtained for the purposes of this investigation and highlights of the results are included in this paper. The USM3D solutions were compared to equivalent solutions at select flow conditions from a real gas Navier- Stokes solver (Loci-CHEM) and a structured overset grid Navier-Stokes solver (OVERFLOW).

Pao, S. Paul↗
NASA NTRS2014

Aeroelastic Analyses of the SemiSpan SuperSonic Transport (S4T) Wind Tunnel Model at Mach 0.95

Detailed aeroelastic analyses of the SemiSpan SuperSonic Transport (S4T) wind tunnel model at Mach 0.95 with a 1.75deg fixed angle of attack are presented. First, a numerical procedure using the Computational Fluids Laboratory 3-Dimensional (CFL3D) Version 6.4 flow solver is investigated. The mesh update method for structured multi-block grids was successfully applied to the Navier-Stokes simulations. Second, the steady aerodynamic analyses with a rigid structure of the S4T wind tunnel model are reviewed in transonic flow. Third, the static analyses were performed for both the Euler and Navier-Stokes equations. Both the Euler and Navier-Stokes equations predicted a significant increase of lift forces, compared to the results from the rigid structure of the S4T wind-tunnel model, over various dynamic pressures. Finally, dynamic aeroelastic analyses were performed to investigate the flutter condition of the S4T wind tunnel model at the transonic Mach number. The condition of flutter was observed at a dynamic pressure of approximately 75.0-psf for the Navier-Stokes simulations. However, it was observed that the flutter condition occurred a dynamic pressure of approximately 47.27-psf for the Euler simulations. Also, the computational efficiency of the aeroelastic analyses for the S4T wind tunnel model has been assessed.

Hur, Jiyoung↗
DOE OSTI2020

Phase field benchmark problems targeting fluid flow and electrochemistry

Here, we continue our development of phase field model benchmark problems with the addition of a third set, complimenting our previously developed problems for diffusion, precipitation, dendritic growth and linear elasticity. These benchmark problems are being jointly developed by the Center for Hierarchical Materials Design (CHiMaD) and the National Institute of Standards and Technology (NIST) along with input from the phase field community. The first problem in this third set targets Stokes flow, with a particular emphasis on flow around an obstruction placed asymmetrically in the domain. While Stokes flow is not traditionally in the canon of phase field problems, it is a class of problems gaining importance in areas such as filtration and water purification. The second problem deals with coupled Cahn-Hilliard diffusion and electrostatic forces, which is an important area in energy storage and battery sciences. We present our own solutions and discuss sources of numerical errors for the Stokes problem as well as simple checks to avoid fundamental issues in the coupled diffusion-electrostatics problem. The latter problem contains some subtleties that we expand on in an Appendix.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On exact-WKB analysis, resurgent structure, and quantization conditions

There are two well-known approaches to studying nonperturbative aspects of quantum mechanical systems: saddle point analysis of the partition functions in Euclidean path integral formulation and the exact-WKB analysis based on the wave functions in the Schrödinger equation. In this work, based on the quantization conditions obtained from the exact-WKB method, we determine the relations between the two formalism and in particular show how the two Stokes phenomena are connected to each other: the Stokes phenomenon leading to the ambiguous contribution of different sectors of the path integral formulation corresponds to the change of the “topology” of the Stoke curves in the exact-WKB analysis. We also clarify the equivalence of different quantization conditions including Bohr-Sommerfeld, path integral and Gutzwiller’s ones. In particular, by reorganizing the exact quantization condition, we improve Gutzwiller’s analysis in a crucial way by bion contributions (incorporating complex periodic paths) and turn it into an exact result. Furthermore, we argue the novel meaning of quasi-moduli integral and provide a relation between the Maslov index and the intersection number of Lefschetz thimbles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. Additionally, the scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Dual-solver hybrid computational approaches for design and analysis of vertical lift vehicles

Abstract The cost of Reynolds-Averaged Navier-Stokes simulations can be restrictive to implement in aeromechanics design and analysis of vertical lift configurations given the cost to resolve the flow on a mesh sufficient to provide accurate aerodynamic and structural loads. Dual-solver hybrid methods have been developed that resolve the configuration and the near field with the Reynolds-Averaged Navier-Stokes solvers, while the wake is resolved with vorticity-preserving methods that are more cost-effective. These dual-solver approaches can be integrated into an organisation’s workflow to bridge the gap between lower-fidelity methods and the expensive Reynolds-Averaged Navier-Stokes when there are complex physics present. This paper provides an overview of different dual-solver hybrid methods, coupling approaches, and future efforts to expand their capabilities in the areas of novel configurations and operations in constrained and turbulent environments.

Engineering↗