Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Anisotropic flow”

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 109 records · Page 6

Turbulence and surface heat transfer near the stagnation point of a circular cylinder in turbulent flow

A turbulent boundary layer flow analysis of the momentum and thermal flow fields near the forward stagnation point due to a circular cylinder in turbulent cross flow is presented. Turbulence modeling length scale, anisotropic turbulence initial profiles and boundary conditions were identified as functions of the cross flow turbulence intensity and the boundary layer flow far field velocity. These parameters were used in a numerical computational procedure to calculate the mean velocity, mean temperature, and turbulence double correlation profiles within the flow field. The effects of the cross flow turbulence on the stagnation region momentum and thermal flow fields were investigated. This analysis predicted the existing measurements of the stagnation region mean velocity and surface heat transfer rate with cross flow Reynolds number and turbulence intensity less than 250,000 and 0.05, respectively.

Wang, C. R.↗

Coupled Thermo-Hydrological-Mechanical-Chemical Behavior of Anisotropic Granite for Geologic Disposal of High-Level Radioactive Waste: A Core-Scale Laboratory Investigation

The coupled thermo-hydrological-mechanical-chemical (THMC) behavior of rock within an Excavation Damaged Zone (EDZ) is critical for the safety and long-term performance of a geological repository for high-level radioactive wastes. While many laboratory experiments have been conducted to investigate EDZ rocks, the flow and deformation characteristics resulting from anisotropic rock textures and microcrack distribution under triaxial loading and elevated temperatures remain poorly understood. Particularly, cracks at various scales serve as fast paths for fluid flow and solute transport and present as focal points of mechanical weakness, which complicate the coupled THMC processes in anisotropic EDZ rocks and challenge modeling predictions. Here, in this study, a series of core-scale experiments was conducted on three granite samples under repository-relevant conditions. These rock samples were obtained from the Grimsel Underground Research Laboratory (URL), featured by anisotropic minerals (represented by bedding layers) and microcrack distributions and coarse cm-scale grain sizes. During the experiments, samples were subjected to an elevated temperature at 90 °C and different triaxial loading conditions either by radial (normal to bedding layers) or axial (parallel to bedding layers) compaction. Water was injected into the samples, and the rock permeability evolutions and effluent water chemistry were monitored closely. For intact samples, thermal expansion of minerals at 90 °C resulted in a large, 75% irreversible permeability reduction and rock strengthening under radial compaction, while thermal impact was limited to a 15% permeability reduction under axial compaction. In contrast, for a sample containing open cracks, the growth of fractures during the experiment resulted in an abrupt permeability increase and fast failure at 90 °C. The effluent water chemistry indicates much more considerable mineral dissolution from large shear sliding than that in rocks dominated by mechanical compaction. These results helped better understand the coupled THMC processes in anisotropic rocks containing cracks, evaluate the behaviors of EDZ rocks, and predict the long-term evolution of EDZ for the performance of the repositories.

Coupled THMC processes↗

Verification of Viscous Goal-Based Anisotropic Mesh Adaptation

Adaptive unstructured mesh techniques have a limited, but growing impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic mesh adaptation mechanics. Anisotropic metric construction methods are evaluated with analytically defined primal and adjoint fields. This allows the comparison of different metric formulations and different implementations of the same formulation without the complications of a flow and adjoint solution method. Unstructured mesh adaptation tools are verified by comparison on analytic primal and dual field before verification on benchmark aerodynamics cases. The documentation of these verification exercises helps to prepare these goal-based methods for routine use in production simulation workflows.

Mesh adaptation↗

Numerical computation of aerodynamics and heat transfer in a turbine cascade and a turn-around duct using advanced turbulence models

The objective of this research is to develop turbulence models to predict the flow and heat transfer fields dominated by the curvature effect such as those encountered in turbine cascades and turn-around ducts. A Navier-Stokes code has been developed using an explicit Runge-Kutta method with a two layer k-epsilon/ARSM (Algebraic Reynolds Stress Model), Chien's Low Reynolds Number (LRN) k-epsilon model and Coakley's LRN q-omega model. The near wall pressure strain correlation term was included in the ARSM. The formulation is applied to Favre-averaged N-S equations and no thin-layer approximations are made in either the mean flow or turbulence transport equations. Anisotropic scaling of artificial dissipation terms was used. Locally variable timestep was also used to improve convergence. Detailed comparisons were made between computations and data measured in a turbine cascade by Arts et al. at Von Karman Institute. The surface pressure distributions and wake profiles were predicted well by all the models. The blade heat transfer is predicted well by k-epsilon/ARSM model, as well as the k-epsilon model. It's found that the onset of boundary layer transition on both surfaces is highly dependent upon the level of local freestream turbulence intensity, which is strongly influenced by the streamline curvature. Detailed computation of the flow in the turn around duct has been carried out and validated against the data by Monson as well as Sandborn. The computed results at various streamwise locations both on the concave and convex sides are compared with flow and turbulence data including the separation zone on the inner well. The k-epsilon/ARSM model yielded relatively better results than the two-equation turbulence models. A detailed assessment of the turbulence models has been made with regard to their applicability to curved flows.

Lakshminarayana, B.↗

Anisotropy of MHD Turbulence at Low Magnetic Reynolds Number

Turbulent fluctuations in MHD flows are known to become dimensionally anisotropic under the action of a sufficiently strong magnetic field. We consider the technologically relevant case of low magnetic Reynolds number and apply the method of DNS of forced flow in a periodic box to generate velocity fields. The analysis based on different anisotropy characteristics shows that the dimensional anisotropy is virtually scale-independent. We also find that, except for the case of very strong magnetic field, the flow is componentally isotropic. Its kinetic energy is practically uniformly distributed among the velocity components.

Zikanov, O.↗

Tracking Dendritic Growth in Hydrogen-Based Hematite Reduction via Computer Vision

The reduction of hematite to metallic iron using hydrogen (H2) as a reducing agent presents a promising pathway for decarbonizing steel production. In this study, we employ a combination of in situ confocal scanning laser microscopy (CSLM) and advanced computer vision techniques to quantitatively analyze dendritic growth of ferrite during H2-based reduction of iron oxide at high temperatures. A workflow integrating Watershed Image Segmentation (WIS) and Lucas-Kanade Optical Flow (LKOF) is developed to extract both global and local kinetic information from time-resolved micrograph sequences. H2 reduction experiments conducted at 1400 degrees C and 1500 degrees C demonstrate a clear correlation between temperature and reduction rate, as evidenced by accuracy of fitted Johnson-Mehl-Avrami-Kolmogorov (JMAK) parameters. Optical flow analysis further elucidates the anisotropic and branched nature of dendritic growth, providing spatially resolved velocity fields that correlate well with global transformation kinetics. The proposed methodology demonstrates strong agreement with experimental measurements and literature values, offering a robust framework for automated image-based analysis to study kinetics through microstructural evolution in the reduction of iron ore, and likely other reaction-diffusion phenomena.

08 HYDROGEN↗

The structure of sheared turbulence near a plane boundary

An analysis is presented of how a plane boundary affects the structure of turbulence in a sheared free stream. A uniform-shear boundary layer (USBL) is formulated with slip velocity condition at the surface, and inhomogeneous rapid distortion theory is applied. The effects of blocking by the surface on the turbulence structure in USBL is compared with those in the shear-free boundary layer (SFBL). Shear produces highly anisotropic eddies elongated in the flow direction. The vertical velocity variance is reduced with shear at all heights, roughly in proportion to the reduction in the homogeneous value, but the shape of the profile remains unchanged only near the surface. The streamwise integral scales increase with shear, indicating elongation of the streamwise extent of eddies.

Lee, Moon J.↗

Large-scale solar velocity fields

Daily observations of Doppler line shifts made with very low spatial resolution with the Stanford magnetograph were used to study the equatorial rotation rate, limb effect on the disk, and the mean meridional circulation. The equatorial rotation rate compared favorably with Howard's results (1977), and the RMS deviation of the daily measurements about the mean value was 1% of the rate, much smaller than the fluctuations reported by Howard and Harvey (1970). The limb effect on the disk in equatorial scans, a redshift at the center of the disk, and central meridian scans used to search for the component of mean meridional circulation symmetric across the equator were discussed, and it was shown that models of the differential rotation driven by an axisymmetric meridional circulation and an anisotropic eddy viscosity predict an equatorwards flow at the surface. It was concluded that the observation of a poleward-directed meridional circulation lends support to the giant cell convection models of the differential rotation over the anisotropic eddy viscosity models.

Duvall, T. L., Jr.↗

Predicting aerodynamic sound utilizing a two-point, two-time turbulence theory

The feasibility of determining two-point, two-time turbulent velocity correlations based on extending closure models for one-point, one-time turbulence correlations is demonstrated. The procedure used is based on a spatial moment integral formulation of the governing equations using approximate, parameterized trial functions for the two-point, two-time velocity correlations. Solution of the equations is shown to give a set of anisotropic length scales and the separation-time-dependent decorrelation of the ensemble averaged turbulent velocities. The analysis is simplified using the assumption of homogeneous stationary turbulence and a constant shear, unidirectional mean flow. It is shown that the anisotropic behavior of measured turbulence correlations can be characterized through this technique.

Hecht, A. M.↗

Birkeland currents in an anisotropic, magnetostatic plasma

The paper derives an expression for the parallel current density for a plasma characterized by negligible bulk flow (magnetostatic velocity and a two-component (anisotropic) pressure tensor by expanding the equilibrium Vlasov equation for each species in the adiabatic parameter until such point as a nonvanishing moment j-parallel = integral d-cubed vv-parallel f is identified. The result is a nonlocal one: it relates j-parallel at one point s along a field line to j-parallel at another (reference) point s0 plus an integral function of the pressure and magnetic field between them. The equation derived by Vasyliunas (1970) is generalized and extended to a plasma in which the pressure tensor is comprised of two separate elements, P-perpendicular and P-parallel. This equation follows when P-perpendicular is set equal to P-parallel and s and s0 are taken to be at the ionosphere and the equator. These results are compared to others in the literature.

Birmingham, Thomas J.↗

Jet Drift in Heavy Ion Collisions

We introduce a sub-eikonal anisotropic contribution to jetbroadening, “jet drift”, that couples to the flow of the medium, showing that this effect results in a deflection of hard partons, and thus jets, in the direction of the medium flow. Next, we study this effect in a full-fledged hybrid transport simulation of √s = 5.02 TeV PbPb collisions at the LHC, tracking trajectories of hard partons with perturbative energy loss and drift. We show that sub-eikonal anisotropic effects, including flow-mediated jet drift, are sensitive to properties of the medium that traditional eikonal isotropic effects are insensitive to, demonstrating that including these effects leads to modifications to jet and hard particle observables that survive averaging over events. We show that jet drift leads to an enhancement of the elliptic flow (v2) of hard particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Martian global groundwater model

A global groundwater flow model was constructed for Mars to study hydrologic response under a variety of scenarios, improving and extending earlier simple cross sectional models. The model is capable of treating both steady state and transient flow as well as permeability that is anisotropic in the horizontal dimensions. A single near surface confining layer may be included (representing in these simulations a coherent permafrost layer). Furthermore, in unconfined flow, locations of complete saturation and seepage are determined. The flow model assumes that groundwater gradients are sufficiently low that DuPuit conditions are satisfied and the flow component perpendicular to the ground surface is negligible. The flow equations were solved using a finite difference method employing 10 deg spacing of latitude and longitude.

Howard, Alan D.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. Laminar flow over a delta wing and turbulent flow over an ONERA M6 wing are verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and a moment. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. The wing adaptation cases in the current study document a clear improvement to existing grid adaptation procedures. The stage is set for the infusion of verified grid adaptation into production fluid flow simulations.

Park, Michael A.↗

Verification of Unstructured Grid Adaptation Components

Adaptive unstructured grid techniques have made limited impact on production analysis workflows where the control of discretization error is critical to obtaining reliable simulation results. Recent progress has matured a number of independent implementations of flow solvers, error estimation methods, and anisotropic grid adaptation mechanics. Known differences and previously unknown differences in grid adaptation components and their integrated processes are identified here for study. Unstructured grid adaptation tools are verified using analytic functions and the Code Comparison Principle. Three analytic functions with different smoothness properties are adapted to show the impact of smoothness on implementation differences. A scalar advection-diffusion problem with an analytic solution that models a boundary layer is adapted to test individual grid adaptation components. The scalar problems illustrate known differences in a grid adaptation component implementation and a previously unknown interaction between components. Laminar flow over a delta wing is verified with multiple, independent grid adaptation procedures to show consistent convergence to fine-grid forces and pitching moment.

Park, Michael A.↗

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra↗

A close-up view of the S87 molecular outflow

Observations of the S87 star-forming region have been made at optical, far-infrared, centimeter, and millimeter wavelengths in order to explain the origin of the previously discovered supersonic molecular gas in this source. S87/IRS 1 is a massive, pre-main-sequence object still embedded in its parent molecular cloud, but disrupting its surroundings through the action of its powerful stellar wind. The shocked wind gas provides the force required to accelerate the surrounding molecular gas to supersonic velocities. An ordered, large-scale magnetic field is postulated to set up an initially anisotropic pressure distribution which channels the flow into two oppositely directed lobes and provides some further acceleration to the already supersonic molecular gas. Only 5 percent of the molecular gas actually reaches escape velocity.

Barsony, M.↗

Positivity-preserving numerical schemes for multidimensional advection

This report describes the construction of an explicit, single time-step, conservative, finite-volume method for multidimensional advective flow, based on a uniformly third-order polynomial interpolation algorithm (UTOPIA). Particular attention is paid to the problem of flow-to-grid angle-dependent, anisotropic distortion typical of one-dimensional schemes used component-wise. The third-order multidimensional scheme automatically includes certain cross-difference terms that guarantee good isotropy (and stability). However, above first-order, polynomial-based advection schemes do not preserve positivity (the multidimensional analogue of monotonicity). For this reason, a multidimensional generalization of the first author's universal flux-limiter is sought. This is a very challenging problem. A simple flux-limiter can be found; but this introduces strong anisotropic distortion. A more sophisticated technique, limiting part of the flux and then restoring the isotropy-maintaining cross-terms afterwards, gives more satisfactory results. Test cases are confined to two dimensions; three-dimensional extensions are briefly discussed.

Leonard, B. P.↗