SEARCH · Engineering Papers
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.
Background detector rates in the Electron Ion Collider depend in part on the pressure in the interaction region. Materials choice, synchrotron radiation induced desorption, conditioning time and pumping configuration all affect the pressure in the system. Simulations of the region using Synrad and Molflow+ coupled simulation codes will be presented for various configruations and conditioning times.
Abstract. Ice-free land comprises 26 % of the Earth's surface and holds liquid water that delineates ecosystems, affects global geochemical cycling, and modulates sea levels. However, we currently lack the capacity to simulate and predict these terrestrial water changes across the full range of relevant spatial (watershed to global) and temporal (monthly to millennial) scales. To address this knowledge gap, we present the Water Table Model (WTM), which integrates coupled components to compute dynamic lake and groundwater levels. The groundwater component solves the 2D horizontal groundwater flow equation using non-linear equation solvers from the C++ PETSc (Portable, Extensible Toolkit for Scientific Computation) library. The dynamic lake component makes use of the Fill–Spill–Merge (FSM) algorithm to move surface water into lakes, where it may evaporate or affect groundwater flow. In a proof-of-concept application, we demonstrate the continental-scale capabilities of the WTM by simulating the steady-state climate-driven water table for the present day and the Last Glacial Maximum (LGM; 21 000 calendar years before present) across the North American continent. During the LGM, North America stored an additional 14.98 cm of sea-level equivalent (SLE) in lakes and groundwater compared to the climate-driven present-day scenario. We compare the present-day result to other simulations and real-world data. Open-source code for the WTM is available on GitHub and Zenodo.
DayCent data and results for:J. L. Field, T. L. Richard, E. A. Smithwick, H. Cai, M. S. Laser, D. S. LeBauer, S. P. Long, K. Paustian, Z. Qin, J. J. Sheehan, P. Smith, M. Q. Wang, L. R. Lynd, Robust paths to net greenhouse gas mitigation and negative emissions via advanced biofuels. Proceedings of the National Academy of Sciences (2020). https://doi.org/10.1073/pnas.1920877117 This zip file contains a UNIX-format DayCent model executable, input files, automation code, and associated directory structure necessary to re-produce the DayCent analysis underlying the manuscript. The main script 'autodaycent.py' (written for Python 2.7) opens an interactive command line routine that facilitates:* Calibrating the DayCent pine growth model.* Initializing DayCent for a set of case studies sites.* Executing an ensemble of model runs representing case study site reforestation, grassland restoration, or conversion to switchgrass cultivation. * Results analysis & generation of manuscript Fig. 3. Note that the interactive analysis code requires that all input files to be contained in the directory structure as uploaded, without modification. Executable versions of the DayCent model (https://www.nrel.colostate.edu/projects/daycent/) compatible with other operating systems are available upon request.Please send questions/comments to John.L.Field@gmail.com
This data release includes experimental data of potentiometric titration for iron oxides. The data in the provided .csv files is not our own experimental data but have been compiled from the multiple literature sources. The master database is L-SCIE (LLNL Surface Complexation/Ion Exchange) database, and the provided .csv files are extracted data from L-SCIE. The .csv files were obtained by using the Lawrence Livermore National Laboratory Surface Complexation Database Converter (SCDC) code written in the R programming language (free licensing available at https://ipo.llnl.gov/technologies/software/llnl-surface-complexation-database-converter-scdc).The released data was used for developing a comprehensive community data-driven surface complexation modeling (SCM) framework for simulating potentiometric titration of mineral surfaces. Compiled community data for ferrihydrite, goethite, hematite, and magnetite are fit to produce representative protolysis constants that can reproduce potentiometric titration data collected from multiple literature sources.
Stability constraints play an important role in integrated modeling. The global plasma stability and macro-instabilities such as internal kink modes, neoclassical tearing modes, and edge localized perturbations limit the plasma performance and can result in large-scale transient events and plasma disruptions. These macro-instabilities can be also beneficial to the plasma performance. For example, the peeling models often lead to the edge localized modes (ELMs). However, the formation of stationary edge harmonic oscillations (EHOs) due to nonlinear interaction of peeling modes can result in a transition to the Quiescent H-mode (QH-mode) without ELMs. The profiles predicted with transport models need to be refined using the MHD stability constraints. This is important especially for transient and nonlinear stages of discharges such as ramp-up access to hybrid and steady state operation, and L- to H-mode transition as well as transition to the QH-mode. Several stability and MHD codes are already included in the OMFIT framework. These codes include BALOO, ELITE, GATO, M3DC1, MARS, and NIMROD. However, the verifications of stability constraints are currently mostly excluded from transport modeling workflows. The only exception is the EPED module which includes the ELITE predictions to limit the pedestal height. Here, we improved the OMFIT workflow to include stability calculations. Including the stability conditions to the integrated modeling workflow improved the robustness of the predictive modeling discharges. Being implemented in the workflow, the stability conditions can be also used for the stability analysis of experimental data. This improved the physics understanding of various discharge scenarios and can be used in experiment planning.
Explore the source record for details and available documents.
Abstract not provided.
Minimizing and understanding errors is critical for quantum science, both in noisy intermediate scale quantum (NISQ) devices and for the quest towards fault-tolerant quantum computation. Rydberg arrays have emerged as a prominent platform in this context with impressive system sizes and proposals suggesting how error-correction thresholds could be significantly improved by detecting leakage errors with single-atom resolution, a form of erasure error conversion. However, two-qubit entanglement fidelities in Rydberg atom arrays have lagged behind competitors and this type of erasure conversion is yet to be realized for matter-based qubits in general. Here we demonstrate both erasure conversion and high-fidelity Bell state generation using a Rydberg quantum simulator. When excising data with erasure errors observed via fast imaging of alkaline-earth atoms, we achieve a Bell state fidelity of $\ge 0.997{1}_{-13}^{+10}$ , which improves to $\ge 0.998{5}_{-12}^{+7}$ when correcting for remaining state-preparation errors. We further apply erasure conversion in a quantum simulation experiment for quasi-adiabatic preparation of long-range order across a quantum phase transition, and reveal the otherwise hidden impact of these errors on the simulation outcome. Our work demonstrates the capability for Rydberg-based entanglement to reach fidelities in the 0.999 regime, with higher fidelities a question of technical improvements, and shows how erasure conversion can be utilized in NISQ devices. These techniques could be translated directly to quantum-error-correction codes with the addition of long-lived qubits.
Magnetized target fusion approach to inertial confinement fusion involves the formation of strong shocks that travel along a magnetized plasma. Shocks, which play a dominant role in thermalizing the upstream kinetic energy generated in the implosion stage, are seldom free from perturbations, and they wrinkle in response to upstream or downstream disturbances. In Z-pinch experiments, significant plasma instability mitigation was observed with pre-embedded axial magnetic fields. To isolate effects, in this work we theoretically study the impact of perpendicular magnetic fields on the planar shock dynamics for different equations of state. For fast magnetosonic shocks in ideal gases, it was found that the magnetic field amplifies the intensity of the perturbations when $\gamma$ > 2 or it weakens them when $\gamma$ < 2. Weak shocks have been found to be stable regardless of the magnetic plasma intensity and gas compressibility; however, for sufficiently strong shocks the magnetic fields can promote a neutral stability/SAE at the shock if the adiabatic index is higher than 1 + $\sqrt{2}$. Finally, results have been validated with numerical simulations performed with the FLASH code.
The FEMP Cyber Security Arsenal is a family of cyber security tools for the federal facility owners and operators. Using these tools, facility owners can evaluate their overall cybersecurity posture. These tools are web-based front-end tools. This is an upgrade of the tools available from the earlier version in DOE CODE ID: #41457
Abstract Tides are often nonstationary due to nonastronomical influences. Investigating variable tidal properties implies a trade-off between separating adjacent frequencies (using long analysis windows) and resolving their time variations (short analysis windows). Previous continuous wavelet transform (CWT) tidal methods resolved tidal species. Here, we present CWT_Multi, a MATLAB code that 1) uses CWT linearity (via the “response coefficient method”) to implement superresolution, i.e., resolving tidal constituents beyond the Rayleigh criterion; 2) provides a Munk–Hasselmann constituent selection criterion appropriate for superresolution; and 3) introduces an objective, time-variable form of inference (“dynamic inference”) based on time-varying data properties. CWT_Multi resolves tidal species on time scales of days, and multiple constituents per species with fortnightly filters. It outputs astronomical phase lags and admittances, analyzes multiple records, and provides power spectra of the signal(s), residual(s), and reconstruction(s); confidence limits; and signal-to-noise ratios. Artificial data and water levels from the Lower Columbia River Estuary (LCRE) and San Francisco Bay Delta (SFBD) are used to test CWT_Multi and compare it to harmonic analysis programs NS_Tide and UTide. CWT_Multi provides superior reconstruction, detiding, dynamic analysis utility, and time resolution of constituents (but with broader confidence limits). Dynamic inference resolves closely spaced constituents (like K 1 , S 1 , and P 1 ) on fortnightly time scales, quantifying impacts of diel power peaking (with a 24-h period, like S 1 ) on water levels in the LCRE. CWT_Multi also helps quantify the impacts of high flows and a salt barrier closing on tidal properties in the SFBD. On the other hand, CWT_Multi does not excel at prediction, and results depend on analysis details, as for any method applied to nonstationary data. Significance Statement Ocean tides, especially in coastal and estuarine systems, are often nonstationary, in the sense that the mean and standard deviation of tidal properties vary over time, usually in response to some nontidal process. We introduce here a MATLAB code, CWT_Multi, that uses wavelet transforms to resolve both tidal species and constituents on time scales from a few days to months. Our code accommodates multiple scalar time series and has typical tidal analysis features like constituent selection and inference, plus two forms of uncertainty analyses. It is flexible, allowing the user to adapt analysis properties to diverse datasets. CWT_Multi is applicable to many problems involving time-variable tides, including sea level rise, compound flooding, sediment transport, and wetland habitat analyses. Application to vector data is a straightforward extension, but further development of our uncertainty analysis is merited. Because nonstationary tidal analysis is rapidly advancing, we also define the features of a “well-formed” analysis code.
Annual cropping systems are common in the Driftless Region of the U.S. Midwest, but soil degradation is prone to happen in such systems due to the rugged topography of the region. Recent rapid increases in row crop area have been noted in this region, with annual precipitation and hydrologic extremes on the rise in recent decades. The aim of this research to use geospatial datasets and tools in order to assess the regional trends in land use, precipitation, and hydrologic change and quantify the relationship between these environmental trends. Between 2006 and 2017, substantial row crop expansion of 10,000 ha or more was common across HUC 8 (Hydrologic Unit Code 8) watersheds in our study area. Expansion occurred mainly on steeper slopes, converting existing grasslands or alfalfa (Medicago sativa L.) to row crops. Classifying land as planted (in row crops), plantable (in row crops or could be converted), and unplantable (unable to be converted) revealed that Driftless Region watersheds have ~30–50% of plantable land available for future expansion. Annual precipitation was highly variable during this time period but had a general increasing trend. On average, precipitation showed higher correlation to streamflow compared to row crop expansion across 27 USGS river gage drainage basins in our study area. However, when the increase in row crop area was significant and was accompanied by increasing precipitation, stronger correlation between row crop area and annual streamflow was exhibited. This finding suggests that row crop expansion acts to enhance the effects of increasing precipitation on local hydrology.
First patch release of v2.0.0 Changes since v2.0.0 [Important change] Fix ocean threading bug seen in debug cases on Chrysalis with Intel 20.0.4. Was introduced around time of v2.0.0 tag. Does not change v2.0.0 answers on Chrysalis because those didn't use threading or debugging. [EAM] Add semi-lagrangian tracer transport for theta-l (F90 and C++), add new algorithm for finding tropopause, add DSCREAM to allow v2 and SCREAM settings in same code such as adjust_ps [EAM-MMF] 60L default, allow C++ back end of RRTMGP (EAM too). [EAMxx] add nu-top functionality, fix forcing functor, add ttype9 and dcmip2012 tests 2.1, 2.2, and 3 HOMME: remove obsolete remap algs, option to specify dynamics alg indep of tracer, new sponge layer, add imex tests [ELM] Add topography-based subgrid (topounits), add FATES-ELM Nitro., Phos. and CH4 coupling, add land-use ts for NARRM, add lulc for SSP3 RCP7, Fix nutrient fertilization exp test and carbon isotope flux, Fix xactive lnd dry deposition, add lake water storage option, fix plant hydraulics 2d params, fix carbon budget calc, fix soil nutrient conc. bug, fix mosart dam bug, add test for new ELM, MOSART features, fix bug in O3 dry dep stomatal resistances, fix plant hydraulics restart BFB error, update mkmapdata. [MOSART] fix bug for reading the latitude from an unstructured input file, fix oversat in bubble test. [MPAS-ocean] Add CFC11, CFC12 tracers, add 2D spherical transport tests, fix del4 tracer mixing, add MARBL ocean tracer mixing, modify harmonic analysis options, add GPU port of vmix routines, fix calc of ML-averaged BV freq. [MPAS-seaice] Change extents of initial polar disks for oRRS18to6v3 grid, fix ice BGC with MARBL, update spherical test cases, fix DON coupling, Remove Cf from sea ice constants. [MPAS-landice] add CRYO1850-4xCO2 compset [CIME] add GCP, ANL GCE, Spock, Perlmutter, deprecate config_compilers.xml, fix and clean-up cmake macros, fix slurm bindings, refactor CIME internal testing, cleanup SCORPIO perf data, allow position independent compset naming, [also] update v2 benchmarking suite, extend e3sm_prod with throughput and memory checks
This is a composite 3D seismic velocity that was constructed from compiled information from several local studies regarding seismic velocities and structural information. This seismic velocity model is provided in NonLinLoc format (slow_len), which is readily usable in NonLinLoc software. Other model formats and versions of the model can be produced using the Python script provided with this data set. Details on how the model was created and prior velocity and structural information was used is provided in the accompanying documentation.
We describe the Studying Quenching in Intermediate-z Galaxies: Gas, $\scriptsize\mathrm{\rm{angu}}\overrightarrow{L}{\rm{ar}}$ momentum, and Evolution ($\scriptsize\mathrm{{SQuIGG\overrightarrow{L}E}}$) survey of intermediate-redshift post-starburst galaxies. We leverage the large sky coverage of the Sloan Digital Sky Survey to select ~1300 recently quenched galaxies at 0.5 < z ≤ 0.9 based on their unique spectral shapes. These bright, intermediate-redshift galaxies are ideal laboratories to study the physics responsible for the rapid quenching of star formation: they are distant enough to be useful analogs for high-redshift quenching galaxies, but low enough redshift that multiwavelength follow-up observations are feasible with modest telescope investments. We use the Prospector code to infer the stellar population properties and nonparametric star formation histories (SFHs) of all galaxies in the sample. We find that $\scriptsize\mathrm{{SQuIGG\overrightarrow{L}E}}$ galaxies are both very massive (M * ~ 10 11.25 M ⊙ ) and quenched, with inferred star formation rates ≲1 M ⊙ yr -1 , more than an order of magnitude below the star-forming main sequence. The best-fit SFHs confirm that these galaxies recently quenched a major burst of star formation: >75% of $\scriptsize\mathrm{{SQuIGG\overrightarrow{L}E}}$ galaxies formed at least a quarter of their total stellar mass in the recent burst, which ended just ~200 Myr before observation. We find that $\scriptsize\mathrm{{SQuIGG\overrightarrow{L}E}}$ galaxies are on average younger and more burst-dominated than most other z ≲ 1 post-starburst galaxy samples. This large sample of bright post-starburst galaxies at intermediate redshift opens a wide range of studies into the quenching process. In particular, the full $\scriptsize\mathrm{{SQuIGG\overrightarrow{L}E}}$ survey will investigate the molecular gas reservoirs, morphologies, kinematics, resolved stellar populations, active galactic nucleus incidence, and infrared properties of this unique sample of galaxies in order to place definitive constraints on the quenching process.
Transient analyses of heat pipe micro-reactors are necessary to ensure that hypothetical accident scenarios do not comprise reactor safety. Due to its small size and reliance on heat-pipes for cooling, the micro-reactor design introduced in this paper is a tightly coupled system which requires multi-physics tools to accurately model transient events. Idaho National Laboratory's DireWolf code suite based on the MOOSE framework is tailor-built to model heat-pipe reactors. This paper demonstrates DireWolf's ability to simulate the coupled thermal-neutronics transient behavior of a heat-pipe micro-reactor. The transient events presented here include an inadvertent rotation of all control drums simultaneously and a sudden complete rotation of a single control drum. A detailed description of each event is provided along with simulation results, including time dependent power and temperature distributions, and discussion. (authors)
The magnetic pre-sheath (MPS) width, L MPS , is a critical parameter to define the sheath potential, which controls the ion trajectory of low-Z species (D, T, He, and C), as well as the prompt re-deposition of high-Z species. In this work, to determine L MPS , we fabricated micro-trenches (30×30×4 µm) via focused ion beam (FIB) milling on a silicon surface and exposed them to L-mode deuterium plasmas in DIII-D via the Divertor Material Evaluation System (DiMES) removable sample exposure probe. The areal distribution of impurity depositions, mainly consisting of carbon, was measured by energy-dispersive X-ray spectroscopy (EDS) to reveal the deuterium ion shadowing effect on the trench floors. The carbon deposition profiles showed that the erosion was maximized for the azimuthal direction of φ = -40° (referenced to the toroidal magnetic field direction) as well as the polar angle of θ = 80°. A Monte Carlo equation-of-motion model, based on a collisionless MPS, was used to calculate the azimuthal and polar deuterium ion angle distributions (IADs) for a range of L MPS = k × ρ i , where ρ i is the ion gyro radius and k = 0.5-4. Then, gross erosion profiles were calculated by a Monte Carlo micro-patterning and roughness (MPR) code for ion sputtering using as input the calculated azimuthal and polar IADs for each value of k . Good agreement with the experimental C deposition profiles was obtained for the case k = 2.5-3.5. This result is consistent with a previous kinetic modeling prediction of k ~ 3, as well as previous analytical investigations that predicted the L MPS to be several ion gyro radii. A validation of theoretical sheath models supports its applicability to ITER and pilot plant divertors to successfully predict plasma-materials interactions.