Engineering Papers⌕ Search

Engineering topics

Bell, J. B.

Publications and source records attributed to Bell, J. B..

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A new re-redistribution scheme for weighted state redistribution with adaptive mesh refinement

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. Here, we present numerical results that demonstrate that the algorithm is conservative when the coarse-fine interface intersects the embedded boundary. Additionally we compare the numerical solution of the Sod shock tube problem in an inclined cylinder with the analytic solution, and we compare the simulation of a shock hitting a cylindrical obstacle with experimental data. Finally we demonstrate the methodology for simulation of the multicomponent compressible Navier-Stokes equations in a piston-bowl geometry, and discuss the computational efficiency gained by not requiring the entire embedded boundary to be defined at the finest level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeling electrokinetic flows with the discrete ion stochastic continuum overdamped solvent algorithm

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the discrete ion stochastic continuum overdamped solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Herein, several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced-charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a “dry diffusion” approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Discrete ion stochastic continuum overdamped solvent algorithm for modeling electrolytes

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Towards the distributed burning regime in turbulent premixed flames

Three-dimensional numerical simulations of canonical statistically steady, statistically planar turbulent flames have been used in an attempt to produce distributed burning in lean methane and hydrogen flames. Dilatation across the flame means that extremely large Karlovitz numbers are required; even at the extreme levels of turbulence studied (up to a Karlovitz number of 8767) distributed burning was only achieved in the hydrogen case. In this case, turbulence was found to broaden the reaction zone visually by around an order of magnitude, and thermodiffusive effects (typically present for lean hydrogen flames) were not observed. In the preheat zone, the species compositions differ considerably from those of one-dimensional flames based a number of different transport models (mixture averaged, unity Lewis number and a turbulent eddy viscosity model). The behaviour is a characteristic of turbulence dominating non-unity Lewis number species transport, and the distinct limit is again attributed to dilatation and its effect on the turbulence. Peak local reaction rates are found to be lower in the distributed case than in the lower Karlovitz cases but higher than in the laminar flame, which is attributed to effects that arise from the modified fuel-temperature distribution that results from turbulent mixing dominating low Lewis number thermodiffusive effects. Finally, approaches to achieve distributed burning at realisable conditions are discussed; factors that increase the likelihood of realising distributed burning are higher pressure, lower equivalence ratio, higher Lewis number and lower reactant temperature.

42 ENGINEERING↗
DOE OSTI2019

Toward Resolved Simulations of Burning Fronts in Thermonuclear X-ray Bursts

© 2019 Published under licence by IOP Publishing Ltd. We discuss the challenges of modeling X-ray bursts in multi-dimensions, review the different calculations done to date, and discuss our new set of ongoing simulations. We also describe algorithmic improvements that may help in the future to offset some of the expense of these simulations, and describe what may be possible with exascale computing.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Explosion-induced ignition and combustion of acetylene clouds

We investigate the explosion-induced ignition and combustion of an acetylene cloud in a rectangular chamber (10 cm × 10 cm × 39 cm). In the experiments, a 0.2-g PETN charge was located at $x = 9.7\,\hbox {cm}$ and a soap bubble (d=5cm) filled with pure acetylene was located at $x = 27\,\hbox {cm}$ as measured relative to the end wall. Detonation of the charge created a blast wave that crushed the soap bubble—inducing mixing with the air. After 0.55 ms, the mixture ignited, forming at turbulent combustion cloud. The flow was modeled using the compressible Navier–Stokes equations assuming unity Lewis number. Arrhenius-based kinetics were used to model ignition. Adaptive mesh refinement was used to capture turbulent mixing on the grid (the MILES approach of J. Boris). Computed pressures were found to be in agreement with measured pressure histories. Finite-rate kinetics were required to capture the ignition processes over the duration of the experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Meeting the Challenges of Modeling Astrophysical Thermonuclear Explosions: Castro, Maestro, and the AMReX Astrophysics Suite

We describe the AMReX suite of astrophysics codes and their application to modeling problems in stellar astrophysics. Maestro is tuned to efficiently model subsonic convective flows while Castro models the highly compressible flows associated with stellar explosions. Both are built on the block-structured adaptive mesh refinement library AMReX. Together, these codes enable a thorough investigation of stellar phenomena, including Type Ia supernovae and X-ray bursts. We describe these science applications and the approach we are taking to make these codes performant on current and future many-core and GPU-based architectures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Turbulent jets with off-source heating

Motivated by anomalous entrainment behaviour in cumulus clouds, Bhat et al. pioneered a laboratory experiment to study turbulent jets subjected to a volumetric heating away from the momentum source. The study concluded that the use of a constant entrainment coefficient was insufficient for the flow, and that the results did not confirm the analysis of Hunt, which suggested that an increase in relative turbulent transport of streamwise momentum could lead to a decrease in entrainment. The present paper re-evaluates theoretical aspects of both studies, and includes a decomposition of the factors contributing to entrainment. We then use the reworked analysis to examine three-dimensional numerical simulations of turbulent jets with off-source heating. The data are consistent with previous work, but give deeper insight not easily obtainable through experiment. Specifically, direct measurement of flux integrals shows that previous inference from experimental measurements of centreline velocity and profile widths under the assumption of self-similarity can lead to underestimation of the mass flux by over 50 % in some cases. Radial profiles of temperature, radial velocity and turbulent correlations show significant departures from self-similarity. The flux measurements show that there is actually an increase in the entrainment coefficient with heating, and that it is locally enhanced by positive forcing and decreased by an increase in turbulent transport of streamwise momentum, thereby confirming the essence of the original proposal of Hunt.

42 ENGINEERING↗
DOE OSTI2016

LOW MACH NUMBER MODELING OF CONVECTION IN HELIUM SHELLS ON SUB-CHANDRASEKHAR WHITE DWARFS. II. BULK PROPERTIES OF SIMPLE MODELS

The dynamics of helium shell convection driven by nuclear burning establish the conditions for runaway in the sub-Chandrasekhar-mass, double-detonation model for SNe Ia, as well as for a variety of other explosive phenomena. We explore these convection dynamics for a range of white dwarf core and helium shell masses in three dimensions using the low Mach number hydrodynamics code MAESTRO. We present calculations of the bulk properties of this evolution, including time-series evolution of global diagnostics, lateral averages of the 3D state, and the global 3D state. We find a variety of outcomes, including quasi-equilibrium, localized runaway, and convective runaway. Our results suggest that the double-detonation progenitor model is promising and that 3D dynamic convection plays a key role.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Three-dimensional direct numerical simulation of turbulent lean premixed methane combustion with detailed kinetics

The interaction of maintained homogeneous isotropic turbulence with lean premixed methane flames is investigated using direct numerical simulation with detailed chemistry. The conditions are chosen to be close to those found in atmospheric laboratory experiments. As the Karlovitz number is increased from 1 to 36, the preheat zone becomes thickened, while the reaction zone remains largely unaffected. A negative correlation of fuel consumption with mean flame surface curvature is observed. With increasing turbulence intensity, the chemical composition in the preheat zone tends towards that of an idealised unity Lewis number flame, which we argue is the onset of the transition to distributed burning, and the response of the various chemical species is shown to fall into broad classes. Smaller-scale simulations are used to isolate the specific role of species diffusion at high turbulent intensities. Diffusion of atomic hydrogen is shown to be related to the observed curvature correlations, but does not have significant consequential impact on the thickening of the preheat zone. It is also shown that susceptibility of the preheat zone to thickening by turbulence is related to the 'global' Lewis number (the Lewis number of the deficient reactant); higher global Lewis number flames tend to be more prone to thickening.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

COMPARISONS OF TWO- AND THREE-DIMENSIONAL CONVECTION IN TYPE I X-RAY BURSTS

We perform the first detailed three-dimensional simulation of low Mach number convection preceding runaway thermonuclear ignition in a mixed H/He X-ray burst. Our simulations include a moderate-sized, approximate network that captures hydrogen and helium burning up through rp-process breakout. We look at the difference between two- and three-dimensional convective fields, including the details of the turbulent convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

The deflagration stage of Chandrasekhar mass models for type Ia supernovae. I. Early evolution

We present high-resolution, full-star simulations of the post-ignition phase of Type Ia supernovae using the compressible hydrodynamics code Castro. Initial conditions, including the turbulent velocity field and ignition site, are imported directly from a simulation of the last few hours of presupernova convection using a low Mach number code, Maestro. Adaptive mesh refinement allows the initial burning front to be modeled with an effective resolution of 36,864{sup 3} zones (136 m zone{sup –1}). The initial rise and expansion of the deflagration front are tracked until burning reaches the star's edge and the role of the background turbulence on the flame is investigated. The effect of artificially moving the ignition location closer to the star's center is explored. The degree to which turbulence affects the burning front decreases with increasing ignition radius since the buoyancy force is stronger at larger radii. Even central ignition—in the presence of a background convective flow field—is rapidly carried off-center as the flame is carried by the flow field. We compare our results to analytic models for burning thermals, and find that they reproduce the general trends of the bubble's size and mass, but underpredict the amount of buoyant acceleration due to simplifying assumptions of the bubble's properties. Overall, we find that the amount of mass that burns prior to flame break out is small, consistent with a gravitationally confined detonation' occurring at a later epoch, but additional burning will occur following breakout that may modify this conclusion.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

An attempt to comprehend Martian weathering conditions through the analysis of terrestrial palagonite samples

Spectroscopic observations of the Martian surface in the invisible to near infrared (0.4-1.0 micron), coupled with measurements made by Viking, have shown that the surface is composed of a mixture of fine-grained weathered and nonweathered minerals. The majority of the weathered components are thought to be materials like smectite clays, scapolite, or palagonite. Until materials are returned for analysis there are two possible ways of proceeding with an investigation of Martian surface processes: (1) the study of weathering products in meteorites that have a Martian origin (SNC's), and (2) the analysis of certain terrestrial weathering products as analogs to the material found in SNC's, or predicted to be present on the Martian surface. We describe some preliminary measurements of the carbon chemistry of terrestrial palagonite samples that exhibit spectroscopic similarities with the Martian surface. The data should aid the understanding of weathering in SNC's and comparisons between terrestrial palagonites and the Martian surface.

Douglas, C.↗
NASA NTRS1981

Fully implicit shock tracking

A tracking method for solving hyperbolic free boundary problems arising in fluid dynamics is presented. The principal feature of the new approach is that the implicit form is used to simultaneously treat both interior points and boundary conditions. In particular, the location of the free boundary (shock) surface is treated implicitly and coupled with all other unknowns. The method is presented in the context of unsteady one-dimensional flow in a variable area duct with an internal shock wave. The fully implicit method and other strategies for advancing the shock are compared for computing a steady solution via a time asymptotic approach. Issues regarding extension of the method to multiple dimensions are also discussed.

Bell, J. B.↗