Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “THREE-DIMENSIONAL BOUNDARY LAYER”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Mixed finite-difference scheme for analysis of simply supported thick plates.

A mixed finite-difference scheme is presented for the stress and free vibration analysis of simply supported nonhomogeneous and layered orthotropic thick plates. The analytical formulation is based on the linear, three-dimensional theory of orthotropic elasticity and a Fourier approach is used to reduce the governing equations to six first-order ordinary differential equations in the thickness coordinate. The governing equations possess a symmetric coefficient matrix and are free of derivatives of the elastic characteristics of the plate. In the finite difference discretization two interlacing grids are used for the different fundamental unknowns in such a way as to reduce both the local discretization error and the bandwidth of the resulting finite-difference field equations. Numerical studies are presented for the effects of reducing the interior and boundary discretization errors and of mesh refinement on the accuracy and convergence of solutions. It is shown that the proposed scheme, in addition to a number of other advantages, leads to highly accurate results, even when a small number of finite difference intervals is used.

Noor, A. K.↗
NASA NTRS1972

Three dimensional potential and current distributions in a Hall generator with assumed velocity profiles

Three-dimensional potential and current distributions in a Faraday segmented MHD generator operating in the Hall mode are computed. Constant conductivity and a Hall parameter of 1.0 is assumed. The electric fields and currents are assumed to be coperiodic with the electrode structure. The flow is assumed to be fully developed and a family of power-law velocity profiles, ranging from parabolic to turbulent, is used to show the effect of the fullness of the velocity profile. Calculation of the square of the current density shows that nonequilibrium heating is not likely to occur along the boundaries. This seems to discount the idea that the generator insulating walls are regions of high conductivity and are therefore responsible for boundary-layer shorting, unless the shorting is a surface phenomenon on the insulating material.

Stankiewicz, N.↗
NASA NTRS1961

Boundary-Layer Transition on a Group of Blunt Nose Shapes at a Mach Number of 2.20

An investigation has been made to study boundary-layer transition on six axisymmetrical blunt bodies of revolution. Model shapes were selected with respect to the degree of favorable pressure gradient over the model surface. Tests were conducted at a Mach number of 2.20 and over a range of free-stream Reynolds number per foot of about 1.4 x 10(exp 6) to 6.5 x 10(exp 6). The tests were made at an angle of attack of 0 deg. with zero heat transfer. For the hemisphere, the flow remained essentially laminar over the model surface length for the entire pressure range of the tests. For a strong favorable pressure gradient followed by any weak favorable, neutral, or adverse gradient, the tendency was for transition to occur at or immediately behind the shoulder. A single strip of three-dimensional roughness in the region of strong favorable pressure gradient did not fix transition on the models at the roughness location except at the maximum test pressures, whereas a second roughness strip added in a region of neutral or adverse pressure gradient did fix transition. Experimental pressure coefficients agreed closely with modified Newtonian theory except in the shoulder region.

Jackson, Mary W.↗
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
DOE OSTI2026

The Role of Subcloud Mesoscale Convergence in Sculpting Convective Updraft Width and Depth

The initiation of deep moist convection is governed in part by the horizontal width of updrafts near cloud base, which limits the deleterious effects of entrainment-driven dilution on buoyant thermals ascending through the free troposphere. However, the factors controlling cloud-base updraft width, which in turn dictates cloud depth, are not well understood. We track the evolving three-dimensional structure of the mesoscale subcloud forcing for vertical motion and near-cloud thermodynamic ingredients within a high-resolution ensemble of simulations of seven realistic daytime orographic convection initiation events to determine their relative roles in controlling cloud width and depth. Statistical analysis of approximately 5000 cloudy updraft samples indicates that the most important contributors to the width of cloudy updrafts across the ensemble are the depth and magnitude of the subcloud mesoscale ascent. However, the depth achieved by clouds is more consistently predicted by the near-cloud ambient relative humidity within the lower to middle free troposphere and convective available potential energy. Therefore, although the width of cloudy updrafts may be partly set at low levels by the mesoscale vertical mass and moisture flux, the likelihood of deep moist convection is governed by the generation of positive buoyancy within cumulus thermals and entrainment-driven dilution that reduces it. The persistence of the low-level mesoscale vertical forcing locally consolidates and vertically transports boundary layer moisture, helping to reduce updraft dilution. However, these factors vary in relative impacts on cloudy updrafts across individual cases, indicating multiple pathways for deep convection initiation.

Convective storms↗
NASA NTRS1944

The "limiting line" in mixed subsonic and supersonic flow of compressible fluids

It is well known that the vorticity for any fluid element is constant if the fluid is non-viscous and the change of state of the fluid is isentropic. When a solid body is placed in a uniform stream, the flow far ahead of the body is irrotational. Then if the flow is further assumed to be isentropic, the vorticity will be zero over the whole filed of flow. In other words, the flow is irrotational. For such flow over a solid body, it is shown by Theodorsen that the solid body experiences no resistance. If the fluid has a small viscosity, its effect will be limited in the boundary layer over the solid body and the body will have a drag due to the skin friction. This type of essentially isentropic irrotational flow is generally observed for a streamlined body placed in a uniform stream, if the velocity of the stream is kept below the so-called "critical speed." At the critical speed or rather at a certain value of the ratio of the velocity of the undisturbed flow and the corresponding velocity of sound, shock waves appear. This phenomenon is called the "compressibility bubble." Along a shock wave, the change of state of the fluid is no longer isentropic, although still adiabatic. This results in an increase in entropy of the fluid and generally introduces vorticity in an originally irrotational flow. The increase in entropy of the fluid is, of course, the consequence of changing part of the mechanical energy into heat energy. In other words, the part of fluid affected by the shock wave has a reduced mechanical energy. Therefore, with the appearance of shock waves, the wake of the streamline body is very much widened, and the drag increases drastically. Furthermore, the accompanying change in the pressure distribution over the body changes the aerodynamic moment acting on it and in the case of an airfoil decreases the lift force. All these consequences of the breakdown of isentropic irrotational flow are generally undesirable in applied aerodynamics. Its occurrence should be delayed as much as possible by modifying the shape or contour of the body. However, such endeavor will be very much facilitated if the cause or the criterion for the breakdown can be found first.

MACH NUMBER↗
DOE OSTI2025

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Bio-inspired alula-based winglet design for enhanced heat transfer in high temperature fin-and-tube heat exchangers

Fin-and-tube heat exchangers (FTHEs) are widely used for high-temperature flue-gas heat recovery, but their performance is often limited by wake regions and non-uniform fin-surface temperatures. This study proposes and numerically evaluates four bio-inspired longitudinal vortex generator (VG) configurations in a high-temperature FTHE with flue-gas inlet temperature ∼1230 K: double-delta, curved double-delta, alula, and a new curved-alula geometry. The reference fin is not hydraulically plain; it already incorporates leading-edge separation columns and convex protrusions, so the alula-type winglets are assessed as downstream add-ons acting on a strongly disturbed flow. In a second step, perforations (one, two and three circular holes) are introduced into the curved-alula VGs to further tailor the flow field. Three-dimensional simulations with the Shear Stress Transpor (SST) $k - ω$ model, temperature-dependent flue-gas properties and conjugate conduction are carried out for gas-side Reynolds numbers $Re_g ≈ 8.0$ x $10^2 - 3.6$ x $10^3$ (mass flow rates 0.5 – 2.5 g/s), and the designs are compared in terms of surface heat flux, Nusselt number, friction factor and hydrothermal performance factor (HTPF). For this already-promoted fin, the additional downstream winglets provide moderate, incremental hydrothermal gains. At the highest Reynolds number, the best non-perforated design (curved-alula) increases surface heat flux from 1630.9 to 1794.7 kW/m² (∼ 10 % gain) and the Nusselt number from 227.6 to 242.6 (∼ 7 % gain), while the friction factor rises from 0.26 to about 0.30, yielding HTPF values close to unity (∼ 0.9 – 1.0). Introducing circular perforations into the curved-alula winglets acts mainly as a wake-bleeding refinement: the three-hole configuration provides a heat flux of 1824.7 kW/m² and a pressure drop of 127.9 Pa, with HTPF in the range ∼ 1.03 – 1.14 and a small (∼ 1 – 3 %) improvement over the solid curved-alula design. Flow-field analysis shows that the perforated curved-alula VGs shrink tube-wake regions, thin the thermal boundary layer and homogenize the fin-surface temperature (outlet-gas temperature ∼ 510 – 520 K and fin-surface temperature ∼ 420 – 421 K for the three-hole case). An optimal flue-gas mass flow rate of ∼ 1 g/s ($Re_g ≈ 1.5$ x $10^3$) is identified, beyond which additional heat-transfer gains are offset by rapidly increasing pressure losses. Overall, the results highlight that initial fin geometry and VG placement are as important as VG shape: alula-based winglets are expected to yield larger relative gains on simpler flat-fin layouts or when positioned closer to the fin leading edge and tube

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

A North Sea in Situ Evaluation of the Fitch Wind Farm Parameterization Within the Mellor-Yamada-Nakanishi-Niino and 3D Planetary Boundary Layer Schemes

Wind resource assessments and wind power forecasts that account for wind farm wakes are sensitive to the choice of planetary boundary layer (PBL) scheme. This work compares the one-dimensional Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme with a three-dimensional PBL (3DPBL) scheme, evaluating predictions made with both schemes against two sets of North Sea in situ observations of wind farm wakes. The optimal PBL scheme varies based on the observations (FINO1 tower vs. aircraft), the quantity of interest (wind speed vs. turbulence kinetic energy [TKE]), and the error metric (bias, centered root mean square error [cRMSE], R2, and earth mover's distance [EMD]). Whereas 3DPBL wind speeds outperform MYNN wind speeds with respect to the cRMSE at the FINO1 site located at a single point within the turbine rotor layer, 3DPBL TKE bias is larger than MYNN TKE bias when compared to aircraft observations taken 100 m above a wind farm. Wind speeds in the aircraft region are ambiguous with regard to which PBL scheme is optimal. Aircraft MYNN wind speeds outperform 3DPBL wind speeds with respect to R2 and cRMSE but underperform with respect to bias and EMD. Future evaluations across broader temporal and spatial scales may offer further insight into model differences.

17 WIND ENERGY↗
NASA NTRS1975

Numerical computation of two-dimensional viscous blunt body flows with an impinging shock

Two-dimensional, viscous, blunt body flows with an impinging shock wave are computed using a time-dependent, finite-difference method to solve the complete set of Navier-Stokes equations. The bow shock wave is treated as a discontinuity, while all interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer are automatically captured in the solution. Numerical results are presented for cases in which shock waves of different strengths are allowed to impinge on the flow field surrounding a circular cylinder resulting in different shock interference patterns. The two-dimensional results are compared qualitatively with existing three-dimensional experiments.

Tannehil, J. C.↗
NASA NTRS1975

Experiments on supersonic boundary-layer separation in three dimensions

A detailed experimental study of three-dimensional separation in supersonic flow is described. Three-dimensional wedges affixed to flat plates were used to generate the flow fields that were used to examine the relationship of cross flow to the extent of separation. Models with strong transverse pressure gradients resulted in a diminished extent of separation about the model plane of symmetry but demonstrated expansive growth of separation to a limit with distance from the plane of symmetry. A secondary flow region was found embedded in the separated flow.

Bachalo, W. D.↗
NASA NTRS1975

Supercritical wing sections II: A handbook

The numerical aspect of theoretical work on transonics and supercritical wing sections are compiled. A model of the trailing edge is introduced which eliminates the loss of 15 to 20 percent experienced with heavily aft-loaded models, and it is indicated how drag creep can be reduced at off-design conditions. A rotated finite difference scheme is presented which can handle supersonic as well as subsonic free stream Mach numbers and leads to an effective three-dimensional program for the computation of transonic flow past an oblique wing. In the case of two-dimensional flow, the method is extended to take into account the displacement thickness computed by a semiempirical turbulent boundary layer correction. A series of supercritical wing sections is discussed together with comparisons between experimental and theoretical data. Computer programs and a brief manual for their operation are listed. It is shown that the programs furnish a physically adequate computer simulation of the compressible flows that arise in problems of transonic aerodynamics.

Bauer, F.↗
NASA NTRS1975

Laminar separation on a blunted cone at high angles of attack

A new technique is reported for calculating the entire flow field on spherically blunted cones at high angles of attack and high laminar Reynolds numbers. An approximate system of parabolic equations obtained from the steady Navier-Stokes equations by assuming the viscous, streamwise derivative terms are small compared to the viscous normal and circumferential derivatives is the basis of the calculations. These equations are valid for both the inviscid and viscous regions, including the circumferential separation zone that develops on the leeward side at high angles of attack. Two different methods are used to obtain the initial conditions for these equations at the sphere cone tangency plane. For small nose Reynolds numbers, an axisymmetric merged layer solution around a sphere is rotated to provide a three-dimensional initial plane of data. For large nose Reynolds numbers, the nose region is solved using an inviscid, three dimensional time dependent solution combined with a boundary layer solution for the viscous flow. The computed flowfield including the leeward separation region is described and compared with data for a 7 deg half angle cone at 10 deg angle of attack, and a blunt 15 deg half angle cone at 15 deg angle of attack.

Lubard, S. C.↗
NASA NTRS1974

Computation of steady axisymmetric flow using a one-dimensional time-dependent method

An iterative numerical method for computing steady, three dimensional, viscous, compressible flow fields, about aerodynamic bodies was studied. In order to develop the iterative method economically, the primary emphasis was directed towards supersonic, axisymmetric flow. However, the technique readily extends to three spatial dimensions. The viscous flow field about a cone-cylinder-flare body was calculated and compared to existing experimental data. Numerical predictions of the cone boundary layer and the flow field shock wave structure agreed with corresponding measurements. A separation was calculated at the cylinder-flare junction in six iterations; however, the size of the vortex did not correspond to the measured size. It was estimated that fifty iterations would be required to properly define the vortex. It was concluded that the iteration technique is of limited value for plane two dimensional and axisymmetrix flows, but of great value for three-dimensional flows.

Walitt, L.↗
NASA NTRS1976

Numerical solution of supersonic three-dimensional free-mixing flows using the parabolic-elliptic Navier-Stokes equations

A numerical method is presented for solving the parabolic-elliptic Navier-Stokes equations. The solution procedure is applied to three-dimensional supersonic laminar jet flow issuing parallel with a supersonic free stream. A coordinate transformation is introduced which maps the boundaries at infinity into a finite computational domain in order to eliminate difficulties associated with the imposition of free-stream boundary conditions. Results are presented for an approximate circular jet, a square jet, varying aspect ratio rectangular jets, and interacting square jets. The solution behavior varies from axisymmetric to nearly two-dimensional in character. For cases where comparisons of the present results with those obtained from shear layer calculations could be made, agreement was good.

Hirsh, R. S.↗
NASA NTRS1977

A three-dimensional viscous/potential flow interaction analysis method for multi-element wings

An analysis method and computer program were developed for the calculation of the viscosity dependent aerodynamic characteristics of multi-element, finite wings in incompressible flow. A fully-three dimensional potential flow program is used to determine the inviscid pressure distribution about the configuration. The potential flow program uses surface source and vortex singularities to represent the inviscid flow. The method is capable of analysing configurations having at most one slat, a main element, and two slotted flaps. Configurations are limited to full span slats or flaps. The configuration wake is allowed to relax as a force free wake, although roll up is not allowed at this time. Once the inviscid pressure distribution is calculated, a series of boundary layer computations are made along streamwise strips.

Dvorak, F. A.↗