Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stress equilibrium”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Directed aging, memory, and nature’s greed

Disordered materials are often out of equilibrium and evolve very slowly in a rugged and tortuous energy landscape. This slow evolution, referred to as aging, is deemed undesirable as it often leads to material degradation. However, we show that aging also encodes a memory of the stresses imposed during preparation. Because of inhomogeneous local stresses, the material itself decides how to evolve by modifying stressed regions differently from those under less stress. Because material evolution occurs in response to stresses, aging can be “directed” to produce sought-after responses and unusual functionalities that do not inherently exist. Aging obeys a natural “greedy algorithm” as, at each instant, the material simply follows the path of most rapid and accessible relaxation. Our experiments and simulations illustrate directed aging in examples in which the material’s elasticity transforms as desired because of an imposed deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Turbulence: The chief outstanding difficulty of our subject

A review of interesting current topics in turbulence research is decorated with examples of popular fallacies about the behavior of turbulence. Topics include the status of the Law of the Wall, especially in compressible flow; analogies between the effects of Reynolds numbers, pressure gradient, unsteadiness and roughness change; the status of Kolmogorov's universal equilibrium theory and local isotropy of the small eddies; turbulence modelling, with reference to universality, pressure-strain modelling and the dissipation equation; and chaos. Fallacies include the mixing-length concept; the effect of pressure gradient on Reynolds shear stress; the separability of time and space derivatives; models of the dissipation equation; and chaos.

Bradshaw, Peter↗
DOE OSTI2023

Deformation mechanisms in an additively manufactured dual-phase eutectic high-entropy alloy

Nanostructured metals and alloys often exhibit high strengths but at the expense of reduced ductility. Through harnessing the far-from-equilibrium processing conditions of laser powder-bed fusion (L-PBF) additive manufacturing, we develop a dual-phase nanolamellar structure comprised of FCC/L1 2 and BCC/B2 phases in a Ni 40 Co 20 Fe 10 Cr 10 Al 18 W 2 eutectic high-entropy alloy (EHEA), which exhibits a combination of ultrahigh yield strength (>1.4 GPa) and large tensile ductility (~17%). The deformation mechanisms of the additively manufactured EHEA are studied via in-situ synchrotron X-ray diffraction and high-resolution transmission electron microscopy. The high yield strength mainly results from effective blockage of dislocation motion by the high density of lamellar interfaces. The refined nanolamellar structures and low stacking fault energy (SFE) promote stacking fault (SF)-mediated deformation in FCC/L1 2 nanolamellae. The accumulation of abundant dislocations and SFs at lamellar interfaces can effectively elevate local stresses to promote dislocation multiplication and martensitic transformation in BCC/B2 nanolamellae. The cooperative deformation of the dual phases, assisted by the semi-coherent lamellar interfaces, gives rise to the large ductility of the as-printed EHEA. In addition, here we also demonstrate that post-printing heat treatment allows us to tune the non-equilibrium microstructures and deformation mechanisms. After annealing, the significantly reduced SFE and thicknesses of the FCC nanolamellae facilitate the formation of massive SFs. The dissolution of nano-precipitates in the BCC/B2 nanolamellae reduces spatial confinement and further promotes martensitic transformation to enhance work hardening. Our work provides fundamental insights into the rich variety of deformation mechanisms underlying the exceptional mechanical properties of the additively manufactured dual-phase nanolamellar EHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗
NASA NTRS1975

Hot corrosion evaluation of aluminide coated superalloys in support of an ASTM Round Robin program

Commercial aluminized coatings on substrates were hot corroded at 900 C in a 0.3 Mach burner rig with 5 ppm synthetic sea salt and at two cycling frequencies. Extensive post-exposure examinations were conducted on the corroded specimens such as metallography, X-ray diffraction, scanning electron microscopy, microprobe raster scans, and spectrographic analyses. Thermodynamic calculations were made of the equilibrium burner flame composition and the calculations were compared to the experimental findings. It was found that localized spalling of the coatings preceded coating failure. It is suggested that the spalling of the coatings is due to the formation of localized stresses caused by the depletion of chromium and aluminum in the coating or the enrichment of the coating with sulfur. For the materials and test conditions investigated, it was found that coating life is dependent only upon the initial coating thickness and not on the type of aluminized coating, the substrate, or the cycle frequency.

Santoro, G.↗
DOE OSTI2022

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Characteristics of Ultralow Aspect Ratio Toroidal Plasmas and Non-Solenoidal Startup and Edge Stability Studies at Near-Unity Aspect Ratio (Final Report)

This grant supported a long-term research and development program to explore the physics and technology of spherical tokamak (ST) plasma behavior at near-unity aspect ratio ( A ). The PEGASUS Toroidal Experiment is a university-scale ( R ) ~ 0.40 m, α ~ .35 m, I p ≤ 0.3 MA, B T ≤ 0.15 T) experiment that accesses A → 1, to provide large values of the normalized plasma current, I N = I p /α B T , or equivalently the toroidal field utilization factor I p / I TF . The research activities focused on: 1) developing the technique of Local Helicity Injection (LHI) as a viable means of initiating and growing an ST plasma without the use of a central solenoid; 2) exploring the equilibrium and stability properties of spherical tokamak plasmas at extremely low aspect ratio; and 3) obtaining new insights in the properties of the enhanced confinement H-mode regime at A ~ 1. The initial (Phase I) facility employed a unique high-stress central solenoid Ohmic for current drive and heating with simple magnet current waveforms. A subsequent major upgrade (Phase II) offered more flexible control of the plasma evolution along with a powerful LHI system for plasma initiation, heating and sustainment. The LHI technique consists of injection of high-density current streams into vacuum toroidal and vertical magnetic fields. These streams become unstable and support magnetic reconnection to allow the system to relax into a tokamak-like magnetic configuration. This program focused on: 1) development of electron current sources that provide sufficiently high injection currents in the challenging edge plasma region; 2) validation of theoretical limits to achievable I p arising from magnetic helicity and energy conservation plus absolute limits arising from relaxation to a minimal energy states; 3) tests of a power balance model to describe the evolution of the plasma current I p (t); and 4) study of magnetic fluctuations and their role in the I p evolution. This culminated in the production of toroidal currents of I p ≥ 0.2 MA with only ~8 kA of injected current. Operation with Ohmic induction alone in Phase I created plasmas with soft limits on the achievable I p and field utilization ( I p / I TF ~ 1) due to the onset of large low-order (2/1 or 3/2) internal resonant tearing modes. With the expanded capabilities in Phase II operation, methods employed to mitigate those modes and achieve I p / I TF > 1 were: 1) I TF ramp-downs; 2) non-solenoidal startup via LHI; and 3) extremely low I TF Ohmic operations assisted by LHI current injector pre-ionization. Extensive use of LHI-driven startup and drive for non-solenoidal operation produced favorable current profiles to stably operate tokamak plasmas with unprecedentedly high I N = I p / aB T ~ 15 and I p / I TF > 2. With strong reconnection-driven anomalous ion heating, this in turn provided access to tokamak plasmas with B t up to 100%, a world record 2 times higher than that previously achieved in advanced or spherical tokamaks. This high B t regime had a large fraction (~50%) of the plasma volume in a region of an absolute minimum- B configuration. The MHD stability limit in this regime was the ideal external kink mode without a stabilizing wall. These results provided the first demonstration of access to tokamak plasmas with very high B t , high i N , high elongation κ, and low internal inductance ℓ i , and as such served as a proof-ofprinciple for access to this long-sought ST regime. The naturally low B T needed for stability at A ≤ 1.2 reduces the L-H mode power threshold, P LH , so that H-mode plasmas were achieved with Ohmic heating alone. The power threshold for the transition to the H-mode regime was found to be 10× the accepted values from the ITER international scaling, confirming initial trends indicated by experiments on MAST and NSTX, and suggesting that the physics behind this transition is not yet fully understood. The density dependence of the transition and the observed insensitivity to whether the plasma is bounded by a limiter or magnetic divertor are consistent with the recent FM 3 model. Multiple unstable toroidal magnetic modes were measured during ELM (Edge Localized Mode) crashes, with generally lower toroidal mode numbers than seen at higher A , which was attributed to the larger peeling mode instability drive at low A . First-ever edge current profile measurements through an ELM crash showed current-carrying filaments were generated, as seen in ELM simulations and in electromagnetic blob transport theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Analysis and Design of Variable Stiffness Composite Cylinders

An investigation of the possible performance improvements of thin circular cylindrical shells through the use of the variable stiffness concept is presented. The variable stiffness concept implies that the stiffness parameters change spatially throughout the structure. This situation is achieved mainly through the use of curvilinear fibers within a fiber-reinforced composite laminate, though the possibility of thickness variations and discrete stiffening elements is also allowed. These three mechanisms are incorporated into the constitutive laws for thin shells through the use of Classical Lamination Theory. The existence of stiffness variation within the structure warrants a formulation of the static equilibrium equations from the most basic principles. The governing equations include sufficient detail to correctly model several types of nonlinearity, including the formation of a nonlinear shell boundary layer as well as the Brazier effect due to nonlinear bending of long cylinders. Stress analysis and initial buckling estimates are formulated for a general variable stiffness cylinder. Results and comparisons for several simplifications of these highly complex governing equations are presented so that the ensuing numerical solutions are considered reliable and efficient enough for in-depth optimization studies. Four distinct cases of loading and stiffness variation are chosen to investigate possible areas of improvement that the variable stiffness concept may offer over traditional constant stiffness and/or stiffened structures. The initial investigation deals with the simplest solution for cylindrical shells in which all quantities are constant around the circumference of the cylinder. This axisymmetric case includes a stiffness variation exclusively in the axial direction, and the only pertinent loading scenarios include constant loads of axial compression, pressure, and torsion. The results for these cases indicate that little improvement over traditional laminates exists through the use of curvilinear fibers, mainly due to the presence of a weak link area within the stiffness variation that limits the ultimate load that the structure can withstand. Rigorous optimization studies reveal that even though slight increases in the critical loads can be produced for designs with an arbitrary variation of the fiber orientation angle, the improvements are not significant when compared to traditional design techniques that utilize ring stiffeners and frames. The second problem that is studied involves arbitrary loading of a cylinder with a stiffness variation that changes only in the circumferential direction. The end effects of the cylinder are ignored, so that the problem takes the form of an analysis of a cross-section for a short cylinder segment. Various load cases including axial compression, pressure, torsion, bending, and transverse shear forces are investigated. It is found that the most significant improvements in load-carrying capability exist for cases which involve loads that also vary around the circumference of the shell, namely bending and shear forces. The stiffness variation of the optimal designs contribute to the increased performance in two ways: lowering the stresses in the critical areas through redistribution of the stresses; and providing a relatively stiff region that alters the buckling behavior of the structure. These results lead to an in-depth optimization study involving weight optimization of a fuselage structure subjected to typical design constraints. Comparisons of the curvilinear fiber format to traditional stiffened structures constructed of isotropic and composite materials are included. It is found that standard variable stiffness designs are quite comparable in terms of weight and load-carrying capability yet offer the added advantage of tailorability of distinct regions of the structure that experience drastically different loading conditions. The last two problems presented in this work involve the nonlinear phenomenon of long tubes under bending. Though this scenario is not as applicable to fuselage structures as the previous problems, the mechanisms that produce the nonlinear effect are ideally suited to be controlled by the variable stiffness concept. This is due to the fact that the dominating influence for long cylinders under bending is the ovalization of the cross-section, which is governed mainly by the stiffness parameters of the cylindrical shell. Possible improvement of the critical buckling moments for these structures is investigated using either a circumferential or axial stiffness variation. For the circumferential case involving infinite length cylinders, it is found that slight improvements can be observed by designing structures that resist the cross-sectional deformation yet do not detract from the buckling resistance at the critical location. The results also indicate that buckling behavior is extremely dependent on cylinder length. This effect is most easily seen in the solution of finite length cylinders under bending that contain an axial stiffness variation. For these structures, the only mechanism that exhibits improved response are those that effectively shorten the length of the cylinder, thus reducing the cross-sectional deformation due to the forced restraint at the ends. It was found that the use of curvilinear fibers was not able to achieve this effect in sufficient degree to resist the deformation, but that ring stiffeners produced the desired response admirably. Thus, it is shown that the variable stiffness concept is most effective at improving the bending response of long cylinders through the use of a circumferential stiffness variation.

Tatting, Brian F.↗
NASA NTRS1985

Environmental Degradation of Materials: Surface Chemistry Related to Stress Corrosion Cracking

Parallel experiments have been performed in order to develop a comprehensive model for stress cracking (SCC) in structural materials. The central objective is to determine the relationship between the activity and selectivity of the microstructure of structural materials to their dissolution kinetics and experimentally measured SCC kinetics. Zinc was chosen as a prototype metal system. The SCC behavior of two oriented single-crystal disks of zinc in a chromic oxide/sodium sulfate solution (Palmerton solution) were determined. It was found that: (1) the dissolution rate is strongly (hkil)-dependent and proportional to the exposure time in the aggressive environment; and (2) a specific slip system is selectively active to dissolution under applied stress and this slip line controls crack initiation and propagation. As a precursor to potential microgrvity experiments, electrophoretic mobility measurements of zinc particles were obtained in solutions of sodium sulfate (0.0033 M) with concentrations of dissolved oxygen from 2 to 8 ppm. The equilibrium distribution of exposed oriented planes as well as their correlation will determine the particle mobility.

Schwarz, J. A.↗
DOE OSTI2020

Beyond Equilibrium Temperature: How the Atmosphere/Interior Connection Affects the Onset of Methane, Ammonia, and Clouds in Warm Transiting Giant Planets

The atmospheric pressure–temperature profiles for transiting giant planets cross a range of chemical transitions. Here we show that the particular shapes of these irradiated profiles for warm giant planets below ~1300 K lead to striking differences in the behavior of nonequilibrium chemistry compared to brown dwarfs of similar temperatures. Our particular focus is H 2 O, CO, CH 4 , CO 2 , and NH 3 in Jupiter- and Neptune-class planets. We show that the cooling history of a planet, which depends most significantly on planetary mass and age, can have a dominant effect on abundances in the visible atmosphere, often swamping trends one might expect based on T eq alone. The onset of detectable CH 4 in spectra can be delayed to lower T eq for some planets compared to equilibrium, or pushed to higher T eq . The detectability of NH 3 is typically enhanced compared to equilibrium expectations, which is opposite to the brown dwarf case. We find that both CH 4 and NH 3 can become detectable at around the same T eq (at T eq values that vary with mass and metallicity), whereas these “onset” temperatures are widely spaced for brown dwarfs. We suggest observational strategies to search for atmospheric trends and stress that nonequilibrium chemistry and clouds can serve as probes of atmospheric physics. As examples of atmospheric complexity, we assess three Neptune-class planets, GJ 436b, GJ 3470b, and WASP-107, all around T eq = 700 K. Tidal heating due to eccentricity damping in all three planets heats the deep atmosphere by thousands of degrees and may explain the absence of CH 4 in these cool atmospheres. Atmospheric abundances must be interpreted in the context of physical characteristics of the planet.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Experimental investigation of turbulent flow through a circular-to-rectangular transition duct

An incompressible, turbulent, swirl-free flow through a circular-to-rectangular transition duct was studied experimentally. The cross-sectional geometry all along the duct was defined using the equation of a superellipse. The three mean velocity components and the six Reynolds stress components were measured at two axial stations downstream from the transition. It is shown that a secondary flow vortex pair which develops along the duct sidewalls significantly distorts the mean and turbulence fields. At the duct exit, the flow is not in local equilibrium, but recovers to local equilibrium conditions in the rectangular extension duct. Analysis demonstrates that conventional wall functions are not applicable at all streamwise locations in the duct.

Davis, D. O.↗
NASA NTRS1991

Chemical vapor deposition modeling: An assessment of current status

The shortcomings of earlier approaches that assumed thermochemical equilibrium and used chemical vapor deposition (CVD) phase diagrams are pointed out. Significant advancements in predictive capabilities due to recent computational developments, especially those for deposition rates controlled by gas phase mass transport, are demonstrated. The importance of using the proper boundary conditions is stressed, and the availability and reliability of gas phase and surface chemical kinetic information are emphasized as the most limiting factors. Future directions for CVD are proposed on the basis of current needs for efficient and effective progress in CVD process design and optimization.

Gokoglu, Suleyman A.↗
NASA NTRS1945

Stresses in and general instability of monocoque cylinders with cutouts II : calculation of the stresses in a cylinder with a symmetric cutout

A numerical procedure is presented for the calculation of the stresses in a monocoque cylinder with a cutout. In the procedure the structure is broken up into a great many units; the forces in these units corresponding to specified distortions of the units are calculated; a set of linear equations is established expressing the equilibrium conditions of the units in the distorted state; and the simultaneous linear equations are solved. A fully worked out numerical example, corresponding to the application of a pure bending moment, gave results in good agreement with experiments carried out earlier at the Polytechnic Institute of Brooklyn.

Hoff, N J↗
DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Amidination of ligands for chemical and field-effect passivation stabilizes perovskite solar cells

Surface passivation has driven the rapid increase in the power conversion efficiency (PCE) of perovskite solar cells (PSCs). However, state-of-the-art surface passivation techniques rely on ammonium ligands that suffer deprotonation under light and thermal stress. We developed a library of amidinium ligands, of interest for their resonance effect–enhanced N–H bonds that may resist deprotonation, to increase the thermal stability of passivation layers on perovskite surfaces. This strategy resulted in a >10-fold reduction in the ligand deprotonation equilibrium constant and a twofold increase in the maintenance of photoluminescence quantum yield after aging at 85°C under illumination in air. Here, implementing this approach, we achieved a certified quasi–steady-state PCE of 26.3% for inverted PSCs; and we report retention of ≥90% PCE after 1100 hours of continuous 1-sun maximum power point operation at 85°C.

Yang, Yi↗
DOE OSTI2019

Novel ductile wellbore cementitious composite for geologic CO 2 storage

CO 2 leakage through damaged wellbore cement sheaths is a major risk of geologic CO 2 storage (GCS), as conventional wellbore cement is brittle and can be damaged due to acid attack and downhole stresses during CO 2 injection. Here we examine a novel fiber-reinforced engineered cementitious composite (ECC) proposed as a substitute to conventional wellbore cement due to its superior ductility and intrinsic crack width control. ECC and conventional wellbore cement coupons were exposed to water in equilibrium with CO 2 at 50 °C and 10 MPa. The samples were retrieved after several days and their mechanical performance was evaluated using a four-point bending test, microhardness, and compressive strength analyses. Optical microscopy and mercury intrusion porosimetry were used to characterize the progression of the carbonation front and pore structures of the specimens. Control experiments were conducted under the same temperature and pressure conditions but with a N 2 headspace to isolate the impact of CO 2 . It was found that carbonation increased the ultimate flexural strength of ECC but decreased its ductility. However, the ductility of carbonated ECC remained higher than that of conventional wellbore cements that exhibited brittle failure under all test conditions. Additionally, ECC exhibited minimal material loss and continued resistance to deformation in comparison to conventional wellbore cements. This suggests that while the exposure of ECC to CO 2 will alter its mechanical properties, altered ECC will continue to exhibit mechanical performance superior to conventional wellbore cement, and therefore shows promise as a highly durable wellbore cementing material for GCS applications.

58 GEOSCIENCES↗
NASA NTRS1999

Turbulence Model Comparisons for a High-Speed Aircraft

Four turbulence models are described and evaluated for transonic flows over the High-Speed Research/industry baseline configuration known as Reference H by using the thin-layer, upwind, Navier-Stokes solver known as CFL3D. The turbulence models studied are the equilibrium model of Baldwin-Lomax (B-L) with the Degani-Schiff (D-S) modifications, the one-equation Baldwin-Barth (B-B) model, the one-equation Spalart-Allmaras (S-A) model, and Menter's two-equation Shear Stress Transport (SST) model. The flow conditions, which correspond to tests performed in the National Transonic Facility (NTF) at Langley Research Center, are a Mach number of 0.90 and a Reynolds number of 30 x 10 (exp. 6) based on mean aerodynamic chord for angles of attack of 1 deg., 5 deg., and 10 deg. The effects of grid topology and the representation of the actual wind tunnel model geometry are also investigated. Computed forces and surface pressures compare reasonably well with the experimental data for all four turbulence models.

Rivers, Melissa B.↗
NASA NTRSDate not supplied

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗