Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reach code”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-I-31: Differences Observed in Radiation Doses Across 2 Similar CT Scanners From Adult Brain-Neck CT Angiography

Purpose: The aim of this study is to evaluate the difference in radiation doses from adult Brain-Neck CT angiography (CTA) between two CT scanners. Methods: We collected CT dose index data (CTDIvol, DLP) from adult Brain-Neck CTA performed with two CT scanners (Sensation 64 (S64) and Definition AS (AS), Siemens Healthcare) performed at two of our facilities from Jan 1st to Dec 31th, 2014. X-ray dose management software (Radmetrics, Bayer Healthcare) was used to mine these data. All exams were performed with Tube Current Modulation (Care Dose 4D), tube voltage of 120 kVp, quality reference mAs of 300, beam collimation of 64*0.6 mm. The rotation time was set to 0.5 sec for S64 and 1.0 sec for AS. We also scanned an anthropomorphic skull and chest phantom under routine Brain-Neck CTA protocol with the two scanners and extracted the tube current values from the raw projection data. Results: The mean CTDIvol and DLP in Brain-Neck CTA was 72 mGy and 2554 mGy*cm for AS, which was substantially larger than the mean values of 46 mGy and 1699 mGy*cm for S64. The maximum tube current was 583 mA for most cases on the S64 while the maximum was 666 mA for AS even though the rotation time set for AS was 1.0 sec. Measurements obtained with the anthropomorphic phantom showed that the tube current reached 583 mA at the shoulder region for S64 while it reached to 666 mA for AS. Conclusion: The results of this study showed that substantially different CT doses can Result from Brain-Neck CTA protocols even when similar scanners and similar settings are used. Though both scanners have a similar maximum mA rating, differences in mA were observed through the shoulders, resulting in substantially different CTDIvol values.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Equilibrium β-limits in classical stellarators

Here, a numerical investigation is carried out to understand the equilibrium β-limit in a classical stellarator. The stepped-pressure equilibrium code is used in order to assess whether or not magnetic islands and stochastic field-lines can emerge at high β. Two modes of operation are considered: a zero-net-current stellarator and a fixed-iota stellarator. Despite the fact that relaxation is allowed, the former is shown to maintain good flux surfaces up to the equilibrium β-limit predicted by ideal-magnetohydrodynamics (MHD), above which a separatrix forms. The latter, which has no ideal equilibrium β-limit, is shown to develop regions of magnetic islands and chaos at sufficiently high β, thereby providing a ‘non-ideal β-limit’. Perhaps surprisingly, however, the value of β at which the Shafranov shift of the axis reaches a fraction of the minor radius follows in all cases the scaling laws predicted by ideal-MHD. We compare our results to the High-Beta-Stellarator theory of Freidberg and derive a new prediction for the non-ideal equilibrium β-limit above which chaos emerges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Advanced instrumentation and analysis methods for in-pile thermal and nuclear measurements: from out-of-pile studies to irradiation campaigns

Research and development on nuclear fuel behavior under irradiations and accelerated ageing of structure materials is a key issue for sustainable nuclear energy in order to meet specific needs by keeping the best level of safety. A new Material Testing Reactor (MTR), the Jules Horowitz Reactor (JHR) currently under construction in the South of France in the CEA Cadarache research centre will offer a real opportunity to perform R and D programs and hence will crucially contribute to the selection, optimization and qualification of innovative materials and fuels. To perform such programs advanced accurate and innovative experiments, irradiation devices that contain material and fuel samples are required to be set up inside or beside the reactor core. These experiments needs beforehand in situ and on line sophisticated measurements to accurately reach specific and determining parameters such as thermal and fast neutron fluxes, nuclear heating and temperature conditions to precisely monitor and control the conducted assays. Consequently, since 2009 CEA and Aix-Marseille University collaborate in order to design and develop a new multi-sensor device which will be dedicated to measuring profiles of such conditions inside the experimental channels of the JHR. These works are performed in the framework of two complementary joint research programs called MAHRI-BETHY and INCORE. These programs couple experimental studies carried out both out-of nuclear fluxes (in laboratory) and under irradiation conditions (in OSIRIS MTR reactor in France and MARIA MTR reactor in Poland) with numerical works realized by thermal simulations (CAST3M code) and Monte Carlo simulations (MCNP code). These programs deal with three main aims. The first one corresponds to the design and/or the test of new in-pile instrumentation. The second one concerns the development of advanced calibration procedures in particular in the case of one specific sensor: a differential calorimeter used to quantify nuclear heating. The last one consists in the development of accurate measurement and analysis methods. The paper will be dedicated to a complete review of the experimental and numerical works performed since 2009 thanks to two parts. The first part will detail a new thermal approach implemented to improve nuclear heating measurements by radiometric calorimeters. New experimental tools (calorimeter prototypes and set-ups such BETHY Bench) developed to perform preliminary out-of-pile studies under suitable conditions will be presented (temperature and velocity of the external cooling fluid, heat source localization and intensity inside the calorimetric cells). Then the response of two kinds of sensors, their calibrations curves and their thermal behaviors will be compared for various parameters. Finally validated numerical thermal and Monte Carlo works will be discussed to propose new improvements. The second parts of the paper will focus on works realized in order to design, develop and test the first prototype of the multi-sensor device called CARMEN [7-9]. The two mock-ups dedicated respectively to neutron measurements and photon measurements will be detailed. The results obtained during two irradiation campaigns inside the periphery of OSIRIS reactor will be shown. The new analysis method will be discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

OpenMP Target Task: Tasking and Target Offloading on Heterogeneous Systems

This work evaluated the use of OpenMP tasking with target GPU offloading as a potential solution for programming productivity and performance on heterogeneous systems. Also, it is proposed a new OpenMP specification to make the implementation of heterogeneous codes simpler by using OpenMP target task, which integrates both OpenMP tasking and target GPU offloading in a single OpenMP pragma. As a test case, the authors used one of the most popular and widely used Basic Linear Algebra Subprogram Level-3 routines: triangular solver (TRSM). To benefit from the heterogeneity of the current high-performance computing systems, the authors propose a different parallelization of the algorithm by using a nonuniform decomposition of the problem. This work used target GPU offloading inside OpenMP tasks to address the heterogeneity found in the hardware. This new approach can outperform the state-of-the-art algorithms, which use a uniform decomposition of the data, on both the CPU-only and hybrid CPU-GPU systems, reaching speedups of up to one order of magnitude. The performance that this approach achieves is faster than the IBM ESSL math library on CPU and competitive relative to a highly optimized heterogeneous CUDA version. One node of Oak Ridge National Laboratory’s supercomputer, Summit, was used for performance analysis.

Valero Lara, Pedro↗
DOE OSTI2023

Exploring the parameter space of MagLIF implosions using similarity scaling. III. Rise-time scaling

Magnetized liner inertial fusion (MagLIF) is a z-pinch magneto-inertial-fusion concept studied at the Z Pulsed Power Facility of Sandia National Laboratories. Two important metrics characterizing current delivery to a z-pinch load are the peak current and the current-rise time, which is roughly the time interval to reach the peak current. It is known that, when driving a z-pinch load with a longer current-rise time, the performance of the z-pinch decreases. However, a theory to understand and quantify this effect is still lacking. Here, in this paper, we utilize a framework based on similarity scaling to analytically investigate the variations in the performance of MagLIF loads when varying the current-rise time, or equivalently, the implosion timescale. To maintain similarity between the implosions, we provide scaling prescriptions of experimental input parameters defining a MagLIF load and derive the expected scaling laws for stagnation conditions and for various performance metrics. We compare predictions of the theory to 2D numerical simulations using the radiation, magneto-hydrodynamic code hydra. For several metrics, we find acceptable agreement between the theory and simulations. Our results show that the voltage φload near the MagLIF load follows a weak scaling law φload ∝$t_{φ}^{-0.12}$ with respect to the characteristic timescale t φ of the voltage source, instead of the ideal φload ∝$t_{φ}^{-1}$ scaling. This occurs because the imploding height of the MagLIF load must increase to preserve end losses. As a consequence of the longer imploding liners, the required total laser preheat energy and delivered electric energy increase. Overall, this study helps understand the trade-offs of the MagLIF design space when considering future pulsed-power generators with shorter and longer current-rise times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

LINE EMISSION FROM RADIATION-PRESSURIZED H II REGIONS. I. INTERNAL STRUCTURE AND LINE RATIOS

The emission line ratios [O III] {lambda}5007/H{beta} and [N II] {lambda}6584/H{alpha} have been adopted as an empirical way to distinguish between the fundamentally different mechanisms of ionization in emission-line galaxies. However, detailed interpretation of these diagnostics requires calculations of the internal structure of the emitting H II regions, and these calculations depend on the assumptions one makes about the relative importance of radiation pressure and stellar winds. In this paper, we construct a grid of quasi-static H II region models to explore how choices about these parameters alter H II regions' emission line ratios. We find that when radiation pressure is included in our models, H II regions reach a saturation point beyond which further increase in the luminosity of the driving stars does not produce any further increase in effective ionization parameter, and thus does not yield any further alteration in an H II region's line ratio. We also show that if stellar winds are assumed to be strong, the maximum possible ionization parameter is quite low. As a result of this effect, it is inconsistent to simultaneously assume that H II regions are wind-blown bubbles and that they have high ionization parameters; some popular H II region models suffer from this inconsistency. Our work in this paper provides a foundation for a companion paper in which we embed the model grids we compute here within a population synthesis code that enables us to compute the integrated line emission from galactic populations of H II regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Few-shot leaf segmentation

In this work, we use few-shot learning to segment the body and vein architecture of Populus trichocarpa leaves from high-resolution scans obtained in the UC Davis common garden. Leaf and vein segmentation are formulated as separate tasks, in which convolutional neural networks (CNNs) are used to iteratively expand partial segmentations until reaching stopping criteria. Our leaf and vein segmentation approaches use just 50 and 8 manually traced images for training, respectively, and are applied to a set of 2,634 top and bottom leaf scans. We show that both methods achieve high segmentation accuracy, in some cases exceeding even human-level segmentation. The leaf and vein segmentations are subsequently used to extract 68 morphological traits using traditional open-source image processing tools, which are validated using real-world physical measurements. For a biological perspective, we perform a genome-wide association study using the "vein density" trait to discover novel genetic architectures associated with multiple physiological processes relating to leaf development and function. In addition to sharing all of the few-shot learning code, we are releasing all images, manual segmentations, model predictions, 68 extracted leaf phenotypes, and a new set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes.

Lagergren, John [Oak Ridge National Lab. (ORNL), O↗
DOE OSTI2023

Supporting information for Few-Shot Learning Enables Population-Scale Analysis of Leaf Traits in Populus trichocarpa

In this work, we use few-shot learning to segment the body and vein architecture of P. trichocarpa leaves from high-resolution scans obtained in the UC Davis common garden. Leaf and vein segmentation are formulated as separate tasks, in which convolutional neural networks (CNNs) are used to iteratively expand partial segmentations until reaching stopping criteria. Our leaf and vein segmentation approaches use just 50 and 8 manually traced images for training, respectively, and are applied to a set of 2,634 top and bottom leaf scans. We show that both methods achieve high segmentation accuracy, in some cases exceeding even human-level segmentation. The leaf and vein segmentations are subsequently used to extract 68 morphological traits using traditional open-source image processing tools, which are validated using real-world physical measurements. For a biological perspective, we perform a genome-wide association study using the vein density trait to discover novel genetic architectures associated with multiple physiological processes relating to leaf development and function. In addition to sharing all of the few-shot learning code (see https://github.com/jlager/few-shot-leaf-segmentation), we are releasing all images, manual segmentations, model predictions, 68 extracted leaf phenotypes, and a new set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes. The data folder includes all images, ground truth segmentations, predicted segmentations, and extracted leaf traits. All images encode the sample ID in the file name by indicating the treatment, block, row, position, and leaf side, respectively. For example, the file, C_1_1_2_bot.jpeg, indicates the control treatment, block 1, row 1, position 2, and the bottom side of the leaf. Tabulated results include position IDs as well as the corresponding genotype IDs. The images folder includes the 2,906 high-resolution leaf scans taken in the field. The leaf_masks folder includes 50 ground truth segmentations used for training the leaf tracing algorithm. The leaf_preds folder includes the 2,906 predicted segmentations from the leaf tracing algorithm. The vein_masks folder includes 8 ground truth segmentations used for training the vein growing algorithm. The vein_preds folder includes the 1,453 predicted segmentations from the vein growing algorithm. The vein_probs folder includes the 1,453 predicted probability maps from the vein growing algorithm before thresholding. The genomes folder includes the set of SNPs called against the v4 P. trichocarpa genome for 1,419 genotypes with a README file detailing the steps taken. The results folder includes: (i) raw values of the 68 predicted leaf traits in digital_traits.tsv, (ii) manually measured values of petiole length and width in manual_traits.tsv, (iii) thin plate spline (TPS) adjusted values of the vein density trait in vein_density_tps_adj.tsv, (iv) best linear unbiased prediction (BLUP) adjusted values of the vein density trait in vein_density_blups.tsv, and (v) GWAS results for the vein density trait, including chromosome positions and corresponding P values, in gwas_results.csv.

09 BIOMASS FUELS↗
DOE OSTI2016

Conditional Tree Reduction in the ADAPT Framework

Probabilistic Risk Assessment (PRA) uses fault tree/event tree analysis to evaluate the environmental impacts of nuclear power plants due to internal or external initiating events. In a fault tree, basic failure events are assembled using primarily And/Or logic to determine the combinations of failures that lead to failure of the system as a whole. An event tree is forward-facing, and begins with a single initiating event. From there the analysis considers what event may occur next, and branches out among the possible configurations of occurrence and non-occurrence of that event. This branching continues until end states are reached. End states may include a safe and stable configuration of the plant, or any number of differing failure states. Both fault and event trees require basic event probabilities as inputs to provide insight on the likelihood of different outcomes. PRA has been applied to nuclear power plants as an analysis tool since the 1975 Reactor Safety Study. In a traditional PRA, the order of events is prescribed by the analyst and each event is typically a binary aleatory uncertainty: occurrence or non-occurrence. Discrete dynamic event tree (DDET) analysis eliminates this subjectivity in ordering of events by using the output of a dynamic system model (simulator) to inform the branching. Branching conditions are triggered by the existence of a relevant plant state in the code, and therefore only occur as physically appropriate. The DDET approach also allows consideration of both epistemic and aleatory uncertainties on a phenomenologically and stochastically consistent platform. A DDET may be orders of magnitude larger than a traditional event tree for the same initiating event, as more of the uncertainty space is likely to be explored and the tree does not require manual assembly. Because of the size, manual inspection of the entire DDET is often infeasible. This problem may be approached from the front by reducing the number of sampled branching values or by online pruning of the tree, thus reducing the total number of branches. Utilizing post-processing techniques to draw actionable insights via clustering has also been proposed. One post-processing tool to reduce the DDET according to user-defined rules has been recently implemented in the Analysis of Dynamic Accident Progression Trees (ADAPT) DDET driver code. The tool creates a 'slice' of the DDET comprised of the set of branches that meet a set of user-input rules, for immediate visualization or for further calculations conditional on the rules being met. This extension is expected to ease analysis of the results of large and complex DDETs, leading to greater applicability of DDETs in PRA. It should be emphasized that this tool is a visualization aide and does not impact the theoretical basis of the DDET approach. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Extension of the ADAPT Framework for Multiple Simulators

Probabilistic Risk Assessment (PRA) uses fault tree/event tree analysis to evaluate the environmental impacts of nuclear power plants due to internal or external initiating events. In a fault tree, basic failure events are assembled using primarily And/Or logic to determine the combinations of failures that lead to failure of the system as a whole. An event tree is forward-facing, and begins with a single initiating event. From there the analysis considers what event may occur next, and branches out among the possible configurations of occurrence and non-occurrence of that event. This branching continues until end states are reached. End states may include a safe and stable configuration of the plant, or any number of differing failure states. Both fault and event trees require basic event probabilities as inputs to provide insight on the likelihood of different outcomes. PRA has been applied to nuclear power plants as an analysis tool since the 1975 Reactor Safety Study. In a traditional PRA, the order of events is prescribed by the analyst and each event is typically a binary aleatory uncertainty: occurrence or non-occurrence. Discrete dynamic event tree (DDET) analysis eliminates this subjectivity in ordering of events by using the output of a dynamic system model (simulator) to inform the branching. Branching conditions are triggered by the existence of a relevant plant state in the code, and therefore only occur as physically appropriate. The DDET approach also allows consideration of both epistemic and aleatory uncertainties on a phenomenologically and stochastically consistent platform. One limitation of the DDET driver codes developed to date (namely ADAPT, ADS, and MCDET ) is that the tree is driven by a single simulator. The goal of this work is a general multi-simulator driver. This limits the parameter space that may be explored in a single DDET to that which may be covered by a single simulator. In this work, the Analysis of Dynamic Accident Progression Trees (ADAPT) code is extended to support multiple arbitrary simulators. This capability will increase the depth and breadth of phenomena that may be analyzed in a single DDET. For example, a large-scale system code may pass the analysis to a more detailed code if certain conditions are met, in order to capture relevant phenomena that would be lost if only the system code were to be used. In addition, the time required to cover a given uncertainty space may be reduced by dynamically identifying and running only the situations that are appropriate given the plant conditions in each simulator. The use of multiple simulators is expected to further reduce subjectivity in PRA by enabling high fidelity modeling of reality. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Dynamic Importance Measures in the ADAPT Framework

Probabilistic Risk Assessment (PRA) uses fault tree/event tree analysis to evaluate the environmental impacts of nuclear power plants due to internal or external initiating events. In a fault tree, basic failure events are assembled using primarily And/Or logic to determine the combinations of failures that lead to failure of the system as a whole. An event tree is forward-facing, and begins with a single initiating event. From there the analysis considers what event may occur next, and branches out among the possible configurations of occurrence and non-occurrence of that event. This branching continues until user-defined end states are reached. End states may include a safe and stable configuration of the plant, or any number of differing failure states. Both fault and event trees require basic event probabilities as inputs to provide insight on the likelihood of different outcomes. PRA has been applied to nuclear power plants as an analysis tool since the 1975 Reactor Safety Study. In a traditional PRA, the order of events is prescribed by the analyst and each event is typically a binary aleatory uncertainty: occurrence or non-occurrence. Discrete dynamic event tree (DDET) analysis eliminates this subjectivity in ordering of events by using the output of a dynamic system model (simulator) to inform the branching. Branching conditions are triggered by the existence of a relevant plant state in the code, and therefore only occur as physically appropriate. The DDET approach also allows consideration of both epistemic and aleatory uncertainties on a phenomenologically and stochastically consistent platform. One current limitation of DDETs is the assessment of results, which is not as well developed as in traditional PRA. A concept that is used in PRA to assess the significance of a basic event is the concept of importance measures (IMs). In traditional PRA, IMs only consider the likelihood of occurrence of non-occurrence of a basic event. The application of importance measures to DDETs must take into consideration not only the occurrence and non-occurrence of an event, but also uncertain timing or severity. This work describes three general dynamic importance measures (DIMs) as implemented in the Analysis of Dynamic Accident Progression Trees (ADAPT) DDET driver code. The measures account for the change in progression of the DDET resulting from the different values of each uncertain variable under investigation. The DIMs are applied to a sample DDET to demonstrate the type of insights that may be gained from each one. These measures will facilitate comparison of the impact of dynamic branching events on a consequence of interest, which also allow DDETs to be used in prioritizing plant investments. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Minibeam radiation therapy for the management of osteosarcomas: A Monte Carlo study

Purpose: Minibeam radiation therapy (MBRT) exploits the well-established tissue-sparing effect provided by the combination of submillimetric field sizes and a spatial fractionation of the dose. The aim of this work is to evaluate the feasibility and potential therapeutic gain of MBRT, in comparison with conventional radiotherapy, for osteosarcoma treatments. Methods: Monte Carlo simulations (PENELOPE/PENEASY code) were used as a method to study the dose distributions resulting from MBRT irradiations of a rat femur and a realistic human femur phantoms. As a figure of merit, peak and valley doses and peak-to-valley dose ratios (PVDR) were assessed. Conversion of absorbed dose to normalized total dose (NTD) was performed in the human case. Several field sizes and irradiation geometries were evaluated. Results: It is feasible to deliver a uniform dose distribution in the target while the healthy tissue benefits from a spatial fractionation of the dose. Very high PVDR values (⩾20) were achieved in the entrance beam path in the rat case. PVDR values ranged from 2 to 9 in the human phantom. NTD{sub 2.0} of 87 Gy might be reached in the tumor in the human femur while the healthy tissues might receive valley NTD{sub 2.0} lower than 20 Gy. The doses in the tumor and healthy tissues might be significantly higher and lower than the ones commonly delivered used in conventional radiotherapy. Conclusions: The obtained dose distributions indicate that a gain in normal tissue sparing might be expected. This would allow the use of higher (and potentially curative) doses in the tumor. Biological experiments are warranted.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

HPC4Mfg with Samsung: Making semiconductor devices cool through HPC ab initio simulations

For decades, the semiconductor technology has followed the Moore’s law, making newer devices more powerful and energy efficient. Recently, however, it has reached a point where the performance and the energy efficiency of the device do not improve with the shrinking device size. One of the fundamental reasons of this deviation from the past trend is the interconnect resistance, which becomes larger with the shrinking size. The devices size is so small that the quantum mechanical effects can no longer be ignored and the traditional continuum simulation tools such as TCAD become inadequate. In this project, Samsung Semiconductor Inc. and Lawrence Berkeley National Laboratory has collaborated to perform first of kind device-scale ab initio simulations to optimize materials and interconnect morphology to minimize interconnect resistance. We have tested the use of LS3DF method and the PEtot_trans approach on top of the folded spectrum method (FSM) Escan code to calculate the scattering state, and to study various effects influence the interconnect conductivity. We found that, the LS3DF can be used to calculate such metallic system. On the other hand, the use of Escan code to solve the linear equation is not practical due to the slow convergence. We have implemented a Chebyshev filter technique to calculate a few hundred eigen states near the scattering state energy E, then use these eigen states as preconditioner to solve the linear equation. We have used this approach to study the different factors which affect the interconnect conductivity, including the shape, the point defect, the temperature, and the grain boundary.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Variant profiling of evolving prokaryotic populations

Genomic heterogeneity of bacterial species is observed and studied in experimental evolution experiments and clinical diagnostics, and occurs as micro-diversity of natural habitats. The challenge for genome research is to accurately capture this heterogeneity with the currently used short sequencing reads. Recent advances in NGS technologies improved the speed and coverage and thus allowed for deep sequencing of bacterial populations. This facilitates the quantitative assessment of genomic heterogeneity, including low frequency alleles or haplotypes. However, false positive variant predictions due to sequencing errors and mapping artifacts of short reads need to be prevented. We therefore created VarCap, a workflow for the reliable prediction of different types of variants even at low frequencies. In order to predict SNPs, InDels and structural variations, we evaluated the sensitivity and accuracy of different software tools using synthetic read data. The results suggested that the best sensitivity could be reached by a union of different tools, however at the price of increased false positives. We identified possible reasons for false predictions and used this knowledge to improve the accuracy by post-filtering the predicted variants according to properties such as frequency, coverage, genomic environment/localization and co-localization with other variants. We observed that best precision was achieved by using an intersection of at least two tools per variant. This resulted in the reliable prediction of variants above a minimum relative abundance of 2%. VarCap is designed for being routinely used within experimental evolution experiments or for clinical diagnostics. The detected variants are reported as frequencies within a VCF file and as a graphical overview of the distribution of the different variant/allele/haplotype frequencies. The source code of VarCap is available at https://github.com/ma2o/VarCap . In order to provide this workflow to a broad community, we implemeted VarCap on a Galaxy webserver, which is accessible at http://galaxy.csb.univie.ac.at .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Processed Lab Data for Neural Network-Based Shear Stress Level Prediction

Machine learning can be used to predict fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. The files are extracted features and labels from lab data (experiment p4679). The features are extracted with a non-overlapping window from the original acoustic data. The first column is the time of the window. The second and third columns are the mean and the variance of the acoustic data in this window, respectively. The 4th-11th column is the the power spectrum density ranging from low to high frequency. And the last column is the corresponding label (shear stress level). The name of the file means which driving velocity the sequence is generated from. Data were generated from laboratory friction experiments conducted with a biaxial shear apparatus. Experiments were conducted in the double direct shear configuration in which two fault zones are sheared between three rigid forcing blocks. Our samples consisted of two 5-mm-thick layers of simulated fault gouge with a nominal contact area of 10 by 10 cm^2. Gouge material consisted of soda-lime glass beads with initial particle size between 105 and 149 micrometers. Prior to shearing, we impose a constant fault normal stress of 2 MPa using a servo-controlled load-feedback mechanism and allow the sample to compact. Once the sample has reached a constant layer thickness, the central block is driven down at constant rate of 10 micrometers per second. In tandem, we collect an AE signal continuously at 4 MHz from a piezoceramic sensor embedded in a steel forcing block about 22 mm from the gouge layer The data from this experiment can be used with the deep learning algorithm to train it for future fault property prediction.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Magnetohydrodynamic simulations preliminarily predict the habitability and radio emission of TRAPPIST-1e

TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.

Sun: magnetic fields↗
DOE OSTI2013

Industrial and natural sources of gaseous elemental mercury in the Almadén district (Spain): An updated report on this issue after the ceasing of mining and metallurgical activities in 2003 and major land reclamation works

Two events during the last decade had major environmental repercussions in Almadén town (Spain). First it was the ceasing of activities in the mercury mine and metallurgical facilities in 2003, and then the finalization of the restoration works on the main waste dump in 2008. The combination of both events brought about a dramatic drop in the emissions of gaseous elemental mercury (GEM) to the atmosphere. Although no one would now call the Almadén area as ‘mercury-free’, the GEM levels have fallen beneath international reference safety levels for the first time in centuries. This has been a major breakthrough because in less than one decade the site went from GEM levels in the order of “tens of thousands” to mere “tens” nanogram per cubic meter. Although these figures are per se a remarkable achievement, they do not mark the end of the environmental concerns in the Almadén district. Two other sites remain as potential environmental hazards. (1) The Las Cuevas mercury storage complex, a partially restored ex-mining site where liquid mercury is being stored. The MERSADE Project (LIFE—European Union) has tested the Las Cuevas complex as a potential site for the installation of a future European prototype safe deposit of surplus mercury from industrial activities. Despite restoration works carried out in 2004, the Las Cuevas complex can still be regarded as hotspot of mercury contamination, with high concentrations above 800 μg g{sup −1} Hg{sub soil} and 300 ng m{sup −3} Hg{sub gas}. However, as predicted by air contamination modeling using the ISC-AERMOD software, GEM concentrations fade away in a short distance following the formation of a NW–SE oriented narrow plume extending for a few hundred meters from the complex perimeter. (2) Far more dangerous from the human health perspective is the Almadenejos area, hosting the small Almadenejos village, the so-called Cerco de Almadenejos (CDA; an old metallurgical precinct), and the mines of La Nueva Concepción, La Vieja Concepción and El Entredicho. The CDA is an old metallurgical site that operated between 1794 and 1861, leaving behind a legacy of extremely contaminated soils (mean concentration=4220 μg g{sup −1} Hg) and GEM emissions that in summer can reach levels up to 4,000–5,000 ng m{sup −3}. Thus the CDA remains the sole ‘urban’ site in the district surpassing GEM international reference safety levels. In order to prevent these emissions, the CDA requires immediate action regarding restoration works. These could involve the full removal of soils or their permanent capping to create an impermeable barrier.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

THE IMPACT OF TYPE Ia SUPERNOVA EXPLOSIONS ON HELIUM COMPANIONS IN THE CHANDRASEKHAR-MASS EXPLOSION SCENARIO

In the version of the single-degenerate scenario of Type Ia supernovae (SNe Ia) studied here, a carbon-oxygen white dwarf explodes close to the Chandrasekhar limit after accreting material from a non-degenerate helium (He) companion star. In the present study, we employ the STELLAR GADGET code to perform three-dimensional hydrodynamical simulations of the interaction of the SN Ia ejecta with the He companion star taking into account its orbital motion and spin. It is found that only 2%-5% of the initial companion mass is stripped off from the outer layers of He companion stars due to the supernova (SN) impact. The dependence of the unbound mass (or the kick velocity) on the orbital separation can be fitted to a good approximation by a power law for a given companion model. After the SN impact, the outer layers of a He donor star are significantly enriched with heavy elements from the low-expansion-velocity tail of SN Ia ejecta. The total mass of accumulated SN-ejecta material on the companion surface reaches about {approx}> 10{sup -3} M{sub Sun} for different companion models. This enrichment with heavy elements provides a potential way to observationally identify the surviving companion star in SN remnants. Finally, by artificially adjusting the explosion energy of the W7 explosion model, we find that the total accumulation of SN ejecta on the companion surface is also dependent on the explosion energy with a power-law relation to a good approximation.

79 ASTRONOMY AND ASTROPHYSICS↗