Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW MODELS”

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 1,765 records · Page 98

Improved representation of laminar and turbulent sheet flow in subglacial drainage models

Subglacial hydrology models struggle to reproduce seasonal drainage patterns that are consistent with observed subglacial water pressures and surface velocities. We modify the standard sheet-flow parameterization within a coupled sheet–channel subglacial drainage model to smoothly transition between laminar and turbulent flow based on the locally computed Reynolds number in a physically consistent way (the ‘transition’ model). We compare the transition model to standard laminar and turbulent models to assess the role of the sheet-flow parameterization in reconciling observed and modelled water pressures under idealized and realistic forcing. Relative to the turbulent model, the laminar and transition models improve seasonal simulations by increasing winter water pressure and producing a more prominent late-summer water pressure minimum. In contrast to the laminar model, the transition model remains consistent with its own internal assumptions across all flow regimes. Based on the internal consistency of the transition model and its improved performance relative to the standard turbulent model, we recommend its use for transient simulations of subglacial drainage.

54 ENVIRONMENTAL SCIENCES↗

Validation of a Node-Centered Wall Function Model for the Unstructured Flow Code FUN3D

In this paper, the implementation of two wall function models in the Reynolds averaged Navier-Stokes (RANS) computational uid dynamics (CFD) code FUN3D is described. FUN3D is a node centered method for solving the three-dimensional Navier-Stokes equations on unstructured computational grids. The first wall function model, based on the work of Knopp et al., is used in conjunction with the one-equation turbulence model of Spalart-Allmaras. The second wall function model, also based on the work of Knopp, is used in conjunction with the two-equation k-! turbulence model of Menter. The wall function models compute the wall momentum and energy flux, which are used to weakly enforce the wall velocity and pressure flux boundary conditions in the mean flow momentum and energy equations. These wall conditions are implemented in an implicit form where the contribution of the wall function model to the Jacobian are also included. The boundary conditions of the turbulence transport equations are enforced explicitly (strongly) on all solid boundaries. The use of the wall function models is demonstrated on four test cases: a at plate boundary layer, a subsonic di user, a 2D airfoil, and a 3D semi-span wing. Where possible, different near-wall viscous spacing tactics are examined. Iterative residual convergence was obtained in most cases. Solution results are compared with theoretical and experimental data for several variations of grid spacing. In general, very good comparisons with data were achieved.

Carlson, Jan-Renee↗

Application of a new K-tau model to near wall turbulent flows

A recently developed K-tau model for near wall turbulent flows is applied to a variety of test cases. The turbulent flows considered include the incompressible flat plate boundary layer with adverse pressure gradients, incompressible flow past a backward facing step, and the supersonic flat plate boundary layer at zero pressure gradient. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.

Abid, R.↗

Coarse Grid Modeling of Turbine Film Cooling Flows Using Volumetric Source Terms

The recent trend in numerical modeling of turbine film cooling flows has been toward higher fidelity grids and more complex geometries. This trend has been enabled by the rapid increase in computing power available to researchers. However, the turbine design community requires fast turnaround time in its design computations, rendering these comprehensive simulations ineffective in the design cycle. The present study describes a methodology for implementing a volumetric source term distribution in a coarse grid calculation that can model the small-scale and three-dimensional effects present in turbine film cooling flows. This model could be implemented in turbine design codes or in multistage turbomachinery codes such as APNASA, where the computational grid size may be larger than the film hole size. Detailed computations of a single row of 35 deg round holes on a flat plate have been obtained for blowing ratios of 0.5, 0.8, and 1.0, and density ratios of 1.0 and 2.0 using a multiblock grid system to resolve the flows on both sides of the plate as well as inside the hole itself. These detailed flow fields were spatially averaged to generate a field of volumetric source terms for each conservative flow variable. Solutions were also obtained using three coarse grids having streamwise and spanwise grid spacings of 3d, 1d, and d/3. These coarse grid solutions used the integrated hole exit mass, momentum, energy, and turbulence quantities from the detailed solutions as volumetric source terms. It is shown that a uniform source term addition over a distance from the wall on the order of the hole diameter is able to predict adiabatic film effectiveness better than a near-wall source term model, while strictly enforcing correct values of integrated boundary layer quantities.

Heidmann, James D.↗

Scaling and modeling of turbulent suspension flows

Scaling factors determining various aspects of particle-fluid interactions and the development of physical models to predict gas-solid turbulent suspension flow fields are discussed based on two-fluid, continua formulation. The modes of particle-fluid interactions are discussed based on the length and time scale ratio, which depends on the properties of the particles and the characteristics of the flow turbulence. For particle size smaller than or comparable with the Kolmogorov length scale and concentration low enough for neglecting direct particle-particle interaction, scaling rules can be established in various parameter ranges. The various particle-fluid interactions give rise to additional mechanisms which affect the fluid mechanics of the conveying gas phase. These extra mechanisms are incorporated into a turbulence modeling method based on the scaling rules. A multiple-scale two-phase turbulence model is developed, which gives reasonable predictions for dilute suspension flow. Much work still needs to be done to account for the poly-dispersed effects and the extension to dense suspension flows.

Chen, C. P.↗

Adaptation of multidimensional group particle tracking and particle wall-boundary condition model to the FDNS code

A particulate two-phase flow CFD model was developed based on the FDNS code which is a pressure based predictor plus multi-corrector Navier-Stokes flow solver. Turbulence models with compressibility correction and the wall function models were employed as submodels. A finite-rate chemistry model was used for reacting flow simulation. For particulate two-phase flow simulations, a Eulerian-Lagrangian solution method using an efficient implicit particle trajectory integration scheme was developed in this study. Effects of particle-gas reaction and particle size change to agglomeration or fragmentation were not considered in this investigation. At the onset of the present study, a two-dimensional version of FDNS which had been modified to treat Lagrangian tracking of particles (FDNS-2DEL) had already been written and was operational. The FDNS-2DEL code was too slow for practical use, mainly because it had not been written in a form amenable to vectorization on the Cray, nor was the full three-dimensional form of FDNS utilized. The specific objective of this study was to reorder to calculations into long single arrays for automatic vectorization on the Cray and to implement the full three-dimensional version of FDNS to produce the FDNS-3DEL code. Since the FDNS-2DEL code was slow, a very limited number of test cases had been run with it. This study was also intended to increase the number of cases simulated to verify and improve, as necessary, the particle tracking methodology coded in FDNS.

Chen, Y. S.↗

LES-Based Analysis of Film Cooling in a High-Pressure Turbine Rotor Including Effects of Purge and Tip Flow

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Secondary Flows↗

Wind Tunnel Testing of AFC over a Deflected Aileron on the High-Lift Common Research Model

Active flow control (AFC) using discrete fluidic actuators distributed along the span, just upstream of the deflected aileron of the 10% scale high-lift version of the Common Research Model (CRM-HL), was evaluated during a wind tunnel test in the NASA Langley 14- by 22-Foot Subsonic Tunnel. For this set of experiments, a new outboard section was fabricated incorporating a deflectable aileron. Aileron deflection angles of 0, 7.5, 16, and 25 were investigated. This experimental investigation is in response to a recent study by Boeing potential to use AFC with a drooped aileron during takeoff to improve lift-to-drag ratio (L/D) by as much as 5%, depending on AFC mass flow rate. AFC is used at deflection angles above the nominal deflection of 7.5 to control the resulting flow separation that occurs and potentially improve L/D. Improvements in aircraft low-speed L/D can affect vehicle range and/or payload. Tuft flow visualization data, steady and unsteady surface pressure data, and force and moment data are used to compare the flowfields with and without AFC. Parameters varied include actuator momentum coefficient and aileron deflection angle. The surface pressure and tuft flow visualization results indicate that without AFC, the flow over a majority of the aileron upper surface is separated for aileron deflection angles larger than 7.5. When AFC is applied, the flow is reattached to the aileron. Force and moment results show that a local increase in lift leads to an increase in L/D of at least 3.5% using mass flow rates thought to be available from air sources onboard a commercial transport.

AFC↗

Wind Tunnel Testing of AFC over a Deflected Aileron on the High-Lift Common Research Model

Active flow control (AFC) using discrete fluidic actuators distributed along the span, just upstream of the deflected aileron of the 10% scale high-lift version of the Common Research Model (CRM-HL), was evaluated during a wind tunnel test in the NASA Langley 14- by 22-Foot Subsonic Tunnel. For this set of experiments, a new outboard section was fabricated incorporating a deflectable aileron. Aileron deflection angles of 0, 7.5, 16, and 25 were investigated. This experimental investigation is in response to a recent study by Boeing potential to use AFC with a drooped aileron during takeoff to improve lift-to-drag ratio (L/D) by as much as 5%, depending on AFC mass flow rate. AFC is used at deflection angles above the nominal deflection of 7.5 to control the resulting flow separation that occurs and potentially improve L/D. Improvements in aircraft low-speed L/D can affect vehicle range and/or payload. Tuft flow visualization data, steady and unsteady surface pressure data, and force and moment data are used to compare the flowfields with and without AFC. Parameters varied include actuator momentum coefficient and aileron deflection angle. The surface pressure and tuft flow visualization results indicate that without AFC, the flow over a majority of the aileron upper surface is separated for aileron deflection angles larger than 7.5. When AFC is applied, the flow is reattached to the aileron. Force and moment results show that a local increase in lift leads to an increase in L/D of at least 3.5% using mass flow rates thought to be available from air sources onboard a commercial transport.

AFC↗

Measurement of Surface Deformation Related to the December 2018 Mt. Etna Eruption Using Time-Series Interferometry and Magma Modeling for Hazard Zone Mapping

Mount Etna has erupted several times since it was first formed. Recently, Mount Etna began erupting again over 24–27 December 2018. Because it erupts frequently, Mount Etna should be observed on a frequent basis. From June 2018 to October 2019, 34 and 56 interferometric synthetic aperture radar (InSAR) images were acquired from the ascending and descending tracks of the Sentinel-1 satellite, respectively. We employed the Stanford Method for Persistent Scatterers (StaMPS) and a refined small baseline subset (SBAS) InSAR method to produce a surface deformation time-series map. In the time-series analysis, the phase signal remained unaltered with time. The Okada model was then applied to the result to generate a modeled interferogram, and the Q-LavHA program was run to generate a lava flow prediction model. A direct comparison of the results showed that Persistent Scatterers Interferometry (PSI)-StaMPS and the refined SBAS technique were comparable in terms of the displacement pattern, with slightly different velocity values obtained for individual points. In particular, a velocity range of −25 to 21 cm/yr was obtained from PSI-StaMPS, whereas a range of −30 to 25 cm/yr was obtained from the refined SBAS method. Upon computation of the vertical and east-west displacement components based on ascending and descending track data using both methods, deformation velocities of 51.5 and 52.5 cm/yr in the westerly direction on the western flank of Mount Etna were obtained from PSI-StaMPS and the refined SBAS method, respectively, whereas on the eastern flank, deformation toward the east was estimated to occur at a velocity of 50.1 or 54.2 cm/yr, respectively. PSIStaMPS estimated a vertical deformation velocity of −5.3 to 18.3 cm/yr, whereas the refined SBAS method produced a velocity range of approximately −7 to 19 cm/yr. The interferogram obtained via Okada modeling showed two fault sources in the 2018 Mount Etna eruption and a total volume change of approximately 12.39 × 10 6 m 3 . From the modeling results, a lava flow prediction model was generated using the Q-LavHA program. The approaches described in this study can be used by government officials, authorities, and other decision-makers to monitor and assess the risk of volcanic activity in the region.

Suci Ramayanti↗

Field Measurements of the 1983 Royal Gardens Lava Flows, Kilauea Volcano, and 1984 Mauna Loa Lava Flow, Hawaii

Theoretical models used in the remote determination of lava flow rheology and compositions rely on estimates of such geometric and flow parameters as volume flow rates, levee heights, and channel dimensions, as well as morphologic and structural patterns on the flow surfaces. Quantitative measures of these variables are difficult to obtain, even under optimum conditions. Detailed topographic profiles across several Hawaiian lava flows that were carefully monitored by the U.S. Geological Survey during their emplacement in 1983 were surveyed in order to test various flow emplacement models. Twenty two accurate channel cross sections were constructed by combining these profiles with digitized pre-flow topographic measurements. Levee heights, shear zone widths, and flow depths could then be read directly from the cross sections and input into the models. The profiles were also compared with ones constructed for some Martian lava flows.

Fink, J.↗

On turbulence modeling for unsteady transonic flows

A detailed examination of the turbulent field in an unsteady transonic flow undergoing shock-induced separation is conducted. Ensemble-averaged mean and fluctuating velocities, obtained from conditionally sampled laser velocimeter data, are described and analyzed to assess the applicability of modeling concepts usually employed in steady-flow problems. Some comparisons with computations employing the Reynolds averaged Navier-Stokes equations with a mixing length turbulence model are then presented to illustrate the status of current predictive capabilities. The results appear to imply that turbulence models developed for steady flows apply and that the model need not reflect all the fine details of the turbulent structure but rather account in an approximate way for the production and destruction of the turbulence.

Marvin, J. G.↗

Digital Distributed Parameter Model for Analysis of Unsteady Flow in Liquid-Filled Lines

A distributed parameter model for computing unsteady flow conditions in liquid systems is presented. The model is adapted for use in conjunction with a digital computer. The analytical method employed is essentially a solution synthesized from the effects of incremental step pressure pulses. The pressure pulses are generated because of incremental flow-rate changes that originate in a fluid system from a variety of sources, including the mechanical motion of the system structure. The pressure pulses propagate throughout the system at sonic velocity and are partially transmitted and reflected at each discontinuity. The velocity change caused by each pressure pulse is obtained from the relation for the characteristic acoustic impedance. Pressure and velocity time histories at any point in the system are obtained by a timewise summation of the contributions of the incremental pressure pulses passing that point. The analysis is presented in a form general enough to be applied to a variety of liquid-filled fluid systems. To illustrate the application of the method to a specific system, the response of a long, straight propellant line to a sinusoidal inlet flow and pressure perturbation is computed. Both a constant-cross-section line and a tapered line are analyzed in the example.

UNSTEADY FLOW↗

A preliminary compressible second-order closure model for high speed flows

A preliminary version of a compressible second-order closure model that was developed in connection with the National Aero-Space Plane Project is presented. The model requires the solution of transport equations for the Favre-averaged Reynolds stress tensor and dissipation rate. Gradient transport hypotheses are used for the Reynolds heat flux, mass flux, and turbulent diffusion terms. Some brief remarks are made about the direction of future research to generalize the model.

Speziale, Charles G.↗

Comparison of turbulence models for powered-lift flow fields

The development of a practical turbulence model for the complex flows with strong streamline curvature associated with powered-lift aircraft was investigated. The approach was to modify existing models that were compatible with the widely used kappa-epsilon model already available in an existing 3D Navier-Stokes code. The state-of-the-art of current turbulence models was examined. Several models that provide curvature corrections to the kappa-epsilon model were numerically investigated. An algebraic Reynolds stress model, which provides additional physics, was also examined. Test cases including an impinging jet, an impinging jet in a crossflow, and two adjacent impinging jets with the resultant upwash fountain were calculated. The first two test cases could be modeled reasonsbly well with some of the models. The upwash fountain could not be modeled with sufficient accuracy, which is consistent with the findings of other investigators.

Roberts, Donald W.↗

Development of a Computational Fluid Dynamics Model for a High-Speed Centrifugal Compressor

Computational fluid dynamics (CFD) has become widely used in the design and analysis of turbomachinery components such as centrifugal compressors. However, CFD is only a limited representation of experimental cases and struggles to model complex flows or can lack small model details to increase computation speed. To make advancements in compressor technology, designers need tools like CFD that can help them predict flow behavior in new designs. The High Efficiency Centrifugal Compressor (HECC) was designed by United Technologies Research Center (UTRC) for NASA to investigate the difficulties behind improving centrifugal compressor technology and to provide an open case to the turbomachinery research community. CFD analysis was completed by UTRC to aid in the design process. The pre-test study significantly overpredicted the compressor’s performance, and while post-test analysis was more representative, its results have yet to be publicly available. CFD technology has also since improved in recent years. In this study, a high-fidelity computational model for the vaneless configuration of HECC operating at its design point was developed and compared to an equivalent experimental case. This model incorporated GPU versions of mesh generation and solver codes from AeroDynamic Solutions, Inc. The development of such a model is significant since there has yet to be a published numerical simulation of HECC’s vaneless configuration. Additionally, the use of GPU codes decreases computation times for the model, which has allowed for the inclusion of impeller blade fillets. Initial results showed a general overprediction of compressor performance by the model. Predictions were more accurate for the impeller compared to predictions for the full compressor stage. Spanwise analyses showed that the model tended to follow flow behavior patterns exhibited in the experiment. Beyond this study, this model will be iterated upon and experimentally validated to further examine HECC performance.

centrifugal compressor↗