Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence optimization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The right conditions for high-precision dynamic temperature and heat capacity measurement via pyrometry and conductivity

The pursuit of accurate bulk temperature T under extreme conditions has been a long-standing goal of the high pressure science community, complicated by a lack of data to inform models. To reach these extremely high-pressure, high-temperature (high P − T) conditions, a combination of dynamic and heated static experiments (e.g., diamond or gem anvil cel experiments) are used. For example, in a diamond anvil cell (DAC) experiment, a sample placed in the DAC is first pressurized. Following pressurization, the sample T is increased either by heating the entire DAC (usually using resistive heating, and limited to ∼1000K) or by applying intense laser power to the sample surfaces. In a dynamic experiment, the process of pressurizing the sample also heats it. In the case of shock physics experiments, such heating is substantial, easily reaching thousands of Kelvin; in our work we have seen T ∼17000K. Most methods of measuring temperature at ambient are not compatible with experiments under these high-pressure, high-temperature conditions: thermocouples break, melt, or have conductivity properties that differ from ambient where they are calibrated; thermometers would melt; both are too slow. As a result most methods are based on non-contact techniques such as x-ray diffraction broadening, neutron scattering, or optical methods. Of these, optical methods using the visible and near-infrared region of the spectrum are the most commonly used as the sources and detectors are readily available. In the case of optical methods the optical depth, and therefore the measurement location, is limited to the surface. When a window or anvil material is used, heat flows from the sample into the window/anvil. Likewise, if the sample undergoes a change in thermodynamic state, such as expansion upon release, different T may be expected. As a result, the surface or apparent temperature T app measurement will differ from the bulk or interior temperature that is desired. This surface measurement must be related to the bulk measurement using thermal transport models and material models. While it is tempting to conclude that one should just use x-ray methods that directly probe the interior, even these methods have been shown to depend on thermal transport and material models. Regardless of the method used to create the high P − T condition, therefore, we must understand the role of thermal transport and material models upon our interpretation of the T measurement, as well as the errors and uncertainties associated with the choice of models used in the analysis. This is a substantial area of research and this paper is by no means a complete survey of the relevant sources of uncertainty. For example, we have yet to begin to address alternate transport models in a detailed manner (e.g., Tan-Ahrens), or the many models that use additional layers to approximate melting, turbulence, or epitaxial phenomena). Likewise, we have not explored the impact upon uncertainty of thermal models that use temperature-dependent thermal transport coefficients, or the wide range of material models that can be applied. Instead, this paper focuses on using one simple model, the Urtiew-Grover model, to understand the sources of error in T measurement so that we may identify how best to focus future research efforts to return the best improvements and avoid working on over-optimizing a single type of measurement. To this end, we work through some of the best and worst case scenarios for T measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantifying and attributing time step sensitivities in present-day climate simulations conducted with EAMv1

Abstract. This study assesses the relative importance of time integration error in present-day climate simulations conducted with the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) at 1∘ horizontal resolution. We show that a factor-of-6 reduction of time step size in all major parts of the model leads to significant changes in the long-term mean climate. Examples of changes in 10-year mean zonal averages include the following: up to 0.5 K of warming in the lower troposphere and cooling in the tropical and subtropical upper troposphere, 1 %–10 % decreases in relative humidity throughout the troposphere, and 10 %–20 % decreases in cloud fraction in the upper troposphere and decreases exceeding 20 % in the subtropical lower troposphere. In terms of the 10-year mean geographical distribution, systematic decreases of 20 %–50 % are seen in total cloud cover and cloud radiative effects in the subtropics. These changes imply that the reduction of temporal truncation errors leads to a notable although unsurprising degradation of agreement between the simulated and observed present-day climate; to regain optimal climate fidelity in the absence of those truncation errors, the model would require retuning. A coarse-grained attribution of the time step sensitivities is carried out by shortening time steps used in various components of EAM or by revising the numerical coupling between some processes. Our analysis leads to the finding that the marked decreases in the subtropical low-cloud fraction and total cloud radiative effect are caused not by the step size used for the collectively subcycled turbulence, shallow convection, and stratiform cloud macrophysics and microphysics parameterizations but rather by the step sizes used outside those subcycles. Further analysis suggests that the coupling frequency between the subcycles and the rest of EAM significantly affects the subtropical marine stratocumulus decks, while deep convection has significant impacts on trade cumulus. The step size of the cloud macrophysics and microphysics subcycle itself appears to have a primary impact on cloud fraction in the upper troposphere and also in the midlatitude near-surface layers. Impacts of step sizes used by the dynamical core and the radiation parameterization appear to be relatively small. These results provide useful clues for future studies aiming at understanding and addressing the root causes of sensitivities to time step sizes and process coupling frequencies in EAM. While this study focuses on EAMv1 and the conclusions are likely model-specific, the presented experimentation strategy has general value for weather and climate model development, as the methodology can help researchers identify and understand sources of time integration error in sophisticated multi-component models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1972

Some effects of swirl on turbulent mixing and combustion

A general formulation of some effects of swirl on turbulent mixing is given. The basis for the analysis is that momentum transport is enhanced by turbulence resulting from rotational instability of the fluid field. An appropriate form for the turbulent eddy viscosity is obtained by mixing length type arguments. The result takes the form of a corrective factor that is a function of the swirl and acts to increase the eddy viscosity. The factor is based upon the initial mixing conditions implying that the rotational turbulence decays in a manner similar to that of free shear turbulence. Existing experimental data for free jet combustion are adequately matched by using the modifying factor to relate the effects of swirl on eddy viscosity. The model is extended and applied to the supersonic combustion of a ring jet of hydrogen injected into a constant area annular air stream. The computations demonstrate that swirling the flow could: (1) reduce the burning length by one half, (2) result in more uniform burning across the annulus width, and (3) open the possibility of optimization of the combustion characteristics by locating the fuel jet between the inner wall and center of the annulus width.

Rubel, A.↗
NASA NTRS1994

Merging Intelligent Systems Technologies with CFD Analysis Strategies: Prototype Development at NASA Ames

A collaborative team of researchers from fields of Computational Fluid Dynamics (CFD), fluid physics, computer architectures, and computer science and knowledge engineering have begun work on a prototype system that addresses several of industry's concerns in using NASA-developed CFD codes as part of the design cycle. A major problem exists in the application of CFD technologies within the aeronautics design cycle due primarily to misunderstandings in the ranges of applicability of the various solver codes or turbulence models. Features that arise during the CFD solution process need to be discriminated and recognized as actual flow features with physical support in the geometry and flow conditions of the problem being solved, or as numerical or non-physical errors arising from mis-application of solver code and its parameters, gridding strategies, or discretization. interpolations. The fundamental concept is to develop an intelligent computational system that can accept the engineer's definition of the problem and construct an optimal CFD solution. To do this requires capturing both the knowledge of how to apply the various CFD tools and how to adapt the application of those tools to flow structures as they evolve during the flow simulation. Embedded within this adaptive system approach is the additional desire to automatically identify and quantify the quality of resolution of the pertinent flow structures, be they genuine or error-induced, and then to adjust the solution strategy accordingly. This paper discusses the status of that prototyping effort.

Thompson, David E.↗
NASA NTRS2004

Review of Research into the Concept of the Microblowing Technique for Turbulent Skin Friction Reduction

A new technology for reducing turbulent skin friction, called the Microblowing Technique (MBT), is presented. Results from proof-of-concept experiments show that this technology could potentially reduce turbulent skin friction by more than 50% of the skin friction of a solid flat plate for subsonic and supersonic flow conditions. The primary purpose of this review paper is to provide readers with information on the turbulent skin friction reduction obtained from many experiments using the MBT. Although the MBT has a penalty for obtaining the microblowing air associated with it, some combinations of the MBT with suction boundary layer control methods are an attractive alternative for a real application. Several computational simulations to understand the flow physics of the MBT are also included. More experiments and computational fluid dynamics (CFD) computations are needed for the understanding of the unsteady flow nature of the MBT and the optimization of this new technology.

Source record↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Gyrokinetic profile prediction and validation of a negative triangularity plasma in ASDEX Upgrade

In this work, gyrokinetic simulations are performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73–93) for a negative triangularity H-mode plasma in ASDEX Upgrade, and compared with experimental measurements. The PORTALS framework (Rodriguez-Fernandez et al 2024 Nucl. Fusion 64 076034) is used to accelerate the prediction of kinetic profiles for this plasma, using surrogate modeling and Bayesian optimization. Ion heat flux, electron heat flux, and electron particle flux are simultaneously matched across the simulated radial regime of the plasma (normalized radius $r/a = 0.35-0.90$), and the resulting ion temperature, electron temperature, and electron density profiles match well with the experimental profile data within this radial range. A synthetic Correlation Electron Cyclotron Emission diagnostic is applied to find well-matched electron temperature fluctuation properties between simulation and experiment. The flux-matched profiles provide a basis for investigation of the turbulence nature across the plasma radius, revealing the dominance of Trapped Electron Mode turbulence at $r/a = 0.35$, the dominance of Ion Temperature Gradient turbulence at $r/a = 0.55$, 0.75, and 0.83, and an instability boundary at $r/a = 0.90$.

gyrokinetic simulation↗
NASA NTRS2000

Use Computer-Aided Tools to Parallelize Large CFD Applications

Porting applications to high performance parallel computers is always a challenging task. It is time consuming and costly. With rapid progressing in hardware architectures and increasing complexity of real applications in recent years, the problem becomes even more sever. Today, scalability and high performance are mostly involving handwritten parallel programs using message-passing libraries (e.g. MPI). However, this process is very difficult and often error-prone. The recent reemergence of shared memory parallel (SMP) architectures, such as the cache coherent Non-Uniform Memory Access (ccNUMA) architecture used in the SGI Origin 2000, show good prospects for scaling beyond hundreds of processors. Programming on an SMP is simplified by working in a globally accessible address space. The user can supply compiler directives, such as OpenMP, to parallelize the code. As an industry standard for portable implementation of parallel programs for SMPs, OpenMP is a set of compiler directives and callable runtime library routines that extend Fortran, C and C++ to express shared memory parallelism. It promises an incremental path for parallel conversion of existing software, as well as scalability and performance for a complete rewrite or an entirely new development. Perhaps the main disadvantage of programming with directives is that inserted directives may not necessarily enhance performance. In the worst cases, it can create erroneous results. While vendors have provided tools to perform error-checking and profiling, automation in directive insertion is very limited and often failed on large programs, primarily due to the lack of a thorough enough data dependence analysis. To overcome the deficiency, we have developed a toolkit, CAPO, to automatically insert OpenMP directives in Fortran programs and apply certain degrees of optimization. CAPO is aimed at taking advantage of detailed inter-procedural dependence analysis provided by CAPTools, developed by the University of Greenwich, to reduce potential errors made by users. Earlier tests on NAS Benchmarks and ARC3D have demonstrated good success of this tool. In this study, we have applied CAPO to parallelize three large applications in the area of computational fluid dynamics (CFD): OVERFLOW, TLNS3D and INS3D. These codes are widely used for solving Navier-Stokes equations with complicated boundary conditions and turbulence model in multiple zones. Each one comprises of from 50K to 1,00k lines of FORTRAN77. As an example, CAPO took 77 hours to complete the data dependence analysis of OVERFLOW on a workstation (SGI, 175MHz, R10K processor). A fair amount of effort was spent on correcting false dependencies due to lack of necessary knowledge during the analysis. Even so, CAPO provides an easy way for user to interact with the parallelization process. The OpenMP version was generated within a day after the analysis was completed. Due to sequential algorithms involved, code sections in TLNS3D and INS3D need to be restructured by hand to produce more efficient parallel codes. An included figure shows preliminary test results of the generated OVERFLOW with several test cases in single zone. The MPI data points for the small test case were taken from a handcoded MPI version. As we can see, CAPO's version has achieved 18 fold speed up on 32 nodes of the SGI O2K. For the small test case, it outperformed the MPI version. These results are very encouraging, but further work is needed. For example, although CAPO attempts to place directives on the outer- most parallel loops in an interprocedural framework, it does not insert directives based on the best manual strategy. In particular, it lacks the support of parallelization at the multi-zone level. Future work will emphasize on the development of methodology to work in a multi-zone level and with a hybrid approach. Development of tools to perform more complicated code transformation is also needed.

Jin, H.↗
NASA NTRS2019

Optical Array Versus Monolithic Telescope Ground Station Cost Assessment

The purpose of this study is to outline the design of an optimal array of optical telescopes to emulate performance of a monolithic 12 m telescope in support of deep-space communications. In this case, optimal means minimizing the initial capital investment and operational cost while maintaining performance requirements of the deep-space link. The design is approached from a practical, engineering perspective. Pulse position modulation (PPM) signal formatting and photon counting detectors are assumed at each telescope in the array. That is, the telescopes function as so-called light buckets, so direct detection (as opposed to coherent reception) of the received signals is assumed, and there is no intention to consider active compensation for atmospheric turbulence-induced phase fluctuations. A parametric analysis among aperture size, detector size, and primary mirror surface quality, in the context of field-of-view expansion, is presented to minimize the cost function.

Romanofsky, Robert R.↗
NASA NTRS1972

An analytical study of reduced-gravity liquid reorientation using a simplified marker and cell technique

A computer program called HOPI was developed to predict reorientation flow dynamics, wherein liquids move from one end of a closed, partially filled, rigid container to the other end under the influence of container acceleration. The program uses the simplified marker and cell numerical technique and, using explicit finite-differencing, solves the Navier-Stokes equations for an incompressible viscous fluid. The effects of turbulence are also simulated in the program. HOPI can consider curved as well as straight walled boundaries. Both free-surface and confined flows can be calculated. The program was used to simulate five liquid reorientation cases. Three of these cases simulated actual NASA LeRC drop tower test conditions while two cases simulated full-scale Centaur tank conditions. It was concluded that while HOPI can be used to analytically determine the fluid motion in a typical settling problem, there is a current need to optimize HOPI. This includes both reducing the computer usage time and also reducing the core storage required for a given size problem.

Betts, W. S., Jr.↗
NASA NTRS2019

A Knowledge-Based Optimization Method for Aerodynamic Design

A new aerodynamic design method, CODISC, has been developed that combines a legacy knowledge-based design method, CDISC, with a simple optimization module known as SOUP. The primary goal of this new design system is to improve the performance gains obtained using CDISC without adding significant computational time. An additional objective of this approach is to reduce the need for a priori knowledge of good initial input variable values, as well as for subsequent manual revisions of those values as the design progresses. Several test cases illustrate the development of the process to date and some of the options available at transonic and supersonic speeds for turbulent flow designs. The test cases generally start from good baseline configurations and, in all cases, were able to improve the performance. Several new guidelines for good initial values for the design variables, as well as new design rules within CDISC itself, were developed from these cases.

SOUP↗
NASA NTRS2014

Assessment of Reduced-Kinetics Mechanisms for Combustion of Jet Fuel in CFD Applications

A computational effort was undertaken to analyze the details of fluid flow in Lean-Direct Injection (LDI) combustors for next-generation LDI design. The National Combustor Code (NCC) was used to perform reacting flow computations on single-element LDI injector configurations. The feasibility of using a reduced chemical-kinetics approach, which optimizes the reaction rates and species to model the emissions characteristics typical of lean-burning gas-turbine combustors, was assessed. The assessments were performed with Reynolds- Averaged Navier-Stokes (RANS) and Time-Filtered Navier Stokes (TFNS) time-integration, with a Lagrangian spray model with the NCC code. The NCC predictions for EINOx and combustor exit temperature were compared with experimental data for two different single-element LDI injector configurations, with 60deg and 45deg axially swept swirler vanes. The effects of turbulence-chemistry interaction on the predicted flow in a typical LDI combustor were studied with detailed comparisons of NCC TFNS with experimental data.

Ajmani, Kumud↗
DOE OSTI2025

Optimize the Autonomous Sensor Fish Device for Understanding Interactions of Aquatic Animals with Marine and Hydrokinetic and Hydro Systems (CRADA 476) Final Report

This project at the Pacific Northwest National Laboratory (PNNL) enhances our understanding of how aquatic life interacts with marine renewable energy (MRE) and hydropower systems through the development and testing of advanced sensor technology. The Marine Sensor Fish (MSF) and its smaller variants, the Sensor Fish Mini (SF Mini) and Flexible Sensor Fish (FSF), measure critical environmental stressors, such as shear forces, pressure changes, and strike impacts, that are essential for assessing the risks posed by energy devices, including tidal turbines, hydropower installations, and other MRE infrastructure. Rigorous field testing has shown that these tools are effective in capturing detailed data on turbulence, pressure variations, strike impacts, and other stressors near these devices, providing valuable insights into the conditions faced by aquatic species. Economically, these tools offer renewable energy developers a practical solution for streamlining environmental assessments and ensuring regulatory compliance, ultimately reducing costs associated with evaluating and mitigating ecological impacts. This research also benefits the public by promoting the growth of renewable energy sources that protect aquatic ecosystems, advancing both sustainable energy production and environmental stewardship.

13 HYDRO ENERGY↗
NASA NTRS1999

Viscous Analysis of Ref. H Based Wing/Bodies

Computations have been performed on the baseline Reference H wing/body configuration, as well as the Wing 704 configuration, an optimized wing and fuselage combination derived from Ref. H through automated optimization. The parabolized Navier-Stokes solver UPS was employed with viscous terms in two directions in an effort to understand the source and level of potential viscous/inviscid interactions. The paper briefly describes the UPS code and the grids used to obtain the solutions before the discussion of results. Results of these computations indicate that viscous/inviscid interaction can contribute increments to both the pressure- and friction-related drag. Computations were performed for wind tunnel conditions-1.675% scale models at a Reynolds number of 4 million per foot. Turbulent flow results were obtained using the Baldwin-Lomax algebraic turbulence model and were compared with laminar flow results. The laminar flow fields were used to obtain upper bounds on potential interaction effects.

Lawrence, Scott L.↗
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↗
NASA NTRS1972

Wall temperature distribution calculation for a rocket nozzle contour

The JANNAF turbulent boundary layer (TBL) computer program, applicable to rocket nozzles, requires a wall temperature distribution among other input parameters to determine boundary layer behavior, heat transfer, and performance degradation. The inclusion of a complete regenerative cooling cycle model with associate geometry, material and fluid property data provides a capability to internally calculate wall temperature profiles on the hot gas and coolant flow-side, as well as the coolant flow bulk temperature variation. Besides the regular heat transfer and performance degradation calculations, the new concept can be used to optimize the cooling cycle, flow requirements, and cooling jacket geometry.

Omori, S.↗
NASA NTRS1992

Fluid flow in a jet and the Ca XIX line profiles observed during solar flares

The present study considers two types of fluid models as methods for reproducing the blueshifts and line broadening observed in soft X-ray lines during the rise of solar flares. The fluid models representing the laminar flow of material in a jet and through a pipe were used to derive the velocity at each location in the flow. The optimal values for the free parameters in the jet model were within the ranges expected for solar flare conditions. The jet was much smaller than typical flare volumes estimated from spatial images, but the resulting densities were consistent with the upper limit of densities derived using density-sensitive line ratios. This supports previous reports of small filling factors and is interpreted as evidence for a chromospheric origin for the flow. Stability analysis indicates that, although the laminar jet can reproduce the observations, instabilities may cause the flow to become turbulent, and that this laminar treatment is not strictly valid.

Bornmann, P. L.↗