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.

SourceResearch connectionFrom yearThrough year

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

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Low-order preconditioning of the Stokes equations

A well-known strategy for building effective preconditioners for higher-order discretizations of some PDEs, such as Poisson's equation, is to leverage effective preconditioners for their low-order analogs. In this work, we show that high-quality preconditioners can also be derived for the Taylor–Hood discretization of the Stokes equations in much the same manner. In particular, we investigate the use of geometric multigrid based on the Q 1 iso Q 2 /Q 1 discretization of the Stokes operator as a preconditioner for the Q 2 /Q 1 discretization of the Stokes system. We utilize local Fourier analysis to optimize the damping parameters for Vanka and Braess–Sarazin relaxation schemes and to achieve robust convergence. Furthermore, these results are then verified and compared against the measured multigrid performance. While geometric multigrid can be applied directly to the Q 2 /Q 1 system, our ultimate motivation is to apply algebraic multigrid within solvers for Q 2 /Q 1 systems via the Q 1 iso Q 2 /Q 1 discretization, which will be considered in a companion paper.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

Scramjet mixing establishment times for a pulse facility

A numerical simulation of the temporally developing flow through a generic scramjet combustor duct is presented for stagnation conditions typical of flight at Mach 13 as produced by a shock tunnel pulse facility. The particular focus is to examine the start up transients and to determine the time required for certain flow parameters to become established. The calculations were made with a Navier-Stokes solver SPARK with temporally relaxing inflow conditions derived from operation of the T4 shock tunnel at the University of Queensland in Australia. Calculations at nominal steady inflow conditions were made for comparison. The generic combustor geometry includes the injection of hydrogen fuel from the base of a centrally located strut. In both cases, the flow was assumed laminar and fuel combustion was not included. The establishment process is presented for viscous parameters in the boundary layer and for parameters related to the fuel mixing.

Rogers, R. Clayton↗

NASA NTRS2007

Ring-Down Spectroscopy for Characterizing a CW Raman Laser

.A relatively simple technique for characterizing an all-resonant intracavity continuous-wave (CW) solid-state Raman laser involves the use of ring-down spectroscopy. As used here, characterizing signifies determining such parameters as threshold pump power, Raman gain, conversion efficiency, and quality factors (Q values) of the pump and Stokes cavity modes. Heretofore, in order to characterize resonant-cavity-based Raman lasers, it has usually been necessary to manipulate the frequencies and power levels of pump lasers and, in each case, to take several sets of measurements. In cases involving ultra-high-Q resonators, it also has been desirable to lock pump lasers to resonator modes to ensure the quality of measurement data. Simpler techniques could be useful. In the present ring-down spectroscopic technique, one infers the parameters of interest from the decay of the laser out of its steady state. This technique does not require changing the power or frequency of the pump laser or locking the pump laser to the resonator mode. The technique is based on a theoretical analysis of what happens when the pump laser is abruptly switched off after the Raman generation reaches the steady state. The analysis starts with differential equations for the evolution of the amplitudes of the pump and Stokes electric fields, leading to solutions for the power levels of the pump and Stokes fields as functions of time and of the aforementioned parameters. Among other things, these solutions show how the ring-down time depends, to some extent, on the electromagnetic energy accumulated in the cavity. The solutions are readily converted to relatively simple equations for the parameters as functions of quantities that can be determined from measurements of the time-dependent power levels. For example, the steady-state intracavity conversion efficiency is given by G1/G2 1 and the threshold power is given by Pin(G2/G1)2, where Pin is the steady-state input pump power immediately prior to abrupt switch-off, G1 is the initial rate of decay of the pump field, and G2 is the final rate of decay of the pump field. Hence, it is possible to determine all the parameters from a single ring-down scan, provided that the measurements taken in that scan are sufficiently accurate and complete.

Matsko, Andrey↗

NASA NTRS2014

On Parametric Sensitivity of Reynolds-Averaged Navier-Stokes SST Turbulence Model: 2D Hypersonic Shock-Wave Boundary Layer Interactions

Examined is sensitivity of separation extent, wall pressure and heating to variation of primary input flow parameters, such as Mach and Reynolds numbers and shock strength, for 2D and Axisymmetric Hypersonic Shock Wave Turbulent Boundary Layer interactions obtained by Navier-Stokes methods using the SST turbulence model. Baseline parametric sensitivity response is provided in part by comparison with vetted experiments, and in part through updated correlations based on free interaction theory concepts. A recent database compilation of hypersonic 2D shock-wave/turbulent boundary layer experiments extensively used in a prior related uncertainty analysis provides the foundation for this updated correlation approach, as well as for more conventional validation. The primary CFD method for this work is DPLR, one of NASA's real-gas aerothermodynamic production RANS codes. Comparisons are also made with CFL3D, one of NASA's mature perfect-gas RANS codes. Deficiencies in predicted separation response of RANS/SST solutions to parametric variations of test conditions are summarized, along with recommendations as to future turbulence approach.

Hypersonic↗

DOE OSTI2023

Anharmonic Exciton‐Phonon Coupling in Metal‐Organic Chalcogenides Hybrid Quantum Wells

Abstract In contrast to inorganic quantum wells, hybrid quantum wells (HQWs) based on metal‐organic semiconductors are characterized by relatively soft lattices, in which excitonic states can strongly couple to lattice phonons. Therefore, understanding the lattice's impact on exciton dynamics is essential for harnessing the optoelectronic potential of HQWs. Beyond 2D metal halide perovskites, layered metal‐organic chalcogenides (MOCs), which are an air‐stable, underexplored material class hosting room‐temperature excitons, can be exploited as photodetectors, light emitting devices, and ultrafast photoswitches. Here, the role of phonons in the optical transitions of the prototypical MOC [AgSePh] ∞ is elucidated. Impulsive stimulated Raman scattering (ISRS) allows the detection of coherent exciton oscillations driven by Fröhlich interaction with low‐energy optical phonons. Steady state absorption and Raman spectroscopies reveal a strong exciton‐phonon coupling (Huang‐Rhys parameter ≈1.7) and its anharmonicity, manifested as a nontrivial temperature‐dependent Stokes shift. The ab initio calculations support these observations, hinting at an anharmonic behavior of the low‐energy phonons <200 cm −1 . These results untangle complex exciton‐phonon interactions in MOCs, establishing an ideal testbed for room‐temperature many‐body phenomena.

36 MATERIALS SCIENCE↗

NASA NTRS1990

Progress of research on water vapor lidar

Research is described on several aspects of stimulated Raman scattering (SRS) of 532 nm laser light in H2, D2, and CH4. The goals of this work are to develop a more thorough understanding of nonlinear processes involving the Raman effect and four-wave mixing, and to find the best way to generate radiation at several wavelengths simultaneously, for lidar applications. Issues addressed are conversion efficiency, optimization of operating conditions (gas pressure, confocal parameter, etc.) and the distribution of output pulse energy over three Stokes components, the first anti-Stokes component, and the zeroth order (pump) wavelength. The described research and results constitute another step in the development of SRS applications for NASA's atmospheric lidar program.

Wilkerson, T. D.↗

NASA NTRS1991

The dynamics of dense and dilute clusters of drops evaporating in large, coherent vortices

The behavior of evaporating clusters of drops embedded into large, coherent vortices is described using a formulation which is valid for both dense and dilute clusters. Drops and gas interact both dynamically and thermodynamically. Dynamic coupling occurs through a force on the drops due to drag resulting from a slip velocity between the two phases. The net interaction force on the gas with drops is due to a source thrust from evaporation plus drag on each drop. The drag coefficient accounts for blowing from the drop surface. Thermodynamic coupling is a result of drop heating and evaporation. Limitations due to drop proximity on heating and evaporation are taken into account. The vortical motion of the drops in the cluster results in the formation of a core region devoid of drops at the center of the vortex, and a shell region containing the drops and surrounding the inner core. Results are presented showing the dependence of the evaporation time, the final to initial volume ratio and the final to initial shell thickness ratio upon the initial air/fuel mass ratio and as a function of the initial tangential velocities, upon the initial Stokes number, initial drop radius and initial outer cluster radius. Differences in behavior between and control parameters of dense and dilute clusters are pointed out by these new results. It is found that for dense clusters the final to initial volume ratio and final to initial shell thickness scale with the initial Stokes number, a new result which must be validated experimentally.

Bellan, J.↗

DOE OSTI2021

Probabilistic learning on manifolds constrained by nonlinear partial differential equations for small datasets

A novel extension of the Probabilistic Learning on Manifolds (PLoM) is presented. It makes it possible to synthesize solutions to a wide range of nonlinear stochastic boundary value problems described by partial differential equations (PDEs) for which a stochastic computational model (SCM) is available and which depend on a vector-valued random control parameter. The cost of a single numerical evaluation of this SCM is assumed to be such that only a limited number of points can be computed for constructing the training dataset (small data). Each point of the training dataset is made up of realizations from a vector-valued stochastic process (the stochastic solution) and the associated random control parameter on which it depends. The presented PLoM constrained by PDE allows for generating a large number of learned realizations of the stochastic process and its corresponding random control parameter. These learned realizations are generated so as to minimize the vector-valued random residual of the PDE in the mean-square sense. Appropriate novel methods are developed to solve this challenging problem. Three applications are presented. The first one is a simple uncertain nonlinear dynamical system with a nonstationary stochastic excitation. The second one concerns the 2D nonlinear unsteady Navier–Stokes equations for incompressible flows in which the Reynolds number is the random control parameter. Here, the last one deals with the nonlinear dynamics of a 3D elastic structure with uncertainties. The results obtained make it possible to validate the PLoM constrained by stochastic PDE but also provide further validation of the PLoM without constraint.

Machine learning↗

NASA NTRS1993

A comparative study of Full Navier-Stokes and Reduced Navier-Stokes analyses for separating flows within a diffusing inlet S-duct

A three-dimensional implicit Full Navier-Stokes (FNS) analysis and a 3D Reduced Navier Stokes (RNS) initial value space marching solution technique has been applied to a class of separated flow problems within a diffusing S-duct configuration characterized by vortex-liftoff. Both the FNS and the RNS solution technique were able to capture the overall flow physics of vortex lift-off, and gave remarkably similar results which agreed reasonably well with the experimental measured averaged performance parameters of engine face total pressure recovery and distortion. However, the Full Navier-Stokes and Reduced Navier-Stokes also consistently predicted separation further downstream in the M2129 inlet S-duct than was indicated by experimental data, thus compensating errors were present in the two Navier-Stokes analyses. The difficulties encountered in the Navier-Stokes separations analyses of the M2129 inlet S-duct center primarily on turbulence model issues, and these focused on two distinct but different phenomena, namely, (1) characterization of low skin friction adverse pressure gradient flows, and (2) description of the near wall behavior of flows characterized by vortex lift-off.

Anderson, B. H.↗

DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1987

NASA Ames Research Center's parabolized Navier-Stokes code - A critical evaluation of heat-transfer predictions

The viscous supersonic/hypersonic flow over a biconic configuration was numerically simulated using the NASA Ames Research Center's Parabolized Navier-Stokes (PNS) code, and results obtained for the effects of various computational parameters on the heat transfer are compared with experimental results. The PNS code is found to provide accurate results for heat transfer parameters using low smoothing coefficient values, assuming that the radial spacing is at some small value and that the body is simple in cross section. For the cases considered, the nominal value of the spatial step size was 0.05, and the grid density was 45 points in the meridional direction and 30 points in the radial direction.

Chaussee, Denny S.↗

NASA NTRS1992

Design of a wing shape for study of hypersonic crossflow transition in flight

Computational fluid dynamics methods were used in the design of a wing shape for study of hypersonic crossflow transition in flight. The flight experiment is to be performed on the delta wing of the first stage of a Pegasus launch vehicle as a piggy-back experiment to support boundary-layer stability code development and validation. The design goal is to obtain crossflow-induced transition at 20-40 percent of the chord for a flight Mach number of approximately six. The present paper describes the design and analysis process utilized to obtain desired glove shape. A variety of schemes were used in the design, ranging from simple empirical crossflow correlations to three-dimensional Navier-Stokes codes in conjunction with linear stability/N-factor computations. The sensitivity to various parameters, such as trajectory variations, allowable wing thickness, leading-edge radius and surface temperature, is also discussed.

Godil, A.↗

DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗

NASA NTRS1994

The Stokes line width and uncertainty relations

For a function given by contour integral the two types (conventions) of asymptotic representations are considered: the usual representation by asymptotic series in inverse powers of large parameters and the special division of contour integral in contributions of high and low saddle points. It is shown that the width of the recessive term formation zone (Stokes strip) in the second convention is determined by uncertainty relation and is much less than the zone width in the first convention. The reasons of such a difference is clarified. The results of the work are useful for understanding of formation region of the exponentially small process arising on the background of the strong one.

Nikishov, A. I.↗

NASA NTRS1993

Sensitivity derivatives for advanced CFD algorithm and viscous modelling parameters via automatic differentiation

The computational technique of automatic differentiation (AD) is applied to a three-dimensional thin-layer Navier-Stokes multigrid flow solver to assess the feasibility and computational impact of obtaining exact sensitivity derivatives typical of those needed for sensitivity analyses. Calculations are performed for an ONERA M6 wing in transonic flow with both the Baldwin-Lomax and Johnson-King turbulence models. The wing lift, drag, and pitching moment coefficients are differentiated with respect to two different groups of input parameters. The first group consists of the second- and fourth-order damping coefficients of the computational algorithm, whereas the second group consists of two parameters in the viscous turbulent flow physics modelling. Results obtained via AD are compared, for both accuracy and computational efficiency with the results obtained with divided differences (DD). The AD results are accurate, extremely simple to obtain, and show significant computational advantage over those obtained by DD for some cases.

Green, Lawrence L.↗

NASA NTRS1975

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

NASA NTRS1986

Effects of droplet interactions on droplet transport at intermediate Reynolds numbers

Effects of droplet interactions on drag, evaporation, and combustion of a planar droplet array, oriented perpendicular to the approaching flow, are studied numerically. The three-dimensional Navier-Stokes equations, with variable thermophysical properties, are solved using finite-difference techniques. Parameters investigated include the droplet spacing, droplet Reynolds number, approaching stream oxygen concentration, and fuel type. Results are obtained for the Reynolds number range of 5 to 100, droplet spacing from 2 to 24 diameters, oxygen concentrations of 0.1 and 0.2, and methanol and n-butanol fuels. The calculations show that the gasification rates of interacting droplets decrease as the droplet spacings decrease. The reduction in gasification rates is significant only at small spacings and low Reynolds numbers. For the present array orientation, the effects of interactions on the gasification rates diminish rapidly for Reynolds numbers greater than 10 and spacings greater than 6 droplet diameters. The effects of adjacent droplets on drag are shown to be small.

Shuen, J. S.↗

NASA NTRS1987

Effects of droplet interactions on droplet transport at intermediate Reynolds numbers

Effects of droplet interactions on drag, evaporation, and combustion of a planar droplet array, oriented perpendicular to the approaching flow, are studied numerically. The three-dimensional Navier-Stokes equations, with variable thermophysical properties, are solved using finite-difference techniques. Parameters investigated include the droplet spacing, droplet Reynolds number, approaching stream oxygen concentration, and fuel type. Results are obtained for the Reynolds number range of 5 to 100, droplet spacings from 2 to 24 diameters, oxygen concentrations of 0.1 and 0.2, and methanol and n-butanol fuels. The calculations show that the gasification rates of interacting droplets decrease as the droplet spacings decrease. The reduction in gasification rates is significant only at small spacings and low Reynolds numbers. For the present array orientation, the effects of interactions on the gasification rates diminish rapidly for Reynolds numbers greater than 10 and spacings greater than 6 droplet diameters. The effects of adjacent droplets on drag are shown to be small.

Shuen, Jian-Shun↗