Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “friction experiments”

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.

823 records · Page 46

Particle Segregation in Collisional Shearing Flows

The size segregation of flowing or shaken grains is a commonly observed phenomenon in industrial processes and in nature. In systems that do not involve much agitation of the grains, several mechanisms that involve gravity have been identified as leading to such segregation. In highly agitated flows, there is a mechanism independent of gravity that is available to drive separation of different grains. This is associated with spatial gradients in the energy of their velocity fluctuations. Because collisions between grains inevitably dissipate energy, collisional granular shear flows are usually of limited extent in the direction transverse to the flow. One consequence of this is that shear flows are strongly influenced by their boundaries. Because grains, on average, slip relative to boundaries, a bumpy or frictional boundary can convert slip energy into fluctuation energy. However, because each collision between a grain and the boundary dissipates fluctuation energy, there is a competition between production and dissipation. In principle, it is possible to design the geometry of the boundary - for example, the size and spacing of the bumps - so that the boundary either produces or dissipates fluctuation energy. This permits the control of the component of the spatial gradient of the fluctuation energy that is normal to the boundary. The gradients in fluctuation energy established by such boundaries may be exploited to drive the separation by size or other properties in a binary mixture of spherical grains. Microgravity makes the visual observations possible by permitting us to employ moderate rates of shear. On earth, the effects of gravity can be minimized by shearing so rapidly that the particle pressure overwhelms gravity. However, in this event, separation takes place too rapidly for visual observation, buoyancy and/or condensation associated with the centripetal acceleration must be accounted for, and the particles can be severely damaged. Because, in the absence of gravity, the only available time scale is proportional to the speed of the moving boundary, this speed can be made arbitrarily slow to permit observations and to avoid particle damage, without altering the phenomenon under study. The primary goal of this research is to carry out a physical experiment in which particle segregation is induced and maintained in a collisional flow of a binary mixture of two different types of spheres.

Jenkins, J. T.↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

Simulated Service and Stress Corrosion Cracking Testing for Friction Stir Welded Spun Form Domes

Damage tolerance testing development was required to help qualify a new spin forming dome fabrication process for the Ares 1 program at Marshall Space Flight Center (MSFC). One challenge of the testing was due to the compound curvature of the dome. The testing was developed on a sub-scale dome with a diameter of approximately 40 inches. The simulated service testing performed was based on the EQTP1102 Rev L 2195 Aluminum Lot Acceptance Simulated Service Test and Analysis Procedure generated by Lockheed Martin for the Space Shuttle External Fuel Tank. This testing is performed on a specimen with an induced flaw of elliptical shape generated by Electrical Discharge Machining (EDM) and subsequent fatigue cycling for crack propagation to a predetermined length and depth. The specimen is then loaded in tension at a constant rate of displacement at room temperature until fracture occurs while recording load and strain. An identical specimen with a similar flaw is then proof tested at room temperature to imminent failure based on the critical offset strain achieved by the previous fracture test. If the specimen survives the proof, it is then subjected to cryogenic cycling with loads that are a percentage of the proof load performed at room temperature. If all cryogenic cycles are successful, the specimen is loaded in tension to failure at the end of the test. This standard was generated for flat plate, so a method of translating this to a specimen of compound curvature was required. This was accomplished by fabricating a fixture that maintained the curvature of the specimen rigidly with the exception of approximately one-half inch in the center of the specimen containing the induced flaw. This in conjunction with placing the center of the specimen in the center of the load train allowed for successful testing with a minimal amount of bending introduced into the system. Stress corrosion cracking (SCC) tests were performed using the typical double beam assembly and with 4-point loaded specimens under alternate immersion conditions in a 3.5% NaCl environment for 90 days. In addition, experiments were conducted to determine the threshold stress intensity factor for SCC (K1SCC) of Al-Li 2195 which to our knowledge has not been determined previously. The successful simulated service and stress corrosion testing helped to provide confidence to continue to Ares 1 scale dome fabrication.

Stewart, Thomas J.↗

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗

Design and Analysis of High Heat Flux Plasma-Facing Components for NSTX Upgrade

The National Spherical Tokamak Experiment (NSTX) upgrade (NSTX-U) requirements lead to enhanced heat loads on plasma-facing components (PFCs) especially in the divertor regions, where normal heat flux density can reach up to 8 MW/m 2 . For these high heat flux (HHF) regions, the PFCs were redesigned, to use castellations which relieved the thermal stresses associated with high incident heat flux. Another design feature of HHF PFCs is the absence of front face mounting holes which create localized areas of high temperature and thermal stress concentrations. Optimized shaping of the front face of the HHF PFCs eliminates regions where the front face is perpendicular to the heat flux direction during normal helicity operation and spreads the heat load. A unique mechanism was designed to connect PFC tiles to the NSTX-U center stack casing using locking features accessible from the low heat flux regions. Isotropic graphite was selected as the HHF PFC tile material. Initial tile design parameters were accessed using analytical expressions for pulsed heat flux loading. A working prototype of the locking mechanism was created during the initial stages of the design to prove the concept performance. This article presents an overview of the divertor HHF PFC design and the results of the thermal and structural analyses performed using ANSYS software. The results of the analyses cover normal operating conditions and disruptions which impose electromagnetic (EM) loads from eddy and halo currents. The 3-D, transient, nonlinear analyses took into account the temperature-dependent properties of the materials, friction interfaces between the parts, and variable electric properties of the parts and interfaces. The results confirm that the HHF PFC tiles remain within the allowable limits for the loads defined by the NSTX-U Recovery Project. In addition, tolerance stack up analyses were performed to ensure tile performance in the worst possible assembly configuration. Finally, the design process was completed successfully, and the NSTX HHF PFC tiles are currently in the production phase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spaceflight Medical Evacuation Risk Assessment Principles - A Qualitative Investigation

BACKGROUND Future human space exploration beyond Low Earth Orbit (LEO) will require innovative solutions in many areas, primarily those that provide direct medical support to crew on long-duration missions and optimize their health and performance. The associated challenges therein will be numerous, including but not necessarily limited to extended one-way or "asynchronous" communication delays, minimal to non-existent resupply, and a prolonged transit time to "definitive care" ranging from 3-days to 9-months. At the same time, long-duration exploration spacecraft and crew will face restrictions on mass, power, volume, and data far more significant than that seen in current LEO settings. Given these limitations, medical risk assessment is of primary importance, especially evaluating the implications of a medical evacuation of an ill or injured crewmember. Such evacuations are complicated, potentially dangerous, and well may be impossible in certain phases of the mission. Regardless, such issues must be weighed against the risks of the injured crew remaining aboard a spacecraft with limited medical resources. OBJECTIVE This qualitative study drew from the experiences of subject matter experts (SMEs) in spaceflight and appropriate analog environments (i.e., military, disaster, and extreme environment fields) to identify unique principles common amongst medical evacuation considerations helpful in informing future risk assessment tools. Appropriate analog environments included austere operational settings where multiple factors (weather, logistics/limited resupply, denied/extreme environments, and patient condition) resulted in a limited ability to provide definitive local medical care. The fundamental principles in question revolved around scenarios where evacuation became a complicating yet necessary consideration and where life-threatening medical concerns had to be weighed against critical mission objective(s). The primary authors collected semi-structured data gathered through in-depth interviews with 16 subject matter experts (SMEs). Interview questions investigated how these SMEs consider and weigh the attendant risks present in medical evacuation scenarios. Among the critical questions posed were those that sought to understand how the SME balanced the challenges and requirements of medical evacuation (or keeping an injured patient "on-site" aka: "prolonged field care") against the evacuation operation's risks on impacting overarching mission success. The team analyzed interview transcripts for common themes and principles using the qualitative methods of thematic analysis based on consensus, co-occurrence, and comparison. As a result, nine primary risk consideration themes and nine contributing factor themes emerged, all of which will ideally inform future medical evacuation decision-making tools and operational decision-making for exploration class missions. Specific aims for this study included: 1. Identification of common principles used to assess risks and benefits of medical evacuations in extreme environments 2. Identification of common points of friction or complication and challenges in extreme environment evacuations

A T Almand↗

Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer

Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.

42 ENGINEERING↗

Multifunctional, High-Temperature Nanocomposites

In experiments conducted as part of a continuing effort to incorporate multifunctionality into advanced composite materials, blends of multi-walled carbon nanotubes and a resin denoted gPETI-330 h (wherein gPETI h is an abbreviation for gphenylethynyl-terminated imide h) were prepared, characterized, and fabricated into moldings. PETI-330 was selected as the matrix resin in these experiments because of its low melt viscosity (<10 poise at a temperature of 280 C), excellent melt stability (lifetime >2 hours at 280 C), and high temperature performance (>1,000 hours at 288 C). The multi-walled carbon nanotubes (MWCNTs), obtained from the University of Kentucky, were selected because of their electrical and thermal conductivity and their small diameters. The purpose of these experiments was to determine the combination of thermal, electrical, and mechanical properties achievable while still maintaining melt processability. The PETI-330/MWCNT mixtures were prepared at concentrations ranging from 3 to 25 weight-percent of MWCNTs by dry mixing of the constituents in a ball mill using zirconia beads. The resulting powders were characterized for degree of mixing and thermal and rheological properties. The neat resin was found to have melt viscosity between 5 and 10 poise. At 280 C and a fixed strain rate, the viscosity was found to increase with time. At this temperature, the phenylethynyl groups do not readily react and so no significant curing of the resin occurred. For MWCNT-filled samples, melt viscosity was reasonably steady at 280 C and was greater in samples containing greater proportions of MWCNTs. The melt viscosity for 20 weightpercent of MWCNTs was found to be .28,000 poise, which is lower than the initial estimated allowable maximum value of 60,000 poise for injection molding. Hence, MWCNT loadings of as much as 20 percent were deemed to be suitable compositions for scale-up. High-resolution scanning electron microscopy (HRSEM) showed the MWCNTs to be well dispersed in the polymer matrices, while high-resolution transmission electron microscopy shows splits in the walls of the MWCNTs but no catastrophic breakage of tubes. To further assess processing characteristics prior to scale-up, samples containing 10, 15, and 20 weight-percent of MWCNTs were processed through a laboratory melting extruder. HRSEM of the extruded fibers shows significant alignment of MWCNTs in the flow direction (see figure). For the samples containing 20 weight-percent of MWCNTs, difficulties were encountered during feeding, and the temperature of a rotor in the extruder rose to 245 C because of buildup of frictional heat; this indicates that materials of this type having MWCNT concentrations .20 weight- percent may not be melt-processable. On the basis of the results from the foregoing characterizations, samples containing 10, 15, and 20 weight-percent of MWCNTs were scaled up to masses of .300 g and used to make specimens having dimensions of 10.2 by 15.2 by 0.32 cm. These specimens were molded by (1) injecting the mixtures, at temperatures between 260 and 280 C, into a tool made of the low-thermal-expansion alloy InvarR and then (2) curing for 1 hour at 371 C. The tool was designed to impart shear during the injection process in an attempt to achieve some alignment of the MWCNTs in the flow direction.

Connell, John W.↗

RELAP5-3D Modeling of High Temperature Test Facility (HTTF) Test PG-26

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors with prismatic fuel and reflector blocks. Several tests have been completed, and more are still planned to at the HTTF, including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This report analyses test PG-26, a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC. PG-26 has been performed at the HTTF between May 30 and June 30, 2019. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heaters. The DCC transient was initiated during the 50 th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test. Version 4.4.2ie of the RELAP5-3D computer code has been used to model the HTTF PG-26 test, and results have been compared to available high-quality measured data. The model used in this study is the quality-controlled HTTF RELAP5-3D model (HTTF base 2018-04-19 QA), originally developed by P. Bayless. The report includes RELAP5-3D results of the “base calculations” as well as some sensitivities to important uncertain model inputs, such as primary helium mass flow rate, core ceramic thermal properties, as well as heat evacuation and loop friction models. Using the base RELAP5-3D model predicts a countercurrent helium flow in the hot duct observed at the beginning of the DCC, but instead of going into a molecular diffusion mode, the model predicts the onset of natural convection. Increasing friction in the core and hot duct prevents the natural convection from happening in some of the simulations. Although some temperatures are well predicted (and even overpredicted), the general tendency is to underpredict the ceramic and helium temperatures and heat removal rates during the DCC, resulting in many of the assessment findings being in minimal or insufficient agreement with the data. It is worth noting that the described discrepancies between measured data and RELAP5-3D predictions are not RELAP5-3D code limitations. More so, they reflect limitations in boundary condition and thermal property knowledge. While the RELAP5-3D calculations of the test provide some insights into what happens during the transient, and point to missing or potentially uncertain data to which the experimenters can direct their attention, the principal conclusion is that the PG-26 test data are insufficient for a system code assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High Temperature Mechanical Characterization of Ceramic Matrix Composites

A high temperature mechanical characterization laboratory has been assembled at NASA Lewis Research Center. One contribution of this work is to test ceramic matrix composite specimens in tension in environmental extremes. Two high temperature tensile testing systems were assembled. The systems were assembled based on the performance and experience of other laboratories and meeting projected service conditions for the materials in question. The systems use frames with an electric actuator and a center screw. A PC based data acquisition and analysis system is used to collect and analyze the data. Mechanical extensometers are used to measure specimen strain. Thermocouples, placed near the specimen, are used to measure the specimen gage section temperature. The system for testing in air has a resistance element furnace with molybdenum disilicide elements and pneumatic grips with water cooling attached to hydraulic alignment devices. The system for testing in an inert gas has a graphite resistance element furnace in a chamber with rigidly mounted, water cooled, hydraulically actuated grips. Unidirectional SiC fiber reinforced reaction bonded Si3N4 and triaxially woven, two dimensional, SiC fiber reinforced enhanced SiC composites were tested in unidirectional tension. Theories for predicting the Young's modulus, modulus near the ultimate strength, first matrix cracking stress, and ultimate strength were applied and evaluated for suitability in predicting the mechanical behavior of SiC/RBSN and enhanced SiC/SiC composites. The SiC/RBSN composite exhibited pseudo tough behavior (increased area under the stress/strain curve) from 22 C to 1500 C. The rule of mixtures provides a good estimate of the Young's modulus of the SiC/RBSN composite using the constituent properties from room temperature to 1440 C for short term static tensile tests in air or nitrogen. The rule of mixtures significantly overestimates the secondary modulus near the ultimate strength. The ACK theory provides the best approximation of the first matrix cracking stress when residual stresses are ignored. The theory of Cao and Thouless, based on Weibull statistics, gave the best prediction for the composite ultimate strength. The enhanced SiC/SiC composite exhibited nonlinear stress/strain behavior from 24 C to 1370 C in air with increased ultimate strain when compared to monolithic SiC. The theory of Yang and Chou with the assumption of a frictional fiber/matrix interface provided the best estimate of the Young's modulus. The theory of Cao and Thouless gave the best estimate for the ultimate strength.

Gyekenyesi, John Z.↗

Coupled Investigation of Fracture Permeability Impact on Reservoir Stress and Seismic Slip Behavior (Final Technical Report)

Enhanced Geothermal Systems (EGS) produce clean energy by circulating fluid through hot rock deep underground and bringing that heat to the surface to generate electricity. For this process to work reliably, fluids must be able to move efficiently through networks of natural or engineered fractures in the rock. Enhancing and maintaining subsurface permeability over time is essential for sustainable energy production. However, fluid injection changes the underground temperature, pressure, rock stress, and chemistry, which can alter permeability and sometimes trigger earthquakes. Predicting these interconnected processes remains a key challenge. To address this, we combined high-temperature laboratory experiments with high-fidelity simulations to better understand how fractures in geothermal reservoirs evolve over time. Our experiments measured how fractures respond to stress, slip, slip rate, and chemical reactions under geothermal conditions. These data were integrated into coupled thermal-hydrological-mechanical-chemical and earthquake (THMC+E) models tailored to the Utah FORGE site. The validated modeling framework improves predictions of reservoir performance and seismic response and helps guide operational decisions. This work reduces technical risk and strengthens the scientific foundation needed to make geothermal energy a reliable and scalable clean energy resource.

15 GEOTHERMAL ENERGY↗

Friction Stir Welding of Curved Plates

Friction stir welding (FSW) is a remarkable technology for making butt and lap joints in aluminum alloys. The process operates by passing a rotating tool between two closely butted plates. This process generates heat and the heated material is stirred from both sides of the plates to generate a high quality weld. Application of this technique has a very broad field for NASA. In particular, NASA is interested in using this welding process to manufacture tanks and curved elements. Therefore, this research has been oriented to the study the FSW of curved plates. The study has covered a number of topics that are important in the model development and to uncover the physical process involve in the welding itself. The materials used for the experimental welds were as close to each other as we could possibly find, aluminum 5454-0 and 5456-0 with properties listed at http://matweb.com. The application of FSW to curved plates needs to consider the behavior that we observed in this study. There is going to be larger force in the normal direction (Fz) as the curvature of the plate increases. A particular model needs to be derived for each material and thickness. A more complete study should also include parameters such as spin rate, tool velocity, and power used. The force in the direction of motion (Fx) needs to be reconsidered to make sure of its variability with respect to other parameters such as velocity, thickness, etc. It seems like the curvature does not play a role in this case. Variations in temperature were found with respect to the curvature. However, these changes seem to be smaller than the effect on Fz. The temperatures were all below the melting point. We understand now that the process of FSW produces a three dimensional flow of material that takes place during the weld. This flow needs to be study in a more detailed way to see in which directions the flow of material is stronger. It could be possible to model the flow using a 2-dimensional model in the particular directions where the flow moves faster. More experimental information is required to enrich the knowledge about FSW, and from this point, derive useful mathematical formulas to optimize the process and the design of the machines that will perform it. More experiments and experimental equipment are required to uncover the mathematics of the process.

Sanchez, Nestor↗

Modeling the Effects of a Backward-Facing Step on Boundary-Layer Transition

We model transition to turbulence in a two-dimensional boundary layer downstream of a backward-facing step (BFS) along a flat plate. With the goal of evaluating the available engineering models for predicting the effects of step excrescences on the transition characteristics, two separate methodologies are used to monitor the streamwise shift in the transition onset location as the step height and the flow speed are varied across the range of a previously reported experiment involving step-height-to-local-displacement-thickness ratios of 0 < h/δ* < 1.6. Unlike the variable N -factor method from the previous literature, both of these methods are general in scope and do not involve any empirical correlations that are specific to step excrescences. The first of these techniques involves an N -factor method that directly accounts for the change in boundary-layer instability characteristics due to the step. Stability computations using the harmonic linearized Navier-Stokes equations (HLNSE), which fully account for the nonparallel-mean-flow effects close to the BFS, indicate that the measured transition locations at nearly all test conditions ( h/δ* < 1.3) correlate well with a computed N -factor of N tr = 7.6, demonstrating a successful stability-based transition criterion related to step excrescences. Linear stability theory, which does not account for nonparallel effects, demonstrates reasonable agreement with the HLNSE results, yielding good predictions for the overall trends, but predicts a somewhat earlier onset of transition than HLNSE. The other methodology used in this work involves transport-equation-based transition models. We first show that the Langtry-Menter y - Re θt transition model cannot accurately predict the location of transition onset for moderate BFS heights because it is unable to accurately account for the flow history effects. Along with the Langtry-Menter transition model, we also show the amplification factor transport model does not produce accurate transition locations for subsonic flow over steps even though it accounts for some flow history effects.

Transition↗