Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical mixing”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Transient modeling/analysis of hyperbolic heat conduction problems employing mixed implicit-explicit alpha method

This paper describes the evaluation of mixed implicit-explicit finite element formulations for hyperbolic heat conduction problems involving non-Fourier effects. In particular, mixed implicit-explicit formulations employing the alpha method proposed by Hughes et al. (1987, 1990) are described for the numerical simulation of hyperbolic heat conduction models, which involves time-dependent relaxation effects. Existing analytical approaches for modeling/analysis of such models involve complex mathematical formulations for obtaining closed-form solutions, while in certain numerical formulations the difficulties include severe oscillatory solution behavior (which often disguises the true response) in the vicinity of the thermal disturbances, which propagate with finite velocities. In view of these factors, the alpha method is evaluated to assess the control of the amount of numerical dissipation for predicting the transient propagating thermal disturbances. Numerical test models are presented, and pertinent conclusions are drawn for the mixed-time integration simulation of hyperbolic heat conduction models involving non-Fourier effects.

Tamma, Kumar K.↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
DOE OSTI2021

3d N=4 Bootstrap and Mirror Symmetry

We investigate the non-BPS realm of 3d {N} = 4 N = 4 superconformal field theory by uniting the non-perturbative methods of the conformal bootstrap and supersymmetric localization, and utilizing special features of 3d {N} = 4 N = 4 theories such as mirror symmetry and a protected sector described by topological quantum mechanics (TQM). Supersymmetric localization allows for the exact determination of the conformal and flavor central charges, and the latter can be fed into the mini-bootstrap of the TQM to solve for a subset of the OPE data. We examine the implications of the Z_2 Z 2 mirror action for the SCFT single- and mixed-branch crossing equations for the moment map operators, and apply numerical bootstrap to obtain universal constraints on OPE data for given flavor symmetry groups. A key ingredient in applying the bootstrap analysis is the determination of the mixed-branch superconformal blocks. Among other results, we show that the simplest known self-mirror theory with SU(2) \times SU(2) S U ( 2 ) × S U ( 2 ) flavor symmetry saturates our bootstrap bounds, which allows us to extract the non-BPS data and examine the self-mirror Z_2 Z 2 symmetry thereof.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Simulating combustion of a seven-component surrogate for a gasoline/ethanol blend including soot formation and comparison with experiments

Combustion of a seven-component surrogate for a research grade 87 octane gasoline mixed with 10% ethanol is investigated experimentally and numerically from the perspective of an isolated droplet burning under conditions that promote one-dimensional gas transport. The numerical analysis included a kinetic mechanism comprised of 398 species and 24,814 reactions and a soot model that accounted for nucleation, surface growth, coalescence/aggregation of soot particles, and luminous flame radiation. Measurements of droplet and flame diameters were made for an initial droplet diameter (D o ) of approximately 0.63 mm. The simulations agreed well with the measurements including the location of the soot shell. Preferential vaporization was revealed by simulations of the liquid concentrations in the droplet. Predicted peak soot volume fractions coincided with temperatures between 1300 K and 1400 K as a soot inception temperature. Simulations were also carried out for D o between 0.25 mm and 5 mm to explore the effect of radiation and D o on burning. Below 0.25 mm radiation was negligible and burning rates and flame temperatures converged to a single value. Increasing D o up to 1.8 mm lowered the burning rate with luminous radiation having a strong effect. When radiation was entirely removed from the model the burning rate was nearly constant. Above D o = 2 mm droplets extinguished almost immediately after ignition. The flame temperature decreased with increasing D o while it increased when radiation was omitted. As a result, the simulations show that soot precursors including polyaromatic hydrocarbons were concentrated around the soot shell.

42 ENGINEERING↗
NASA NTRS1982

Features of a reattaching turbulent shear layer subject to an adverse pressure gradient

Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step with superimposed adverse pressure gradient. Mean velocities, Reynolds stresses and triple-products measured by a laser Doppler velocimeter are presented for two cases of adverse pressure gradient. Mixing lengths, eddy viscosities, production, convection, turbulent diffusion, and dissipation terms are extracted from the data. These data are compared with various mixing length and eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and the algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model

Driver, D. M.↗
NASA NTRS1985

On the effect of quarter-point element size on fracture criteria

The effect of the quarter-point element size on the solution parameters of fracture criteria is investigated in view of its effect on the computed stress intensity factors. An analytical relationship between the error in calculating crack propagation increment and stress intensity for single mode problems is obtained while, for mixed mode problems, the investigation is based on numerical experimentation. The singular element size is found to have a substantial effect on both crack propagation angle and increment for the geometries considered.

Yehia, N. A. B.↗
NASA NTRS1986

The interaction of a surface cold front with a prefrontal thermodynamically well-mixed boundary layer

The dynamics of a surface cold front as it proceeds from sea to shore and interacts with a thermodynamically well-mixed boundary layer are investigated by means of numerical simulations using a modified version of the two-dimensional model of Reeder and Smith (1987). The theoretical basis of the model is explained, and the results are presented graphically and characterized in detail. Good agreement is found with observational data on Australian summer cold fronts, and the transverse frontal circulation is shown to be significantly affected by diabatic heating.

Reeder, Michael J.↗
NASA NTRS1990

Numerical simulation of solidification in a horizontal cylindrical annulus charged with an aqueous salt solution

A computational study of solidification of a binary Na2CO3 solution in a horizontal cylindrical annulus is performed using a continuum formulation with a control-volume based, finite-difference scheme. The initial conditions were selected to facilitate the study of counter thermal and solutal convection, accompanied by extensive mushy region growth. Numerical results are compared with experimental data with mixed success. Qualitative agreement is obtained for the overall solidification process and associated physical phenomena. However, the plume thickness calculated for the solutally-driven convective upflow is substantially smaller than the observed value. Evolution of double-diffusive layers is predicted, but over a time scale much smaller than that observed experimentally. Good agreement is obtained between predicted and measured results for solid growth, but the mushy region thickness is significantly overpredicted.

Neilson, D. G.↗
NASA NTRS1990

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

This research is involved with the implementation of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program to extend the present capabilities of this method was initiated for the treatment of chemically reacting flows. In the DNS efforts, the focus is on detailed investigations of the effects of compressibility, heat release, and non-equilibrium kinetics modelings in high speed reacting flows. Emphasis was on the simulations of simple flows, namely homogeneous compressible flows, and temporally developing high speed mixing layers.

Givi, Peyman↗
NASA NTRS1991

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

This research is involved with the implementations of advanced computational schemes based on large eddy simulations (LES) and direct numerical simulations (DNS) to study the phenomenon of mixing and its coupling with chemical reactions in compressible turbulent flows. In the efforts related to LES, a research program was initiated to extend the present capabilities of this method for the treatment of chemically reacting flows, whereas in the DNS efforts, focus was on detailed investigations of the effects of compressibility, heat release, and nonequilibrium kinetics modeling in high speed reacting flows. The efforts to date were primarily focussed on simulations of simple flows, namely, homogeneous compressible flows and temporally developing hign speed mixing layers. A summary of the accomplishments is provided.

Givi, Peyman↗
NASA NTRS1993

Time accurate simulations of compressible shear flows

The objectives of this research are to employ direct numerical simulation (DNS) to study the phenomenon of mixing (or lack thereof) in compressible free shear flows and to suggest new means of enhancing mixing in such flows. The shear flow configurations under investigation are those of parallel mixing layers and planar jets under both non-reacting and reacting nonpremixed conditions. During the three-years of this research program, several important issues regarding mixing and chemical reactions in compressible shear flows were investigated.

Givi, Peyman↗
NASA NTRS2007

Time-Accurate Numerical Simulations of Synthetic Jet Quiescent Air

The unsteady evolution of three-dimensional synthetic jet into quiescent air is studied by time-accurate numerical simulations using a second-order accurate mixed explicit-implicit fractional step scheme on Cartesian grids. Both two-dimensional and three-dimensional calculations of synthetic jet are carried out at a Reynolds number (based on average velocity during the discharge phase of the cycle V(sub j), and jet width d) of 750 and Stokes number of 17.02. The results obtained are assessed against PIV and hotwire measurements provided for the NASA LaRC workshop on CFD validation of synthetic jets.

Rupesh, K-A. B.↗
NASA NTRS2014

The Prediction and Analysis of Jet Flows and Scattered Turbulent Mixing Noise about Flight Vehicle Airframes

Jet flows interacting with nearby surfaces exhibit a complex behavior in which acoustic and aerodynamic characteristics are altered. The physical understanding and prediction of these characteristics are essential to designing future low noise aircraft. A new approach is created for predicting scattered jet mixing noise that utilizes an acoustic analogy and steady Reynolds-averaged Navier-Stokes solutions. A tailored Green's function accounts for the propagation of mixing noise about the airframe and is calculated numerically using a newly developed ray tracing method. The steady aerodynamic statistics, associated unsteady sound source, and acoustic intensity are examined as jet conditions are varied about a large flat plate. A non-dimensional number is proposed to estimate the effect of the aerodynamic noise source relative to jet operating condition and airframe position.The steady Reynolds-averaged Navier-Stokes solutions, acoustic analogy, tailored Green's function, non-dimensional number, and predicted noise are validated with a wide variety of measurements. The combination of the developed theory, ray tracing method, and careful implementation in a stand-alone computer program result in an approach that is more first principles oriented than alternatives, computationally efficient, and captures the relevant physics of fluid-structure interaction.

Miller, Steven A. E.↗
NASA NTRS2018

NASA Observations and Modeling During ICE-POP

Recap: NASA-Specific Objectives for ICE-POP: Provide real-time observational and NWP data in support of ICE-POP, participate in significant international science effort; GPM (Global Precipitation Measurement) Ground Validation and NASA Weather Program -Direct/physical validation of active/passive satellite-based snowfall retrieval algorithms over coastline and mountains; melting layer interaction with terrain -Physics of snow, coupling to snow water equivalent rate and satellite remote sensor retrieval algorithm assumptions - -Size distributions, types/habit, water equivalent, profiles -NU-WRF (NASA-Unified Weather Research and Forecasting) Model plus Observational analyses: Movement toward “level IV products” leverage intensive and multi-faceted NWP (Numerical Weather Prediction) component -Model precipitation processes (liquid, mixed phase and frozen); Build model testing database for further active/passive remote sensing algorithm development (e.g., satellite data simulators) -"Integrated" validation of products in operational context.

Precipitation Science↗