Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “backbone”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Resonant X-ray Diffraction Reveals the Location of Counterions in Doped Organic Mixed Ionic Conductors

Organic mixed ionic–electronic conductors (OMIECs) have the potential to enable diverse new technologies, ranging from novel in situ biosensors to flexible energy storage devices and neuromorphic computing platforms. However, their complex behavior in functional films involving electrolyte-induced swelling, ion ingress, and electrochemical doping inhibits rational material design. Of critical importance is an understanding of the specific location of the ions in the volumetrically doped material, yet this information is not readily available. Here, in this report, we present the use of grazing- incidence resonant X-ray diffraction (RXRD, also known as anomalous diffraction) at S and Cl K-edges to determine the structure of a doped, prototypical, semicrystalline polymer OMIEC based on oligo(ethylene glycol) substitution of regioregular polythiophene. The RXRD measurement provides two key insights. Quantitative analysis of the RXRD allows the determination of the position of the ion relative to the polymer backbone in the crystalline regions. We find that the anion is relatively distant from the backbone, nearer to the lamella mid-plane naively in conflict with expected Coulombic attraction between the ion and the doped polymer polaron. Comparison of RXRD to Cl – fluorescence (total Cl – ) allows determination of the relative order of doping between the crystalline and amorphous regions. We find preferential doping of the crystalline regions. Both insights, the preferential doping of crystals at low potential and the specific location of the counterion with respect to the polymer backbone, are critical to developing a microscopic understanding of transport in OMIECs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Quantum-Based Approach to Predict Primary Radiation Damage in Polymeric Networks

Initial atomistic-level radiation damage in chemically reactive materials is thought to induce reaction cascades that can result in undesirable degradation of macroscale properties. Ensembles of quantum-based molecular dynamics (QMD) simulations can accurately predict these cascades, but extracting chemical insights from the many underlying trajectories is a labor-intensive process that can require substantial a priori intuition. We develop here a general and automated graph-based approach to extract all chemically distinct structures sampled in QMD simulations and apply our approach to predict primary radiation damage of polydimethylsiloxane (PDMS), the main constituent of silicones. A postprocessing protocol is developed to identify underlying polymer backbone structures as connected components in QMD trajectories. Furthermore, these backbones form a repository of radiation-damaged structures. A scheme for extracting and updating a library of isomorphically distinct structures is proposed to identify the spanning set and aid chemical interpretation of the repository. The analyses are applied to ensembles of cascade QMD simulations in which the four element types in PDMS are selectively excited in primary knock-on atom events. Our approach reveals a much higher degree of combinatorial complexity in this system than was inferred through radiolysis experiments. Probabilities are extracted for radiation-induced network changes including formation of branch points, carbon linkages, cycles, bond scissions, and carbon uptake into the Si–O siloxane backbone network. The general analysis framework presented here is readily extendable to modeling chemical degradation of other polymers and molecular materials and provides a basis for future quantum-informed multiscale modeling of radiation damage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Folding Coarse-Grained Oligomer Models with PyRosetta

Non-biological foldamers are a promising class of macromolecules that share similarities to classical biopolymers such as proteins and nucleic acids. Currently, designing novel foldamers is a non-trivial process, often involving many iterations of trial synthesis and characterization until folded structures are observed. In this work, we aim to tackle these foldamer design challenges using computational modeling techniques. We developed CG PyRosetta, an extension to the popular protein folding python package, PyRosetta, which introduces coarse-grained (CG) residues into PyRosetta, enabling the folding of toy CG foldamer models. Although these models are simplified, they can help explore overarching physical hypotheses about how oligomers can form. Through systematic variation of CG parameters in these models, we can investigate various folding hypotheses at the CG scale to inform the design process of new foldamer chemistries. In this study, we demonstrate CG PyRosetta’s ability to identify minimum energy structures with a diverse structural search over a range of simple models, as well as two hypothesis-driven parameter scans investigating the effects of side-chain size and internal backbone angle on secondary structures. We are able to identify several types of secondary structures from single- and double-helices to sheet-like and knot-like structures. Here, we show how side-chain size and backbone bond angle both play an important role in the structure and energetics of these toy models. Optimal side-chain sizes promote favorable packing of side chains, while specific backbone bond angles influence the specific helix type found in folded structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Percolated Ionic Aggregate Morphologies and Decoupled Ion Transport in Precise Sulfonated Polymers Synthesized by Ring-Opening Metathesis Polymerization

Here, we describe a set of precise single-ion conducting polymers that form self-assembled percolated ionic aggregates in glassy polymer matrices and have decoupled transport of metal cations. These precise single-ion conductors (SICs), synthesized by a scalable ring-opening metathesis polymerization, consist of a polyethylene backbone with a sulfonated phenyl group pendant on every fifth carbon and are fully neutralized by a counterion X + (Li + , Na + , or Cs + ). Experimental X-ray scattering measurements and fully atomistic molecular dynamics (MD) simulations are in good agreement. The MD simulations show that the ionic groups nanophase separate from the polymer backbone to form percolating ionic aggregates. Using graph theory, we find that within the Li + - and Na + -neutralized polymers the percolated aggregates exhibit planar and ribbon-like configurations at intermediate length scales, while the percolated aggregates within the Cs + -neutralized polymers are more isotropic. Electrical impedance spectroscopy measurements show that the ionic conductivities exhibit Arrhenius behavior, with conductivities of 10 –7 to 10 –6 S/cm at 180 °C. In the MD simulations, the cations move between sulfonate groups in the percolated aggregates, larger ions travel further, and overall cations travel further than the polymer backbones, indicating a decoupled ion-transport mechanism. Thus, the percolated ionic aggregates in these polymers can serve as pathways to facilitate decoupled ion motion through a glassy polymer matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Work of Mechanical Degradation in Elongating Polymer Melts

Molecular dynamics simulations are used to study the mechanical degradation of well-entangled polymer melts during uniaxial extensional flow. Simulations measure the transient rise in extensional stresses and relate them to the molecular alignment and scission of chain backbones. Intermolecular entanglements couple chain scission in space and time, making degradation sensitive to deformation history and strain rate in ways not displayed by dilute polymer solutions. The rate of chain scission is nonmonotonic and peaks at strains corresponding to the maximum extensibility of entanglement segments but prior to the full extension of chain backbones. We measure a specific work per scission event w* and decompose it into separate contributions associated with chain alignment, chemical bond breaking, and scission-induced plasticity. We find chain scission in melts requires activating plastic dissipation that is multiple orders of magnitude larger than the chemical work required to break a covalent backbone bond. Our findings underscore the critical need to consider bulk polymer mechanics and rheology in designing efficient mechanical degradation and mechanochemical processes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Precise Linker Length and Dynamic Bond Exchange Control Penetrant Diffusion in Dense Vitrimers

Polymer networks with dynamic covalent bonds have been investigated for their self-healing ability, recyclability, and potential as more sustainable materials. Recent results have indicated that in some cases, bond exchange can enhance the transport of penetrants in dense networks, pointing to their potential for separations of membranes. Here, imine dynamic bonds in ethylene oxide (EO) networks with precise linker lengths were synthesized to investigate the transport of N,N′-bis(2,5-di-tert-butylphenyl)-3,4,9,10-perylenedicarboximide (BTBP), a large, anisotropic dye molecule. Networks with mesh sizes smaller than, comparable to, and greater than the size of the penetrant axes were investigated to probe the effects of bond exchange and network confinement on transport. Mesh sizes, which ranged from 0.5 to 1.62 nm, were determined from shear rheology, glass transitions by calorimetry, and probe diffusion coefficients by fluorescence recovery after photobleaching. Permanent networks with identical EO chain lengths were prepared as control samples, and up to a 3 orders of magnitude increase in diffusion coefficient is observed in the dynamic systems for short linkers containing 13 backbone atoms. The longest linkers with 71 backbone atoms show no difference between the permanent and dynamic networks. Linkers shorter than 11 backbone atoms, corresponding to a mesh size smaller than the penetrant small axis, diffusion is no longer observable on the experimental time scale, indicating a sharp cutoff attributed to the precise linkers and narrow mesh size distribution. The dynamic imine exchange time scales were compared to the diffusive hopping times of penetrants and indicate that exchange can occur during a diffusive displacement. Furthermore, these findings provide insights into the factors affecting penetrant transport in dense polymers and inspire the development of next-generation selective polymer membranes.

Diffusion↗
DOE OSTI2025

NEXAFS Spectroscopy of P3HT and PBTTT at the Sulfur K-Edge

The sulfur K-edge near-edge X-ray absorption fine-structure (NEXAFS) spectra of the common conjugated polymers P3HT and PBTTT are studied from both experimental and theoretical perspectives. Experimental angle-resolved spectra are measured to characterize both the dominant peaks and the dichroism of the polymers. First-principles calculations using the density functional theory-based many-body X-ray absorption spectroscopy (MBXAS) method are performed for the two polymers as well as for the thiophene and thienothiophene units that make up the conjugated backbones of these polymers. Through this combined approach, we are able to confidently assign the observed peaks to specific molecular orbitals and identify the orientation of their transition dipole moments (TDMs) with respect to the coordinate frame of the polymer backbone. Here, in particular, we are able to establish the character and orthogonal nature of the three main low-energy peaks at: (i) 2473.5 eV, 1s → (S–C)­π* with TDM along the π-stacking direction; (ii) 2474.1 eV, 1s → (S–C)­σ* with TDM along the backbone; and (iii) 2475.4 eV, 1s → (S–C)­σ* with TDM perpendicular to the first two. By performing both gas-phase and solid-state simulations, and with reference to the NEXAFS spectra of thiophene and thienothiophene building blocks, the influences of polymerization and molecular packing are also explored.

Chantler, Paul Alexander [Monash University, Clayt↗
DOE OSTI2025

Single-Molecule Conductance through Hybrid Radially and Linearly π-Conjugated Macromolecules Reveals an Unusual Intramolecular π-Interaction

We describe the design, synthesis, and single-molecule junction conductance of π-electron molecules bearing both radial and linear π-conjugation pathways, whereby cycloparaphenylene (CPP) radial cores are π-extended linearly with aryl alkyne substituents as models for previously reported CPP-arylene ethynylene conjugated polymers. Although radially and linearly conjugated molecules have been studied previously in isolation as junction-bridging molecular electronic units, this is the first study to examine molecules where both topologies are operative. Our results reveal that the presence of radial CPP components within the junction-spanning pathway leads to a reduction in the conductance of the backbone compared to model linear phenyl substituents. Through tight-binding and DFT-based calculations, we attribute this conductance change to intramolecular van der Waals (vdW) interactions between the CPP ring and the junction-spanning arylene-ethynylene molecular backbone. These interactions induce changes in the dihedral angles of the backbone, leading to a reduced overlap of π orbitals within the molecular junction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A Unique Case of the “Goldilocks Rule” in Solid-State Electrolytes: Two Are Good, Four Are Too Many

We report the syntheses of two new series of methacrylate monomers with different backbones: ureidopyrimidinone (PU) and boron-substituted urea pyrimidine (U), which enhance both the mechanical and electrochemical properties of the solid-state electrolyte (SSE) while improving the cycle life of lithium iron phosphate (LiFePO 4 , LFP) cathodes. The PU backbone is characterized by four hydrogen bonds (H-bonds), while the U backbone bears only two. Importantly, our research reveals that two H-bonds in these monomers are optimal; in contrast, four are excessive. The exceptional mechanical properties and processability of the SSE with the U series additives, resulting from the optimal H-bonds, were unexpectedly achieved. This leads to the establishment of a “Goldilocks rule” for additive design. The key strategies include: 1) reducing hydrogen-bonding (H-bonding) sites by changing pyrimidinone to pyrimidine and 2) shifting from intermolecular to intramolecular H-bonding and π−π bonding. Furthermore, this reduction in H-bonding also offers significant advantages in processability. The advancement can be extended to electrode fabrication, making the manufacturing of all-solid-state batteries more practical and efficient.

hydrogen bonding↗
DOE OSTI2020

Sulfurized Polyacrylonitrile for High-Performance Lithium–Sulfur Batteries: In-Depth Computational Approach Revealing Multiple Sulfur’s Reduction Pathways and Hidden Li + Storage Mechanisms for Extra Discharge Capacity

Like no other sulfur host material, polyacrylonitrile-derived sulfurized carbon (SPAN) promises improved electrochemical performance for lithium–sulfur batteries, based on its compatibility with carbonate solvents and ability to prevent self-discharge and shuttle effect. However, a complete understanding of the SPAN’s lithiation mechanism is still missing because its structural features vary widely with synthesis conditions, and its electrochemical performance deviates from elemental sulfur. This study continues our research on the elucidation of the SPAN’s structural characteristics and lithiation mechanisms via computational approaches. Here, our models reproduce most experimental data regarding carbon’s graphitization level and conjugated ordering, sulfur–carbon covalent bonding, sulfur loading, and elemental composition. Our simulations emulate the discharge voltage observed in experiments for the first discharge, which reveals that sulfur follows multiple reduction pathways based on its interaction with the carbon backbone. Sulfur reduction takes place above 1.0 V vs Li/Li + mostly in the SPAN-like material, with no long-chain lithium polysulfide formation. Below 1.0 V vs Li/Li + , the backbone’s electrochemical activity occurs via multiple C–Li and N–Li interactions, mostly with edge carbon atoms and pyridinic nitrogen. Moreover, we identify Li + binding sites throughout the graphitized backbone that might lead to prohibited energy costs for Li + deintercalation, which may explain the irreversible capacity loss between the first and second discharges. This work improves understanding of lithiation mechanisms in sulfurized carbon, which is useful for rationally designing SPAN synthesis pathways tailored to increase sulfur loading and enhanced electrochemical performance.

25 ENERGY STORAGE↗
DOE OSTI2021

A mathematical analysis of carbon fixing materials that grow, reinforce, and self-heal from atmospheric carbon dioxide

Carbon fixing materials are a new class of self-healing, self-reinforcing materials we have introduced that utilize ambient CO 2 to chemically add to an ever extending carbon backbone. This class of materials can utilize biological or non-biological photocatalysts and support a wide range of potential backbone chemistries. However, there is no analysis to date that describes their fundamental limits in terms of chemical kinetics and mass transfer. In this computational study, we employ a reaction engineering, and materials science analysis to answer basic questions about the maximum growth rate, photocatalytic requirements and limits of applicable materials. Our proposed mathematical framework envelops three main functions required for carbon fixing materials: (1) adsorption of CO 2 from air, (2) photocatalytic reduction of CO 2 into selected monomers, and (3) polymerization of CO 2 -derived monomers. First, by performing a Damköhler number analysis, we derive criteria for the cross over from kinetic control to mass transfer limited growth, setting upper limits on performance of a potential photocatalyst. Next, we analyze photocatalytic reduction of CO 2 to single carbon products, using known catalytic pathways and kinetic data. We identify formaldehyde as a C1 intermediate having unique potential for incorporation into the material backbone of carbon fixing materials. As an example, we find that a diamond-shaped reaction network graph for CO 2 reduction, passing through CO and HCOOH intermediates, accurately describes kinetic data for cobalt-promoted TiO 2 nanoparticles at room temperature and 1 atm CO 2 . Finally, as an applied case study, we consider and analyze the example of a carbon fixing poly(oxymethylene) system with embedded catalyst promoting the photocatalytic reduction of CO 2 to formaldehyde. The latter then trimerizes to a trioxane monomer which subsequently polymerizes to polyoxymethylene. This reaction engineering analysis introduces benchmarks for carbon fixing materials with respect to achievable rates of photocatalysis, adsorption, and polymerization. These results should prove valuable for the design, evaluation, and benchmarking of this emergent and new class of environmentally sustainable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Large magnetic anisotropy of a decorated spin-chain system K 2 Co 3 (MoO 4 ) 3 (OH) 2

The magnetic structure of K 2 Co 3 (MoO 4 ) 3 (OH) 2 is studied in detail. The material has a half-sawtooth one-dimensional (1-D) structure containing two unique Co 2+ ions, one in the chain backbone and one on the apex of the sawtooth creating a series of isosceles triangles along the b-axis. These triangles can be a source of magnetic frustration. The ability to grow large single crystals enables detailed magnetic measurements with the crystals oriented in a magnetic field along the respective axes. It has a Curie–Weiss temperature θ CW of 5.3(2) K with an effective magnetic moment of 4.8(3)μ B /Co. The material is highly anisotropic with a sharp antiferromagnetic ordering transition at 7 K with a metamagnetic transition at 2 kOe. Neutron diffraction was used to determine the magnetic structure and revealed a magnetic structure with canted spins along the backbone of the chain while spins along the sawtooth caps maintained a colinear orientation, arranging antiferromagnetically relative to the backbone spins. In conclusion, the parallel chains arrange antiferromagnetically relative to each other along the c-axis and ferromagnetically along the a-axis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Manipulation of intramolecular hydrogen bonds in conjugated pseudoladder polymer for semiconductivity and solution-processability

The conformational coplanarity and local rigidity of π-conjugated backbones are critical for the semiconducting performance of organic electronic materials. The conformational coplanarity and local rigidity of π-conjugated backbones are critical for the semiconducting performance of organic electronic materials. While fusing the aromatic system into a ladder-type structure effectively enhances these properties, it also often results in poor solution processability and hence limits their transition to device application. To address this challenge, an intramolecularly hydrogen-bonded pseudoladder polymer (HPLP) system based on alternating hydrogen bond donating benzobisimidazole (BBI) and hydrogen bond accepting benzodifuran (BDF) units, is designed and synthesized. A Boc-protected precursor of HPLP allows for feasible solution processing of the polymer into thin films. Subsequently, in situ thermal Boc-deprotection generates the HPLP polymer, in which intramolecular hydrogen bonds form between each pair of neighboring repeating units, inducing coplanarity and rigidity throughout the entire backbone. This process is accompanied by a significant red-shift of the absorption spectrum, reduced bandgap, and enhanced rigidity, as confirmed by NMR, UV-Vis, and density functional theory analyses. HPLP films exhibit a three-order-of-magnitude enhancement in charge carrier mobility compared to the Boc-protected precursor and demonstrate excellent solvent resistance in organic thin-film transistors.

Miranda, Octavio↗
DOE OSTI2024

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Science & Technology - Other Topics↗
DOE OSTI2021

Dependence of thick filament structure in relaxed mammalian skeletal muscle on temperature and interfilament spacing

Contraction of skeletal muscle is regulated by structural changes in both actin-containing thin filaments and myosin-containing thick filaments, but myosin-based regulation is unlikely to be preserved after thick filament isolation, and its structural basis remains poorly characterized. Here, we describe the periodic features of the thick filament structure in situ by high-resolution small-angle x-ray diffraction and interference. We used both relaxed demembranated fibers and resting intact muscle preparations to assess whether thick filament regulation is preserved in demembranated fibers, which have been widely used for previous studies. We show that the thick filaments in both preparations exhibit two closely spaced axial periodicities, 43.1 nm and 45.5 nm, at near-physiological temperature. The shorter periodicity matches that of the myosin helix, and x-ray interference between the two arrays of myosin in the bipolar filament shows that all zones of the filament follow this periodicity. The 45.5-nm repeat has no helical component and originates from myosin layers closer to the filament midpoint associated with the titin super-repeat in that region. Cooling relaxed or resting muscle, which partially mimics the effects of calcium activation on thick filament structure, disrupts the helical order of the myosin motors, and they move out from the filament backbone. Compression of the filament lattice of demembranated fibers by 5% Dextran, which restores interfilament spacing to that in intact muscle, stabilizes the higher-temperature structure. The axial periodicity of the filament backbone increases on cooling, but in lattice-compressed fibers the periodicity of the myosin heads does not follow the extension of the backbone. Thick filament structure in lattice-compressed demembranated fibers at near-physiological temperature is similar to that in intact resting muscle, suggesting that the native structure of the thick filament is largely preserved after demembranation in these conditions, although not in the conditions used for most previous studies with this preparation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Self-assembled complex micelle phase stability in ABA-type triblock copolymers and related core-shell bottlebrushes

We report the synthesis and temperature-dependent morphologies of a series of polylactide--poly (-decalactone)--polylactide (LDL) triblock copolymers with M n = 16.0 – 18.1 kg / mol and volume fractions f L = 0.27 – 0.31 and associated core-shell bottlebrush (csBB) polymers, which derive from enchaining LDL triblocks through a polymerizable midchain functionality. While the LDL triblocks form micellar Frank-Kasper A15 and phases due to the conformational asymmetry of this monomer pair, the csBB morphologies sensitively depend on the backbone degree of polymerization ( N bb ). At low N bb values, micellar Frank-Kasper phases with the brush backbone situated in the matrix domain are stable, albeit with a modest reduction in the mean interfacial curvature evidenced by a to A15 order-to-order transition. However, larger N bb values drive csBBs to form hexagonally packed cylinders phases. This N bb -dependent phase behavior is rationalized in terms of a star-to-bottlebrush transition. At low N bb values, the csBBs are akin to star polymers with pointlike junctions that can support complex micelle packings. As N bb increases, the csBBs adopt cylindrical molecular geometries with extended backbones situated in the matrix domain that prefer hexagonally packed cylinders morphologies. Published by the American Physical Society 2024

Materials Science↗
DOE OSTI2020

A Method to Represent a Well in a Three‐Dimensional Discrete Fracture Network Model

Abstract In discrete fracture network (DFN) modeling, fractures are randomly generated and placed in the model domain. The rock matrix is considered impermeable. Small fractures and isolated fractures are often ignored to reduce computational expense. As a result, the rock matrix between fractures could be large and intersections may not be found between a well introduced in the model and the hydraulically connected fracture networks (fracture backbones). To overcome this issue, this study developed a method to conceptualize a well in a three‐dimensional (3D) DFN using two orthogonal rectangular fractures oriented along the well's axis. Six parameters were introduced to parameterize the well screen and skin zone, and to control the connectivity between the well and the fracture backbones. The two orthogonal fractures were discretized using a high‐resolution mesh to improve the quality of flow and transport simulations around and along the well. The method was successfully implemented within dfnWorks 2.0 (Hyman et al. 2015) to incorporate a well in a 3D DFN and to track particles leaving an injection well and migrating to a pumping well. Verification of the method against MODFLOW/MODPATH found a perfect match in simulated hydraulic head and particle tracking. Using three examples, the study showed that the method ensured the connectivity between wells and fracture backbones, and honored the physical processes of flow and transport along and around wells in DFNs. Recommendations are given for estimating the values of the six introduced well parameters in a real‐world case study.

Pham, Hai (ORCID:0000000310164394)↗