Search NASA⌕ Search

SEARCH · Search NASA

Results for “High velocity change”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Acoustic and optical phonon frequencies and acoustic phonon velocities in Si-doped AlN thin films

We report the results of the study of the acoustic and optical phonons in Si-doped AlN thin films grown by metal–organic chemical vapor deposition on sapphire substrates. The Brillouin–Mandelstam and Raman light scattering spectroscopies were used to measure the acoustic and optical phonon frequencies close to the Brillouin zone center. The optical phonon frequencies reveal non-monotonic changes, reflective of the variations in the thin film strain and dislocation densities with the addition of Si dopant atoms. The acoustic phonon velocity decreases monotonically with increasing Si dopant concentration, reducing by ∼300 m/s at the doping level of 3 × 1019 cm−3. The knowledge of the acoustic phonon velocities can be used for the optimization of the ultra-wide bandgap semiconductor heterostructures and for minimizing the thermal boundary resistance of high-power devices.

Physics↗
DOE OSTI2021

Suppression of electroconvective and morphological instabilities by an imposed cross flow of the electrolyte

Electroconvective (EC) instability and its influence on surface morphological perturbations (Morph) are important in many applications, including electrodialysis, batteries, and fuel cells. In this work, we study the effects of a two-dimensional channel flow on the EC and Morph instabilities using two approaches. In the bulk analysis, we derive the asymptotic solutions for small and large wave numbers by neglecting the space charge layer and imposing a second kind electroosmosis slip velocity boundary condition on the electroneutral bulk region. In the full analysis, the instability of the entire region of the liquid electrolyte is analyzed using the ultraspherical spectral method. Both studies show that the flow significantly affects the EC instability. The imposed flow distorts the concentration field, causing a sheltering effect which hinders the ion transport from low- to high-concentration regions, and therefore suppresses the EC instability below a certain wave number. In combination with the viscous stabilization of the high wave number modes by the space charge layer, a sufficiently strong imposed flow can fully suppress the EC instability. The increment in the critical voltage for the instability onset is roughly proportional to the square root of the product of imposed velocity and double layer thickness. The imposed flow has a smaller effect on the Morph instability, except that it may remove the morphological modes resulting from the EC instability and thereby change the wave number of the most unstable mode.

42 ENGINEERING↗
DOE OSTI2020

Scaling Law for Cloud-Rise Velocity vs. Scaled Height of Burst (SHOB)

The shockwave produced by a nuclear detonation leaves a low-density region in the vicinity of air zero (i.e., the point of detonation). As the shockwave degenerates into an acoustic wave, the buoyant force takes over and causes the low-density region to begin to rise at speeds that can reach as high as ~200-300 mph. The objective of this project was to find out if this initial cloud rise followed any scaling law, by stimulating a cloud rise for a 1 kiloton (kt) detonation and determining if a correlation existed between the scaled-height-of-burst (SHOB) and the scaled-cloud-rise velocity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Signal Differentiation of Moving Magnetic Nanoparticles for Enhanced Biodetection and Diagnostics

Magnetic nanoparticles are extensively utilized as markers/signal labelling in various biomedical applications. Detecting and distinguishing magnetic signals from similarly sized moving magnetic nanoparticles in microfluidic systems is crucial yet challenging for biosensing. In this study, we have developed an original method to detect and differentiate magnetic signals from moving superparamagnetic (SPM) and ferrimagnetic (FM) nanoparticles of comparable sizes. Our approach utilizes a highly sensitive magnetic-coil-based sensor that harnesses the combined effects of giant magnetoimpedance (GMI) and an LC-resonance circuit, offering performance superior to that of conventional GMI sensors. Iron oxide nanoparticles, which have similar particle sizes but differing coercivities (zero for SPM and non-zero for FM) or similar zero coercivities but differing particle sizes, flow through the magnetic coil at controlled velocities. Their distinct effects are analyzed through changes in the complex impedance of the sensing system. Our findings provide a unique pathway for utilizing SPM and FM nanoparticles as innovative magnetic markers to identify specific biological entities, thereby expanding their potential applications.

Hwang, Kee Young↗
DOE OSTI2022

An improved pressure drop correlation for modeling localized effects in a pebble bed reactor

Advances in the development of pebble bed reactors (PBRs) has created a desire for accurate and cost-effective simulation tools for design scoping studies and safety analysis. The current state-of-the-art for these simulations is the use of porous media models, although these models rely on correlations to capture the effects of flow features that are not explicitly modeled. One of the areas where correlation accuracy is currently lacking is in the near-wall region of the bed. In this region, the presence of the wall causes the pebbles to pack more orderly, drastically changing the geometry and flow behavior in this region. This work presents a new generalized pressure drop correlation for PBRs based on the KTA equation. A high-to-low methodology is applied, where large eddy simulation (LES) is performed on two beds of 1568 and 1700 pebbles to generate a high-fidelity dataset. The flow fields are then averaged in time and separated into concentric rings of 0.05 D peb width. Average porosity, velocity, and pressure drop are extracted for each ring and the friction and form losses are calculated. The Reynolds number range for this study is 625–10,000, and thus the form losses are dominant over the friction losses. The form losses across the rings are investigated, and a correction term for the form loss calculation is determined and applied to the KTA equation to drastically improve the capability of modeling localized porosity effects in a porous media code. Finally, the improved correlation reduces near-wall velocity prediction error from over 50% with the KTA correlation to around 5%. Agreement in pressure drop prediction between LES and porous media simulations is also improved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Quantifying Turbine-Level Risk to Golden Eagles Using a High-Fidelity Updraft Model and a Stochastic Behavioral Model

To minimize the effects of wind farms on Golden Eagle (Aquila chrysaetos) populations while enabling sustainable development of renewable energy resources, it is important to understand how eagles interact with atmospheric flows, terrain features, and anthropogenic structures. Models that predict migratory flight paths provide one tool that helps us grasp how the location of wind farms may influence interactions and impacts on migrating Golden Eagles. The current state-of-the-art in predicting migratory flight paths uses a deterministic fluid-flow analogy to predict eagle trajectory using only an orographic updraft potential computed from topographical features. This model does not take into account variables, such as thermal updrafts and time varying atmospheric conditions that are known to influence migratory behavior. In this work, we improve on the model with the objective of developing tools that advance our understanding of how atmospheric flows and terrain features affect migratory eagle behavior and their interactions with wind farms. Specifically, we 1) incorporate both orographic and thermal updraft information in simulating eagle flight paths; 2) incorporate stochasticity into eagle travel patterns to better capture the influence of exogenous factors on, and the inherent stochasticity of eagle behavior; 3) consider spatio-temporal atmospheric data at wind-farm-scale when computing updraft potential; and 4) account for how atmospheric conditions and the direction of migration change seasonally and how these changes affect eagle migratory flight behavior. We tested the model using a 50km by 50km region with 50 m resolution in the western United States. We simulated 900 independent, probabilistic eagle tracks during southerly and northerly migration, assuming eagles solely rely on orographic updrafts. The preliminary results indicate that the inclusion of finer resolution atmospheric data allows for the inclusion of realistic conditions that an eagle experiences. The stochasticity in eagle tracks provides a platform to include uncertainty in eagle decision making and help produce robust eagle presence maps. We will deploy updraft and downdraft velocities computed using a high-fidelity, wind farm scale, computational fluid dynamics solver under development at National Renewable Energy Laboratory. This work is a first step in the development of a predictive and generalizable eagle behavior model at the wind farm scale that does not rely on empirical data collection. Although the current model is intended for migratory eagles, we will extend and refine this model to inform the development of additional behavioral modes, including resident eagle behavior. This modeling approach improves our ability to understand eagle use of the landscape at a fine scale, and it is our hope that this work will ultimately help advance strategies that minimize the impact of wind development on Golden Eagle populations.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2021

Resolving the Circumstellar Environment of the Galactic B[e] Supergiant Star MWC 137.II. Nebular Kinematics and Stellar Variability

The Galactic B[e] supergiant MWC 137 is surrounded by a large-scale optical nebula. To shed light on the physical conditions and kinematics of the nebula, we analyze the optical forbidden emission lines [N ii] λ λ 6548,6583 and [S ii] λ λ 6716,6731 in long-slit spectra taken with ALFOSC at the Nordic Optical Telescope. The radial velocities display a complex behavior but, in general, the northern nebular features are predominantly approaching while the southern ones are mostly receding. The electron density shows strong variations across the nebula with values spreading from about zero to ~800 cm{sup -3}. Higher densities are found closer to MWC 137 and in regions of intense emission, whereas in regions with high radial velocities the density decreases significantly. We also observe the entire nebula in the two [S ii] lines with the scanning Fabry–Perot interferometer attached to the 6 m telescope of the Special Astrophysical Observatory. These data reveal a new bow-shaped feature at P.A. = 225°–245° and a distance 80" from MWC 137. A new Hα image has been taken with the Danish 1.54 m telescope on La Silla. No expansion or changes in the nebular morphology appear within 18.1 yr. We derive a mass of 37{sub -5}{sup +9} M{sub ⊙} and an age of 4.7 ± 0.8 Myr for MWC 137. Furthermore, we detect a period of 1.93 d in the time series photometry collected with the TESS satellite, which could suggest stellar pulsations. Other, low-frequency variability is seen as well. Whether these signals are caused by internal gravity waves in the early-type star or by variability in the wind and circumstellar matter currently cannot be distinguished.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

A detailed study of the barium central star of the planetary nebula Abell 70

ABSTRACT We present a detailed study of the barium star at the heart of the planetary nebula Abell 70. Time-series photometry obtained over a period of more than 10 yr demonstrates that the barium-contaminated companion is a rapid rotator with temporal variability due to spots. The amplitude and phasing of the photometric variability change abruptly; however, there is no evidence for a change in the rotation period (P = 2.06 d) over the course of the observations. The co-addition of 17 high-resolution spectra obtained with Ultraviolet and Visual Échelle Spectrograph mounted on the Very Large Telescope allows us to measure the physical and chemical properties of the companion, confirming it to be a chromospherically active, late G-type sub-giant with more than +1 dex of barium enhancement. We find no evidence of radial velocity variability in the spectra, obtained over the course of approximately 130 d with a single additional point some 8 yr later, with the radial velocities of all epochs approximately −10 km s −1 from the previously measured systemic velocity of the nebula. This is perhaps indicative that the binary has a relatively long period (P ≳ 2 yr) and high eccentricity (e ≳ 0.3), and that all the observations were taken around radial velocity minimum. However, unless the binary orbital plane is not aligned with the waist of the nebula or the systemic velocity of the binary is not equal to the literature value for the nebula, this would imply an unfeasibly large mass for the nebular progenitor.

Jones, David (ORCID:0000000339475946)↗
DOE OSTI2020

Development of aerosol activation in the double-moment Unified Model and evaluation with CLARIFY measurements

Representing the number and mass of cloud and aerosol particles independently in a climate, weather prediction or air quality model is important in order to simulate aerosol direct and indirect effects on radiation balance. Here we introduce the first configuration of the UK Met Office Unified Model in which both cloud and aerosol particles have “double-moment” representations with prognostic number and mass. The GLObal Model of Aerosol Processes (GLOMAP) aerosol microphysics scheme, already used in the Hadley Centre Global Environmental Model version 3 (HadGEM3) climate configuration, is coupled to the Cloud AeroSol Interacting Microphysics (CASIM) cloud microphysics scheme. We demonstrate the performance of the new configuration in high-resolution simulations of a case study defined from the CLARIFY aircraft campaign in 2017 near Ascension Island in the tropical southern Atlantic. We improve the physical basis of the activation scheme by representing the effect of existing cloud droplets on the activation of new aerosol, and we also discuss the effect of unresolved vertical velocities. We show that neglect of these two competing effects in previous studies led to compensating errors but realistic droplet concentrations. While these changes lead only to a modest improvement in model performance, they reinforce our confidence in the ability of the model microphysics code to simulate the aerosol–cloud microphysical interactions it was designed to represent. Capturing these interactions accurately is critical to simulating aerosol effects on climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Homologous mutations in human β, embryonic, and perinatal muscle myosins have divergent effects on molecular power generation

Mutations at a highly conserved homologous residue in three closely related muscle myosins cause three distinct diseases involving muscle defects: R671C in β-cardiac myosin causes hypertrophic cardiomyopathy, R672C and R672H in embryonic skeletal myosin cause Freeman–Sheldon syndrome, and R674Q in perinatal skeletal myosin causes trismus-pseudocamptodactyly syndrome. It is not known whether their effects at the molecular level are similar to one another or correlate with disease phenotype and severity. To this end, we investigated the effects of the homologous mutations on key factors of molecular power production using recombinantly expressed human β, embryonic, and perinatal myosin subfragment-1. We found large effects in the developmental myosins but minimal effects in β myosin, and magnitude of changes correlated partially with clinical severity. The mutations in the developmental myosins dramatically decreased the step size and load-sensitive actin-detachment rate of single molecules measured by optical tweezers, in addition to decreasing overall enzymatic (ATPase) cycle rate. In contrast, the only measured effect of R671C in β myosin was a larger step size. Our measurements of step size and bound times predicted velocities consistent with those measured in an in vitro motility assay. Finally, molecular dynamics simulations predicted that the arginine to cysteine mutation in embryonic, but not β, myosin may reduce pre-powerstroke lever arm priming and ADP pocket opening, providing a possible structural mechanism consistent with the experimental observations. This paper presents direct comparisons of homologous mutations in several different myosin isoforms, whose divergent functional effects are a testament to myosin’s highly allosteric nature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Thermal management of a high temperature sodium sulphur battery stack

The sodium sulfur battery is an advanced secondary battery with high potential for grid-level storage due to their high energy density, low cost of the reactants, and high open-circuit voltage. However, as the operating temperature of the battery is high (about 300 °C), effective thermal management is required to prevent thermal runaway under high current density operation. Here, to develop efficient thermal management strategies, a detailed, thermo-electrochemical, non-isothermal, distributed dynamic model of a sodium-sulfur battery stack is developed. Models of three thermal management strategies are developed and analyzed in this work: active cooling, passive cooling, and hybrid cooling. The active cooling strategy uses air as the cooling medium whereas, the passive cooling strategy uses a phase change material. In the hybrid cooling strategy, both active and passive cooling strategies are used. Due to the high operating temperature of the sodium sulfur batteries, the rejected heat can be utilized effectively and consequently a high variability in the heat rejection rate may not be acceptable. Performance of these strategies is analyzed under high current density operation by evaluating the variability in the cell temperature during charge/discharge cycles, temperature difference between the cells based on their locations, and variability in the heat rejection rate. The phase change-based hybrid thermal management strategy with an embedded controller maintained the temperature variation within ±/- 2.8 °C. Considering the tradeoff between the capital cost for PCM, fan power requirement, and variability in the heat rejection rate, the optimal quantity of the PCM and the maximum air velocity for the maximum current density of 260 mA/cm 2 (0.55 C-rate) were found be 0.74 g/Wh and 3 m/s, respectively.

25 ENERGY STORAGE↗
DOE OSTI2023

High-resolution simulations of a turbulent boundary layer impacting two obstacles in tandem

High-fidelity large-eddy simulations of the flow around two rectangular obstacles are carried out at a Reynolds number of 10 000 based on the freestream velocity and the obstacle height. The incoming flow is a developed turbulent boundary layer. Mean-velocity components, turbulence fluctuations, and the terms of the turbulent-kinetic-energy budget are analyzed for three flow regimes: skimming flow, wake interference, and isolated roughness. Three regions are identified where the flow undergoes the most significant changes: the first obstacle’s wake, the region in front of the second obstacle, and the region around the second obstacle. In the skimming-flow case, turbulence activity in the cavity between the obstacles is limited and mainly occurs in a small region in front of the second obstacle. In the wake-interference case, there is a strong interaction between the freestream flow that penetrates the cavity and the wake of the first obstacle. Furthermore, this interaction results in more intense turbulent fluctuations between the obstacles. In the isolated-roughness case, the wake of the first obstacle is in good agreement with that of an isolated obstacle. Separation bubbles with strong turbulent fluctuations appear around the second obstacle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. Here, two-level GDSW (Generalized Dryja–Smith–Widlund) type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the MPI-parallel implementation of multi-level Schwarz preconditioners provided by the package FROSch (Fast and Robust Schwarz)from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To our knowledge, this is the first time two-level Schwarz preconditioners are applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The pre-conditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the sub-domains. Several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored as well. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as non uniform meshes for the Greenland ice sheet are considered. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32 K processor cores (8 K MPI-ranks and 4 OpenMP threads) and 566 M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4 K processor cores (and MPI-ranks) and 68 M degrees of freedom for the coupled problem.

58 GEOSCIENCES↗
DOE OSTI2021

Dynamic Strength of Copper at High Pressures Using Pressure Shear Plate Experiments

Measurement of strength at high pressures is critical to understanding the behavior of materials subjected to extreme loading conditions. Recent developments have extended the pressure shear plate impact (PSPI) technique to the high-pressure regime. Such experiments provide a unique opportunity to extract the complete stress-strain behavior of materials at relatively high pressures. The modified technique includes a new fiber-optic heterodyne transverse velocity interferometer system and new analysis methods to account for the inelastic response of the anvil plates. In this study, PSPI experiments are conducted on oxygen-free high conductivity (OFHC) copper at pressures ranging from 10 to 43 GPa and at strain rates of ~10 5 s -1 . Complete stress-strain curves of copper are obtained using a hybrid methodology that relies on numerical simulations and particle velocity records gathered in these experiments. Experiments are conducted at similar pressures and strain rates using two types of anvils that undergo different levels of plastic strains: tool steel and tungsten carbide, to check the robustness of the hybrid method. The behavior of copper at these high pressures and strain rates is compared with previous literature data, and the effect of pressure and strain rate on the strength behavior of copper are discussed. Here, it was observed that the rate of change in yield strength change in pressure was 2.5 times the rate of change of shear modulus with pressure. The strength models that predict the response of copper at these extreme conditions are also examined and validated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Orographic Cirrus and Its Radiative Forcing in NCAR CAM6

Cirrus clouds play an important role in the Earth's radiative energy budget, thereby affecting the climate state and climate change. Orographic gravity wave (OGW)-induced sub-grid scale vertical velocity (i.e., cooling rate) is not resolved by large-scale models and its impact on ice formation in cirrus clouds is not well quantified. In this study, one sub-grid scale OGW scheme (e.g., McFarlane) is used in the Community Atmosphere Model version 6 (CAM6) to generate vertical velocity variance (σw) for cirrus formation. Results from the default model and simulations with the OGW-induced σw are evaluated against the DOE ARM Small Particles in Cirrus (SPARTICUS) campaign observations. The OGW based on the McFarlane scheme increases the sub-grid scale σw over mountains compared to the default model and improves the model agreement with the SPARTICUS observations. Larger σ w due to OGWs can trigger more frequent homogeneous nucleation in orographic cirrus and generates a higher number concentration of ice crystals observed during the SPARTICUS campaign. Moreover, our evaluation of the model simulations against satellite observations indicates that the McFarlane scheme generates high in-cloud ice number concentrations (>200 L -1 ) in the upper troposphere over mountains and high plateaus at mid- and high-latitudes of the winter hemisphere as shown in the observations. More ice crystals with smaller sizes absorb more infrared radiation (+0.523 ± 0.125 W m -2 ). The net radiative cloud forcing change at the top of the atmosphere is +0.330 W m -2 due to the orographic cirrus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Development of a Testbed for Pipeline Flushing - 20425

Pipelines which carry high-level waste within the DOE complex should be properly flushed in the cases where stationary or moving beds of solid sediment occur, or lines are prone to hydrogen gas buildup. Present guidelines establish a minimum for the amount of water (flush volume) and flush velocity values used for post-transfer flushing operations to achieve satisfactory cleaning of pipelines. However, further studies are needed to find optimal operating velocity modes/values in order to minimizing the flush volume and consequent downstream waste. These efforts are significantly helpful to DOE waste remediation sites by preserving tank storage, preventing additional waste processing, and minimizing dilution and changes in waste chemistry. An experimental test loop was recently developed for study of non-Newtonian slurry flushing at the Florida International University (FIU). This loop was designed to create sediment beds of various materials and bed heights and investigate parameters that effect the efficiency of flushing operations. The objective is to find flush velocity values/modes which lead to satisfactory cleaning of transport lines with a minimum amount of water usage. A 165 ft experimental test loop made of 3-inch carbon steel pipes was constructed and used in initial testing campaigns. This pipe loop was equipped with real-time monitoring instruments to record variation of pressure, density, mass flow rate, and temperature during testing. Post-flush in-situ no-flow evaluations were performed using ultrasonic, endoscopy, and visual inspection. Additional post-flush evaluations were conducted using filtration and density monitoring of residuals within a designated circulations loop. These provisions enable comparison of pipeline cleanness before the start and after termination of flushing operations for efficiency evaluation. This paper presents efforts associated with flushing of kaolin-water mixtures at different concentrations and with various flush volume and flush velocity values/modes. In initial testing, efforts focused on creating repeatable sediment beds inside the pipeline in both fully-flooded (no pre-flush drainage) and gravity-drained (pre-flush drained) systems. Bed characterization (i.e., measurement of sediment height and solids concentration) was conducted to ensure about consistency of initial conditions between tests. Preliminary results for flushing study of 10% vol. kaolin-water mixtures starting from a fully-flooded initial condition showed that flushing with 1-line volume (65 gallons) could remove up to 72 percent of solids from the pipeline. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗