Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mass movements”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hurricanes alter 10 Be concentrations in tropical river sediment but do not change regional erosion rate estimates

Tropical islands, including many in island arcs, are subjected to recurring disturbances from extreme storms such as tropical cyclones. To test whether such storms influence cosmogenic nuclide concentrations such that they do not reflect long-term rates of erosion, we measured meteoric and in situ 10Be in river sediment samples from Dominica, an andesitic island in the Caribbean, before and after category five Hurricane Maria (2017). Populations of before- and after-storm concentrations are statistically indistinguishable (n=7 pairs for in situ 10 Be, n=11 pairs for meteoric 10 Be). 10 Be concentrations vary from -138% to +73% within before-after sample pairs relative to the mean of the pair. These new data suggest that the effects of extreme storms on the depth and amount of near-surface erosion on Dominica vary spatially. Our data support calculations of Niemi, et al. (2005) and Yanites, et al. (2009) suggesting that basin-by-basin comparisons of erosion rates based on cosmogenic nuclides should be approached with caution in small (<~100 km 2 ) watersheds affected by mass movements and extreme storms. Erosion rates determined from in situ 10Be in Dominica (geometric mean=0.102 mm y-1, n=12) are low compared to similarly steep and wet areas globally and correlate positively with the spatial density of mass movements.

58 GEOSCIENCES↗
DOE OSTI2022

Narrow-Band Least-Squares Infrasound Array Processing

Infrasound data from arrays can be used to detect, locate, and quantify a variety of natural and anthropogenic sources from local to remote distances. However, many array processing methods use a single broad frequency range to process the data, which can lead to signals of interest being missed due to the choice of frequency limits or simultaneous clutter sources. In this work, we introduce a new open-source Python code that processes infrasound array data in multiple sequential narrow frequency bands using the least-squares approach. We test our algorithm on a few examples of natural sources (volcanic eruptions, mass movements, and bolides) for a variety of array configurations. Our method reduces the need to choose frequency limits for processing, which may result in missed signals, and it is parallelized to decrease the computational burden. Improvements of our narrow-band least-squares algorithm over broad-band least-squares processing include the ability to distinguish between multiple simultaneous sources if distinct in their frequency content (e.g., microbarom or surf vs. volcanic eruption), the ability to track changes in frequency content of a signal through time, and a decreased need to fine-tune frequency limits for processing. We incorporate a measure of planarity of the wavefield across the array (sigma tau, στ) as well as the ability to utilize the robust least trimmed squares algorithm to improve signal processing and insight into array performance. Our implementation allows for more detailed characterization of infrasound signals recorded at arrays that can improve monitoring and enhance research capabilities.

58 GEOSCIENCES↗
DOE OSTI2023

The age and potential causes of the giant Green Lake Landslide, Fiordland, New Zealand

Abstract Landslide deposits preserved in the geological record afford opportunities to better inform hillslope and seismic hazard and risk models, particularly in regions where observational records are short. In the Southern Alps of New Zealand, small coseismic landslides are frequent, but the geological record preserves several instances of more substantial (> 1 km 3 ) but infrequent mass failures. With an estimated volume of 27 km 3 , the giant Green Lake Landslide represents one of the largest subaerial landslides on Earth. Previous work has suggested this deep-seated mass movement was most likely triggered by high-intensity seismic shaking, but that local structural weakness and/or glacial debuttressing may help to explain the anomalously large failure volume. Resolving the potential contribution of the latter is important given predictions concerning alpine deglaciation in the coming decades to centuries. Key to resolution are secure chronologies of landslide emplacement and past glacier change. Here we present in situ cosmogenic 10 Be exposure ages from the Green Lake Landslide that suggest an emplacement age of 15.5 ± 0.7 ka. Recent work shows that glacial retreat in the region was underway by 19 ka, indicating that the Green Lake Landslide was emplaced 3–4 kyr after the onset of glacier retreat. Given the geometry of the former confining valley glacier, we expect that the deglaciation-landslide age gap is closer to the upper end of this estimate. If correct, this conclusion places greater weight on the roles of local geological structure and/or a great earthquake as factors contributing to the exceptionally large volume of this event.

58 GEOSCIENCES↗
DOE OSTI2022

Effect of grain dissolution on sloping ground

The static and dynamic stability of natural or constructed slopes can be affected by dissolution or dissolution-like phenomena. Their underlying mechanisms, however, remain unclear. New experimental results and discrete element simulations provide particle-level and macroscale information on the consequences of mineral dissolution on slope behavior. At the microscale, load-carrying grain arches develop around dissolving particles, the porosity increases, and contact force chains evolve to form a honeycomb topology. At the macroscale, while vertical settlements are the prevailing deformation pattern, lateral granular movements that create mass wasting are prominent in sloping ground, even under the quasi-static granular loss. Horizontal grain displacement is maximum at the surface and decreases linearly with the distance from the slope surface to become zero at the bottom boundaries, much like vertical granular displacement along the depth. Sediments with smaller friction angles and steeper slopes experience greater displacement, both vertically and horizontally. Slopes become flatter after dissolution, with the reduction in slope angle directly related to the loss in ground elevation, ΔH/H o . Yet, because of the porous fabric that results from dissolution, vertical shortening is less than the upper bound, estimated from the loss in the solid mass fraction, ΔH/H o ≈SF. Under water-saturated conditions, the post-dissolution fabric may lead to sudden undrained shear and slope slide.

42 ENGINEERING↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of a Particle Flow Control Mechanism

The next generation of Concentrated Solar Power (CSP) systems are utilizing solid particles as the heat transfer medium. The state-of-the-art particle flow control mechanism currently utilized is the slide-gate mechanism which functions by linearly actuating a gate across a tear-drop shaped opening. The sensitivity, which is a measure of the change in mass flow rate per unit movement of the slide-gate, is 0.354 g/s-mm[1]. The slide-gate’s main shortcomings include particle lodging and particle leakage between the plates, and difficulty actuating under the pressure from the particles. The goal of a new device is to minimize or eliminate these flaws. The new proposed device functions similar to a chuck mechanism in a drill where a rotation of an outer ring produces a change in outlet area through moving “jaws” diagonally. This mechanism aims to solve the issues of the slide-gate by having the jaws on the outside of the particle flow volume to reduce the probability of particle lodging, and the jaw tips are designed to overlap as they actuate to reduce particle leakage. Additionally, the vertical orientation of the device as well as the inclined movement of the jaws allows for easier actuation of the device under pressure. Experiments involving 3D printed prototypes have been carried out, producing a max mass flowrate of about 50 g/s and a sensitivity curve with the equation 𝑦 (𝑔/𝑠) = 0.0851𝑥^2(𝑔/𝑠−𝑚𝑚2) +3.6993𝑥(𝑔/𝑠−𝑚𝑚). These experiments validated the device’s capability to address the shortcomings of the slide-gate. Currently, an aluminum model of the device is being tested at room temperature, and a stainless-steel model is being developed to be tested at 700oC, which is the expected operating temperature of the CSP system. Furthermore, simulations are on-going to predict the mass flowrate of different device configurations at room and high temperatures.

14 SOLAR ENERGY↗
DOE OSTI2022

Monitoring Carbon Storage Sites with Time-Lapse Gravity Surveys

Monitoring technologies are deployed on carbon storage sites to demonstrate that the injected CO2 is contained within the storage complex during and long after the active life of the project and that the CO2 behavior is conforming to expectations. Detection of gravity variations accompanying mass redistributions in the subsurface caused by fluid movements in reservoirs provides a unique means for monitoring the dynamics of a carbon sequestration site. Time-lapse gravity surveys have a long history of effectively monitoring temporal density changes in the subsurface. However, monitoring of carbon storage sites with time-lapse gravity surveys using surface or borehole measurements remains limited, although it gained more attention during the last decade. With the recent progresses in instrument capabilities, and the increased number of fluid storage sites (e.g., wastewater, CO2, natural gas, etc.), the gravity method (with its unique ability to remotely estimate both the mass of an injected fluid and its location) has a promising future in monitoring. This chapter provides background information regarding the gravity method and discusses its modeling and application to various carbon storage sites.

Carbon Sequestration, time-lapse gravity, monitori↗
DOE OSTI2026

A worldwide climatology of extreme air masses

Extreme temperature events are among the most damaging weather phenomena. In a warming world, more heat extremes and fewer cold extremes are expected in most regions in the future, a trade-off that warrants further understanding of such events. Here, we track and analyze large, persistent areas of hot and cold extreme temperatures in parallel, relative to the location and time of year, to quantify overall regional exposure to extreme temperatures. To accomplish this, we compare the frequencies, movements, trends, and sources/sinks in each type of extreme air mass, calling them extreme cold or extreme hot air masses (ECAMs/EHAMs). For most land regions, ECAMs occur more often than EHAMs, and ECAMs are more common in each hemisphere’s winter, when cold-air advection is strongest and most widespread. Average movement of ECAMs has a stronger equatorward component in winter than in summer, while movement of EHAMs is eastward all year, with less meridional movement than ECAMs. EHAMs have become more common almost everywhere on land, and the reverse is true for ECAMs, with the strongest trends in the northern hemisphere occurring in autumn, especially in the Arctic. The number of EHAMs around the world is increasing at a higher rate (+ 2.32 events/year) than ECAM numbers are decreasing (-1.27 events/year), causing a net increase in extremes, especially in some midlatitude regions. These results help to document the different processes driving hot and cold extremes, and thus the asymmetric and regionally varying trends in frequency for extreme air masses.

Ryan, James M↗
DOE OSTI2023

Particulate mass migration and mixing in cylindrically contained explosions

Abstract To explore particulate movement near the plasma of chemical explosions, rugged tracer particles were placed within and on the exterior of metal charges and electrically detonated. The particles were collected on/in the porous walls of plastic cylinders at diameters that correlated to the plasma width during different phases of the explosion. The particles’ positions were determined by Boolean logic analysis of their luminescent intensity. The cylinders which caught particles from the initial phases of the explosion retained placement information, while wider cylinders showed uniform mixing. These results/analysis methodology can help improve the understanding of particulate mixing in harsh environments. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2026

In situ quantification of fracture slip induced by hydraulic injections in a deep borehole: A comparison of two different borehole techniques

In situ measurements of fracture deformation during fluid injection are rare, yet essential for understanding the mechanical response of fractured rock. In this study, we evaluate the reliability of two methods by comparing their slip vector estimates: high-resolution borehole acoustic televiewer images captured before and after injection tests, and displacement data from a three-component borehole deformation probe recorded during the injections. Acoustic televiewer images capture only final in-plane displacement, whereas three-component borehole deformation measurements provide full 3D, transient fracture movement. Four injection tests in a fractured granitic rock mass along an inclined borehole at the Bedretto Underground Laboratory (Switzerland), beneath more than 1100 m of overburden, were analyzed. The two methods yielded consistent kinematics and comparable slip magnitudes, typically in the range of 0.2 – 0.6 mm. Angular differences between estimated slip directions ranged from 0.2° to 28°. For two test intervals with pre-existing fractures optimally oriented for slip, the measured slip vectors agreed with forward-modeled slip vectors based on the regional in situ stress field, supporting the Wallace–Bott hypothesis. Measurable shear displacement was identified in two test intervals containing fractures oriented nearly perpendicular to the least principal stress, including a newly formed hydraulic fracture. This demonstrates that even fractures unfavorably oriented for slip can undergo shear due to local stress perturbations around an inclined borehole.

Acoustic televiewer↗
DOE OSTI2026

Queen bees offload pesticide burden to eggs when social buffering is overwhelmed

Honey bee colonies pollinate about one-third of the world’s food crops, and their rapid decline directly threatens agricultural productivity and ecosystem stability. Understanding how colony-level social defenses influence pesticide fate and the circumstances under which they fail is therefore a crucial question in pollinator biology. We used biological accelerator mass spectrometry (BioAMS), a sensitive radiotracer technique, to track the movement of a model pesticide through a small honey bee colony under laboratory conditions. We tested the hypothesis that social buffering protects honey bees from toxic accumulation and that this protection can be overcome, leading to maternal offloading of the pesticide to developing eggs. Consistent with this hypothesis, our results identified three key mechanisms governing chemical movement within a social insect colony: (1) worker bees initially decrease dietary pesticide levels by 95% through diet filtering and deposition in honeycombs, though this declines to 86% by day 10; (2) queen bees maintain markedly lower pesticide levels than workers but, over time, they accumulate the pesticide in their ovaries and transfer it into developing eggs, revealing a previously undocumented protective mechanism in reproductive individuals; and (3) the presence of a queen bee shifts colony-wide chemical distribution by concentrating worker exposure and increasing pesticide deposition in wax. Our findings show that honey bee colonies function as integrated detoxification networks, in which chemical fate depends on complex social behaviors and caste-specific physiology. When social buffering is overwhelmed, reproductive queens may survive by transferring their chemical burden to their offspring.

Biological and medical sciences↗
DOE OSTI2023

CHALLENGES IN THE DEVELOPMENT OF THE ELECTROREFINING PROCESS AT Y-12

In order to ensure future capabilities of the Y-12 site as older buildings retire, a new process was developed, the Metal Purification Process. The purpose of this new process is to provide a simpler and more efficient uranium metal purification that the previous multi-staged, complex chemical processes. The basis for this technology currently exists in the US DOE complex, but it has not been utilized on a large scale for uranium at higher enrichments. The system requires larger-geometry vessels and furnaces in order to meet the through-put requirements. The crux of the process is the electrorefining cell, consisting of several concentric liners and an unfavorable geometry cylindrical crucible. Within the crucible, a molten Li-CL-KCL-UCL3 salt electrolyte is present. ‘Dirty’ metal is placed into a loading basket and loaded into an electrorefining cell. Here, the metal reacts into the molten salt within the cell to form UCL3, while ‘clean’ metal simultaneously plates out from the salt electrolyte in the form of metal dendrites. These dendrites are then collected and removed from the cell. After removal from the cell, the dendrites are taken to a furnace to remove adhered salts. Then, the salt-free dendrites are moved to a separate furnace for consolidation into the final product form. The development of the Criticality Safety Evaluation for this complex system brought with it many NCS Challenges and lessons learned. These challenges include: design decisions for ensuring subcriticalty in the electrorefining cell during abnormal conditions, the interface of the main glovebox system with auxiliary systems such as the purification system and designated storage, implementation of mass tracking, and concerns from production and operation regarding the movement of material within the system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The 300 Marines: characterizing the US Marines with perfect scores on their physical and combat fitness tests

Few US Marines earn perfect 300 scores on both their Physical Fitness Test (PFT) and Combat Fitness Test (CFT). The number 300 invokes the legendary 300 Spartans that fought at the Battle of Thermopylae, which inspired high physical fitness capabilities for elite ground forces ever since. Purpose: Determine distinguishing characteristics of the “300 Marines” (perfect PFT and CFT scores) that may provide insights into the physical and physiological requirements associated with this capability. These tests have been refined over time to reflect physical capabilities associated with Marine Corps basic rifleman performance. Materials and methods: Data were analyzed from US Marines, including 497 women (age, 29 ± 7 years; height 1.63 ± 0.07 m; body mass, 67.4 ± 8.4 kg) and 1,224 men (30 ± 8 years; 1.77 ± 0.07 m; 86.1 ± 11.1 kg). Marines were grouped by whether they earned perfect 300 scores on both the PFT and CFT (300 Marines) or not. We analyzed group differences in individual fitness test events and body composition (dual-energy x-ray absorptiometry). Results: Only 2.5% (n = 43) of this sample earned perfect PFT and CFT scores (n = 21 women; n = 22 men). Compared to sex-matched peers, 300 Marines performed more pull-ups, with faster three-mile run, maneuver-under-fire, and movement-to-contact times (each p < 0.001); 300 Marines of both sexes had lower fat mass, body mass index, and percent body fat (each p < 0.001). The lower percent body fat was explained by greater lean mass (p = 0.041) but similar body mass (p = 0.085) in women, whereas men had similar lean mass (p = 0.618), but lower total body mass (p = 0.025). Conclusion: Marines earning perfect PFT and CFT scores are most distinguished from their peers by their maneuverability, suggesting speed and agility capabilities. While both sexes had considerably lower percent body fat than their peers, 300 Marine women were relatively more muscular while men were lighter.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Large eddy simulation of particle dispersion analysis in a pressurized reactive fluid–particle system

Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new combustion technologies. Because of high pressure, the volumetric fraction of particles in a pressurized oxy-combustor could be higher than that in the conventional combustor, consequently, the particle dynamics in the pressurized vessel need to be evaluated accordingly. Specifically, the present computational work analyzes the pattern of the particle dispersion by extracting the particle mass flow rate and the particle size distribution in the pressurized combustor. A large eddy simulation (LES) with the environment pressure 15-bar and 100 kWth thermal input of fuel particle has been performed. The particle mass flow rates sampled in several cross-section planes strongly fluctuate, indicating the impacts of the flame as well as the gas-phase flow on the particle movement. Animations of the particle-gas flow show that the flow recirculation plays a major role in producing the pulsing phenomenon of the particle mass flow rate.

Li, Lei↗
DOE OSTI2025

Atmospheric wind energization of ocean weather

Ocean weather comprises vortical and straining mesoscale motions, which play fundamentally different roles in the ocean circulation and climate system. Vorticity determines the movement of major ocean currents and gyres. Strain contributes to frontogenesis and the deformation of water masses, driving much of the mixing and vertical transport in the upper ocean. While recent studies have shown that interactions with the atmosphere damp the ocean’s mesoscale vortices O(100) km in size, the effect of winds on straining motions remains unexplored. Here, we derive a theory for wind work on the ocean’s vorticity and strain. Using satellite and model data, we discover that wind damps strain and vorticity at an equal rate globally, and unveil striking asymmetries based on their polarity. Subtropical winds damp oceanic cyclones and energize anticyclones outside strong current regions, while subpolar winds have the opposite effect. A similar pattern emerges for oceanic strain, where subtropical convergent flow is damped along the west-equatorward east-poleward direction and energized along the east-equatorward west-poleward direction. These findings reveal energy pathways through which the atmosphere shapes ocean weather.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Center-of-Mass Corrections in Associated Particle Imaging

Associated particle imaging (API) utilizes the inelastic scattering of neutrons produced in deuterium-tritium (DT) fusion reactions to obtain 3-D isotopic distributions within an object. The locations of the inelastic scattering centers are calculated by measuring the arrival time and position of the associated alpha particle produced in the fusion reactions, and the arrival time of the prompt gamma created in the neutron scattering event. While the neutron and its associated particle move in opposite directions in the center-of-mass (COM) system, in the laboratory system the angle is slightly less than 180°, and the COM movement must be taken into account in the reconstruction of the scattering location. Furthermore, the fusion reactions are produced by ions of different momenta, and thus the COM velocity varies, resulting in an uncertainty in the reconstructed positions. In this article, we analyze the COM corrections to this reconstruction by simulating the energy loss of beam ions in the target material and identifying sources of uncertainty in these corrections. Here we show that an average COM velocity calculated using the ion beam direction and energy can be used in the reconstruction and discuss errors as a function of ion beam energy, composition, and alpha detection location. When accounting for the COM effect, the mean of the reconstructed locations can be considered a correctable systematic error leading to a shift/tilt in the reconstruction. However, the distribution of reconstructed locations also have a spread that will introduce an error in the reconstruction that cannot be corrected. In this article, we will use the known stopping powers of ions in materials and reaction cross sections to examine the reconstruction uncertainties. We also discuss the impact of this effect on our API system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Promoting electrochemical rates by concurrent ionic-electronic conductivity enhancement in high mass loading cathode electrode

Enhancing the fast charging capacity of thick electrodes with high mass loading is imperative in expediting the widespread adoption of electric vehicles. Nonetheless, the insufficient charge transfer kinetics of thick electrodes hinder the movement of effective electrons and ions, hence diminishing capacity at high current rates. In this work, we applied sustainable and biodegradable cellulose nanocrystals (CNCs) as electrode additives. It is the first time to simultaneously improve the electronic conductivity by optimizing the carbon dispersion and establishing electron transfer networks, as well as boosting the ionic conductivity of electrodes by shortening the ion transfer pathway. Specifically, the LiNi 0.6 Mn 0.2 Co 0.2 O 2 electrodes incorporating 1% dual functional CNCs additive exhibit improved effective electrical conductivity from 0.11 to 0.16 S/m and risen effective ionic conductivity from 0.36 to 0.62 S/m, in comparison to counterpart electrodes without CNCs. Therefore, the 1% CNC electrode with a high mass loading of 27.0 mg/cm 2 delivers a discharge capacity of 128 mAh/g at 1 C, which is superior to that of the CNC-free electrodes (95 mAh/g). In short, this study presents a novel environmentally friendly, economically viable, and dual-functional electrode additive that enhances both electronic and ionic conductivities with the aim of facilitating the widespread adoption of fast-charging high mass loading electrodes.

25 ENERGY STORAGE↗
DOE OSTI2023

Investigating Inherent Numerical Stabilization for the Moist, Compressible, Non‐Hydrostatic Euler Equations on Collocated Grids

This study investigates inherent numerical dissipation due to upwind fluxes and reconstruction strategies for collocated Finite-Volume integration of the Euler equations. Idealized supercell simulations are used without any explicit dissipation. Flux terms are split into: mass flux, pressure, and advected quantities. They are computed with the following upwind strategies: central, advectively upwind, and acoustically upwind. This is performed for third and ninth-order-accurate reconstructions with and without Weighted Essentially Non-Oscillatory limiting. Acoustic-only upwinding for pressure and mass flux terms and advective-only upwinding for advected quantities is the most flexible simplification found. It reduces data movement and computations. Assuming a constant speed of sound in acoustic upwinding gives similar results to using the true speed of sound. Dissipation from upwind adapts automatically to grid spacing, time step, reconstruction accuracy, and flow smoothness. While stability is maintained even at 21st-order spatial accuracy, there is a limit to the spatial order of accuracy for which upwinding alone can create a realizable solution in the conditions of this study. Convex combinations of upwind and central solutions for flux terms also reduced dissipation, but as the central proportion grows, solutions become physically unrealizable. The range of length scales of the kinetic energy spectra can be extended along k −5/3 to smaller spatial scales by reducing dissipation either with higher-order reconstructions or using convex combinations of upwind and central fluxes. However, not all extensions of the length scale range along k −5/3 exhibit physically realizable solutions, even though the spectra appear to be physical.

54 ENVIRONMENTAL SCIENCES↗