Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Imperfect 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 55 records · Page 3

Giant Planets in Reflected Light: What Science Can We Expect?

Interpreting the reflection spectra of cool giant planets will be a challenge. Spectra of such worlds are expected to be primarily shaped by scattering from clouds and hazes and punctuated by absorption bands of methane, water, and ammonia. While the warmest giants may be cloudless, their atmospheres will almost certainly sport substantial photochemical hazes. Furthermore the masses of most direct imaging targets will be constrained by radial velocity observations, their radii, and thus atmospheric gravity, will be imperfectly known. The uncertainty in planet radius and gravity will compound with uncertain aerosol properties to make estimation of key absorber abundances difficult. To address such concerns our group is developing atmospheric retrieval tools to constrain quantities of interest, particular gas mixing ratios. We have applied our Markov Chain Monte Carlo methods to simulated data of the quality expected from the WFIRST CGI instrument and found that given sufficiently high SNR data we can confidentially identify and constrain the abundance of methane, cloud top pressures, gravity, and the star-planet-observer phase angle. In my presentation I will explain the expected characteristics of cool extrasolar giant planet reflection spectra, discuss these and other challenges in their interpretation, and summarize the science results we can expect from direct imaging observations.

methane↗

Compositions of Spherules and Rock Surfaces at Meridiani

The Alpha Particle X-ray Spectrometers (APXS) on the Mars Exploration Rovers (MER) have proven extremely valuable for analyzing rocks and soils on the surface of Mars. The precision of their compositional measurements has been shown to be phenomenal through analyses of the compositionally very uniform Meridiani soils. Through combined use of the rock abrasion tool (RAT) and the analytical instruments on the in-situ deployment device (IDD), analyses of the interiors of fine-grained and texturally uniform rocks with surfaces ground flat have been made under conditions that are nearly ideal for this mode of analysis. The APXS has also been used frequently to analyze materials whose characteristics, surface morphologies, and sample-detector geometries are less than ideal, but the analyses of which are nonetheless very useful for understanding the makeup of the target materials. Such targets include undisturbed rocks with irregular and sometimes coated surfaces and mixed targets such as soils that include fine-grained components as well as coarse grains and pieces of rocks. Such target materials include the well known hematite-rich concretions, referred to as blueberries because of their multispectral color, size, and mode of occurrence. In addition to non-ideal target geometry, such mixed materials also present a heterogeneous target in terms of density. These irregularities violate the assumptions commonly associated with analyses done in the laboratory to achieve the highest possible accuracy. Here we acknowledge the irregularities and we examine the inferences drawn from specific chemical trends obtained on imperfect targets in light of one of the potential pitfalls of natural materials on the surface of Mars, namely thin dust coatings.

Mittlefehldt, David W.↗

Implementation of a Binary Neural Network on a Passive Array of Magnetic Tunnel Junctions

The increasing scale of neural networks and their growing application space have produced demand for more energy- and memory-efficient artificial-intelligence-specific hardware. Avenues to mitigate the main issue, the von Neumann bottleneck, include in-memory and near-memory architectures, as well as algorithmic approaches. In this report we leverage the low-power and the inherently binary operation of magnetic tunnel junctions (MTJs) to demonstrate neural network hardware inference based on passive arrays of MTJs. In general, transferring a trained network model to hardware for inference is confronted by degradation in performance due to device-to-device variations, write errors, parasitic resistance, and nonidealities in the substrate. To quantify the effect of these hardware realities, we benchmark 300 unique weight matrix solutions of a two-layer perceptron to classify the Wine dataset for both classification accuracy and write fidelity. Despite device imperfections, we achieve software-equivalent accuracy of up to 95.3% with proper tuning of network parameters in 15 x 15 MTJ arrays having a range of device sizes. The success of this tuning process shows that new metrics are needed to characterize the performance and quality of networks reproduced in mixed signal hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multiple-scales model of the shock-cell structure of imperfectly expanded supersonic jets

The present investigation is concerned with the development of an analytical model of the quasi-periodic shock-cell structure of an imperfectly expanded supersonic jet. The investigation represents a part of a program to develop a mathematical theory of broadband shock-associated noise of supersonic jets. Tam and Tanna (1982) have suggested that this type of noise is generated by the weak interaction between the quasi-periodic shock cells and the downstream-propagating large turbulence structures in the mixing layer of the jet. In the model developed in this paper, the effect of turbulence in the mixing layer of the jet is simulated by the addition of turbulent eddy-viscosity terms to the momentum equation. Attention is given to the mean-flow profile and the numerical solution, and a comparison of the numerical results with experimental data.

Tam, C. K. W.↗

Gamma-ray spectrometry analysis methods for radioisotope dissolution and mixing for nuclear forensics applications

We report the production of surrogate nuclear explosive debris involves the dissolution of irradiated fissile material, elemental isolation, and recombination based on required fission product ratios. Gamma-ray spectrometry is used following each step to assess the radioactive fission product inventory. A high-purity germanium detector is used to measure each sample, and the detector efficiencies are determined prior to sample measurements. Through radiochemical techniques, elements are isolated and recombined into samples with specific fission product ratios. As a result, a matrix-based approach can account for imperfect chemical separations to calculate the proper mixture volumes to obtain intended volatile-to-refractory fission product ratios, or R-values.

07 ISOTOPE AND RADIATION SOURCES↗

Towards a Regolith Maturity Index for Howardites

The Dawn spacecraft has just arrived at asteroid 4 Vesta, parent of the howardite, eucrite and diogenite (HED) meteorites [1], to begin a yearlong surface study from orbit [2]. As Dawn will view a debris-covered surface, understanding the formation and mixing processes for the debris layer will strongly aid surface data interpretations. Howardites are polymict breccias mainly composed of clasts derived from basaltic (eucritic) and orthopy-roxenitic (diogenitic) parent materials [3]. Some howardites are poorly reworked (fragmental howardites) whilst others have been extensively gardened in an active regolith (regolithic howardites) [4]. The latter may represent an ancient, well-mixed regolith, whilst the former may be from more recent ejecta deposits [4]. Due to environmental differences, regolith development on Vesta differs in detail from that on the Moon [4-6]. We have been developing petrological criteria to apply to howardite thin sections to determine their relative regolithic maturity, which we are fine-tuning with comparison to noble gas data [7, 8]. Whilst we previously emphasized the abundance of reworked clasts (fragmental and impact-melt breccia clasts), this is an imperfect criterion: one howardite with abundant re-worked clasts (EET 99408) shows no evidence of solar wind Ne (SW-Ne), yet, two of our alleged fragmental howardites have clear SW-Ne signatures (LEW 85313, MET 00423) [7, 8]. We are now investigating the diversity in minor and trace element contents of low-Ca pyroxene clasts in howardites as a measure of regolith grade, and will begin analyses of such grains within reworked clasts. Our hypothesis is that regolithic howardites (and the breccia clasts they contain) will show greater diversity because they sampled more diverse diogenitic plutons than fragmental howardites, which formed from ejecta from only a few impacts [e.g. 4]. Our initial LA-ICP-MS work showed ranges in trace element diversity in low-Ca pyroxenes (estimated from the standard error of the mean of analyses), where those howardites considered of medium to high regolithic grade showed greater diversity [9]. Our EMPA results (from a larger howardite suite) show an overall greater diversity in our putative medium to high regolithic grade howardites, though there are exceptions. The greatest diversity is found for paired howardites GRO 95574 and GRO 95581, which were not considered regolithic in our initial study. We will continue investigating avenues to determine regolith maturity in thin section, factoring in bulk rock compositional data, and will coordinate these studies with noble gas results.

Mittlefehldt, David W.↗

Measurements of the Flow in the Vicinity of an Additively Manufactured Turbine Leading-Edge Using X-Ray Particle Tracking Velocimetry

X-ray particle tracking velocimetry (PTV) is performed, for the first time, to measure the velocity field inside a leading-edge of a turbine blade made by laser-additive-manufacturing (LAM) process. The traditional showerhead holes were replaced by a porous matrix in the leading-edge. The flow through such a leading-edge piece cannot be faithfully recreated by traditional prototype testing methods due to the surface roughness and imperfections caused by LAM process. Hence, direct measurement is the only option. However, it is difficult to measure flow inside such pieces with traditional velocimetry measurements due to the existence of metallic walls. Moreover, small internal size and high flow speeds call for a measurement technique with high spatial and temporal resolutions. To address these issues, we performed time-resolved X-ray PTV using the Advanced Photon Source (APS) synchrotron facility at the Argonne National Laboratory (ANL). A hydraulic system was constructed to run water, mixed with seeding particles, through the leading-edge piece. A high-speed camera captured the images of the seeding particles, which were later processed to create particle tracks. The time-averaged velocity field showed distinct pairs of vortices located in front of the porous outlet inside the leading-edge piece. The inlet channel showed reversed flow due to partial obstruction by the porous inlet of the test piece. Such knowledge of the flow field inside a leading-edge of a turbine blade will help us to design better cooling paths leading to higher cooling efficiency and increased life-span of a turbine blade.

42 ENGINEERING↗

Comparison of wind-farm control strategies under realistic offshore wind conditions: wake quantities of interest

Wind-farm control strategies aim to increase the efficiency, and therefore lower the levelized cost of energy, of wind farms. This is done by using turbine settings such as the yaw angle, blade pitch angles, or generator torque to manipulate the wake that negatively affects downstream turbines in the farm. Two inherently different wind-farm control methods have been identified in the literature: wake steering (WS) and active wake mixing (AWM). As one of two companion papers focused on understanding practical aspects of these two wind-farm control strategies using large-eddy simulation (LES), we below analyze the wake quantities of interest for a single wind turbine performing WS and AWM, while the companion article (Frederik et al., 2025) focuses on turbine quantities of interest including power and structural loads for the same computational setup and also includes two-turbine arrays with full and partial wake overlap. The simulations, which are based in the LES solver AMR-Wind, are tailored to have inflow conditions representative of measurements from a site off the East Coast of the US, including with strong veer and low turbulence. The turbine, which is modeled in OpenFAST and coupled to the LES, is the IEA 15 MW, an open-source offshore design. After presenting an overview of the wake recovery for the different wake-control cases, the analysis probes the fluid-dynamic causes for the different performance of the arrays reported in the companion article by examining control volumes around the wakes and the budget of the mean-flow kinetic energy (MKE) within these volumes. In the high veer environment considered, the MKE recovery is dominated by mean convection, and this is shown to especially benefit the WS strategy when a neighboring turbine is directly downstream: there is ≈65 % more available power for a downstream turbine than in the baseline case, and this power is gained primarily through mean convection on the left-tip and top-tip faces of the control volume. However, the case with imperfect knowledge of the exact wind direction favors the pulse-type AWM strategy, largely because of ≈9 % increased turbulent entrainment from aloft versus the baseline that could be related to an apparent resistance to skewing in the pulsed wake. The general reduced effectiveness of helix-type and other individual-pitch-based AWM strategies for inflow with high veer and low turbulence as reported in the companion paper is due, in part, to low magnitudes of phase-averaged turbulent entrainment. Two main findings of this study are thus that veer has a significant impact on the effectiveness of different wake-control strategies and that pulse-type AWM may be a useful strategy when the objective is power maximization in realistic, offshore flow environments with imperfect knowledge of the exact wake overlap position on the downstream turbine.

17 WIND ENERGY↗

Multi-Stage Modeling With Recourse Decisions for Solving Stochastic Complementarity Problems With an Application in Energy

This paper presents a multi-stage model with recourse decisions for solving complementarity problems in a competitive electricity market under uncertainty, while also considering renewable energy technologies and battery storage utilization. The model is based on a Nash-Cournot formulation of imperfect competition among power producers. We analyze the value of variable renewable energy (VRE) and battery storage under different uncertainties, such as demand level and VRE availability. To illustrate the proposed model, we apply it to three- bus five-player model and analyze different cases varying costs, including a user-optimal perspective (with market power) and a system-optimal perspective (with central planning). We also consider the potential for congestion in the system by restricting the transmission capacity between a single interface that connects two buses. Our findings show that increasing the battery storage capacity results in a decrease in the need for perfect information about future uncertainties. Additionally, as the model allows for more uncertainty, it becomes more apparent that the stochastic mixed complementarity problem (MCP) has an advantage over a deterministic equivalent. We propose the use of the Value of the Stochastic Equilibrium Solution (VSES) as a quality metric to compare the stochastic MCP with its deterministic equivalent. Overall, expanding battery storage capacity can lower the maximum, mean, and variance values of delivered prices, but there are diminishing returns to this approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, John P.↗

Surface flaw reliability analysis of ceramic components with the SCARE finite element postprocessor program

The SCARE (Structural Ceramics Analysis and Reliability Evaluation) computer program on statistical fast fracture reliability analysis with quadratic elements for volume distributed imperfections is enhanced to include the use of linear finite elements and the capability of designing against concurrent surface flaw induced ceramic component failure. The SCARE code is presently coupled as a postprocessor to the MSC/NASTRAN general purpose, finite element analysis program. The improved version now includes the Weibull and Batdorf statistical failure theories for both surface and volume flaw based reliability analysis. The program uses the two-parameter Weibull fracture strength cumulative failure probability distribution model with the principle of independent action for poly-axial stress states, and Batdorf's shear-sensitive as well as shear-insensitive statistical theories. The shear-sensitive surface crack configurations include the Griffith crack and Griffith notch geometries, using the total critical coplanar strain energy release rate criterion to predict mixed-mode fracture. Weibull material parameters based on both surface and volume flaw induced fracture can also be calculated from modulus of rupture bar tests, using the least squares method with known specimen geometry and grouped fracture data. The statistical fast fracture theories for surface flaw induced failure, along with selected input and output formats and options, are summarized. An example problem to demonstrate various features of the program is included.

Gyekenyesi, J. P.↗

Linear Quantum Coupler for Clean Bosonic Control

Quantum computing with superconducting circuits relies on high-fidelity driven nonlinear processes. An ideal quantum nonlinearity would selectively activate desired coherent processes at high strength, without activating parasitic mixing products or introducing additional decoherence. The wide bandwidth of the Josephson nonlinearity makes this difficult, with undesired drive-induced transitions and decoherence limiting qubit readout, gates, couplers, and amplifiers. Significant strides have recently been made into building better “quantum mixers,” with promise being shown by Kerr-free three-wave mixers that suppress driven frequency shifts, and balanced quantum mixers that explicitly forbid a significant fraction of parasitic processes. We propose a novel mixer that combines both these strengths, with engineered selection rules that make it essentially linear (not just Kerr-free) when idle, and activate clean parametric processes even when driven at high strength. Furthermore, its ideal Hamiltonian is simple to analyze analytically, and we show that this ideal behavior is first-order insensitive to dominant experimental imperfections. We expect this mixer to allow significant advances in high-𝑄 control, readout, and amplification.

Quantum circuits↗

Improving the Optical Quality Factor of the WGM Resonator

Resonators usually are characterized with two partially dependent values: finesse (F) and quality factor (Q). The finesse of an empty Fabry-Perot (FP) resonator is defined solely by the quality of its mirrors and is calculated as F=piR(exp 1/2)/(1-R). The maximum up-to-date value of reflectivity R approximately equal to 1 - 1.6 x 10(exp -6) is achieved with dielectric mirrors. An FP resonator made with the mirrors has finesse F=1.9 x 10(exp 6). Further practical increase of the finesse of FP resonators is problematic because of the absorption and the scattering of light in the mirror material through fundamental limit on the reflection losses given by the internal material losses and by thermodynamic density fluctuations on the order of parts in 109. The quality factor of a resonator depends on both its finesse and its geometrical size. A one-dimensional FP resonator has Q=2 F L/lambda, where L is the distance between the mirrors and lambda is the wavelength. It is easy to see that the quality factor of the resonator is unlimited because L is unlimited. F and Q are equally important. In some cases, finesse is technically more valuable than the quality factor. For instance, buildup of the optical power inside the resonator, as well as the Purcell factor, is proportional to finesse. Sometimes, however, the quality factor is more valuable. For example, inverse threshold power of intracavity hyperparametric oscillation is proportional to Q(exp 2) and efficiency of parametric frequency mixing is proportional to Q(exp 3). Therefore, it is important to know both the maximally achievable finesse and quality factor values of a resonator. Whispering gallery mode (WGM) resonators are capable of achieving larger finesse compared to FP resonators. For instance, fused silica resonators with finesse 2.3 x 10(exp 6) and 2.8 x 10(exp 6) have been demonstrated. Crystalline WGM resonators reveal even larger finesse values, F=6.3 x 10(exp 6), because of low attenuation of light in the transparent optical crystals. The larger values of F and Q result in the enhancement of various nonlinear processes. Low-threshold Raman lasing, optomechanical oscillations, frequency doubling, and hyperparametric oscillations based on these resonators have been recently demonstrated. Theory predicts a possibility of nearly 10(exp 14) room-temperature optical Q-factors of optical crystalline WGM resonators, which correspond to finesse levels higher than 10(exp 9). Experiments have shown numbers a thousand times lower than that. The difference occurs due to media imperfections. To substantially reduce the optical losses caused by the imperfections, a specific, multi-step, asymptotic processing of the resonator is implemented. The technique has been initially developed to reduce microwave absorption in dielectric resonators. One step of the process consists of mechanical polishing performed after high temperature annealing. Several steps repeat one after another to lead to significant reduction in optical attenuation and, as a result, to the increase of Q-factor as well as finesse of the resonator which demonstrates a CaF2 WGM resonator with F greater than 10(exp 7) and Q greater than 10(exp 11).

Savchenkov, Anatoliy↗

Surface and subsurface oceanographic features drive forage fish distributions and aggregations: Implications for prey availability to top predators in the US Northeast Shelf ecosystem

Abstract Forage fishes are a critical food web link in marine ecosystems, aggregating in a hierarchical patch structure over multiple spatial and temporal scales. Surface‐level forage fish aggregations (FFAs) represent a concentrated source of prey available to surface‐ and shallow‐foraging marine predators. Existing survey and analysis methods are often imperfect for studying forage fishes at scales appropriate to foraging predators, making it difficult to quantify predator–prey interactions. In many cases, general distributions of forage fish species are known; however, these may not represent surface‐level prey availability to predators. Likewise, we lack an understanding of the oceanographic drivers of spatial patterns of prey aggregation and availability or forage fish community patterns. Specifically, we applied Bayesian joint species distribution models to bottom trawl survey data to assess species‐ and community‐level forage fish distribution patterns across the US Northeast Continental Shelf (NES) ecosystem. Aerial digital surveys gathered data on surface FFAs at two project sites within the NES, which we used in a spatially explicit hierarchical Bayesian model to estimate the abundance and size of surface FFAs. We used these models to examine the oceanographic drivers of forage fish distributions and aggregations. Our results suggest that, in the NES, regions of high community species richness are spatially consistent with regions of high surface FFA abundance. Bathymetric depth drove both patterns, while subsurface features, such as mixed layer depth, primarily influenced aggregation behavior and surface features, such as sea surface temperature, sub‐mesoscale eddies, and fronts influenced forage fish diversity. In combination, these models help quantify the availability of forage fishes to marine predators and represent a novel application of spatial models to aerial digital survey data.

54 ENVIRONMENTAL SCIENCES↗

$t_1$-Noise eliminated dipolar heteronuclear multiple-quantum coherence solid-state NMR spectroscopy

Heteronuclear correlation (HETCOR) spectroscopy is one of the key tools in the arsenal of the solid-state NMR spectroscopist to probe chemical and spatial proximities between two different nuclei and enhance spectral resolution. Dipolar heteronuclear multiple-quantum coherence (D-HMQC) is a powerful technique that can be potentially utilized to obtain 1 H detected 2D HETCOR solid-state NMR spectra of any NMR active nucleus. A long-standing problem in 1 H detected D-HMQC solid-state NMR experiments is the presence of t 1 -noise which reduces sensitivity and impedes spectral interpretation. In this contribution, we describe novel pulse sequences, termed t 1 -noise eliminated (TONE) D-HMQC, that minimize t 1 -noise and can provide higher sensitivity and resolution than conventional D-HMQC. Monte-Carlo and numerical simulations confirm that t 1 -noise in conventional D-HMQC primarily occurs because random MAS frequency fluctuations cause variations in the NMR signal amplitude from scan to scan, leading to imperfect cancellation of uncorrelated signals by phase cycling. The TONE D-HMQC sequence uses 1 H π-pulses to refocus the evolution of 1 H CSA across each SR4 2 1 recoupling block, improving the stability of the pulse sequence to random MAS frequency fluctuations. The 1 H refocusing pulses also restore the orthogonality of in-phase and anti-phase magnetization for all crystallite orientations at the end of each recoupling block, enabling the use of 90° flip-back or LG spin-lock trim pulses to reduce the intensity of uncorrelated signals. We demonstrate the application of these methods to acquire 1 H detected 2D 1 H{ 35 Cl} and 1 H{ 13 C} HETCOR spectra of histidine·HCl·H 2 O with reduced t 1 -noise. To show generality, we also apply these methods to obtain 2D 1 H{ 17 O} spectra of 20%- 17 O fmoc-alanine and for the first time at natural abundance, 2D 1 H{ 25 Mg} HETCOR spectra of magnesium hydroxide. The TONE D-HMQC sequences are also used to probe 1 H– 25 Mg and 1 H– 27 Al proximities in Mg–Al layered double hydroxides and confirm the even mixing of Mg and Al in these materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Global Ozone Climatology from Ozone Soundings via Trajectory Mapping: A Stratospheric Perspective

This study explores a domain-filling trajectory approach to generate a global ozone climatology from sparse ozonesonde data. Global ozone soundings of 51,898 profiles at 116 stations over 44 years (1965-2008) are used, from which forward and backward trajectories are performed for 4 days, driven by a set of meteorological reanalysis data. Ozone mixing ratios of each sounding from the surface to 26 km altitude are assigned to the entire path along the trajectory. The resulting global ozone climatology is archived monthly for five decades from the 1960s to the 2000s with grids of 5 degree 5 degree 1 km (latitude, longitude, and altitude). It is also archived yearly from 1965 to 2008. This climatology is validated at 20 ozonesonde stations by comparing the actual ozone sounding profile with that found through the trajectories, using the ozone soundings at all the stations except one being tested. The two sets of profiles are in good agreement, both individually with correlation coefficients between 0.975 and 0.998 and root mean square (RMS) differences of 87 to 482 ppbv, and overall with a correlation coefficient of 0.991 and an RMS of 224 ppbv. The ozone climatology is also compared with two sets of satellite data, from the Satellite Aerosol and Gas Experiment (SAGE) and the Optical Spectrography and InfraRed Imager System (OSIRIS). Overall, the ozone climatology compares well with SAGE and OSIRIS data by both seasonal and zonal means. The mean difference is generally under 20 above 15 km. The comparison is better in the northern hemisphere, where there are more ozonesonde stations, than in the southern hemisphere; it is also better in the middle and high latitudes than in the tropics, where assimilated winds are imperfect in some regions. This ozone climatology can capture known features in the stratosphere, as well as seasonal and decadal variations of these features. Furthermore, it provides a wealth of detail about longitudinal variations in the stratosphere such as the spring ozone maximum over the Canadian Arctic. It also covers higher latitudes than current satellite data. The climatology shows clearly the depletion of ozone from the 1970s to the mid 1990s and ozone recovery in the 2000s. When this climatology is used as the upper boundary condition in an Environment Canada operational chemical forecast model, the forecast is improved in the vicinity of the upper tropospherelower stratosphere region. As this ozone climatology is neither dependent on a priori data or photochemical modeling, it provides independent information and insight that can supplement satellite data and model simulations and enhance our understanding of stratospheric ozone.

SAGE↗

Density variance dynamics in disparate shock tubes

This report discusses two experiments which investigate thin layer, heavy curtain fragmentation from the perspective of a Reynolds-Averaged mix model, in drastically disparate experimental regimes. The first, the centimeter/millisecond-scale “Horizontal Shock Tube” (HST) is a compressed-gas piston-driven shock tube experiment. The second, the “Multishock thin layer” (Mshock) experiment performed at the National Ignition Facility, is a micrometer/nanosecond-scale laser-driven shock tube experiment. Both are situations in which a heavy plane layer (a ‘curtain’) is initially suspended in a lighter medium. After being shocked from at least one side, the layer translates while its interfaces evolve due to the excitation of the Richtmyer-Meshkov instability at its surfaces. The evolution of density variance, which initially exists only on the surface of the layer, as it comes to encompass the whole layer interior is used as a description of layer fragmentation and dissolution. These experiments have each been simulated in the Los Alamos National Laboratory multi-physics code xRAGE, which includes fundamental hydrodynamics, extended plasma physics and radiation effects which are important to drive the high-energy density experiment, and the Besnard-Harlow-Rauenzahn (BHR) turbulence model. In each, the principal diagnostic for comparison is an experimental metric for the density (co)variance, b, which tracks the moments of the density field at the curtain interfaces and body. Due to experimental constraints in different regimes (i.e. optical diagnostics can be deployed on conventional shock tubes, while the plasma shock tubes must be imaged by x-rays; interfaces can be imposed to specification on laser-driven experiments, which are stored in the solid phase, while conventional experiments have imperfect control of the flow fields which separate the layer, etc.) the experiments are not designed to be perfect scaled cognates of one another. However, despite the separation of six orders of magnitude of scaling in time, and four in space, we are able to demonstrate that the same turbulence model, operating in the same fashion in the same computer code, is able to reproduce results in each experiment, by tracking evolution due to common relevant physics. Additionally, we will present preliminary work toward density variance comparisons in a single-interface Richtmyer-Meshkov configuration, the conventional fluid “Vertical Shock Tube” (VST) experiment, and the Modal Initial Conditions (ModCons) campaign fielded at the OMEGA-EP laser facility.

42 ENGINEERING↗