Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stokes parameters”

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 217 records · Page 12

Combustion instability analysis

A theory and computer program for combustion instability analysis are presented. The basic theoretical foundation resides in the concept of entropy-controlled energy growth or decay. Third order perturbation expansion is performed on the entropy-controlled acoustic energy equation to obtain the first order integrodifferential equation for the energy growth factor in terms of the linear, second, and third order energy growth parameters. These parameters are calculated from Navier-Stokes solutions with time averages performed on as many Navier-Stokes time steps as required to cover at least one peak wave period. Applications are made for a 1-D Navier-Stokes solution for the Space Shuttle Main Engine (SSME) thrust chamber with cross section area variations taken into account. It is shown that instability occurs when the mean pressure is set at 2000 psi with 30 percent disturbances. Instability also arises when the mean pressure is set at 2935 psi with 20 percent disturbances. The system with mean pressures and disturbances more adverse that these cases were shown to be unstable.

Chung, T. J.↗

Numerical Study of High-Temperature Jet Flow Using RANS/LES and PANS Formulations

Two multi-scale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the Partially Averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for Partial Averaged Navier-Stokes (PANS). It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide valuable tool for accurate jet noise predictions. Solutions from these models are compared to RANS results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid RANS/LES and PANS in simulating such flow phenomena.

Abdol-Hamid, Khaled S.↗

The Primordial Inflation Polarization Explorer (PIPER)

The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne instrument to measure the gravity-wave signature of primordial inflation through its distinctive imprint on the polarization of the cosmic microwave background. PIPER combines cold (1.5 K) optics, 5120 bolometric detectors, and rapid polarization modulation using VPM grids to achieve both high sensitivity and excellent control of systematic errors. A series of flights alternating between northern and southern hemisphere launch sites will produce maps in Stokes I, Q, U, and V parameters at frequencies 200, 270, 350, and 600 GHz (wavelengths 1500, 1100, 850, and 500 microns) covering 85% of the sky. We describe the PIPER instrument and discuss the current status and expected science returns from the mission.

Kogut, Alan J.↗

Inlet Vortex Generator Design for the X-59 Low Boom Flight Demonstrator

An array of vortex generators was designed for the X-59 Low Boom Flight Demonstrator aircraft inlet. Prior to the a High Speed Inlet Distortion Test, concern regarding potentially high measurements of steady state or dynamic distortion prompted an effort to examine the use of passive flow control to reduce distortion levels. Design of experiments methodology was used to develop two matrices of vortex generators placed inside the inlet serpentine duct and upstream of the inlet on the compression surface, varying on several geometric parameters. Reynolds-Averaged Navier-Stokes simulations were performed using the FUN3D flow solver at Mach 1.45 to analyze the impact of each configuration on total pressure recovery, circumferential distortion, and radial distortion. Results from simulations of the baseline aircraft are compared with data acquired in the NASA Glenn Research Center’s 8- by 6-foot Supersonic Wind Tunnel. A vortex generator configuration with vanes placed upstream of the inlet terminal shock was shown to significantly increase total pressure recovery and decrease circumferential distortion with no penalty to radial distortion levels. Inlet mass flow sweeps at Mach 1.45, 0.6, and 0.3 indicated no detrimental effects from the vanes at off-design conditions.

Vortex generators↗

Evaluation of turbulence induced noise in coherent anti-Stokes Raman scattering

The effect of turbulence in a transonic wind tunnel on coherent anti-Stokes Raman scattering is considered. The driving pump and Stokes waves are taken to be coaxially propagating Gaussian beam waves which are focused on the Raman active medium through the turbulent boundary layer of the flow tube. The random index of refraction variations in the layer are modeled as phase perturbations of the driving waves which cause a reduction of the mean on-axis field and an increase in the mean diameter of the beams. Effective Gaussian beam parameters are developed and the radiated anti-Stokes power calculated as a function of the phase screen parameters. A significant reduction in signal strength occurs for realistic estimates of the phase screen parameter appropriate to a confined transonic flow. A method for estimating the signal degradation which could be applied to other experimental situations is presented.

Elliott, R. A.↗

Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law

Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities ( η ). Constraining geofluid η at relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values of η . When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrect η . We also suggest the uncertainty in η is best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.

Geochemistry & Geophysics↗

Polarization effects on hard target calibration of lidar systems

The theory of hard target calibration of lidar backscatter data, including laboratory measurements of the pertinent target reflectance parameters, is extended to include the effects of polarization of the transmitted and received laser radiation. The bidirectional reflectance-distribution function model of reflectance is expanded to a 4 x 4 matrix allowing Mueller matrix and Stokes vector calculus to be employed. Target reflectance parameters for calibration of lidar backscatter data are derived for various lidar system polarization configurations from integrating sphere and monostatic reflectometer measurements. It is found that correct modeling of polarization effects is mandatory for accurate calibration of hard target reflectance parameters and, therefore, for accurate calibration of lidar backscatter data.

Kavaya, Michael J.↗

Target Calibration: Retro-Reflection Mechanisms at 2 Micron

With the growing use of 2-micrometer eye-safe solid-state laser transmitters in coherent Doppler lidar applications, there is increased interest in the reflectance properties of hard target calibration materials at this wavelength. One of the factors to be considered is the polarization characteristic of the transmitter/receiver. The two fundamental lidar systems are either based upon the backscattering of linearly polarized or circularly polarized light. Consequently, the response of calibration materials to the particular state of polarization of the incident radiation is an important consideration. Ideally the calibration materials should have a similar reflectance response to that of the remote scattering medium that is being observed. Also the material should be most efficiently adapted to the polarization characteristics of the lidar system. There are two measurement parameters that are useful in presenting the characteristics of polarization of reflecting materials. The first is the linear polarization ratio: mu (sub l) =I (sub ol)/I (sub sl), where ol denotes opposite linear and sl denotes same linear. The second is circular polarization ratio: mu (sub c) = I (sub sc)/I (sub oc) where ol denotes opposite circular and sc denotes same circular. Notice that the ratios are inverted with respect to same polarization. It is the relative values of these polarization ratios that are used to characterize the retro-reflectance mechanism. Measurement of these parameters is conveniently obtained by using the Stokes vectors corresponding to the polarized backscatter at and near the retro-angle.

Haner, David A.↗

Upper semicontinuous global attractors for viscous flow

A particular form of regularization of the Navier-Stokes equations is studied. It is shown that as the regularization parameter goes to zero, the global attractor for the regularized system converges to the global attractor for the conventional Navier-Stokes system.

Ou, Yuh-Roung↗

Validating CAR - A comparison study of experimentally-derived and computer-generated reach envelopes

In the present investigation, Crewstation Assessment of Reach (CAR) results in the form of male hand reach envelopes were generated and compared with an anthropometric survey performed by Kennedy (1978) to determine the extent of the validity of the CAR model with respect to experimentally-derived anthropometric data. The CAR-generated reach envelopes extensively matched the Kennedy envelopes. The match was particularly good in the areas to the front and side from which the reach originated. Attention is given to the crewstation model, the operator sample population, the CAR analysis, aspects of validation methodology, and the modeling of experimental parameters.

Harris, R.↗

Continuation in a parameter - Experience with viscous and free surface flows

The results of modifications in continuation methods applied to obtain solutions to the Navier-Stokes systems of equations for incompressible, two-dimensional, steady flows are reported. It is shown that parameter continuation permits prediction of accurate, initial estimates for iterative processing of nonlinear finite difference and finite element equations of motions. The new parameter steps are derived from values of the preceding parameter steps. The accuracy of the estimates is ensured with appropriate choices of the step size. The continuation predictor/iterative corrector is demonstrated to trace the branches of parameter space along which steady flow states are found, and techniques are available for tracing multiply branching paths. The techniques are applied to solving the Navier-Stokes equations for flow through a rotating square channel, the formation of a falling liquid curtain, and gyrostatic equilibria of rotating cylindrical drops.

Kheshgi, H. S.↗

Numerical and analytical study of fluid dynamic forces in seals and bearings

Study results for rotor/bearing/seal system models are compared with both a numerical model based on a transformed, conservative version of the three-dimensional Navier-Stokes equation, and an analytical model based on 'lumped' fluid parameters. Rotor-destabilizing factors are related to flow field rotation; it is shown that these factors can be reduced through a decrease in the fluid average circumferential velocity. A transformation of the Navier-Stokes equations to those for a rotating observer and a numerical code with a nonorthogonal body-fitted grid are used to carry out numerical experiments demonstrating the character of the complex flow field.

Tam, L. T.↗

Modification of the PARC Navier-Stokes Code to predict rocket engine nozzle performance

The PARC2D Navier-Stokes Code was modified to compute the performance parameters for rocket engine nozzles. The perfect gas code was applied to the SSME engine nozzle for inviscid, laminar and turbulent flow. Inviscid computations compare well with Rocketdyne computations. Performance degradation due to the boundary layers is very reasonable. Application of the code to nontraditional nozzle geometries and to low Reynolds nozzles is demonstrated. Modification of the code for equilibrium H2/O2 chemistry is described. Thermodynamic and equilibrium constants are determined from statistical mechanics and the transport properties from exact kinetic theory, using collison integrals determined from appropriate intermolecular potentials. The equilibrium code was used to compute the SSME flowfield. Modifications of the flowfield due to the change of composition are described.

Collins, Frank G.↗

CFD modeling of turbulent flows around the SSME main injector assembly using porosity formulation

Hot gas turbulent flow distribution around the main injector assembly of the Space Shuttle Main Engine (SSME) and liquid oxygen (LOX) flow distributions through the LOX posts have great effect on the combustion phenomenon inside the main combustion chamber. An advanced computational fluid dynamics (CFD) analysis will help to provide more accurate and efficient characterization of this type of flow field. In order to design a CFD model to be an effective engineering analysis tool with good computational turn-around time and still maintain good accuracy in describing the flow features, the concept of porosity is employed to describe the effects of blockage and drag force due to the presence of the LOX posts in the turbulent flow field around the main injector assembly of the SSME. A validated non-isotropic porosity model is developed and incorporated into an existing Navier-Stokes flow solver (FDNS). Volume and surface porosity parameters, which are based on the configurations of local LOX post clustering, are introduced into the governing equations.

Cheng, Gary C.↗

Numerical Study of Flow Past a Circular Cylinder Using RANS, Hybrid RANS/LES and PANS Formulations

Two multiscale type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k(sub epsilon)) model with a RANS/LES transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) model, where the unresolved kinetic energy parameter (f(sub k)) is allowed to vary as a function of grid spacing and the turbulence length scale. Solutions from these models are compared to RANS results and experimental data for a stationary and rotating cylinder. The parameter f(sub k) varies between zero and one and has the characteristic to be equal to one in the viscous sub layer, and when the RANS turbulent viscosity becomes smaller than the LES viscosity. The formulation, usage methodology, and validation example are presented to demonstrate the enhancement of PAB3D's time-accurate and turbulence modeling capabilities. The models are compared to RANS results and experimental data for turbulent separated flows (TS) and laminar separated flows (LS) around stationary and rotating cylinders. For a stationary cylinder, the TS case is accurately simulated using the general two-equation k(sub epsilon) turbulence model (eddy-viscosity model). PAB3D accurately predicts the drag coefficient (CD), lift coefficient (CL) and the Strouhal number (St). The LS case was a challenge for the RANS computation with an eddy-viscosity turbulence model. The RANS/LES and PANS performed well and showed marked improvements over the RANS solution. The modified PANS model was the most accurate. For the rotating cylinder, the spin ratio varied from zero to one, and the PANS results were in good agreement with published experimental data. RANS/LES and PANS capture both temporal and spatial fluctuations and produce large-scale structures that do not occur in the RANS simulation. The current results show promise for the capability of PANS in simulating unsteady and complex flow phenomena.

Elmiligui, Alaa A.↗

Reduced-Order Modeling of a Heaving Airfoil

A reduced-order model of a flapping airfoil is developed using Proper Orthogonal Decomposition (POD). The proper basis functions, developed from snapshots of full Navier-Stokes simulations, are used for a Galerkin projection of the governing equations. The resulting coupled, nonlinear ordinary di.erential equations have a low dimension because the first few basis members capture most of the energy of the flow. The reduced-order model is used to simulate heaving motions that are both similar to and different from the motion(s) used to generate the basis functions, and the errors in the model are quantified. Several methods are used to generate mode sets that can be used over a range of heaving parameters, including snapshots from one, two, and multiple Navier-Stokes simulations. As snapshots from additional simulations are added to the decomposition, the mode sets become richer and can simulate a wider range of parameter space, at some computational cost. Whereas the POD method is fully applicable in three dimensions, the simulation technique based on a body-fixed and body-fitted grid suffers large overhead when extended to three dimensions. To reduce the overhead, an embedding technique is discussed which embeds the solid wing into a fixed Cartesian grid. The wing, which can now have multiple pieces and also be flexible, is represented by a distribution of body forces. This distribution is determined to give exactly the flow around a flapping wing.

Haj-Hariri, H.↗

Interfacial dynamics of two liquids under an oscillating gravitational field

The evolution of two miscible liquids meeting at an initially sharp interface inside a cavity under microgravity g-jitter conditions is studied numerically. The response of the interface and kinematics of the flowfield to various g-jitter accelerations and aspect ratio variations is shown. The interface region acts like a vortex source sheet, and it can be unstable to Kelvin-Helmholtz and Rayleigh-Taylor instabilities. The vortices produced along the interface can serve as a stirring mechanism to promote local mixing. Below the critical Stokes-Reynolds number, the destabilization of the interface results in deformation into wavy structures. In some parameter regions, these structures oscillate in time; in others, they are quasi-steady. Above the critical Stokes-Reynolds number, 'chaotic' instability results, and the interface breaks into concentration pockets. The morphology of the initial breakup is similar to observed Rayleigh-Taylor instability. Subsequent mixing of the two fluids after the breakup of the interface is then very rapid.

Duval, Walter M. B.↗

A Numerical Testbed for Remote Sensing of Aerosols, and its Demonstration for Evaluating Retrieval Synergy from a Geostationary Satellite Constellation of GEO-CAPE and GOES-R

We present a numerical testbed for remote sensing of aerosols, together with a demonstration for evaluating retrieval synergy from a geostationary satellite constellation. The testbed combines inverse (optimal-estimation) software with a forward model containing linearized code for computing particle scattering (for both spherical and non-spherical particles), a kernel-based (land and ocean) surface bi-directional reflectance facility, and a linearized radiative transfer model for polarized radiance. Calculation of gas absorption spectra uses the HITRAN (HIgh-resolution TRANsmission molecular absorption) database of spectroscopic line parameters and other trace species cross-sections. The outputs of the testbed include not only the Stokes 4-vector elements and their sensitivities (Jacobians) with respect to the aerosol single scattering and physical parameters (such as size and shape parameters, refractive index, and plume height), but also DFS (Degree of Freedom for Signal) values for retrieval of these parameters. This testbed can be used as a tool to provide an objective assessment of aerosol information content that can be retrieved for any constellation of (planned or real) satellite sensors and for any combination of algorithm design factors (in terms of wavelengths, viewing angles, radiance and/or polarization to be measured or used). We summarize the components of the testbed, including the derivation and validation of analytical formulae for Jacobian calculations. Benchmark calculations from the forward model are documented. In the context of NASA's Decadal Survey Mission GEOCAPE (GEOstationary Coastal and Air Pollution Events), we demonstrate the use of the testbed to conduct a feasibility study of using polarization measurements in and around the O2 A band for the retrieval of aerosol height information from space, as well as an to assess potential improvement in the retrieval of aerosol fine and coarse mode aerosol optical depth (AOD) through the synergic use of two future geostationary satellites, GOES-R (Geostationary Operational Environmental Satellite R-series) and TEMPO (Tropospheric Emissions: Monitoring of Pollution). Strong synergy between GEOS-R and TEMPO are found especially in their characterization of surface bi-directional reflectance, and thereby, can potentially improve the AOD retrieval to the accuracy required by GEO-CAPE.

Aerosols↗