SEARCH · Engineering Papers

Results for “kappa”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1974

Interplanetary diffusion coefficients for cosmic rays

Information on the cosmic-ray diffusion coefficient, kappa, derived from near-earth observations of the solar modulation of galactic electron fluxes and from the near-earth power spectra of the interplanetary magnetic field, has been used to study the heliocentric radial dependence of kappa, and to derive limits on the spatial extent of the solar modulation region. Representing kappa, as a separable function of radius r and rigidity, and assumming kappa(r) proportional to r to the n-th power, we can place a limit on the power law exponent, n not greater than 1.2. The distance of the modulation boundary is a function of n, and, e.g., for n = 0, falls into the range of 6-25 AU.

Cummings, A. C.↗

NASA NTRS1976

Consequences of using nonlinear particle trajectories to compute spatial diffusion coefficients

The propagation of charged particles through interstellar and interplanetary space has often been described as a random process in which the particles are scattered by ambient electromagnetic turbulence. In general, this changes both the magnitude and direction of the particles' momentum. Some situations for which scattering in direction (pitch angle) is of primary interest were studied. A perturbed orbit, resonant scattering theory for pitch-angle diffusion in magnetostatic turbulence was slightly generalized and then utilized to compute the diffusion coefficient for spatial propagation parallel to the mean magnetic field, Kappa. All divergences inherent in the quasilinear formalism when the power spectrum of the fluctuation field falls off as K to the minus Q power (Q less than 2) were removed. Various methods of computing Kappa were compared and limits on the validity of the theory discussed. For Q less than 1 or 2, the various methods give roughly comparable values of Kappa, but use of perturbed orbits systematically results in a somewhat smaller Kappa than can be obtained from quasilinear theory.

Goldstein, M. L.↗

NASA NTRS1984

Asymmetric instabilities in buoyancy-driven flow in a tall vertical annulus

Linear stability of the one-dimensional flow between infinite vertical coaxial cylinders induced by heating the inner cylinder is considered for various ratios (kappa) of the inner radius to the outer radius, and for Prandtl numbers (P) appropriate to air and water. For air with P = 0.71 the least stable disturbance is nonaxisymmetric for kappa less than 0.44 and is axisymmetric for kappa greater than 0.44, and in either case the instability is due to the action of the shear forces. For P = 3.5, the situation is similar, except that the asymmetric shear mode is superseded by an axisymmetric instability driven by buoyancy forces for kappa = 0.03-0.16. Wave speeds, wavenumbers, and critical Grashof numbers for these cases and for the case of zero Prandtl number are given.

Mcfadden, G. B.↗

NASA NTRS1990

Effects of dynamic long-period ocean tides on changes in earth's rotation rate

As a generalization of the zonal response coefficient first introduced by Agnew and Farrell (1978), the zonal response function kappa of the solid earth-ocean system is defined as the ratio, in the frequency domain, of the tidal change in earth's rotation rate to the tide-generating potential. Amplitudes and phases of kappa for the monthly, fortnightly, and nine-day lunar tides are estimated from 2 1/2 years of VLBI UT1 observations, corrected for atmospheric angular momentum effects using NMC wind and pressure series. Using the dynamic ocean tide model of Dickman (1988, 1989), amplitudes and phases of kappa for an elastic earth-ocean system are predicted. The predictions confirm earlier results which found that dynamic effects of the longer-period ocean tides reduce the amplitude of kappa by about 1 percent.

Nam, Young↗

NASA NTRS1993

The Th and U abundances in chondritic meteorites

We present new analyses of Th-232/U-238 in CI and CM meteorites. The relative abundance of these nuclides is important in estimates of the age of r-process elements. The cosmochronology based upon the Th-232/U-238 ratio (kappa) depends on the precise determinations of these two different elements in meteorites and on the production ratios. Both parameters are subject to substantial errors. Recent recalculations of this chronology have used selected values from compilations but do not adequately address the errors in terms of a reliable data base. Morgan and Lovering provided extensive neutron activation analyses for ordinary chondrites which yield an average kappa of 3.6 +/- 0.4. Their work on carbonaceous chondrites showed a wide range in kappa from 2 to 6. More recent investigations by isotopic dilution have established the following: (1) highly variable kappa from 2.7 to 11 in Allende Ca-Al-rich inclusions and a value of 3.6 in the Orgueil CI1 chondrite; (2) a range from 2.71 to 6.63 for 7 L-type chondrites and a range from 2.7 to 4.4 for 6 L, H, and LL chondrites. A further investigation of this subject matter is presented.

Chen, J. H.↗

NASA NTRS1994

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗

NASA NTRS1995

Particle chaos and pitch angle scattering

Pitch angle scattering is a factor that helps determine the dawn-to-dusk current, controls particle energization, and it has also been used as a remote probe of the current sheet structure. Previous studies have interpreted their results under the exception that randomization will be greatest when the ratio of the two timescales of motion (gyration parallel to and perpendicular to the current sheet) is closet to one. Recently, the average expotential divergence rate (AEDR) has been calculated for particle motion in a hyperbolic current sheet (Chen, 1992). It is claimed that this AEDR measures the degree of chaos and therefore may be thought to measure the randomization. In contrast to previous expectations, the AEDR is not maximized when Kappa is approximately equal to 1 but instead increases with decreasing Kappa. Also contrary to previous expectations, the AEDR is dependent upon the parameter b(sub z). In response to the challenge to previous expectations that has been raised by this calculation of the AEDR, we have investigated the dependence of a measure of particle pitch angle scattering on both the parameters Kappa and b(sub z). We find that, as was previously expected, particle pitch angle scattering is maximized near Kappa = 1 provided that Kappa/b(sub z) greater than 1. In the opposite regime, Kappa/b(sub z) less than 1, we find that particle pitch angle scattering is still largest when the two timescales are equal, but the ratio of the timescales is proportional to b(sub z). In this second regime, particle pitch angle scattering is not due to randomization, but is instead due to a systematic pitch angle change. This result shows that particle pitch angle scattering need not be due to randomization and indicates how a measure of pitch angle scattering can exhibit a different behavior than a measure of chaos.

Burkhart, G. R.↗

NASA NTRS1997

The Spectral and Statistical Properties of Turbulence Generated by a Vortex/Blade-Tip Interaction

The perpendicular interaction of a streamwise vortex with the tip of a lifting blade was studied in incompressible flow to provide information useful to the accurate prediction of helicopter rotor noise and the understanding of vortex dominated turbulent flows. The vortex passed 0.3 chord lengths to the suction side of the blade tip, providing a weak interaction. Single and two-point turbulence measurements were made using sub-miniature four sensor hot-wire probes 15 chord lengths downstream of the blade trailing edge; revealing the mean velocity and Reynolds stress tensor distributions of the turbulence, as well as its spanwise length scales as a function of frequency. The single point measurements show the flow downstream of the blade to be dominated by the interaction of the original tip vortex and the vortex shed by the blade. These vortices rotate about each other under their mutual induction, winding up the turbulent wakes of the blades. This interaction between the vortices appears to be the source of new turbulence in their cores and in the region between them. This turbulence appears to be responsible for some decay in the core of the original vortex, not seen when the blade is removed. The region between the vortices is not only a region of comparatively large stresses, but also one of intense turbulence production. Velocity autospectra measured near its center suggests the presence quasi-periodic large eddies with axes roughly parallel to a line joining the vortex cores. Detailed two-point measurements were made on a series of spanwise cuts through the flow so as to reveal the turbulence scales as they would be seen along the span of an intersecting airfoil. The measurements were made over a range of probe separations that enabled them to be analyzed not only in terms of coherence and phase spectra but also in terms of wave-number frequency (kappa-omega) spectra, computed by transforming the measured cross-spectra with respect to the spanwise separation of the probes. These data clearly show the influence of the coherent eddies in the spiral wake and the turbulent region between the cores. These eddies produce distinct peaks in the upwash velocity kappa-omega spectra, and strong anisotropy manifested both in the decay of the kappa-omega spectrum at larger wave-numbers and in differences between the kappa-omega spectra of different components. None of these features are represented in the von Karman spectrum for isotropic turbulence that is often used in broadband noise computations. Wave-number frequency spectra measured in the cores appear to show some evidence that the turbulence outside sets tip core waves, as has previously been hypothesized. These spectra also provide for the first time a truly objective method for distinguishing velocity fluctuations produced by core wandering from other motions.

Devenport, William J.↗

NASA NTRS2013

Characterization of the Gaseous Companion k Andromedae B* New Keck and LBTI High-contrast Observations

Context. We previously reported the direct detection of a low mass companion at a projected separation of 55+/-2 astronomical units around the B9 type star kappa Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. Aims. We present new angular differential imaging (ADI) images of the system at 2.146 (K(sub s)), 3.776 (L'), 4.052 (NB 4.05) and 4.78 micrometers (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We aim to determine the near-infrared spectral energy distribution (SED) of the companion and use it to characterize the object. Methods. We used analysis methods adapted to ADI to extract the companion flux. We compared the photometry of the object to reference young/old objects and to a set of seven PHOENIX-based atmospheric models of cool objects accounting for the formation of dust. We used evolutionary models to derive mass estimates considering a wide range of plausible initial conditions. Finally, we used dedicated formation models to discuss the possible origin of the companion. Results. We derive a more accurate J = 15.86 +/- 0.21, H = 14.95 +/- 0.13, K(sub s) = 14.32 +/- 0.09 mag for kappa And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at K(sub s) and L' band. We derive NB 4.05 = 13.0 +/- 0.2 and M' = 13.3 +/- 0.3 mag and estimate Log(base 10)(L/solar luminosity) = −3.76 +/- 0.06. Atmospheric models yield T(sub eff) = 1900(+100/−200) K. They do not set constrains on the surface gravity. "Hot-start" evolutionary models predict masses of 14(+25/−2) Jupiter mass based on the luminosity and temperature estimates, and considering a conservative age range for the system (30(+120/−10) million years). "warm-start" evolutionary tracks constrain the mass to M greater than or equal to 11 Jupiter mass. Conclusions. The mass of kappa Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of kappa and b.

spectral energy distribution (SED)↗

NASA NTRS2018

Magnetospheric Multiscale Observations of Turbulent Magnetic and Electron Velocity Fluctuations in Earth's Magnetosheath Downstream of a Quasi-Parallel Bow Shock

We present statistical single-spacecraft observations of magnetic and electron velocity fluctuations in Earth's magnetosheath, likely in the vicinity of the magnetopause, downstream of a bow shock immersed in quasi-parallel interplanetary magnetic field conditions, a situation conducive to plasma turbulence in the downstream flow. These fluctuations exhibit scale-dependent behavior, wherein histograms of their Partial Variance of Increments (PVIB or PVIV(sub e)) demonstrate highly non-Gaussian forms at small scales and are reasonably well-described by kappa distributions, albeit with fitted values of the kappa parameter only slightly larger than 1.5, exemplifying their power law nature at large values of PVI. At larger scales, the PVI histograms lose their non-Gaussian nature and are well described by both Gaussian and kappa distributions with large values of the kappa parameter. The PVI histograms furthermore exhibit kurtosis that increases with decreasing scale, a characteristic that is much more prominent in the magnetic fluctuations than in the electron velocity fluctuations. This feature that is not yet explained. In both cases, the results are characteristic of turbulent intermittency.

Pollock, C. J.↗

DOE OSTI2024

Non-Thermal Solar Wind Electron Velocity Distribution Function

The quiet-time solar wind electrons feature non-thermal characteristics when viewed from the perspective of their velocity distribution functions. They typically have an appearance of being composed of a denser thermal “core” population plus a tenuous energetic “halo” population. At first, such a feature was empirically fitted with the kappa velocity space distribution function, but ever since the ground-breaking work by Tsallis, the space physics community has embraced the potential implication of the kappa distribution as reflecting the non-extensive nature of the space plasma. From the viewpoint of microscopic plasma theory, the formation of the non-thermal electron velocity distribution function can be interpreted in terms of the plasma being in a state of turbulent quasi-equilibrium. Such a finding brings forth the possible existence of a profound inter-relationship between the non-extensive statistical state and the turbulent quasi-equilibrium state. The present paper further develops the idea of solar wind electrons being in the turbulent equilibrium, but, unlike the previous model, which involves the electrostatic turbulence near the plasma oscillation frequency (i.e., Langmuir turbulence), the present paper considers the impact of transverse electromagnetic turbulence, particularly, the turbulence in the whistler-mode frequency range. It is found that the coupling of spontaneously emitted thermal fluctuations and the background turbulence leads to the formation of a non-thermal electron velocity distribution function of the type observed in the solar wind during quiet times. This demonstrates that the whistler-range turbulence represents an alternative mechanism for producing the kappa-like non-thermal distribution, especially close to the Sun and in the near-Earth space environment.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Dark Energy Survey Year 3 results: Imprints of cosmic voids and superclusters in the Planck CMB lensing map

ABSTRACT The CMB lensing signal from cosmic voids and superclusters probes the growth of structure in the low-redshift cosmic web. In this analysis, we cross-correlated the Planck CMB lensing map with voids detected in the Dark Energy Survey Year 3 (Y3) data set (∼5000 deg2), expanding on previous measurements that used Y1 catalogues (∼1300 deg2). Given the increased statistical power compared to Y1 data, we report a 6.6σ detection of negative CMB convergence (κ) imprints using approximately 3600 voids detected from a redMaGiC luminous red galaxy sample. However, the measured signal is lower than expected from the MICE N-body simulation that is based on the ΛCDM model (parameters Ωm = 0.25, σ8 = 0.8), and the discrepancy is associated mostly with the void centre region. Considering the full void lensing profile, we fit an amplitude $A_{\kappa }=\kappa _{{\rm DES}}/\kappa _{{\rm MICE}}$ to a simulation-based template with fixed shape and found a moderate 2σ deviation in the signal with Aκ ≈ 0.79 ± 0.12. We also examined the WebSky simulation that is based on a Planck 2018 ΛCDM cosmology, but the results were even less consistent given the slightly higher matter density fluctuations than in MICE. We then identified superclusters in the DES and the MICE catalogues, and detected their imprints at the 8.4σ level; again with a lower-than-expected Aκ = 0.84 ± 0.10 amplitude. The combination of voids and superclusters yields a 10.3σ detection with an Aκ = 0.82 ± 0.08 constraint on the CMB lensing amplitude, thus the overall signal is 2.3σ weaker than expected from MICE.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Magnon heat transport in a two-dimensional Mott insulator

Whether or not anomalies in the thermal conductivity in insulating cuprates can be attributed to antiferromagnetic order and magnons in a 2D Mott insulator remains an intriguing open question. To shed light on this issue, we investigate the thermal conductivity $\kappa$ and its relationship with the specific heat $c_v$ in the half-filled 2D single-band Hubbard model, using the numerically exact determinant quantum Monte Carlo algorithm and maximum entropy analytic continuation. At low temperatures where the charge degrees of freedom are gapped-out and $c_v$ exhibits a clear magnon peak, we observe that thermal conductivity $\kappa$ also tends to form a peak at similar temperatures. Reducing temperature further produces a sharp upturn in $\kappa$, associated with an increasing mean-free path. Here we identify this as the high-temperature side of the anomalous peak in insulating cuprates, where the mean-free path eventually is cut-off by other scattering effects, including phonons, disorder, and physical size. Different scattering effects in our model are identified and analyzed in the thermal diffusivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗

DOE OSTI2021

Measurement of the Higgs boson production rate in association with top quarks in final states with electrons, muons, and hadronically decaying tau leptons at $\sqrt{s} =$ 13 TeV

The rate for Higgs (${\mathrm{H}} $) bosons production in association with either one (${\mathrm{t}} {\mathrm{H}} $) or two (${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $) top quarks is measured in final states containing multiple electrons, muons, or tau leptons decaying to hadrons and a neutrino, using proton–proton collisions recorded at a center-of-mass energy of $13\,\text {Te}\text {V} $ by the CMS experiment. The analyzed data correspond to an integrated luminosity of 137$\,\text {fb}^{-1}$. The analysis is aimed at events that contain ${\mathrm{H}} \rightarrow {\mathrm{W}} {\mathrm{W}} $, ${\mathrm{H}} \rightarrow {\tau } {\tau } $, or ${\mathrm{H}} \rightarrow {\mathrm{Z}} {\mathrm{Z}} $ decays and each of the top quark(s) decays either to lepton+jets or all-jet channels. Sensitivity to signal is maximized by including ten signatures in the analysis, depending on the lepton multiplicity. The separation among ${\mathrm{t}} {\mathrm{H}} $, ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and the backgrounds is enhanced through machine-learning techniques and matrix-element methods. The measured production rates for the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ and ${\mathrm{t}} {\mathrm{H}} $ signals correspond to $0.92 \pm 0.19\,\text {(stat)} ^{+0.17}_{-0.13}\,\text {(syst)} $ and $5.7 \pm 2.7\,\text {(stat)} \pm 3.0\,\text {(syst)} $ of their respective standard model (SM) expectations. The corresponding observed (expected) significance amounts to 4.7 (5.2) standard deviations for ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $, and to 1.4 (0.3) for ${\mathrm{t}} {\mathrm{H}} $ production. Assuming that the Higgs boson coupling to the tau lepton is equal in strength to its expectation in the SM, the coupling $y_{{\mathrm{t}}}$ of the Higgs boson to the top quark divided by its SM expectation, $\kappa _{{\mathrm{t}}}=y_{{\mathrm{t}}}/y_{{\mathrm{t}}}^{\mathrm {SM}}$, is constrained to be within $-0.9< \kappa _{{\mathrm{t}}}< -0.7$ or $0.7< \kappa _{{\mathrm{t}}}< 1.1$, at 95% confidence level. This result is the most sensitive measurement of the ${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $ production rate to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1979

Fermi liquid viscosity in a finite geometry

Forced flow of a Fermi liquid is studied for a cell geometry consisting of two planes with a separation on the order of the mean free path. An approximate transport equation is used to derive an integral equation for the velocity profile, which is solved numerically. Results for the total flux through the cell, which determines the dissipation, are given as a function of the Knudsen number kappa (ratio of cell thickness to mean free path). Effects of specular reflection at the boundaries are considered. It is found that the dissipation has a minimum at kappa of 1/2, and behaves linearly for kappa not less than 3. Implications for present experimentation are discussed.

Jaffe, J. E.↗

NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗

NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme for the 3-D Navier-Stokes equations and a two equation model of turbulence

An LU implicit multigrid algorithm is developed to calculate 3-D compressible viscous flows. This scheme solves the full 3-D Reynolds-Averaged Navier-Stokes equation with a two-equation kappa-epsilon model of turbulence. The flow equations are integrated by an efficient, diagonally inverted, LU implicit multigrid scheme while the kappa-epsilon equations are solved, uncoupled from the flow equations, by a block LU implicit algorithm. The flow equations are solved within the framework of the multigrid method using a four-grid level W-cycle, while the kappa-epsilon equations are iterated only on the finest grid. This treatment of the Reynolds-Averaged Navier-Stokes equations proves to be an efficient method for calculating 3-D compressible viscous flows.

Yokota, Jeffrey W.↗