SEARCH · Engineering Papers

Results for “data replication”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Metabolomic, photoprotective, and photosynthetic acclimatory responses to post‐flowering drought in sorghum

Abstract Climate change is globally affecting rainfall patterns, necessitating the improvement of drought tolerance in crops. Sorghum bicolor is a relatively drought‐tolerant cereal. Functional stay‐green sorghum genotypes can maintain green leaf area and efficient grain filling during terminal post‐flowering water deprivation, a period of ~10 weeks. To obtain molecular insights into these characteristics, two drought‐tolerant genotypes, BTx642 and RTx430, were grown in replicated control and terminal post‐flowering drought field plots in California's Central Valley. Photosynthetic, photoprotective, and water dynamics traits were quantified and correlated with metabolomic data collected from leaves, stems, and roots at multiple timepoints during control and drought conditions. Physiological and metabolomic data were then compared to longitudinal RNA sequencing data collected from these two genotypes. The unique metabolic and transcriptomic response to post‐flowering drought in sorghum supports a role for the metabolite galactinol in controlling photosynthetic activity through regulating stomatal closure in post‐flowering drought. Additionally, in the functional stay‐green genotype BTx642, photoprotective responses were specifically induced in post‐flowering drought, supporting a role for photoprotection in the molecular response associated with the functional stay‐green trait. From these insights, new pathways are identified that can be targeted to maximize yields under growth conditions with limited water.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

HopBox: An image analysis pipeline to characterize hop cone morphology

Abstract Hop cone morphology can influence picking and drying ability, and color can impact consumer preference and may be indicative of quality. However, these characteristics are not generally evaluated in hop breeding programs due to the tedious nature of trait quantification and the extensive variation among cones within a genotype. We developed the HopBox, which is a simply constructed light box with a camera mount, and a publicly available image processing pipeline that identifies hop cones within color‐corrected images, reads a QR code within the image, and outputs data on hop cone length, width, area, perimeter, openness, weight, color, and density. The trained model was applied to images of 500 cones each from 15 replicated advanced hop genotypes from the USDA‐ARS breeding program in Prosser, Washington. Analysis of variance revealed significant ( p < 0.001) differences between genotypes for all traits measured, enabling breeders to discriminate between genotypes for selection purposes. Broad sense heritability for all traits ranged from 0.23 to 0.59. A random sampling of hop cones from the complete dataset revealed that imaging only 5–10 cones adequately captured genotypic variation and provided acceptable rank correlations ( r s > 0.75); however, increasing the sample size to 30 provided optimal precision. Instructions for constructing a HopBox and the code for the analysis pipeline are publicly available online and have wide applicability for hop breeding and research.

Altendorf, Kayla R.↗

DOE OSTI2020

Liquid–Solid Mixtures of Ga Metal Infused with Cu Microparticles and Nanoparticles for Microscale and Nanoscale Patterning of Solid Metals at Room Temperature

This paper demonstrates and characterizes a simple ink for nanopatterning of solid metallic structures under ambient conditions by taking advantage of the low melting point of gallium and its affinity to form intermetallics with other metals. Bare copper micro- and nanoparticles readily mix with liquid gallium near room temperature to form a paste that ultimately amalgamates into a dense solid when mixed at the appropriate concentration ratios. The paste has sufficient shelf life that can be extruded from a nozzle at modest temperatures to form solid three-dimensional (3D) shapes. Additionally, the paste can be molded at room temperature to replicate feature sizes ranging from a few millimeters down to hundreds of nanometers. In situ X-ray diffraction (XRD) and thermo-mechanical analysis (TMA) data show gallium and copper readily interdiffuse to form the thermodynamically expected intermetallic phase. We describe the capabilities and limitations of a simple way to pattern solid metals in an additive fashion (syringe-based extrusion) and with high resolution (molding) at or near room temperature. Here, the use of a paste that solidifies provides a novel route for 3D printing of solid metals at ambient temperatures as well as the creation of micro- and nanostructured metallic surfaces that may be useful for optics, non-wetting surfaces, or electronic microcomponents.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Discovering useful genetic variation in the seed parent gene pool for sorghum improvement

Multi-parent populations contain valuable genetic material for dissecting complex, quantitative traits and provide a unique opportunity to capture multi-allelic variation compared to the biparental populations. A multi-parent advanced generation inter-cross (MAGIC) B-line (MBL) population composed of 708 F 6 recombinant inbred lines (RILs), was recently developed from four diverse founders. These selected founders strategically represented the four most prevalent botanical races (kafir, guinea, durra, and caudatum) to capture a significant source of genetic variation to study the quantitative traits in grain sorghum [Sorghum bicolor (L.) Moench]. MBL was phenotyped at two field locations for seven yield-influencing traits: panicle type (PT), days to anthesis (DTA), plant height (PH), grain yield (GY), 1000-grain weight (TGW), tiller number per meter (TN) and yield per panicle (YPP). High phenotypic variation was observed for all the quantitative traits, with broad-sense heritabilities ranging from 0.34 (TN) to 0.84 (PH). The entire population was genotyped using Diversity Arrays Technology (DArTseq), and 8,800 single nucleotide polymorphisms (SNPs) were generated. A set of polymorphic, quality-filtered markers (3,751 SNPs) and phenotypic data were used for genome-wide association studies (GWAS). We identified 52 marker-trait associations (MTAs) for the seven traits using BLUPs generated from replicated plots in two locations. We also identified desirable allelic combinations based on the plant height loci (Dw1, Dw2, and Dw3), which influences yield related traits. Additionally, two novel MTAs were identified each on Chr1 and Chr7 for yield traits independent of dwarfing genes. We further performed a multi-variate adaptive shrinkage analysis and 15 MTAs with pleiotropic effect were identified. The five best performing MBL progenies were selected carrying desirable allelic combinations. Since the MBL population was designed to capture significant diversity for maintainer line (B-line) accessions, these progenies can serve as valuable resources to develop superior sorghum hybrids after validation of their general combining abilities via crossing with elite pollinators. Further, newly identified desirable allelic combinations can be used to enrich the maintainer germplasm lines through marker-assisted backcross breeding.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

A Physics-Based Data-Driven Approach for Modeling of Environmental Degradation in Elastomers

Abstract Elastomers are now commonly used in a number of industries, including aerospace, structure, transportation, shipbuilding, and automotive, due to their excellent workability, formability, and flexibility. During their activity, elastomers are subjected to harsh environmental conditions, which decreases their resilience. False predictions made early in their lives can have major financial and environmental implications. Elastomers’ performance and properties, such as strength, durability, and density, are influenced by chemical changes in these materials, known as degradation, which occurs over time. This process can alter the morphology of a polymer matrix as well as cause chain scission and cross-linking, resulting in different behaviors than that of the unaged material. To demonstrate the effect of thermaloxidative aging on the mechanical behavior of elastomers, several experimental and theoretical models have been proposed. In view of the large volume of experimental data available on micro-structural evolution in the course of aging, we propose a physics-based data-driven approach to overcome the shortcomings of both phenomenological and micro-mechanical models. This work presents a novel thermodynamically consistent, multiagent machine-learned model for predicting the constitutive behavior of cross-linked elastomers during environmental aging, such as thermo-oxidative and hydrolytic aging for various states of deformation. Single mechanism degradation changes the polymer matrix over time where it is causing chain scission, reduction of cross-links, and morphology change. To capture the idealized Mullins effect and permanent set due to the effect of single aging mechanisms on nonlinear mechanical responses of elastomers, we propose a data-driven model for simulating inelastic elements in a polymer matrix. By using a sequential order reduction, we were able to reduce the 3D stress-strain tensor mapping problem to a small number of super-constrained 1D mapping problems. To systematically classify such mapping problems into a few categories, an assembly of multiple replicated conditional neural network learning agents (L-agents) is used based on our recent work. Each category is represented by a different type of agent. The effect of deformation history, aging time, and aging temperature is captured by this model. The model is validated using a broad collection of data, ranging from our experimental results to data from the literature. In addition, thermodynamic consistency and frame independence are investigated. The most significant achievements of this model are its precision, simplicity, and prediction of inelasticity under various states of deformation. The model’s accuracy and simplicity make it a good option for commercial and industrial applications. Conveniently, due to the model modular nature, it can be expanded in the future to include viscoelasticity and non-isotropic formation for better precision.

Ghaderi, Aref↗

DOE OSTI2021

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Single- and multi-component chiral supraparticles as modular enantioselective catalysts

Nanoscale biological assemblies exemplified by exosomes, endosomes and capsids, play crucial roles in all living systems. Supraparticles (SP) from inorganic nanoparticles (NPs) replicate structural characteristics of these bioassemblies, but it is unknown whether they can mimic their biochemical functions. Here, we show that chiral ZnS NPs self-assemble into 70–100 nm SPs that display sub-nanoscale porosity associated with interstitial spaces between constituent NPs. Similarly to photosynthetic bacterial organelles, these SPs can serve as photocatalysts, enantioselectively converting L- or D-tyrosine (Tyr) into dityrosine (diTyr). Experimental data and molecular dynamic simulations indicate that the chiral bias of the photocatalytic reaction is associated with the chiral environment of interstitial spaces and preferential partitioning of enantiomers into SPs, which can be further enhanced by co-assembling ZnS with Au NPs. Besides replicating a specific function of biological nanoassemblies, these findings establish a path to enantioselective oxidative coupling of phenols for biomedical and other needs.

Science & Technology - Other Topics↗

DOE OSTI2019

Posterior samples of the parameters of binary black holes from Advanced LIGO, Virgo’s second observing run

This paper presents a parameter estimation analysis of the seven binary black hole mergers—GW170104, GW170608, GW170729, GW170809, GW170814, GW170818, and GW170823—detected during the second observing run of the Advanced LIGO and Virgo observatories using the gravitational-wave open data. We describe the methodology for parameter estimation of compact binaries using gravitational-wave data, and we present the posterior distributions of the inferred astrophysical parameters. We release our samples of the posterior probability density function with tutorials on using and replicating our results presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

EMT Simulation of Large PV Plant & Power Grid for Disturbance Analysis

One of the challenges faced during the analysis of various disturbances with inverter-based resources is the ability to replicate the same disturbance in simulations. The main issue arises from lack of electromagnetic transient (EMT) models and adequate fidelity of models. To overcome this problem, in this paper, EMT simulation of power grid and one of the affected photovoltaic (PV) plants is performed to replicate the observations during the Angeles Forest disturbance in 2018. The processes to (i) develop the EMT model of the power grid from traditional transient stability (TS) data, and (ii) develop high-fidelity EMT model of PV plant from data collected on a PV plant are described in the paper. The requirements in the high-fidelity EMT model of the PV plant are also described. Finally, simulation of the Angeles Forest disturbance is performed and similar simulation results (as compared to the observed disturbances) are observed.

Debnath, Suman↗

DOE OSTI2014

SU-E-J-65: Evaluation of a Radiation-Induced Cell Proliferation Probability Formula Using Monte Carlo Simulation

Purpose: To evaluate the analytic formula of the cell death probability after single fraction dose. Methods: Cancer cells endlessly divide, but radiation causes the cancer cells to die. Not all cells die right away after irradiation. Instead, they continue dividing for next few cell cycles before they stop dividing and die. At the end of every cell cycle, the cell decides if it undertakes the mitotic process with a certain probability, Pdiv, which is altered by the radiation. Previously, by using a simple analytic model of radiobiology experiments, we obtained a formula of Pdeath (= 1 − Pdiv). A question is if the proposed probability can reproduce the well-known survival data of the LQ model. In this study, we evaluated the formula by doing a Monte Carlo simulation of the cell proliferation process. Starting with Ns seed cells, the cell proliferation process was simulated for N generations or until all cells die. We counted the number of living cells at the end. Assuming that the cell colony survived when more than Nc cells were still alive, the surviving fraction S was estimated. We compared the S vs. dose, or S-D curve, with the LQ model. Results: The results indicated that our formula does not reproduce the experimentally observed S-D curve without selecting appropriate α and α/β. With parameter optimization, there was a fair agreement between the MC result and the LQ curve of dose lower than 20Gy. However, the survival fraction of MC decreased much faster in comparison to the LQ data for doses higher than 20 Gy. Conclusion: This study showed that the previously derived probability of cell death per cell cycle is not sufficiently accurate to replicate common radiobiological experiments. The formula must be modified by considering its cell cycle dependence and some other unknown effects.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

SPRUCE Deep Peat Microbial Diversity, CO2 and CH4 Production in Response to Nutrient, Temperature, and pH Treatments during Incubation Studies

This dataset contains the results of a microcosm incubation study on deep peat collected from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experimental site within the Marcell Experimental Forest in northern Minnesota, USA in September 2014. Microcosms were monitored for carbon dioxide (CO2) and methane (CH4) production, and microbial community dynamics were assessed using quantitative polymerase chain reaction (qPCR) and amplicon sequencing. The experiment was designed with a full factorial design with elevated temperature, nitrogen (N), (P), and pH treatments was used with samples from each transect serving replicates. In all, 96 microcosms were constructed to account for the 16 treatment combinations (Nitrogen x Phosphorous x pH x temperature), 2 time points, and 3 replicates. Temperature treatments were 6 degrees C, to mimic the SPRUCE ambient plot temperatures, and 15 degrees C to mimic the SPRUCE +9 degrees C treatment. These data show that one year of deep peat warming, microbial community structure and abundance of bacterial, archaeal, fungal, and methanogenic populations showed strong vertical stratification across the peat depth profile yet no clear response to the temperature treatments (Kluber et al, 2020). In an effort to identify factors that may be limiting decomposition and microbial community change in deep peat, we conducted an anerobic microcosm incubation of deep peat (150-200 cm depth). This dataset contains one data file in comma-separate values (*.csv) format. Additional metadata are provided: a data dictionary and a file-level metadata file in comma-separate values (.csv) format and a user guide in PDF (*.pdf) format.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Lasers for the observation of multiple order nuclear reactions

Nuclear reaction rates become nonlinear with respect to flux (cm −2 s −1 ) in extreme environments such as those found during stellar nucleosynthesis and terrestrial nuclear detonations. To observe these effects directly in the laboratory, extremely high particle fluences (cm −2 ) are necessary but not sufficient. Reactor-based neutron sources, such as the Institut Laue-Langevin’s high-flux neutron reactor, were previously the closest to meeting this challenge, albeit over ∼hour time scales. In ultra-high flux environments, where multiple reactions occur on picosecond time scales, nuclei are unable to return to their ground states between reactions; consequently, reactions take place on excited nuclei. To accurately model high-flux environments, data on the cross-sections of excited nuclear states are required, which differ significantly from those of ground states due to spin/parity effects. In order to replicate these effects in the laboratory, short high-fluence pulses on the order of the lifetime of a typical nuclear excited state (generally ≲1 ns) are required. Particle beams generated by high-intensity lasers are uniquely positioned to meet this need with the potential to produce fluences of 10 17 protons/cm 2 and 10 22 neutrons/cm 2 over a few pico-seconds or less. In addition to providing a quantitative analysis of the rates of multiple rapid reactions in general, the present work examines a number of laser-based experiments that could be conducted in the near future to observe multiple rapid reactions for laboratory-based astrophysics and the measurement of exotic cross-sections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Status of Research Activities on LVDTs: Modeling and Simulation Development and Performance Characterization of Commercial LVDTs

FY-22 research focused on two major areas. The first area was modeling and simulation with linked testing that will advance the Linear Variable Differential Transformers (LVDT) technology for use by stakeholders requiring LVDTs in upcoming irradiation tests. The second area focused on testing LVDTs purchased from U.S. suppliers that were investigated in FY-21. These cross-cutting development activities will ensure stakeholders have the current state-of-the-art LVDT-based technologies for deployment in future irradiation tests. Two stand-alone reports are included in the appendixes. The purpose of this report is to provide a convenient body of work to access the reports and to evaluate the data and conclusions provided. Appendix A, “Testing of the RDP Translation Transducer,” contains a complete evaluation of an RDP Electrosense linear transducer that has the capability of replicating the form, fit, and function of LVDTs currently being used at Idaho National Laboratory (INL). Appendix B is an extract from the report, “Advanced Sensors and Instrumentation Program at Boise State University (BSU), August 2022, Linear Variable Differential Transformers.” This report contains modeling and testing results conducted at BSU and INL.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Cooperative Automated Cohort Driving on Connected Infrastructure, Arterial Roadways, and Highways: Final Project Demonstration and System-of-Systems Model Correlation

This project seeks to synergize vehicle automated driving and connectivity data to improve mobility and energy efficiency of groups of mixed vehicles operating in close proximity (vehicle cohort) on various infrastructure. A custom cellular communication network links vehicles operating as a cohort with infrastructure to a centralized system-of-systems digital twin with an AI-based optimal behavior planner. The data contained in this set are from final testing and technology demonstrations to U.S. Department of Energy staff at the American Center for Mobility. The data contain single-lane, single-light scenarios; multi-lane, multi-light arterial scenarios; and limited-access highway scenarios. All test cases were derived from simulations and replicated on the test track. The project employed two and four light-duty vehicles with connectivity and drive automation for the testing. The baseline scenario without connectivity was run under the control of the system-of-systems centralized planner but operating each vehicle with an intelligent driver model controlling the velocity, lane utilization, and vehicle gap. This was to ensure the highest compatibility with the simulation in terms of dynamic behavior. The connected cohort case utilized AI optimization to perform coordinated and cooperative control for energy, as well as safe, comfortable behavior for the cohort. The dataset is appropriately named with unconnected and connected designations, with comparisons sharing the same run index number. The included PowerPoint and PDF files describe the test setup and provide an overview of results from the project. ![image](de-EE0009209_March_2023_Data_Arterial_Scenario_Results.png)

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2015

Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg2 CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters Ω m , σ 8 , and w, and replicating the galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of ≤ 5%, and compute the likelihood in the three-dimensional parameter space (Ω m , σ 8 , w) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain w without external data. When the power spectrum and peak counts are combined, the area of the error “banana” in the (Ω m , σ 8 ) plane reduces by a factor of ≈ two, compared to using the power spectrum alone. For a flat Λ cold dark matter model, combining both statistics, we obtain the constraint σ 8 (Ω m /0.27)0.63 = 0.85 +0.03 -0.03 .

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations

Abstract Motivation Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data. Results We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization. Availability and implementation It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Du, Lin (ORCID:0009000289869812)↗

DOE OSTI2016

Backbone chemical shift assignments and secondary structure analysis of the U1 protein from the Bas-Congo virus

The Bas-Congo virus (BASV) is the first rhabdovirus associated with a human outbreak of acute hemorrhagic fever. The single-stranded, negative-sense RNA genome of BASV contains the five core genes present in all rhabdoviral genomes plus an additional three genes, annotated U1, U2, and U3, with weak (<21%) sequence similarity only to a handful of genes observed in a few other rhabdoviral genomes. The function of the rhabdoviral U proteins is unknown, but, they are hypothesized to play a role in viral infection or replication. To better understand this unique family of proteins, a construct containing residues 27–203 of the 216-residue U1 protein (BASV-U1*) was prepared. By collecting data in 0.5 M urea it was possible to eliminate transient association enough to enable the assignment of most of the observable 1 H N , 1 Hα, 15 N, 13 Cα, 13 Cβ, and 13 C´ chemical shifts for BASV-U1* that will provide a foundation to study its solution properties. Here, the analyses of these chemical shifts along with 15 N-edited NOESY data enabled the identification of the elements of secondary structure present in BASV-U1*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

An updated Lagrangian LBM–DEM–FEM coupling model for dual-permeability fissured porous media with embedded discontinuities

Many engineering applications and geological processes involve embedded discontinuities in porous media across multiple length scales (e.g. rock joints, grain boundaries, deformation bands and faults). Understanding the multiscale path-dependent hydro-mechanical responses of these interfaces across length scales is of ultimate importance for applications such as CO 2 sequestration, hydraulic fracture and earthquake rupture dynamics. While there exist mathematical frameworks such as extended finite element and assumed strain to replicate the kinematics of the interfaces, modeling the cyclic hydro-mechanical constitutive responses of the interfaces remains a difficult task. Here this paper presents a semi-data-driven multiscale approach that obtains both the traction-separation law and the aperture-porosity–permeability relation from micro-mechanical simulations performed on representative elementary volumes in the finite deformation range. To speed up multiscale simulations, the incremental constitutive updates of the mechanical responses are obtained from discrete element simulations at the representative elementary volume whereas the hydraulic responses are generated from a neural network trained with data from lattice Boltzmann simulations. These responses are then linked to a macroscopic dual-permeability model. This approach allows one to bypass the need of deriving multi-physical phenomenological laws for complex loading paths. More importantly, it enables the capturing of the evolving anisotropy of the permeabilities of the macro- and micro-pores. A set of numerical experiments are used to demonstrate the robustness of the proposed model.

42 ENGINEERING↗