Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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,495 records · Page 83

Collisionless shocks in the heliosphere: Reviews of current research

The present conference on heliospheric collisionless shocks considers such macrostructure-, microstructure-, and particle acceleration-related topics as interplanetary shock phenomena near and within 1 AU, as well as beyond, planetary bow shocks, shock formation and evolution in the solar atmosphere, MHD and gasdynamic theories for planetary bow waves, and subcritical collisionless shock waves. Also discussed are ion reflection, gyration, and dissipation at supercritical shocks, the numerical simulation of quasi-perpendicular collisionless shocks, electron distributions near collisionless shocks, the microtheory of collisionless shock current layers, plasma waves and instabilities, the electron foreshock, upstream suprathermal ions, and both diffusive and shock drift acceleration.

Tsurutani, Bruce T.↗

A two-dimensional numerical simulation of a supersonic, chemically reacting mixing layer

Research has been undertaken to achieve an improved understanding of physical phenomena present when a supersonic flow undergoes chemical reaction. A detailed understanding of supersonic reacting flows is necessary to successfully develop advanced propulsion systems now planned for use late in this century and beyond. In order to explore such flows, a study was begun to create appropriate physical models for describing supersonic combustion, and to develop accurate and efficient numerical techniques for solving the governing equations that result from these models. From this work, two computer programs were written to study reacting flows. Both programs were constructed to consider the multicomponent diffusion and convection of important chemical species, the finite rate reaction of these species, and the resulting interaction of the fluid mechanics and the chemistry. The first program employed a finite difference scheme for integrating the governing equations, whereas the second used a hybrid Chebyshev pseudospectral technique for improved accuracy.

Drummond, J. Philip↗

Tidal evolution of globular clusters. I - Method

Tidal evolution of globular clusters is regulated by both Galactic tidal effects and internal relaxation processes. In order to investigate the tidal evolution of globular clusters, a numerical scheme which utilizes a Fokker-Planck approach as well as direct numerical integration of the restricted three-body problem is developed. In the inner regions of the cluster, stellar orbits are mapped with the cluster's gravitational potential and orbit-averaged diffusion coefficients. In the outer regions, the Galactic tidal field is explicitly included in the direct orbital integration. This method is presented here with some tests on King-Michie models.

Oh, K. S.↗

Role of acoustics in flame/vortex interactions

The role of acoustics in flame/vortex interactions is examined via asymptotic analysis and numerical simulation. The model consists of a one-step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes which are allowed to mix and react by convection and diffusion in the presence of an acoustic field or a time-varying pressure field of small amplitude. The main emphasis is on the influence of the acoustics on the ignition time and flame structure as a function of vortex Reynolds number and initial temperature differences of the reactants.

Jackson, T. L.↗

Time-dependent cosmic ray shocks with injection - A progress report

The study presents preliminary results of a two-fluid hydrodynamical simulation of cosmic-ray-modified shock structure by extending the steady-state calculations of Zank et al. (1991) to the time-dependent case. The present simulations model diffusive shock acceleration of the energetic particle component at planar shocks. A new model for the injection of cosmic ray particles is presented. The quadratic artificial viscosity algorithm is found to be numerically unstable and therefore unusable for strong cosmic ray shocks. There is visible ringing and overshoot in the various profiles. A large enhancement in the fluid density appears downstream of the subshock. This 'overcompression' results from the upstream fluid being compressed in the precursor, then traversing a strong shock with compression ratio near the adiabatic limit of four.

Donohue, D. J.↗

Diffusion Limited Supercritical Water Oxidation (SCWO) in Microgravity Environments

Tests designed to quantify the gravitational effects on thermal mixing and reactant injection in a Supercritical Water Oxidation (SCWO) reactor have recently been performed in the Zero Gravity Facility (ZGF) at NASA s Glenn Research Center. An artificial waste stream, comprising aqueous mixtures of methanol, was pressurized to approximately 250 atm and then heated to 450 C. After uniform temperatures in the reactor were verified, a controlled injection of air was initiated through a specially designed injector to simulate diffusion limited reactions typical in most continuous flow reactors. Results from a thermal mapping of the reaction zone in both 1-g and 0-g environments are compared. Additionally, results of a numerical model of the test configuration are presented to illustrate first order effects on reactant mixing and thermal transport in the absence of gravity.

Hicks, M. C.↗

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗

Implementation of Charged Particle Behavior in Discrete Element Method (DEM) Simulations

Lunar landers will agitate the surface of the Moon with an exhaust plume during descent which will, due to the particulate nature of the lunar regolith, loosen and eject grains from the surface. This ejection is not only coupled with the charged plume gas, but also results in significant particle-particle interactions. Settling of these grains after plume effects have subsided takes much longer than expected in a ballistic trajectory. The prevailing hypothesis attributes this behavior to the accumulated charge on the particles. We are thus developing a discrete element method (DEM) approach to explore these charged particle interactions on the lunar surface. The Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Improved for General Granular and Granular Heat Transfer Simulations (LIGGGHTS) software package provides a DEM modeling framework for granular interactions. It includes many complexities such as non-spherical particle shapes, cohesion and frictional forces, and heat transfer, but has no provision for inter-particle electrostatic forces and charge transfer that are important to examine in the lunar environment. In this work, a standard Coulomb potential and a Yukawa potential are integrated into the LIGGGHTS framework to provide a basis for particle-particle electrostatic interactions, as well as a gravitational potential to enable inter-grain gravitational interactions. A preliminary approach to charge transfer between grains incorporating properties such as work function and electrical conductivity to the library of available material characteristics will be presented. Several scenarios have been simulated that include charged particle interactions within a diffuse granular gas, settling of charged grains into a regolith bed, sliding of granular material along an incline, and vibration of settled grains to produce a behavior similar to fluidization. There are numerous challenges to incorporate realistic interactions between complex lunar particles. Currently, grains are modeled to behave as if the entirety of the charge acts at the center of mass, such as conductors with spherical symmetry and insulators with homogeneously distributed charge. We are developing improvements that will include the use of non-spherical particle geometries, as well as reasonable approximations of insulating/dielectric materials that have non-uniform charge distributions. The cases simulated thus far will be examined in a relevant environment within a vacuum chamber to validate the simulations. These simulations will be bounded by experiments utilizing high-speed camera observations of the motion for validation. The grains in the experiment will exchange charge during their motion and this can be quantified by collection within a charge measurement device such as a Faraday cup. Such a device may be modeled within the software by defining an integration region and computing the contained charge as a function of simulation time, allowing for side-by-side comparison of simulated and measured bulk charging results. Any differences will be reconciled by updating the mathematical mechanisms described within the simulation suite. Successfully combining results from experiments within a relevant environment into the LIGGGHTS framework will improve modeling of the charged grain dynamics experienced on the Moon to provide insights into dust behavior for future lunar exploration missions.

Electrostatics↗

On the importance of the heat and mass transfer resistances in internally-cooled liquid desiccant dehumidifiers and regenerators

Liquid desiccant heat and mass exchangers are a promising technology for efficient humidity control in buildings. Many researchers have investigated these exchangers, often using numerical models to predict their performance. However, there is a lack of information in the literature on the magnitude of the heat and mass transfer resistances, both for the dehumidifier (which absorbs moisture from the air) and the regenerator (which heats the liquid desiccant to re-concentrate it). This article focuses on internally-cooled, 3-fluid exchangers in a parallel plate geometry. Water heats or cools a desiccant across a plate, and the desiccant absorbs or releases water into an airstream through a membrane. A sensitivity analysis was used to estimate the importance of each of the heat and mass transfer resistances (air, membrane, desiccant, plate, water), and how it changes with different design geometries. The results show that, for most designs, the latent and sensible heat transfer of the dehumidifier is dominated by the air mass transfer resistance and air heat transfer resistance, respectively. The air mass transfer resistance is also important for the regenerator, but much less so; the change in the desiccant equilibrium humidity ratio due to a change in either temperature or desiccant mass fraction is much higher at the regenerator's higher temperatures. This increases the importance of (1) getting heat from the water to the desiccant/membrane interface, and (2) diffusing salt ions quickly away from the desiccant/membrane interface. The membrane heat transfer and water heat transfer resistances were found to be the least important. These results can help inform decisions about what simplifying assumptions to make in numerical models, and can also help in designing these exchangers by understanding which resistances are most important.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The implications of the COBE diffuse microwave radiation results for cosmic strings

We compare the anisotropies in the cosmic microwave background radiation measured by the COBE experiment to those predicted by cosmic string theories. We use an analytic model for the Delta T/T power spectrum that is based on our previous numerical simulations of strings, under the assumption that cosmic strings are the sole source of the measured anisotropy. This implies a value for the string mass per unit length of 1.5 +/- 0.5 x 10 exp -6 C-squared/G. This is within the range of values required for cosmic strings to successfully seed the formation of large-scale structures in the universe. These results clearly encourage further studies of Delta T/T and large-scale structure in the cosmic string model.

Bennett, David P.↗

Effects of Fe atoms on hardening of a nickel matrix: Nanoindentation experiments and atom-scale numerical modeling

Significant hardening effect due to Fe concentrations in Ni-based alloys with face-centered-cubic structure has been studied by using a combined experimental and atomistic-based computational approach via nanoindentation tests. The obtained experimental load–displacement data for the [0 0 1] crystal orientation reached a qualitative good agreement with molecular dynamics simulations results, leading to strong evidence that the main strengthening factors are associated to sluggish dislocation diffusion, reduced defect sizes and the nucleation of tetrahedral stacking faults. Here, interstitial type prismatic dislocation loops mainly formed by $\frac{1}{6}$$\langle$1 1 2$\rangle$ Shockley dislocations are nucleated during the loading process, where their interaction leads to the formation of pyramidal shaped stacking fault which are mainly created by $\frac{1}{3}$$\langle$1 0 0$\rangle$ Hirth dislocations lines. Observing both types of defects coexisting in the same plastic deformation zone by both approaches. Reported mechanical data, measured experimentally and interpreted numerically, are also in accordance with microstructural SEM and TEM investigations.

36 MATERIALS SCIENCE↗

Planck 2018 results XII. Galactic astrophysics using polarized dust emission

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle ψ, and dispersion function of polarization angles S . The joint distribution (one-point statistics) of p and NH confirms that the mean and maximum polarization fractions decrease with increasing NH. The uncertainty on the maximum observed polarization fraction, pmax = 22.0 -1.4 +3.5% at 353 GHz and 80' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and S found earlier is seen to be present at high latitudes. This follows the S ∝ p -1 relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, ψ, and S for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S × p, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between NH = 10 20 cm -2 and NH = 2 x 10 22 cm -2 , out of the Galactic plane, this product S × p only decreases by about 25%. Because S is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of S × p with the dust temperature T d , whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion S tend, on the contrary, to have colder dust than lines of sight with low p and high S . We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R P/p and R S/V , which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p V /E(B - V), of 13% at high Galactic latitude, compatible with the polarization fraction p ≈ 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

79 ASTRONOMY AND ASTROPHYSICS↗

NOAA Air Resources Laboratory Atmospheric Turbulence and Diffusion Division Contribution to LAFE Field Campaign Report

The Land-Atmosphere Feedback Experiment (LAFE) was conducted between 1 and 31 August, 2017 at the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) observatory. The main scientific goal of LAFE was to investigate feedbacks occurring between different land surface types and the overlying atmosphere to improve turbulence parameterizations in numerical weather prediction models. LAFE collaborators included scientists from the following institutions: Cleveland State University Cooperative Institute for Meteorological Satellite Studies (CIMSS) Cooperative Institute for Research in Environmental Sciences (CIRES) Leibniz University Institute of Meteorology and Climatology Marquette University NASA Goddard Space Flight Center National Center for Atmospheric Research (NCAR) National Severe Storms Laboratory (NSSL) NOAA Earth System Research Laboratory (ESRL) Pacific Northwest National Laboratories (PNNL) University of Hohenheim Institute of Physics and Meteorology University of Tennessee Space Institute (UTSI). In addition to helping plan the scientific objectives of LAFE, collaborators deployed several lidars (i.e., temperature, water, vapor, and wind) and remote-sensing instruments (e.g., microwave radiometers, atmospheric emitted radiance interferometers [AERIs], etc.) to complement the suite of instrument platforms already present at the SGP site. The measurement platforms, as well as scientific objectives of LAFE, are summarized in Wulfmeyer et al. (2018). To complement these platforms, the NOAA Atmospheric Turbulence and Diffusion Division (ATDD) installed micrometeorological towers to continuously measure fluxes of momentum, sensible heat, latent heat, and CO 2 at 2 m and 10 m above ground level (AGL) at three sites during LAFE. Tower measurements also included temperature, humidity, wind speed and direction, pressure, rainfall, soil temperature, and soil moisture. The three micrometeorological towers were installed between 9 and 13 July, 2017 and were removed between 12 and 14 September, 2017. During the experiment, the footprint of Tower 1 included a mixture of soybean and grassland that transitioned to mature soybean crop during the second half of the experiment. Tower 2’s footprint encompassed a field of grazed pasture, and Tower 3’s footprint was a mature soybean crop. The meteorological and flux measurements from these towers were/are being used to: Increase the spatial density of measurement platforms that were already available at the ARM facility and that were deployed during LAFE Improve the characterization of surface fluxes around the SGP site Explore the validity of applying Monin-Obukhov Similarity Theory (MOST) over different land surface types and improve knowledge of flux-gradient relationships Provide initial conditions and be used for verification for large-eddy simulations (LES) Evaluate numerical weather prediction models, e.g., the High-Resolution Rapid Refresh (HRRR), to identify and correct for biases within these models Evaluate the role of low-level advection at each site, which would be present in the flux-divergence data sets obtained from the three towers.

Southern Great Plains, unmanned aerial systems, su↗

Neutron diffusion calculation in heterogeneous geometry based on local/global iteration using proper orthogonal decomposition

This study newly proposes a heterogeneous core calculation method based on local/global iteration using proper orthogonal decomposition (POD). By using the singular value decomposition (SVD) and the low-rank approximation, appropriate POD bases for expanding the neutron flux can be obtained from snapshot data of the neutron flux obtained by fine mesh calculations. By projection using the POD bases, the dimension of the target equation (e.g., discretized neutron diffusion equation) can be dramatically reduced. In the proposed method, POD is effectively applied to each single assembly calculation (local calculation). Furthermore, using the local/global iteration, the effective neutron multiplication factor and the neutron flux distribution in the whole core geometry can be obtained by combining the numerical results of the local calculation for each fuel assembly and the global calculation for the whole core. As a feasibility study, the proposed method is applied to a one-dimensional heterogeneous core analysis, and the accuracy is investigated by changing the total number of POD bases. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Laminar flamelet modeling of turbulent diffusion flames

In modeling turbulent combustion, decoupling the chemistry from the turbulence is of great practical significance. In cases in which the equilibrium chemistry model breaks down, laminar flamelet modeling (LFM) is a promising approach to decoupling. Here, the validity of this approach is investigated using direct numerical simulation of a simple chemical reaction in two-dimensional turbulence.

Mell, W. E.↗

Multiscale models and stochastic simulation methods for computing rare but key binding events in cell biology

The main difficulty in simulating diffusion processes at a molecular level in cell microdomains is due to the multiple scales involving nano- to micrometers. Few to many particles have to be simulated and simultaneously tracked while there are exploring a large portion of the space for binding small targets, such as buffers or active sites. Bridging the small and large spatial scales is achieved by rare events representing Brownian particles finding small targets and characterized by long-time distribution. These rare events are the bottleneck of numerical simulations. A naive stochastic simulation requires running many Brownian particles together, which is computationally greedy and inefficient. Solving the associated partial differential equations is also difficult due to the time dependent boundary conditions, narrow passages and mixed boundary conditions at small windows. We present here two reduced modeling approaches for a fast computation of diffusing fluxes in microdomains. The first approach is based on a Markov mass-action law equations coupled to a Markov chain. The second is a Gillespie's method based on the narrow escape theory for coarse-graining the geometry of the domain into Poissonian rates. The main application concerns diffusion in cellular biology, where we compute as an example the distribution of arrival times of calcium ions to small hidden targets to trigger vesicular release.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations for solvent-based carbon capture. Part 2: Chemical absorption across a wetted wall column: Original Research Article: Hierarchical calibration and validation framework of bench-scale computational fluid dynamics simulations

Part 1 of this paper presents a numerical model for non-reactive physical mass transfer across a wetted wall column (WWC). In Part 2, we improved the existing computational fluid dynamics (CFD) model to simulate chemical absorption occurring in a WWC as a bench-scale study of solvent-based carbon dioxide (CO 2 ) capture. In this study, to generate data for WWC model validation, CO 2 mass transfer across a monoethanolamine (MEA) solvent was first measured on a WWC experimental apparatus. The numerical model developed in this work can account for both chemical absorption and desorption of CO 2 in MEA. In addition, the overall mass transfer coefficient predicted using traditional/empirical correlations is conducted and compared with CFD prediction results for both steady and wavy falling films. A Bayesian statistical calibration algorithm is adopted to calibrate the reaction rate constants in chemical absorption/desorption of CO 2 across a falling film of MEA. The posterior distributions of the two transport properties, i.e., Henry's constant and gas diffusivity in the non-reacting nitrous oxide (N 2 O)/MEA system obtained from Part 1 of this study, serves as priors for the calibration of CO 2 reaction rate constants after using the N 2 O/CO 2 analogy method. Finally, the calibrated model can be used to predict the CO 2 mass transfer in a WWC for a wider range of operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗