Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADVECTION”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Limitation of parallel flow in double diffusive convection: Two- and three-dimensional transitions in a horizontal porous domain

Two-dimensional (2D) and three-dimensional (3D) numerical simulations of double diffusion natural convection in an elongated enclosure filled with a binary fluid saturating a porous medium are carried out in the present work. The Boussinesq approximation is made in the formulation of the problem, and Neumann boundary conditions for temperature and concentration are adopted, respectively, on vertical and horizontal walls of the cavity. The used numerical method is based on the control volume approach, with the third order quadratic upstream interpolation scheme in approximating the advection terms. A semi implicit method algorithm is used to handle the velocity-pressure coupling. To avoid the excessively high computer time inherent to the solution of 3D natural convection problems, full approximation storage with full multigrid method is used to solve the problem. A wide range of the controlling parameters (Rayleigh-Darcy number Ra, lateral aspect ratio Ay, Lewis number Le, and the buoyancy ration N) is investigated. We clearly show that increasing the depth of the cavity (i.e., the lateral aspect ratio) has an important effect on the flow patterns. The 2D perfect parallel flows obtained for small lateral aspect ratio are drastically destabilized by increasing the cavity lateral dimension. This yields a 3D fluid motion with a much more complex flow pattern and the usually considered 2D parallel flow model cannot be applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

On the nonlinear stability of a quasi-two-dimensional drift kinetic model for ion temperature gradient turbulence

We study a quasi-two-dimensional electrostatic drift kinetic system as a model for near-marginal ion temperature gradient driven turbulence. A proof is given for the nonlinear stability of this system under conditions of linear stability. This proof is achieved using a transformation that diagonalizes the linear dynamics and also commutes with nonlinear E × B advection. For the case when linear instability is present, a corollary is found that forbids nonlinear energy transfer between appropriately defined sets of stable and unstable modes. It is speculated that this may explain the preservation of linear eigenmodes in nonlinear gyrokinetic simulations. Based on this property, a dimensionally reduced (∞×∞→1) system is derived that may be useful for understanding dynamics around the critical gradient of Dimits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Transformed Fourier and Fick equations for the control of heat and mass diffusion

We review recent advances in the control of diffusion processes in thermodynamics and life sciences through geometric transforms in the Fourier and Fick equations, which govern heat and mass diffusion, respectively. We propose to further encompass transport properties in the transformed equations, whereby the temperature is governed by a three-dimensional, time-dependent, anisotropic heterogeneous convection-diffusion equation, which is a parabolic partial differential equation combining the diffusion equation and the advection equation. We perform two dimensional finite element computations for cloaks, concentrators and rotators of a complex shape in the transient regime. We precise that in contrast to invisibility cloaks for waves, the temperature (or mass concentration) inside a diffusion cloak crucially depends upon time, its distance from the source, and the diffusivity of the invisibility region. However, heat (or mass) diffusion outside cloaks, concentrators and rotators is unaffected by their presence, whatever their shape or position. Finally, we propose simplified designs of layered cylindrical and spherical diffusion cloaks that might foster experimental efforts in thermal and biochemical metamaterials.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Direct experimental visualization of the global Hamiltonian progression of two-dimensional Lagrangian flow topologies from integrable to chaotic state

Countless theoretical/numerical studies on transport and mixing in two-dimensional (2D) unsteady flows lean on the assumption that Hamiltonian mechanisms govern the Lagrangian dynamics of passive tracers. However, experimental studies specifically investigating said mechanisms are rare. Moreover, they typically concern local behavior in specific states (usually far away from the integrable state) and generally expose this indirectly by dye visualization. Laboratory experiments explicitly addressing the global Hamiltonian progression of the Lagrangian flow topology entirely from integrable to chaotic state, i.e., the fundamental route to efficient transport by chaotic advection, appear non-existent. This motivates our study on experimental visualization of this progression by direct measurement of Poincaré sections of passive tracer particles in a representative 2D time-periodic flow. This admits (i) accurate replication of the experimental initial conditions, facilitating true one-to-one comparison of simulated and measured behavior, and (ii) direct experimental investigation of the ensuing Lagrangian dynamics. The analysis reveals a close agreement between computations and observations and thus experimentally validates the full global Hamiltonian progression at a great level of detail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The appearance, motion, and disappearance of three-dimensional magnetic null points

While theoretical models and simulations of magnetic reconnection often assume symmetry such that the magnetic null point when present is co-located with a flow stagnation point, the introduction of asymmetry typically leads to non-ideal flows across the null point. To understand this behavior, we present exact expressions for the motion of three-dimensional linear null points. The most general expression shows that linear null points move in the direction along which the magnetic field and its time derivative are antiparallel. Null point motion in resistive magnetohydrodynamics results from advection by the bulk plasma flow and resistive diffusion of the magnetic field, which allows non-ideal flows across topological boundaries. Null point motion is described intrinsically by parameters evaluated locally; however, global dynamics help set the local conditions at the null point. During a bifurcation of a degenerate null point into a null-null pair or the reverse, the instantaneous velocity of separation or convergence of the null-null pair will typically be infinite along the null space of the Jacobian matrix of the magnetic field, but with finite components in the directions orthogonal to the null space. Not all bifurcating null-null pairs are connected by a separator. Furthermore, except under special circumstances, there will not exist a straight line separator connecting a bifurcating null-null pair. The motion of separators cannot be described using solely local parameters because the identification of a particular field line as a separator may change as a result of non-ideal behavior elsewhere along the field line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Particle acceleration and reconnection in the solar wind

An emerging paradigm for the dissipation of magnetic turbulence in the supersonic solar wind is via localized quasi-2D small-scale magnetic island reconnection processes. An advection-diffusion transport equation for a nearly isotropic particle distribution describes particle transport and energization in a region of interacting magnetic islands [1; 2]. The dominant charged particle energization processes are 1) the electric field induced by quasi-2D magnetic island merging, and 2) magnetic island contraction. The acceleration of charged particles in a “sea of magnetic islands” in a super-Alfvénic flow, and the energization of particles by combined diffusive shock acceleration (DSA) and downstream magnetic island reconnection processes are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Verification of BOUT++ by the method of manufactured solutions

BOUT++ is a software package designed for solving plasma fluid models. It has been used to simulate a wide range of plasma phenomena ranging from linear stability analysis to 3D plasma turbulence and is capable of simulating a wide range of drift-reduced plasma fluid and gyro-fluid models. A verification exercise has been performed as part of a EUROfusion Enabling Research project, to rigorously test the correctness of the algorithms implemented in BOUT++, by testing order-of-accuracy convergence rates using the Method of Manufactured Solutions (MMS). We present tests of individual components including time-integration and advection schemes, non-orthogonal toroidal field-aligned coordinate systems and the shifted metric procedure which is used to handle highly sheared grids. The flux coordinate independent approach to differencing along magnetic field-lines has been implemented in BOUT++ and is here verified using the MMS in a sheared slab configuration. Finally, we show tests of three complete models: 2-field Hasegawa-Wakatani in 2D slab, 3-field reduced magnetohydrodynamics (MHD) in 3D field-aligned toroidal coordinates, and 5-field reduced MHD in slab geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A comparison of Eulerian and semi-Lagrangian approaches for modeling stream water quality

This paper describes and compares some of the advantages and limitations of Eulerian and Lagrangian approaches to water quality modeling and introduces a mixed Eulerian-Lagrangian (or semi-Lagrangian) methodology that captures the strengths of both approaches. The semi-Lagrangian modeling approach is applied to advection-dominated rivers, and flexibly ensures unconditional stability for all time step durations and grid segmentations. The semi-Lagrangian modeling approach is demonstrated by applying it to estimate the dissolved oxygen concentrations in the Sava River in Slovenia, focusing on aspects of the methodology and findings that would be of broad interest to managers of water quality in fluvial water bodies. Results of comparisons of the semi-Lagrangian model with the Eulerian-based QUAL2K model in steady and non-steady scenarios demonstrate that while both models are fully capable of producing satisfactory results when optimally configured, the semi-Lagrangian approach offers accuracy and stability without sensitivity to the interaction of time step size and computational grid segmentation scheme.

Sava River↗
DOE OSTI2013

THE TRANSITION LAYER BETWEEN THE THIN ACCRETION DISK AND THE HOT CORONA: CONSTRAINTS FROM THE THERMAL X-RAY LINE EMISSION

The observed continuum spectra in the low-luminosity active galactic nuclei, M81, can be reproduced by the advection-dominated accretion flow (ADAF)+thin disk/corona accretion flow model. We find that the best fit to the spectral energy distribution of M81 requires accretion rate m-dot =3 Multiplication-Sign 10{sup -4} for R {sub tr} = 60 R {sub S}. The observed H-like and He-like iron lines may probably be emitted from the transition layer between the vertically connected hot corona and cold thin disk. We model a power-law function of height z for the electron temperature in the transition layer, and calculate the thermal H-like and He-like iron line emission from M81. Combining the observations with our results, we find that the thickness of the transition layer is {approx}14%-23% of the thickness of the corona at the same radius. We discuss that the observed stronger H-like line emission than He-like line emission cannot be explained with the pure ADAF model. Our results provide useful constraints on the theoretical model for disk-corona connection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

EFFECTS OF RESISTIVITY ON MAGNETIZED CORE-COLLAPSE SUPERNOVAE

We studied the role of turbulent resistivity in the core-collapse of a strongly magnetized massive star, carrying out two-dimensional resistive-MHD simulations. Three cases with different initial strengths of magnetic field and rotation are investigated: (1) a strongly magnetized rotating core, (2) a moderately magnetized rotating core, and (3) a very strongly magnetized non-rotating core. In each case, one ideal-MHD model and two resistive-MHD models are computed. As a result of these computations, each model shows an eruption of matter assisted by magnetic acceleration (and also by centrifugal acceleration in the rotating cases). We found that resistivity attenuates the explosion in cases 1 and 2, while it enhances the explosion in case 3. We also found that in the rotating cases, the main mechanisms for the amplification of a magnetic field in the post-bounce phase are an outward advection of the magnetic field and a twisting of poloidal magnetic field lines by differential rotation, which are somewhat dampened down with the presence of resistivity. Although magnetorotational instability seems to occur in the rotating models, it plays only a minor role in magnetic field amplification. Another impact of resistivity is that on the aspect ratio. In the rotating cases, a large aspect ratio of the ejected matter, >2.5, attained in an ideal-MHD model is reduced to some extent in a resistive model. These results indicate that resistivity possibly plays an important role in the dynamics of strongly magnetized supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

PRE-SUPERNOVA EVOLUTION OF ROTATING SOLAR METALLICITY STARS IN THE MASS RANGE 13-120 M {sub Sun} AND THEIR EXPLOSIVE YIELDS

We present the first set of a new generation of models of massive stars with a solar composition extending between 13 and 120 M {sub Sun }, computed with and without the effects of rotation. We included two instabilities induced by rotation: the meridional circulation and the shear instability. We implemented two alternative schemes to treat the transport of the angular momentum: the advection-diffusion formalism and the simpler purely diffusive one. The full evolution from the pre-main sequence up to the pre-supernova stage is followed in detail with a very extended nuclear network. The explosive yields are provided for a variety of possible mass cuts and are available at the Web site http://www.iasf-roma.inaf.it/orfeo/public{sub h}tml. We find that both the He and the CO core masses are larger than those of their non-rotating counterparts. Also the C abundance left by the He burning is lower than in the non-rotating case, especially for stars with an initial mass of 13-25 M {sub Sun }, and this affects the final mass-radius relation, basically the final binding energy, at the pre-supernova stage. The elemental yields produced by a generation of stars rotating initially at 300 km s{sup -1} do not change substantially with respect to those produced by a generation of non-rotating massive stars, the main differences being a slight overproduction of the weak s-component and a larger production of F. Since rotation also affects the mass-loss rate, either directly or indirectly, we find substantial differences in the lifetimes as O-type and Wolf-Rayet subtypes between the rotating and non-rotating models. The maximum mass exploding as Type IIP supernova ranges between 15 and 20 M {sub Sun} in both sets of models (this value depends basically on the larger mass-loss rates in the red supergiant phase due to the inclusion of the dust-driven wind). This limiting value is in remarkably good agreement with current estimates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

WINDS, CLUMPS, AND INTERACTING COSMIC RAYS IN M82

We construct a family of models for the evolution of energetic particles in the starburst galaxy M82 and compare them to observations to test the calorimeter assumption that all cosmic ray energy is radiated in the starburst region. Assuming constant cosmic ray acceleration efficiency with Milky Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron populations as a function of energy. Cosmic rays are injected with Galactic energy distributions and electron-to-proton ratio via Type II supernovae at the observed rate of 0.07 yr{sup -1}. From the cosmic ray spectra, we predict the radio synchrotron and {gamma}-ray spectra. To more accurately model the radio spectrum, we incorporate a multiphase interstellar medium in the starburst region of M82. Our model interstellar medium is highly fragmented with compact dense molecular clouds and dense photoionized gas, both embedded in a hot, low density medium in overall pressure equilibrium. The spectra predicted by this one-zone model are compared to the observed radio and {gamma}-ray spectra of M82. {chi}{sup 2} tests are used with radio and {gamma}-ray observations and a range of model predictions to find the best-fit parameters. The best-fit model yields constraints on key parameters in the starburst zone of M82, including a magnetic field strength of {approx}250 {mu}G and a wind advection speed in the range of 300-700 km s{sup -1}. We find that M82 is a good electron calorimeter but not an ideal cosmic-ray proton calorimeter and discuss the implications of our results for the astrophysics of the far-infrared-radio correlation in starburst galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE GLOBAL SOLAR MAGNETIC FIELD-IDENTIFICATION OF TRAVELING, LONG-LIVED RIPPLES

We have examined the global structure of the solar magnetic field using data from the Fe I spectral line at 5250.2 A obtained at the 150 foot tower telescope at the Mt. Wilson Observatory. For each point on the solar surface, we find the value of the magnetic field in the meridional plane, B{sub m} , by averaging over all available observations using a cosine weighting method. We have revised our cosine weighting method by now taking into account more fully the highest latitude geometry. We use the annual variation in the latitude of the disk center, b{sub 0}, to deduce the tilt angle of the field relative to the local vertical so that we can find the radial component of the field, B{sub r} , from B{sub m} . We find this tilt angle to be small except for a near-polar zone where a tilt-angle model can reduce the annual variation. The reduced annual variation in the deduced B{sub r} allows us to study dB{sub r} /dt and associated deviations in B{sub r} from a smoothed B{sub r} with a smoothing width of 2.5 yr. These functions make evident the presence of small amplitude (3-5 G) but spatially coherent ripples with a semi-regular periodicity of one to three years. At any given time, the half-wavelength (peak to trough) is between 15 Degree-Sign and 30 Degree-Sign of latitude. These patterns are ubiquitous and in many cases drift from near the equator to the poles over a time period of roughly two years. The drift rate pattern is not compatible with simple advection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ACCRETION DISK TEMPERATURES OF QSOs: CONSTRAINTS FROM THE EMISSION LINES

We compare QSO emission-line spectra to predictions based on theoretical ionizing continua of accretion disks. The observed line intensities do not show the expected trend of higher ionization with theoretical accretion disk temperature as predicted from the black hole mass and accretion rate. Consistent with earlier studies, this suggests that the inner disk does not reach temperatures as high as expected from standard disk theory. Modified radial temperature profiles, taking account of winds or advection in the inner disk, achieve better agreement with observation. The emission lines of radio-detected and radio-undetected sources show different trends as a function of the theoretically predicted disk temperature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A PHYSICAL LINK BETWEEN JET FORMATION AND HOT PLASMA IN ACTIVE GALACTIC NUCLEI

Recent observations suggest that in black hole X-ray binaries jet/outflow formation is related to the hot plasma in the vicinity of the black hole, either in the form of an advection-dominated accretion flow at low accretion rates or in a disk corona at high accretion rates. We test the viability of this scenario for supermassive black holes using two samples of active galactic nuclei distinguished by the presence (radio-strong) and absence (radio-weak) of well-collimated, relativistic jets. Each is centered on a narrow range of black hole mass but spans a very broad range of Eddington ratios, effectively simulating in a statistical manner the behavior of a single black hole evolving across a wide spread in accretion states. Unlike the relationship between the radio and optical luminosity, which shows an abrupt break between high- and low-luminosity sources at an Eddington ratio of {approx}1%, the radio emission-a measure of the jet power-varies continuously with the hard X-ray (2-10 keV) luminosity, roughly as L{sub R} {proportional_to} L{sub X}{sup 0.6-0.75}. This relation, which holds for both radio-weak and radio-strong active galaxies, is similar to the one seen in X-ray binaries. Jet/outflow formation appears to be closely linked to the conditions that give rise to the hot, optically thin coronal emission associated with accretion flows, both in the regime of low and high accretion rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

A NUMERICAL TREATMENT OF ANISOTROPIC RADIATION FIELDS COUPLED WITH RELATIVISTIC RESISTIVE MAGNETOFLUIDS

We develop a numerical scheme for solving fully special relativistic, resistive radiation magnetohydrodynamics. Our code guarantees conservation of total mass, momentum, and energy. The radiation energy density and the radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation field, in contrast to the Eddington approximation, as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly the Harten-Lax-van Leer scheme; the HLLC and HLLD schemes are also solved for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in an optically thick regime or a large magnetic Reynolds number regime. Although we need to invert 4 Multiplication-Sign 4 matrices (for the gas-radiation interaction) and 3 Multiplication-Sign 3 matrices (for the magnetic energy dissipation) at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek-type magnetic reconnection, revealing the reduction of the reconnection rate via radiation drag.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

SWIFT AND RXTE OBSERVATIONS OF THE BLACK HOLE TRANSIENT GX 339-4 IN THE HARD STATE BETWEEN OUTBURSTS

We use simultaneous Swift and RXTE observations of the black hole binary GX 339-4 to measure the inner radius of its accretion disk in the hard state down to 0.4% L{sub Edd} via modeling of the thermal disk emission and the relativistically broadened iron line. For the luminosity range covered in this work, our results rule out a significantly truncated disk at 100-1000 R{sub g} as predicted by the advection-dominated accretion flow paradigm. The measurements depend strongly on the assumed emission geometry, with most results providing no clear picture of radius evolution. If the inclination is constrained to roughly 20 Degree-Sign , however, the measurements based on the thermal disk emission suggest a mildly receding disk at a luminosity of 0.4% L{sub Edd}. The iron abundance varies between {approx}1 and 2 solar abundances, with the i = 20 Degree-Sign results indicating a negative correlation with luminosity, though this is likely due to a change in disk illumination geometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

ATMOSPHERIC HEAT REDISTRIBUTION ON HOT JUPITERS

Infrared light curves of transiting hot Jupiters present a trend in which the atmospheres of the hottest planets are less efficient at redistributing the stellar energy absorbed on their daysides—and thus have a larger day-night temperature contrast—than colder planets. To this day, no predictive atmospheric model has been published that identifies which dynamical mechanisms determine the atmospheric heat redistribution efficiency on tidally locked exoplanets. Here we present a shallow-water model of the atmospheric dynamics on synchronously rotating planets that explains why heat redistribution efficiency drops as stellar insolation rises. Our model shows that planets with weak friction and weak irradiation exhibit a banded zonal flow with minimal day-night temperature differences, while models with strong irradiation and/or strong friction exhibit a day-night flow pattern with order-unity fractional day-night temperature differences. To interpret the model, we develop a scaling theory which shows that the timescale for gravity waves to propagate horizontally over planetary scales, τ{sub wave}, plays a dominant role in controlling the transition from small to large temperature contrasts. This implies that heat redistribution is governed by a wave-like process, similar to the one responsible for the weak temperature gradients in the Earth's tropics. When atmospheric drag can be neglected, the transition from small to large day-night temperature contrasts occurs when τ{sub wave}∼√(τ{sub rad}/Ω), where τ{sub rad} is the radiative relaxation time and Ω is the planetary rotation frequency. Alternatively, this transition criterion can be expressed as τ{sub rad} ∼ τ{sub vert}, where τ{sub vert} is the timescale for a fluid parcel to move vertically over the difference in day-night thickness. These results subsume the more widely used timescale comparison for estimating heat redistribution efficiency between τ{sub rad} and the horizontal day-night advection timescale, τ{sub adv}. Only because τ{sub adv} ∼ τ{sub vert} for hot Jupiters does the commonly assumed timescale comparison between τ{sub rad} and τ{sub adv} yield approximately correct predictions for the heat redistribution efficiency.

79 ASTRONOMY AND ASTROPHYSICS↗