4,030 research outputs found
A non-hybrid method for the PDF equations of turbulent flows on unstructured grids
In probability density function (PDF) methods of turbulent flows, the joint
PDF of several flow variables is computed by numerically integrating a system
of stochastic differential equations for Lagrangian particles. A set of
parallel algorithms is proposed to provide an efficient solution of the PDF
transport equation, modeling the joint PDF of turbulent velocity, frequency and
concentration of a passive scalar in geometrically complex configurations. An
unstructured Eulerian grid is employed to extract Eulerian statistics, to solve
for quantities represented at fixed locations of the domain (e.g. the mean
pressure) and to track particles. All three aspects regarding the grid make use
of the finite element method (FEM) employing the simplest linear FEM shape
functions. To model the small-scale mixing of the transported scalar, the
interaction by exchange with the conditional mean model is adopted. An adaptive
algorithm that computes the velocity-conditioned scalar mean is proposed that
homogenizes the statistical error over the sample space with no assumption on
the shape of the underlying velocity PDF. Compared to other hybrid
particle-in-cell approaches for the PDF equations, the current methodology is
consistent without the need for consistency conditions. The algorithm is tested
by computing the dispersion of passive scalars released from concentrated
sources in two different turbulent flows: the fully developed turbulent channel
flow and a street canyon (or cavity) flow. Algorithmic details on estimating
conditional and unconditional statistics, particle tracking and particle-number
control are presented in detail. Relevant aspects of performance and
parallelism on cache-based shared memory machines are discussed.Comment: Accepted in Journal of Computational Physics, Feb. 20, 200
A numerical stabilization framework for viscoelastic fluid flow using the finite volume method on general unstructured meshes
A robust finite volume method for viscoelastic flow analysis on general
unstructured meshes is developed. It is built upon a general-purpose
stabilization framework for high Weissenberg number flows. The numerical
framework provides full combinatorial flexibility between different kinds of
rheological models on the one hand, and effective stabilization methods on the
other hand. A special emphasis is put on the velocity-stress-coupling on
co-located computational grids. Using special face interpolation techniques, a
semi-implicit stress interpolation correction is proposed to correct the
cell-face interpolation of the stress in the divergence operator of the
momentum balance. Investigating the entry-flow problem of the 4:1 contraction
benchmark, we demonstrate that the numerical methods are robust over a wide
range of Weissenberg numbers and significantly alleviate the high Weissenberg
number problem. The accuracy of the results is evaluated in a detailed mesh
convergence study
Derivation and evaluation of an approximate analysis for three-dimensional viscous subsonic flow with large secondary velocities
An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls
INITIAL MICROSEISMIC RECORDINGS AT THE ONSET OF UNCONVENTIONAL HYDROCARBON DEVELOPMENT IN THE ROME TROUGH, EASTERN KENTUCKY
The Cambrian Rogersville Shale is a part of a hydrocarbon system in the Rome Trough of eastern Kentucky and West Virginia that can only be produced unconventionally. In Kentucky, the Rogersville Shale ranges in depth from ~1,800 to ~3,700 m below the surface with the crystalline basement ~1,000 m lower than the formation’s base. Baseline Rome Trough microseismicity data were collected, focusing on wastewater injection wells and recently completed and planned unconventional hydrocarbon test wells in the Rogersville Shale, using thirteen broadband seismic stations installed between June, 2015 and June, 2016 and existing University of Kentucky and central and eastern United States network stations. In addition, the network’s minimum detection threshold, the magnitude at which the theoretical signal exceeds the noise by a factor of 3 between 1 and 20 Hz for at least 4 stations, was estimated for the project area. Thirty-eight local and regional events were located and magnitudes were calculated for each event. No events were proximal to operating disposal or hydrocarbon test wells, nor did any occur in the eastern Kentucky’s Rome Trough. The minimum detection threshold varies between 0.4 and 0.7 Mw from 0000-1100 UTC and 0.6 to 0.9 Mw from 1100-2300 UTC
A General Return-Mapping Framework for Fractional Visco-Elasto-Plasticity
We develop a fractional return-mapping framework for power-law
visco-elasto-plasticity. In our approach, the fractional viscoelasticity is
accounted through canonical combinations of Scott-Blair elements to construct a
series of well-known fractional linear viscoelastic models, such as
Kelvin-Voigt, Maxwell, Kelvin-Zener and Poynting-Thomson. We also consider a
fractional quasi-linear version of Fung's model to account for stress/strain
nonlinearity. The fractional viscoelastic models are combined with a fractional
visco-plastic device, coupled with fractional viscoelastic models involving
serial combinations of Scott-Blair elements. We then develop a general
return-mapping procedure, which is fully implicit for linear viscoelastic
models, and semi-implicit for the quasi-linear case. We find that, in the
correction phase, the discrete stress projection and plastic slip have the same
form for all the considered models, although with different property and
time-step dependent projection terms. A series of numerical experiments is
carried out with analytical and reference solutions to demonstrate the
convergence and computational cost of the proposed framework, which is shown to
be at least first-order accurate for general loading conditions. Our numerical
results demonstrate that the developed framework is more flexible, preserves
the numerical accuracy of existing approaches while being more computationally
tractable in the visco-plastic range due to a reduction of in CPU time.
Our formulation is especially suited for emerging applications of fractional
calculus in bio-tissues that present the hallmark of multiple viscoelastic
power-laws coupled with visco-plasticity
Research reports: 1991 NASA/ASEE Summer Faculty Fellowship Program
The basic objectives of the programs, which are in the 28th year of operation nationally, are: (1) to further the professional knowledge of qualified engineering and science faculty members; (2) to stimulate an exchange of ideas between participants and NASA; (3) to enrich and refresh the research and teaching activities of the participants' institutions; and (4) to contribute to the research objectives of the NASA Centers. The faculty fellows spent 10 weeks at MSFC engaged in a research project compatible with their interests and background and worked in collaboration with a NASA/MSFC colleague. This is a compilation of their research reports for summer 1991
- …