Search NASA⌕ Search

SEARCH · Search NASA

Results for “Conditional diffusion models”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Low Stretch Diffusion Flames Over a Solid Fuel

A unique new way to study low gravity flames in normal gravity has been developed. To study flame structure and extinction characteristics in low stretch environments, a normal gravity low-stretch diffusion flame is generated using a cylindrical PMMA sample of varying large radii. Foutch and T'ien used the radiative loss as well as a densimetric Froude number to characterize the blowoff (small Da) and quenching extinction (large Da) boundaries in stagnation-point diffusion flames under various convective conditions. An important conclusion of this study was that the shape and location of the extinction boundary, as well as a number of important flame characteristics, were almost identical for the buoyant, forced, and mixed convective environments they modeled. This theory indicates it should be possible to understand a material's burning characteristics in the low stretch environment of spacecraft (induced by fans and crew movements) by understanding its burning characteristics in an equivalent Earth-based stretch environment (induced by normal gravity buoyancy). Similarly, the material's burning characteristics in Lunar or Martian stretch environments (induced by partial gravity buoyancy) can be assessed. Equivalent stretch rates can be determined as a function of gravity, imposed flow, and geometry. A generalized expression for stretch rate which captures mixed convection includes both buoyant and forced stretch is defined as a = a(sub f) ((1 + (a(sub b))exp 2/(a(sub b))exp 2))exp 1/2. For purely buoyant flow, the equivalent stretch rate is a(sub b) = [(rho(exp e)-rho(exp *)/rho(sub e)][g/R](exp 1/2). For purely forced flow, the equivalent stretch rate is characterized by either a(sub f)= 2U(sub infinity)/R for a cylinder, or a(sub f)=U(sub jet)/d(sub jet) for a jet impinging on a planar surface. In these experiments, the buoyant stretch is varied through R, the radius of curvature, but the buoyant stretch could also be varied through g, the gravity level. In this way the effect of partial gravity, such as those found on the Moon (1/6 g) or Mars (1/3 g) can be captured in the definition of flame stretch.

Olson, S. L.↗
NASA NTRS2001

ALEXIS Observations of the Monogem Ring Supernova Remnant

The subject grant is for the analysis of ALEXIS observations of the Monogem Ring supernova remnant using the diffuse all-sky maps produced from the ALEXIS all-sky survey. The work is to produce ratio maps of the three energy bands provided by ALEXIS, analyze the ratio data to constrain the intervening neutral hydrogen column density and the temperature and elemental abundances of the X-ray emitting gas, compare the structure to that observed in the ROSAT maps, and incorporate the results into current supenova remnant evolution models. The work outlined above has been significantly delayed since the ALEXIS diffuse all-sky maps took longer to produce than anticipated. Unfortunately, the ALEXIS satellite suffered a failure of the Pegasus launch vehicle which left the satellite in a partially functioning condition. The attitude control system of the spacecraft was unable to operate as planned and this has greatly increased the complexity of the aspect solution. Our colleagues at Los Alamos have made progress in producing these maps and are nearing completion of the final maps. However, the quality of the data have been significantly compromised by the overall lower exposure due to the spacecraft problems and the higher background of the micro-channel plate detectors. We have compared the ALEXIS and ROSAT maps of this region of the sky and there is no obvious signal in the ALEXIS maps of the Monogem Ring. We are now exploring correlation techniques to determine if there is indeed a faint signal in the ALEXIS maps. Although, the project has been a disappointment so far, the data may still provide a valuable lower limit on the neutral hydrogen column density. This is a far cry from our original intentions, but would still be valuable science. Given the large delays in producing the ALEXIS sky maps, this work will continue past the end of the grant period.

Plucinsky, Paul↗
NASA NTRSDate not supplied

NASA HECC Geometry and Performance Review Part 1: Validation of a Computational Model for the Vaneless Diffuser Configuration with As-Manufactured Impeller Geometry

An investigation of the NASA High Efficiency Centrifugal Compressor (HECC) vaneless diffuser configuration was performed. This multipart investigation focused on validating the computational model against experimental data. The validated model was then used to explore the effects of observed geometric differences between the as-manufactured and design-intent impellers. Lastly, the validated model was used for comparison against an experimental tip clearance study. Details regarding each of the parts are detailed below: A multipart investigation of the NASA High Efficiency Centrifugal Compressor (HECC) vaneless diffuser configuration was performed. Part I of the investigation focused on describing the HECC vaneless computational model and its validation against the experimental data published in GT2023-103128. To generate the model, a verification of the manufactured impeller was performed and was confirmed to match the impeller CAD. Section data from the impeller CAD, referred to as As-Manufactured, were used in the model. Steady RANS simulations using a multiblock structured mesh were used to conduct a mesh convergence study. The resultant mesh was then used for simulations at both design-speed and off-design speed conditions. The simulations compared well to data in both 1D and 2D terms. Based on the observed results the model was deemed validated and was used in the latter two parts of the multipart investigation.

centrifugal compressor↗
NASA NTRS1997

Wave Driven Exothermic Heating in the Mesopause Region

A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.

Hickey, Michael P.↗
DOE OSTI2025

On the Efficacy of Repeat Voltage Holds for Conditioning and Calendar Life Testing of Graphite and Silicon Cells

Voltage-hold (V-hold) protocols have shown promise toward calendar lifetime analysis of cells with graphite (Gr) and silicon (Si) anodes. In this work, repeat V-holds are performed on Gr and Si cells paired with lithium iron phosphate cathodes to delineate their beneficial role in formation and conditioning. We find that V-hold at the top of charge supplements constant current cycling in conditioning the cell to higher capacities for both Gr and Si cells after the first V-hold. A reduced order model provides the irreversible capacity proportions of each V-hold. With each repeat V-hold, parasitic loss of lithium to the solid electrolyte interphase (SEI) decreases on both Gr and Si cells. Gr cells show the square-root-of-time capacity loss behavior within 200 h of V-hold, indicative of its fast relaxation and low impact of reference performance test cycles on the SEI growth. Lifetime estimates from repeat V-holds on Gr can reach years. Si exhibits longer transition times from kinetic to diffusion-limited SEI growth, evidenced by the 400 h and 200 h holds showing square-root-of-time and linear behavior, respectively. Lifetime predictions from repeat V-holds on Si only reach 1–2 months, highlighting its limitations. Recommended duration of V-holds for Si cells should be ≥400 h.

25 ENERGY STORAGE↗
NASA NTRS1994

Observations of Methane and Ethylene Diffusion Flames Stabilized Around a Blowing Porous Sphere Under Microgravity Conditions

This paper presents the experimental and theoretical results for expanding methane and ethylene diffusion flames in microgravity. A small porous sphere made from a low-density and low-heat-capacity insulating material was used to uniformly supply fuel at a constant rate to the expanding diffusion flame. A theoretical model which includes soot and gas radiation is formulated but only the problem pertaining to the transient expansion of the flame is solved by assuming constant pressure infinitely fast one-step ideal gas reaction and unity Lewis number. This is a first step toward quantifying the effect of soot and gas radiation on these flames. The theoretically calculated expansion rate is in good agreement with the experimental results. Both experimental and theoretical results show that as the flame radius increases, the flame expansion process becomes diffusion controlled and the flame radius grows as gamma t. Theoretical calculations also show that for a constant fuel mass injection rate a quasi-steady state is developed in the region surrounded by the flame and the mass flow rate at any location inside this region equals the mass injection rate.

Atreya, Arvind↗
DOE OSTI2025

Simulations of the churning mode: Toroidally symmetric plasma convection and turbulence around the X-points in a snowflake divertor

Using a reduced MHD model, extended to include field-aligned thermal conduction, we present numerical simulations of the churning mode (CM): a toroidally symmetric, non-linear plasma vortex in the vicinity of the null points in a snowflake (SF) divertor [D. D. Ryutov et al., Phys. Scr. 89, 088002 (2014)]. Simulations are carried out across a range of inter-null separations, d xx ⁠, and inter-null orientations, θ⁠, primarily in conditions relevant to the MAST-U tokamak. We find that, when d xx is small, the CM induces additional transport across the X-points when β pm ≳8%, where β pm is the ratio of the plasma pressure in the null region to poloidal magnetic pressure at the midplane. This transport also increases approximately linearly as d xx is reduced. A diffusive model of this transport is shown to predict the total transport across the null points, where diffusion coefficients of up to ∼ 10 2 m 2 s −1 centered on a small region around the X-points are used. However, the CM also results in significant changes to the flux surfaces in the null region, which is not captured by this diffusive model. The changes in magnetic geometry mean the fractional exhaust power delivered to each divertor leg is highly sensitive to β pm ⁠, d xx ⁠, and θ⁠. For small values of θ⁠, the CM can induce a change in topology, redirecting exhaust power from a secondary divertor leg on the high field side to one on the low field side. Similar behavior is found in the fraction of exhaust power going to the inner and outer divertor. Such changes in the flux surfaces may not be captured by Grad–Shafranov solvers and so may be a source of error in the magnetic reconstruction of SF experiments. We consistently find that the fractional exhaust power going to a secondary divertor leg on the high field side is small, consistent with SF experiments.

Computer simulation↗
NASA NTRS1973

A preliminary sensitivity analysis of the coupled diffusion and chemistry model

The sensitivity of the coupled chemistry/diffusion model's outputs to a wide range of variation of the model's independent variables has been investigated. It is shown that the efficiency with which the now catalytic cycle destroys ambient O3 is extremely sensitive to the amount of NO emitted and to the relative rates of turbulent diffusion and chemical reactions. For representative conditions in the stratosphere, a tenfold variation of either the turbulence intensity or the reaction rate constant or the source strength can vary the efficiency from 1% to 50%. If the duration of Phase 3 is a significant fraction of the total residence time of the plume, then these efficiency variations can alter O3 depletion rates by more than a factor of two. These results, therefore, point toward those variables which must be accurately defined or measured if one is to adequately predict the effect of SST operations on the ambient inventory of O3 in the lower stratosphere.

Hilst, G. R.↗
NASA NTRS1977

Auroral energy deposition and neutral composition changes observed simultaneously by ESRO 4 and AE-C at different altitudes

Neutral composition data obtained simultaneously from ESRO 4 and AE-C during geomagnetically disturbed conditions at different altitudes (160 and 230 km) are used to investigate the atmospheric response to geomagnetic storms and to infer information regarding the excitation mechanism. The data are compared with a theoretical model that estimates the composition effects in terms of wind induced diffusion. A parametric study was conducted bearing on the influence of energy deposition at different altitudes and with varying latitudinal extent. In one of the observed events the composition effects at 160 km are substantially smaller than at 230 km for which we inferred by comparison with the theory that the energy mainly was deposited at 150 km altitude over a wide latitude range. Another event required energy deposition at somewhat lower altitudes near 120 km with a more localized energy source. Significant variations of the turbopause level were not necessary to explain the observed variations.

Trinks, H.↗
NASA NTRS1981

Fundamentals of dendritic solidification. I - Steady-state tip growth. II - Development of sidebranch structure

Systematic measurements of dendrite tip radius and growth velocity in succinonitrile reveal that consideration of dendrite tip stability should be incorporated into the heat transfer theory to determine the steady-state dendritic growth condition. The dendritic stability criterion measured is 2 alpha d0/VR squared = 0.0195, where V is the dendritic growth velocity, R is the dendritic tip radius, alpha is the liquid thermal diffusivity, and d0 is a capillary length defined in the text. Several dendritic stability models are reviewed and discussed in comparison to the present experimental results.

Huang, S.-C.↗
NASA NTRS1999

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid interface shape and convection patterns during solidification of pure Al and Al-4 wt. % Cu alloy. Steady state calculations were performed for different gravity levels and orientations. Furthermore, transient modeling was undertaken to investigate effect of the solidification velocity. The furnace configuration used in this analysis is the proposed International Space Station Furnace. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. The Solidification model was adopted from previous work and was based on the finite element code FIDAP. Thermocouple data and quenched interface shape from a pure Al sample flown on the Life Sciences and Microgravity Spacelab (LMS) mission, July 1996 was used for model validation. Good agreement was obtained between the predicted interface shape and that measured from the quenched LMS sample. It was found that the imposed temperature boundary condition must be sufficiently smooth for the heat fluxes in the model to be self-consistent. The model predicted that alloy sample Al-Cu is more sensitive to variations in the gravity level compared to pure Al. For Al-Cu alloy, solute diffusivity is approximately 10(exp 4) times smaller than thermal diffusivity of pure Al. Hence, in a microgravity environment the weak convection has no measurable effect on the heat fluxes, but is still strong enough to affect concentration distribution. Since the alloy melting temperature is determined by concentration the interface shape depends on the level of convection. This was found not to be the case for pure metal.

Bune, Andris V.↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2025

Ab Initio Design of High-Entropy Thermal/ Environmental Barrier Coatings

Next generation thermal/environmental barrier coatings (TEBC) require carefully balancing various properties including phase stability, thermal conductivity, coefficient of thermal expansion (CTE), mechanical properties, and resistance against hot corrosion and water vapor recession. This work mainly focuses on rapid design of cost-effective high entropy rare-earth disilicates and aluminum garnets to protect SiC-based ceramic matrix composites and nickel-based superalloys in the hot section of gas turbine engines using density functional theory methods. Our calculations identify several low-cost high entropy TEBC exhibiting ultralow thermal conductivity at 1500 K and desirable CTE while maintaining good mechanical properties, including Er1/2Y3/4Yb3/4Si2O7, Gd1/4Er1/4Y3/4Yb3/4Si2O7, Eu1/4Er1/4Y3/4Yb3/4Si2O7, and (Y1/4Gd1/4Er1/4Yb1/4)3Al5O12. This work also aims to gain fundamental understanding of oxygen diffusion in model disilicates. Minimizing oxidizer (such as water vapor and oxygen) permeability through the EBC layer can significantly decrease the growth rate of thermally grown oxide and extend the service life of the coating system. Oxygen diffusion mechanisms including formation energy of defects under varying oxygen conditions and defect migration energy barriers will be presented.

coefficient of thermal expansion↗
NASA NTRS1982

Importance of heterogeneous processes to tropospheric chemistry - Studies with a one-dimensional model

A one-dimensional, time-dependent model of tropospheric air composition is developed which incorporates several heterogeneous physical and chemical processes. The model includes the interaction of gases, aerosols, and hydrometeors through the physical mechanisms of nucleation, condensation, evaporation, coagulation, coalescence, and deliquescence. Precipitation, sedimentation, and dry deposition act to remove material from the atmosphere, while chemical transformations occur both in the vapor and the condensed phases. The model also incorporates the sources and vertical diffusion of gases and particles, as well as changes in the solar intensity caused by light-scattering from aerosols and clouds. Preliminary simulations made using this model indicate that rainout and washout processes strongly influence the distributions of tropospheric gases and aerosols under certain conditions.

Turco, R. P.↗
NASA NTRS2010

Ares I Reaction Control System Propellant Feedline Decontamination Modeling

The objective of the work presented here is to quantify the effects of purge gas temperature, pressure, and mass flow rate on Hydrazine (Hz) decontamination rates of the Ares I Roll Control System and Reaction Control System. A survey of experts in this field revealed the absence of any decontamination rate prediction models. Three basic decontamination methods were identified for analysis and modeling. These include low pressure eduction, high flow rate purge, and pulse purge. For each method, an approach to predict the Hz mass transfer rate, as a function of system pressure, temperature, and purge gas mass flow rate, is developed based on the applicable physics. The models show that low pressure eduction is two orders of magnitude more effective than the high velocity purge, which in turn is two orders of magnitude more effective than the pure diffusion component of pulse purging of deadheads. Eduction subjects the system to low pressure conditions that promote the extraction of Hz vapors. At 120 F, Hz is saturated at approximately 1 psia. At lower pressures and 120 F, Hz will boil, which is an extremely efficient means to remove liquid Hz. The Hz boiling rate is predicted by equating the rate at which energy is added to the saturated liquid Hz through heaters at the tube outer wall with the energy removed from the liquid through evaporation. Boil-off fluxes were predicted by iterating through the range of local pressures with limits set by the minimum allowed pressure of 0.2 psia and maximum allowed wall temperature of 120 F established by the heaters, which gives a saturation pressure of approximately 1.0 psia. Figure 1 shows the resulting boil-off fluxes as a function of local eduction pressure. As depicted in figure 1, the flux is a strong inverse function of eduction pressure, and that minimizing the eduction pressure maximizes the boil-off flux. Also, higher outer wall temperatures lead to higher boil-off fluxes and allow for boil-off over a greater range of eduction pressures.

Pasch, James J.↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
NASA NTRS1999

Growing wheat in Biosphere 2 under elevated CO2: observations and modeling

Spring wheat (Triticum aestivum L., cv. Yecora Rojo) was grown in the intensive agricultural biome (IAB) of Biosphere 2 during the l995-l996 winter/spring season. Environmental conditions were characterized by a day/night temperature regime of 27/17 degrees C, relative humidity (RH) levels around 45%, mean atmospheric CO2 concentration of 450 ppmv, and natural light conditions with mean intensities about half of outside levels. Weekly samples of above-ground plant matter were collected throughout the growing season and phenological events recorded. A computer model, CERES-Wheat, previously tested under both field and controlled conditions, was used to simulate the observed crop growth and to help in data analysis. We found that CERES-Wheat simulated the data collected at Biosphere 2 to within 10% of observed, thus suggesting that wheat growth inside the IAB was comparable to that documented in other environments. The model predicts phenological stages and final dry matter (DM) production within l0% of the observed data. Measured DM production rates, normalized for light absorbed by the crop. suggested photosynthetic efficiencies intermediate between those observed under optimal field conditions and those recorded in NASA-Controlled Ecological Life-Support Systems (CELSS). We suggest that such a difference can be explained primarily in terms of low light levels inside the IAB, with additional effects due to elevated CO2 concentrations and diffuse light fractions.

Non-NASA Center↗
NASA NTRS2018

Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Comparisons with experimental data showed a sensitivity of the current prediction to the modeling of the sidewalls. This was found to be common among various computational studies in the literature where periodic boundary conditions were used in the spanwise direction, as was the case in the present work. Thus, although the experiment was quasi-two-dimensional, the present simulation was determined to be two-dimensional. Quantities present in the exact equation of the Reynolds-stress transport, i.e., production, molecular diffusion, turbulent transport, pressure diffusion, pressure strain, dissipation, and turbulent mass flux were calculated. Reynolds-stress budgets were compared with past large-eddy simulation and direct numerical simulation datasets in the undisturbed portion of the turbulent boundary layer to validate the current approach. The budgets in SBLI showed the growth in the production term for the primary normal stress and energy transfer mechanism was led by the pressure strain term in the secondary normal stresses. The pressure diffusion term, commonly assumed as negligible by turbulence model developers, was shown to be small but non-zero in the normal stress budgets, however it played a key role in the primary shear stress budget.

shock wave/boundary-layer interactio↗