Search NASA⌕ Search

SEARCH · Search NASA

Results for “patterns in complex systems”

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 1,009 records · Page 56

Predicting impacts of chemicals from organisms to ecosystem service delivery: A case study of insecticide impacts on a freshwater lake

Assessing and managing risks of anthropogenic activities to ecological systems is necessary to ensure sustained delivery of ecosystem services for future generations. Ecological models provide a means of quantitatively linking measured risk assessment endpoints with protection goals, by integrating potential chemical effects with species life history, ecological interactions, environmental drivers and other potential stressors. In this work we demonstrate how an ecosystem modeling approach can be used to quantify insecticide-induced impacts on ecosystem services provided by a lake from toxicity data for organism-level endpoints. We used a publicly available aquatic ecosystem model AQUATOX that integrates environmental fate of chemicals and their impacts on food webs in aquatic environments. By simulating a range of exposure patterns, we illustrated how exposure to a hypothetical insecticide could affect aquatic species populations (e.g., recreational fish abundance) and environmental properties (e.g., water clarity) that would in turn affect delivery of ecosystem services. Different results were observed for different species of fish, thus the decision to manage the use of the insecticide for ecosystem services derived by anglers depends upon the favored species of fish. In our hypothetical shallow reservoir, water clarity was mostly driven by changes in food web dynamics, specifically the presence of zooplankton. In contrast to the complex response by fishing value, water clarity increased with reduced insecticide use, which produced a monotonic increase in value by waders and swimmers. Our research clearly showed the importance of considering nonlinear ecosystem feedbacks where the presence of insecticide changed the modeled food-web dynamics in unexpected ways. Our study introduces one of the main advantages of using ecological models for risk assessment, namely the ability to generalize to meaningful levels of organization and to facilitate quantitative comparisons among alternative scenarios and associated trade-offs among them while explicitly accounting for different groups of beneficiaries.

54 ENVIRONMENTAL SCIENCES↗

Structural and metal-insulator transitions in rhenium-based double perovskites via orbital ordering

Re-based double perovskites (DPs) have garnered substantial attention due to their high Curie temperatures (T C ) and display of complex interplay of structural and metal-insulator transitions (MIT). In this report we systematically study the ground-state electronic and structural properties for a family of Re-based DPs A 2 BReO 6 (A=Sr, Ca and B=Cr, Fe), which are related by a common low-energy Hamiltonian, using density functional theory +U calculations. We show that the on-site interaction U of Re induces orbital ordering (denoted C-OO), with each Re site having an occupied d xy orbital and a C-type alternation among d xz /d yz , resulting in an insulating state consistent with experimentally determined insulators Sr 2 CrReO 6 , Ca 2 CrReO 6 , and Ca 2 FeReO 6 . The threshold value of U Re for orbital ordering is reduced by inducing E g octahedral distortions of the same C-type wavelength (denoted C-OD), which serves as a structural signature of the orbital ordering; octahedral tilting also reduces the threshold. The C-OO and the concomitant C-OD are a spontaneously broken symmetry for the Sr-based materials (i.e., a 0 a 0 c - tilt pattern), while not for the Ca-based systems (i.e., a - a - b + tilt pattern). Spin-orbit coupling does not qualitatively change the physics of the C-OO/C-OD, but can induce relevant quantitative changes. We prove that a single set of U Cr , U Fe , U Re capture the experimentally observed metallic state in Sr 2 FeReO 6 and insulating states in other three systems. We predict that the C-OO is the origin of the insulating state in Sr 2 CrReO 6 , and that the concomitant C-OD may be experimentally observed at sufficiently low temperatures (i.e., space group P4 2 /m) in pure samples. Additionally, given our prescribed values of U, we show that the C-OO induced insulating state in Ca 2 CrReO 6 will survive even if the C-OD amplitude is suppressed (e.g., due to thermal fluctuations). The role of the C-OO/C-OD in the discontinuous, temperature driven MIT in Ca 2 FeReO 6 is discussed.

36 MATERIALS SCIENCE↗

Gamow-shell-model description of Li isotopes and their mirror partners

Background: Weakly bound and unbound nuclei close to particle drip lines are laboratories of new nuclear structure physics at the extremes of neutron/proton excess. The comprehensive description of these systems requires an open quantum system framework that is capable of treating resonant and nonresonant many-body states on equal footing. Purpose: In this work, we develop the complex-energy configuration interaction approach to describe binding energies and spectra of selected 5 ≤ $\textit{A}$ ≤ 11 nuclei. Method: We employ the complex-energy Gamow shell model (GSM) assuming a rigid 4 He core. The effective Hamiltonian, consisting of a core-nucleon Woods-Saxon potential and a simplified version of the Furutani-Horiuchi-Tamagaki interaction with the mass-dependent scaling, is optimized in the $\textit{sp}$ space. To diagonalize the Hamiltonian matrix, we employ the Davidson method and the Density Matrix Renormalization Group technique. Results: We report our optimized GSM Hamiltonian offers a good reproduction of binding energies and spectra with the root-mean-square (rms) deviation from experiment of 160 keV. Since the model performs well when used to predict known excitations that have not been included in the fit, it can serve as a reliable tool to describe poorly known states. A case in point is our prediction for the pair of unbound mirror nuclei 10 Li – 10 N in which a huge Thomas-Ehrman shift dramatically alters the pattern of low-energy excitations. Conclusion: The new model will enable comprehensive studies of structure and reactions aspects of light drip-line nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Additive manufacturing of amorphous metal soft magnetic composites

Soft magnet alloys are used as magnetic cores for electric motors, transformers, wind turbines and other power generation systems. Soft magnetic cores are expensive and time consuming to manufacture in the complex shapes required for next-generation devices using conventional press and sinter powder metallurgy. The objective of this effort is to additively manufacture high performance soft magnets, with reduced cost and reduced material waste and 10x lower energy (core) loss at the high operating frequencies of many electric machines. Laser powder bed fusion additive manufacturing is used as the fabrication method. Electrical steel and amorphous alloys atomized powders are used as a feedstock materials. Magnetic cores are printed in topology optimized structures such as the Hilbert curve because this has been shown to reduce energy losses by minimizing the eddy currents that circulate within the magnet at high frequencies. In our project, we succeed in printing FeSi 3.5wt% and FeSi 6.5wt% electrical steels, and iron-based soft magnetic amorphous alloys in the shape of Hilbert and Peano curve topology optimized structures. We found that the Peano curve has a higher cut-off frequency than the Hilbert curve, and that amorphous alloys have high cut-off frequencies and higher mechanical hardness than electrical steels. Processing conditions such as laser power, scan speed, and hatching pattern were optimized to achieve high density prints, and optimize magnetic performance. We find that the printing of amorphous alloy soft magnetic cores may be technoeconomically feasible for large scale applications such as transformers, for which supply chain issues and the labor costs of manual fabrication of magnetic cores is prohibitive in some cases.

36 MATERIALS SCIENCE↗

Computationally Guided and Experimentally Validated Design of Custom Chelators for Critical Mineral Recovery

Selective, high throughput separation of target critical metals from complex environments such as fly ash leachates and mining process streams presents a significant challenge for economical production. Custom chelators and sorbents are an attractive technology for selective metal extraction, however it can be difficult to predict their performance, and significant experimental efforts are often required to develop chelating technologies. Here, we present a computational strategy focused on modelling chelator-metal binding interactions and benchmark these results versus experimental data. A computational pipeline combining forcefield, semiempirical, and meta-GGA methods with a thermodynamic framework optimized for error cancellation has been developed to predict binding energies of chelator complexes towards critical mineral recovery applications. This approach, originally validated on [2.2.2] cryptates binding mono- and divalent cations, demonstrated robust predictive capabilities with an R2 of 0.850 against experimental aqueous binding energies. The workflow includes metadynamics for exploring high-dimensional potential energy surfaces and a cluster-continuum model for accurate yet computationally efficient solvation modeling. Error cancellation between solvation energies of free and chelator-coordinated ions enables faster convergence, even with finite cluster sizes. Initial studies on the cryptates revealed consistent metal-ligand coordination patterns, with systematic variations influenced by ion size and charge, highlighting key structural features linked to binding selectivity. Further studies of a proprietary chelator have resulted in identification of previously unreported selectivity towards economically significant metals, which in-house experiments have confirmed, demonstrating the feasibility of this approach. By applying this methodology to new chelators targeting critical minerals such as lithium, cobalt, nickel and other strategic metals, we aim to accelerate the discovery of next-generation chelators for efficient recovery, recycling, and separation processes. This computational framework serves as the backbone of a high-throughput design pipeline tailored for sustainable resource utilization and may be applied to a wide range of systems to meet experimental needs.

computational materials↗

Multitemporal and Multiscaled Fractal Analysis of Landsat Satellite Data Using the Image Characterization and Modeling System (ICAMS)

The Image Characterization And Modeling System (ICAMS) is a public domain software package that is designed to provide scientists with innovative spatial analytical tools to visualize, measure, and characterize landscape patterns so that environmental conditions or processes can be assessed and monitored more effectively. In this study ICAMS has been used to evaluate how changes in fractal dimension, as a landscape characterization index, and resolution, are related to differences in Landsat images collected at different dates for the same area. Landsat Thematic Mapper (TM) data obtained in May and August 1993 over a portion of the Great Basin Desert in eastern Nevada were used for analysis. These data represent contrasting periods of peak "green-up" and "dry-down" for the study area. The TM data sets were converted into Normalized Difference Vegetation Index (NDVI) images to expedite analysis of differences in fractal dimension between the two dates. These NDVI images were also resampled to resolutions of 60, 120, 240, 480, and 960 meters from the original 30 meter pixel size, to permit an assessment of how fractal dimension varies with spatial resolution. Tests of fractal dimension for two dates at various pixel resolutions show that the D values in the August image become increasingly more complex as pixel size increases to 480 meters. The D values in the May image show an even more complex relationship to pixel size than that expressed in the August image. Fractal dimension for a difference image computed for the May and August dates increase with pixel size up to a resolution of 120 meters, and then decline with increasing pixel size. This means that the greatest complexity in the difference images occur around a resolution of 120 meters, which is analogous to the operational domain of changes in vegetation and snow cover that constitute differences between the two dates.

Quattrochi, Dale A.↗

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗

In Situ Inference for Earth System Predictability

An understanding of future evolution in precipitation extremes is critical to numerous DOE mission questions. Extreme events are by nature short time-scale events that are difficult to diagnose in available model data. Accurate modeling of extreme events necessarily requires high spatial resolution at the storm scale locally. However, the environment in which storms grow is dependent on global, remote, processes. These complex spatiotemporal relationships are impossible to diagnose at resolutions required to accurately model storms responsible for extreme precipitation. At exascale, climate simulations will produce results at fine enough resolution to investigate these relationships. However, the resulting data from these simulations will be far too large to save for post-simulation analysis. We advocate for fitting statistical models inside the simulations as they run, a context known as in situ, which will facilitate scientific investigations using the full fine-scale data stream. Figure 1 shows an example of the type of model we could consider, a Bayesian hierarchical spatial regression model. Precipitation extremes at each grid cell are modeled using extreme value distributions. Since extremes are rare, fitting models to individual grid cells can result in high variance and poor estimates. Instead, the model can be made more robust by smoothing the parameters of the extreme value model across space. Additionally, the parameters themselves can be functionally linked to other variables elsewhere in the simulation. Thus, we can use the fine-scale data to build more robust models for extremes that link extreme behavior to other climate patterns.

54 ENVIRONMENTAL SCIENCES↗

Comparing multi-source urban flood indicators: satellite, simulation, and citizen-reported data

Urban flooding arises from complex mechanisms, making it challenging to capture accurately with a single detection method. This study evaluates three complementary approaches to detect flooding across three Chicago neighborhoods: (i) Sentinel-1 synthetic aperture radar (SAR), offering weather-independent, high-resolution (10 m) imagery of surface inundation; (ii) the storm water management model (SWMM), simulating combined sewer overflow and drainage performance; and (iii) citizen-generated 311 service requests, capturing observed flooding impacts. By analyzing six storms ranging from severe to mild, we examine how each source uniquely contributes to identifying urban flood events. SAR imagery effectively identifies standing water but can miss brief flooding due to satellite revisit constraints. SWMM provides detailed insights into system-wide drainage behavior yet may underestimate localized street-level flooding. Meanwhile, 311 calls reflect real-world flooding impacts but are vulnerable to underreporting. Statistical overlap analysis highlights chronic flood hotspots repeatedly identified across multiple detection methods, indicating persistent infrastructure and topographic vulnerabilities. Temporal analysis further reveals that while SWMM flooding aligns closely with rainfall peaks, 311 calls typically precede or persist beyond these peaks. Our findings emphasize the value of using satellite observations, hydrological modeling, and resident-reported data in a complementary manner to better interpret patterns in flood timing, severity, and spatial distribution—providing insights that can inform targeted infrastructure improvements and contribute to urban flood resilience planning.

311↗

Simulation of the Multi-Wake Evolution of Two Sandia National Labs/National Rotor Testbed Turbines Operating in a Tandem Layout

The future of wind power systems deployment is in the form of wind farms comprised of scores of such large turbines, most likely at offshore locations. Individual turbines have grown in span from a few tens of meters to today’s large turbines with rotor diameters that dwarf even the largest commercial aircraft. These massive dynamical systems present unique challenges at scales unparalleled in prior applications of wind science research. Fundamental to this effort is the understanding of the wind turbine wake and its evolution. Furthermore, the optimization of the entire wind farm depends on the evolution of the wakes of different turbines and their interactions within the wind farm. In this article, we use the capabilities of the Common ODE Framework (CODEF) model for the analysis of the effects of wake–rotor and wake-to-wake interactions between two turbines situated in a tandem layout fully and partially aligned with the incoming wind. These experiments were conducted in the context of a research project supported by the National Rotor Testbed (NRT) program of Sandia National Labs (SNL). Results are presented for a layout which emulates the turbine interspace and relative turbine emplacement found at SNL’s Scaled Wind Technologies Facility (SWiFT), located in Lubbock, Texas. The evolution of the twin-wake interaction generates a very rich series of secondary transitions in the vortex structure of the combined wake. These ultimately affect the wake’s axial velocity patterns, altering the position, number, intensity, and shape of localized velocity-deficit zones in the wake’s cross-section. This complex distribution of axial velocity patterns has the capacity to substantially affect the power output, peak loads, fatigue damage, and aeroelastic stability of turbines located in subsequent rows downstream on the farm.

Baruah, Apurva (ORCID:0000000252354068)↗

Reactions of Air Transport Flight Crews to Displays of Weather During Simulated Flight

Display of information in the cockpit has long been a challenge for aircraft designers. Given the limited space in which to present information, designers have had to be extremely selective about the types and amount of flight related information to present to pilots. The general goal of cockpit display design and implementation is to ensure that displays present information that is timely, useful, and helpful. This suggests that displays should facilitate the management of perceived workload, and should allow maximal situation awareness. The formatting of current and projected weather displays represents a unique challenge. As technologies have been developed to increase the variety and capabilities of weather information available to flight crews, factors such as conflicting weather representations and increased decision importance have increased the likelihood for errors. However, if formatted optimally, it is possible that next generation weather displays could allow for clearer indications of weather trends such as developing or decaying weather patterns. Important issues to address include the integration of weather information sources, flight crew trust of displayed weather information, and the teamed reactivity of flight crews to displays of weather. Past studies of weather display reactivity and formatting have not adequately addressed these issues; in part because experimental stimuli have not approximated the complexity of modern weather displays, and in part because they have not used realistic experimental tasks or participants. The goal of the research reported here was to investigate the influence of onboard and NEXRAD agreement, range to the simulated potential weather event, and the pilot flying on flight crew deviation decisions, perceived workload, and perceived situation awareness. Fifteen pilot-copilot teams were required to fly a simulated route while reacting to weather events presented in two graphical formats on a separate visual display. Measures of flight crew reactions included performance-based measures such as deviation decision accuracy, and judgment-based measures such as perceived decision confidence, workload, situation awareness, and display trust. Results demonstrated that pilots adopted a conservative reaction strategy, often choosing to deviate from weather rather than ride through it. When onboard and NEXRAD displays did not agree, flight crews reacted in a complex manner, trusting the onboard system more but using the NEXRAD system to augment their situation awareness. Distance to weather reduced situation awareness and heightened workload levels. Overall, flight crews tended to adopt a participative leadership style marked by open communication. These results suggest that future weather displays should exploit the existing benefits of NEXRAD presentation for situation awareness while retaining the display structure and logic inherent in the onboard system.

Bliss, James P.↗

Experimental warming and its legacy effects on root dynamics following two hurricane disturbances in a wet tropical forest

Abstract Tropical forests are expected to experience unprecedented warming and increases in hurricane disturbances in the coming decades; yet, our understanding of how these productive systems, especially their belowground component, will respond to the combined effects of varied environmental changes remains empirically limited. Here we evaluated the responses of root dynamics (production, mortality, and biomass) to soil and understory warming (+4°C) and after two consecutive tropical hurricanes in our in situ warming experiment in a tropical forest of Puerto Rico: Tropical Responses to Altered Climate Experiment (TRACE). We collected minirhizotron images from three warmed plots and three control plots of 12 m 2 . Following Hurricanes Irma and María in September 2017, the infrared heater warming treatment was suspended for repairs, which allowed us to explore potential legacy effects of prior warming on forest recovery. We found that warming significantly reduced root production and root biomass over time. Following hurricane disturbance, both root biomass and production increased substantially across all plots; the root biomass increased 2.8‐fold in controls but only 1.6‐fold in previously warmed plots. This pattern held true for both herbaceous and woody roots, suggesting that the consistent antecedent warming conditions reduced root capacity to recover following hurricane disturbance. Root production and mortality were both related to soil ammonium nitrogen and microbial biomass nitrogen before and after the hurricanes. This experiment has provided an unprecedented look at the complex interactive effects of disturbance and climate change on the root component of a tropical forested ecosystem. A decrease in root production in a warmer world and slower root recovery after a major hurricane disturbance, as observed here, are likely to have longer‐term consequences for tropical forest responses to future global change.

54 ENVIRONMENTAL SCIENCES↗

Vertical Cloud Climatology During TC4 Derived from High-Altitude Aircraft Merged Lidar and Radar Profiles

Aircraft lidar works by shooting laser pulses toward the earth and recording the return time and intensity of any of the light returning to the aircraft after scattering off atmospheric particles and/or the Earth s surface. The scattered light signatures can be analyzed to tell the exact location of cloud and aerosol layers and, with the aid of a few optical assumptions, can be analyzed to retrieve estimates of optical properties such as atmospheric transparency. Radar works in a similar fashion except it sends pulses toward earth at a much larger wavelength than lidar. Radar records the return time and intensity of cloud or rain reflection returning to the aircraft. Lidar can measure scatter from optically thin cirrus and aerosol layers whose particles are too small for the radar to detect. Radar can provide reflection profiles through thick cloud layers of larger particles that lidar cannot penetrate. Only after merging the two instrument products can accurate measurements of the locations of all layers in the full atmospheric column be achieved. Accurate knowledge of the vertical distribution of clouds is important information for understanding the Earth/atmosphere radiative balance and for improving weather/climate forecast models. This paper describes one such merged data set developed from the Tropical Composition, Cloud and Climate Coupling (TC4) experiment based in Costa Rica in July-August 2007 using the nadir viewing Cloud Physics Lidar (CPL) and the Cloud Radar System (CRS) on board the NASA ER-2 aircraft. Statistics were developed concerning cloud probability through the atmospheric column and frequency of the number of cloud layers. These statistics were calculated for the full study area, four sub-regions, and over land compared to over ocean across all available flights. The results are valid for the TC4 experiment only, as preferred cloud patterns took priority during mission planning. The TC4 Study Area was a very cloudy region, with cloudy profiles occurring 94 percent of the time during the ER-2 flights. One to three cloud layers were common, with the average calculated at 2.03 layers per profile. The upper troposphere had a cloud frequency generally over 30%, reaching 42 percent near 13 km during the study. There were regional differences. The Caribbean was much clearer than the Pacific regions. Land had a much higher frequency of high clouds than ocean areas. One region just south and west of Panama had a high probability of clouds below 15 km altitude with the frequency never dropping below 25% and reaching a maximum of 60% at 11-13 km altitude. These cloud statistics will help characterize the cloud volume for TC4 scientists as they try to understand the complexities of the tropical atmosphere.

Hlavka, Dennis↗

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗

Spectra and characteristics of somatic mutations induced by ionizing radiation in hematopoietic stem cells

Spectra and frequencies of spontaneous and X-ray-induced somatic mutations were revealed with mouse long-term hematopoietic stem cells (LT-HSCs) by whole-genome sequencing of clonal cell populations propagated in vitro from single isolated LT-HSCs. SNVs and small indels were the most common types of somatic mutations, and increased up to twofold to threefold by whole-body X-irradiation. Base substitution patterns in the SNVs suggested a role of reactive oxygen species in radiation mutagenesis, and signature analysis of single base substitutions (SBS) revealed a dose-dependent increase of SBS40. Most of spontaneous small deletions were shrinkage of tandem repeats, and X-irradiation specifically induced small deletions out of tandem repeats (non-repeat deletions). Presence of microhomology sequences in non-repeat deletions suggested involvement of microhomology mediated end-joining repair mechanisms as well as nonhomologous end-joining in radiation-induced DNA damages. We also identified multisite mutations and structural variants (SV), i.e., large indels, inversions, reciprocal translocations, and complex variants. The radiation-specificity of each mutation type was evaluated from the spontaneous mutation rate and the per-Gy mutation rate estimated by linear regression, and was highest with non-repeat deletions without microhomology, followed by those with microhomology, SV except retroelement insertions, and multisite mutations; these types were thus revealed as mutational signatures of ionizing radiation. Further analysis of somatic mutations in multiple LT-HSCs indicated that large fractions of postirradiation LT-HSCs originated from single LT-HSCs that survived the irradiation and then expanded in vivo to confer marked clonality to the entire hematopoietic system, with varying clonal expansion and dynamics depending on radiation dose and fractionation.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Trainable Gene Regulation Networks with Applications to Drosophila Pattern Formation

This chapter will very briefly introduce and review some computational experiments in using trainable gene regulation network models to simulate and understand selected episodes in the development of the fruit fly, Drosophila melanogaster. For details the reader is referred to the papers introduced below. It will then introduce a new gene regulation network model which can describe promoter-level substructure in gene regulation. As described in chapter 2, gene regulation may be thought of as a combination of cis-acting regulation by the extended promoter of a gene (including all regulatory sequences) by way of the transcription complex, and of trans-acting regulation by the transcription factor products of other genes. If we simplify the cis-action by using a phenomenological model which can be tuned to data, such as a unit or other small portion of an artificial neural network, then the full transacting interaction between multiple genes during development can be modelled as a larger network which can again be tuned or trained to data. The larger network will in general need to have recurrent (feedback) connections since at least some real gene regulation networks do. This is the basic modeling approach taken, which describes how a set of recurrent neural networks can be used as a modeling language for multiple developmental processes including gene regulation within a single cell, cell-cell communication, and cell division. Such network models have been called "gene circuits", "gene regulation networks", or "genetic regulatory networks", sometimes without distinguishing the models from the actual modeled systems.

Mjolsness, Eric↗

Evaluating Offshore Infrastructure Integrity

Drilling in the offshore environment involves a complex network of infrastructure including pipelines, platforms, rigs, subsea installations, ports, and terminals. Government and industry partners have developed this network over many decades and it remains a critical part of the United States (U.S.) energy portfolio. Many of the major components of this system have been designed with a 20- to 30-year lifespan, yet consistent and growing energy demands support the need to extend the design life of existing infrastructure or repurpose it for secondary needs (i.e. enhanced oil recovery, carbon storage, and new wells). As a result, a growing portion of the offshore infrastructure in the U.S. is approaching or has exceeded its original design life. A critical step in ensuring the continued safe and effective operation of offshore infrastructure is developing a comprehensive understanding of the state of offshore infrastructure and the factors that effect it. The purpose of this project is to assess the current state of existing infrastructure and identify the factors involved in infrastructure degradation through the development and application of big data analytics, machine learning, and advanced spatio-temporal analysis. The project leverages existing data at NETL and combines it with new information on offshore oil and gas structures and the ambient offshore environment in an effort to identify patterns associated with infrastructure integrity. Building on the identified trends and patterns, this project incorporates exploratory analytics and spatial analysis tools in conjunction with machine learning and statistical models to characterize the condition of existing platforms in the offshore environment and predict their risk of failure.

02 PETROLEUM↗