Search NASA⌕ Search

SEARCH · Search NASA

Results for “Transonic”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparing theoretical predictions of radiation-free velocities of edge dislocations to molecular dynamics simulations

Transonic defect motion is of interest for high strain-rate plastic deformation as well as for crack propagation. Ever since Eshelby's 1949 prediction [J. D. Eshelby, Proc. Phys. Soc. A 62, 307 (1949)] in the isotropic limit of a “radiation-free” transonic velocity $v$ RF = $\sqrt{2}$c T , where shock waves are absent, there has been speculation about the significance of radiation-free velocities (if they truly exist) for defect mobility. Here, in this work, we argue that they do not play any significant role in dislocation dynamics in metals, based on comparing theoretical predictions of radiation-free velocities for transonic edge dislocations with molecular dynamics simulations for two face-centered cubic metals: Ag, where theory predicts radiation-free states, and Cu, where it does not.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2024

Mobility of twinning dislocations in copper up to supersonic speeds

Understanding the mobility of twinning dislocations is important for multiscale modeling of crystal plasticity, especially at high strain rates, where such dislocations may reach transonic or supersonic speeds. Here, we used molecular dynamics simulations to investigate the relationship between dislocation velocity and the applied resolved shear stress of an edge twinning dislocation in copper up to supersonic speeds. The twinning dislocation mobility relation is composed of two branches separated by a band of forbidden velocities. The lower velocity branch is limited by the first transverse sound speed ~2000 m/s while the upper branch stretches from ~3500 m/s in the transonic regime to supersonic velocities. Twinning dislocations cannot undergo uniform steady-state motion at velocities within the forbidden band. Our simulation results also reveal that edge twinning dislocation motion in copper is kink-mediated. We discuss the implications of our findings for the motion of twins, twinning dislocations, and twinning dislocation kinks in copper.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Revolutionizing Turbine Cooling with Micro-Architectures Enabled by Direct Metal Laser Sintering

The objective of this research effort was to explore innovative cooling architectures enabled by additive manufacturing techniques for improved turbine cooling performance. The ability to create complex internal geometries was leveraged to better distribute coolant as well as to integrate inherently unstable flow devices to enhance internal and external heat transfer. This was accomplished with a multi-faceted approach including analytical, experimental, and computational components. This final report documents progress during the total 4.25 year effort that was extended (at no additional cost) from the original 3-year cooperative agreement. Prior to starting the project, significant effort was invested in investigating innovative cooling designs from the open literature. When the project started, we quickly down-selected to 4 configurations: impingement jet fluidic oscillators (for internal leading edge cooling), reverse film cooling (for pressure surface), sweeping film cooling jets (for suction surface), and trailing edge slot cooling with microchannels. During the 2nd year, the reverse film cooling design was shown to be sub-optimal and work on that topic was halted. The three remaining technologies were integrated into a large scale nozzle guide vane and installed into a low-speed linear cascade facility for further interrogation. Initial results showed the sweeping film cooling jets (on the suction surface) to be the most promising technology compared to a baseline diffusion shaped film hole (777 design). The benefit was particularly evident at high blowing ratios (>1.5) when the 777 coolant film separated from the downstream surface. Benefits were also evident at elevated freestream turbulence levels. The leading edge cooling with unsteady jets was less effective in terms of peak or average cooling – however it was superior to round hole impingement cooling in terms of spatial uniformity. Finally, the trailing edge cooling design with micro-channels was scrapped in favor of a pinned arrangement with centerbody. During the 3rd year, the same three technologies were integrated into a transonic linear cascade for an assessment of compressibility effects. Again, the sweeping film cooling jets proved superior at high blowing ratios while the unsteady leading edge impingement and trailing edge pin-fins with centerbody designs yielded mixed results. Finally, during the 4th and final year, these 3 technologies were integrated into a direct metal laser sintered (DMLS) nozzle guide vane for testing in a high temperature, transonic annular vane cascade. The facility matches the flow temperature, Mach number, and coolant pressure ratios of an actual gas turbine. Before doing this, an additional series of tests were completed to validate a redesigned trailing edge cooling architecture. This final design included a centerbody with triangular pins between the centerbody and the vane external skin. The majority of these triangular pins are fabricated with a 30% gap to the centerbody – reducing their pressure drop considerably while still providing excellent heat transfer augmentation. This additional test campaign required an additional 3-month extension request to complete testing in the high temperature NGV test facility (Turbine Reacting Flow Rig – TuRFR). Testing of the DMLS vane in TuRFR pitted the 3 innovative cooling technologies against more traditional technologies (shaped 777 film hole, round impingement jet, and full pin arrangement in the trailing edge. The sweeping film cooling jets (on the suction surface) showed approximately 15% improvement in overall effectiveness compared to the 777 film hole while the trailing edge showed up to 20% improvement. The leading edge impingement was disappointing with a substantially lower effectiveness than traditional direct impingement. In summary, the study was successful in demonstrating that DMLS-enabled cooling technologies can yield significant gains in cooling performance.

03 NATURAL GAS↗
DOE OSTI2020

Aerodynamic Data Fusion Toward the Digital Twin Paradigm

This paper considers the fusion of two aerodynamic data sets originating from differing types of physical or computer experiments. This paper specifically addresses the fusion of 1) noisy and in-complete fields from wind-tunnel measurements and 2) deterministic but biased fields from numerical simulations. These two data sources are fused in order to estimate the true field that best matches measured quantities that serve as the ground truth. For example, two sources of pressure fields about an aircraft are fused based on measured forces and moments from a wind-tunnel experiment. A fundamental challenge in this problem is that the true field is unknown and cannot be estimated with 100% certainty. A Bayesian framework is employed to infer the true fields conditioned on measured quantities of interest; essentially a statistical correction to the data is performed. The fused data may then be used to construct more accurate surrogate models suitable for early stages of aerospace design. An extension of the proper orthogonal decomposition with constraints is also introduced to solve the same problem. In this work, both methods are demonstrated on fusing the pressure distributions for flow past the RAE2822 airfoil and the Common Research Model wing at transonic conditions. Comparison of both methods reveals that the Bayesian method is more robust when data are scarce and capable of also accounting for uncertainties in the data. Furthermore, given adequate data, the proper-orthogonal-decomposition-based and Bayesian approaches lead to surprisingly similar results.

42 ENGINEERING↗
DOE OSTI2020

Variable-fidelity multipoint aerodynamic shape optimization with output-based adapted meshes

This work presents a method to control the discretization error in multipoint aerodynamic shape optimization using output-based adapted meshes. The meshes are adapted via adjoint-based error estimates, taking into account both the objective and constraint output errors. A multi-fidelity optimization framework is then developed by taking advantage of the variable fidelity offered by adaptive meshes. The objective functional and its sensitivity at each design point (operating condition) are first evaluated on the same initial coarse mesh, which is then subsequently adapted for each design point individually as the shape optimization proceeds. The effort to set up the optimization is minimal since the initial mesh can be fairly coarse and easy to generate. As the shape approaches the optimal design, the mesh at each design point becomes finer, in regions necessary for that particular operating condition. The multi-fidelity framework is tightly coupled with the objective error estimation to ensure the optimization accuracy at each fidelity. Computational savings arise from a reduction of the mesh size when the design is far from optimal and avoiding an exhaustive search on low-fidelity meshes. The proposed method is demonstrated on multipoint drag minimization problems of a transonic airfoil with lift and area constraints. Improved accuracy and efficiency are shown compared to traditional fixed-fidelity optimization with a fixed computational mesh.

42 ENGINEERING↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Numerical study on spray collapse process of ECN spray G injector under flash boiling conditions

This work features a numerical study of the spray collapse process that occurs under flash-boiling, gasoline-direct-injection conditions. The Engine Combustion Network's Spray G injector was modeled, and iso-octane, n-hexane, and propane were used as fuels with the goal to investigate fuel property effects. The development of individual jets and jet-to-jet interaction were analyzed to reveal the physical mechanisms of spray collapse. The results agreed qualitatively well with the experimental data of propane available from the literature and the predictions revealed the model's ability to capture the near-nozzle shock structures witnessed in the experiments. Simulations showed that the flow development of flashing sprays was strongly tied with the fuel properties. Specifically, under-expanded jets were not observed with iso-octane due to the fuel's low vapor pressure. In the case of n-hexane, the jets were under-expanded, and a low-pressure core was observed inside the jets due to the flow expansion. With propane, the very strong expansion caused the flow to become transonic and generate a unique shock structure near the nozzle exit. The shock structure revealed the formation of primary shock cells within the jets and secondary shock cells induced by jet-to-jet interaction. Based on the results, the spray collapse under flash boiling conditions was mainly attributed to two effects, i.e., the decrease in pressure at the spray center induced by the low-pressure cores inside and between the jets, and the isolation effect due to the jet-to-jet interaction which contributed to maintain the pressure differential between the spray center and the ambient.

42 ENGINEERING↗
DOE OSTI2020

Towards understanding the development and characteristics of under-expanded flash boiling jets

Here in this research, the n-hexane jets discharged from a straight channel were studied numerically to reveal the development and characteristics of under-expanded flashing jets. Firstly, the development of flash boiling jets was investigated and compared to that of two-phase jets without phase transition. It was shown that the two-phase jets without phase transition can be under-expanded only if the flow is choked. In contrast, flash boiling jets can be under-expanded whether or not the choked conditions are verified. Secondly, the factors affecting the flashing jet behaviors were discussed. The relaxation time in phase transition, the expansion ratio, and the injection pressure have been identified as three key factors determining the jet characteristics. With high enough expansion ratio and low enough relaxation time, the expansion was so strong that the flow became transonic and shock structures similar to those of under-expanded gaseous jets appeared near the orifice exit. Finally, two types of choking have been identified for flash boiling jets, i.e., vaporization-induced choking and expansion-induced choking. The former one is mainly ascribed to the vapor formation in the orifice which reduces the local speed of sound, while the latter one is due to the expansion of under-expanded two-phase flow and the resultant flow acceleration.

42 ENGINEERING↗
DOE OSTI2022

CAMERA: A method for cost-aware, adaptive, multifidelity, efficient reliability analysis

Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distributions, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. Here, we propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential acquisition function based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest at the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set and probability estimation on synthetic test functions and two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier-Stokes equations. We show that our method predicts the failure boundary and probability more accurately and at a fraction of the computational cost compared with using just a single expensive high-fidelity model. Finally, we show that our sequential approach is guaranteed to asymptotically converge to the true failure boundary with high probability.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
DOE OSTI2020

Scalability of OpenFOAM Density-Based Solver with Runge–Kutta Temporal Discretization Scheme

Compressible density-based solvers are widely used in OpenFOAM, and the parallel scalability of these solvers is crucial for large-scale simulations. In this paper, we report our experiences with the scalability of OpenFOAM’s native rhoCentralFoam solver, and by making a small number of modifications to it, we show the degree to which the scalability of the solver can be improved. The main modification made is to replace the first-order accurate Euler scheme in rhoCentralFoam with a third-order accurate, four-stage Runge-Kutta or RK4 scheme for the time integration. The scaling test we used is the transonic flow over the ONERA M6 wing. This is a common validation test for compressible flows solvers in aerospace and other engineering applications. Numerical experiments show that our modified solver, referred to as rhoCentralRK4Foam, for the same spatial discretization, achieves as much as a 123.2% improvement in scalability over the rhoCentralFoam solver. As expected, the better time resolution of the Runge–Kutta scheme makes it more suitable for unsteady problems such as the Taylor–Green vortex decay where the new solver showed a 50% decrease in the overall time-to-solution compared to rhoCentralFoam to get to the final solution with the same numerical accuracy. Finally, the improved scalability can be traced to the improvement of the computation to communication ratio obtained by substituting the RK4 scheme in place of the Euler scheme. All numerical tests were conducted on a Cray XC40 parallel system, Theta, at Argonne National Laboratory.

Li, Sibo↗
DOE OSTI2020

SIERRA/Aero User Manual - Version 4.56

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

42 ENGINEERING↗
DOE OSTI2020

SIERRA/Aero Theory Manual - Version 4.56

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

42 ENGINEERING↗
DOE OSTI2020

SIERRA/Aero Theory Manual - Version 4.58

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

42 ENGINEERING↗
DOE OSTI2020

SIERRA/Aero User Manual - Version 4.58

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

42 ENGINEERING↗