Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The design, implementation, and performance of the LZ calibration systems

LUX-ZEPLIN (LZ) is a tonne-scale experiment searching for direct dark matter interactions and other rare events. It is located at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. The core of the LZ detector is a dual-phase xenon time projection chamber (TPC), designed with the primary goal of detecting Weakly Interacting Massive Particles (WIMPs) via their induced low energy nuclear recoils. Surrounding the TPC, two veto detectors immersed in an ultra-pure water tank enable reducing background events to enhance the discovery potential. Intricate calibration systems are purposely designed to precisely understand the responses of these three detector volumes to various types of particle interactions and to demonstrate LZ's ability to discriminate between signals and backgrounds. In this paper, we present a comprehensive discussion of the key features, requirements, and performance of the LZ calibration systems, which play a crucial role in enabling LZ's WIMP-search and its broad science program. The thorough description of these calibration systems, with an emphasis on their novel aspects, is valuable for future calibration efforts in direct dark matter and other rare-event search experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Skyrmion superconductivity: DMRG evidence for a topological route to superconductivity

It was recently suggested that the topology of magic-angle twisted bilayer graphene's (MATBG) flat bands could provide a novel mechanism for superconductivity distinct from both weakly coupled BCS theory and the d-wave phenomenology of the high-T c cuprates. In this work, we examine this possibility using a density matrix renormalization group (DMRG) study of a model which captures the essential features of MATBG's symmetry and topology. Here, using large-scale cylinder-DMRG calculations to obtain the ground state and its excitations as a function of the electron doping, we find clear evidence for superconductivity driven by the binding of electrons into charge-2e skyrmions. Remarkably, this binding is observed even in the regime where the unscreened Coulomb repulsion is by far the largest energy scale, demonstrating the robustness of this topological, all-electronic pairing mechanism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Impact of beam–beam effects on absolute luminosity calibrations at the CERN Large Hadron Collider

At the Large Hadron Collider (LHC), absolute luminosity calibrations obtained by the van der Meer (vdM) method are affected by the mutual electromagnetic interaction of the two beams. The colliding bunches experience relative orbit shifts, as well as optical distortions akin to the dynamic-β effect, that both depend on the transverse beam separation and must therefore be corrected for when deriving the absolute luminosity scale. In the vdM regime, the beam–beam parameter is small enough that the orbit shift can be calculated analytically. The dynamic-β corrections to the luminometer calibrations, however, had until the end of Run 2 been estimated in the linear approximation only. In this report, the influence of beam–beam effects on the vdM-based luminosity scale is quantified, together with the associated systematic uncertainties, by means of simulations that fully take into account the non-linearity of the beam–beam force, as well as the resulting non-Gaussian distortions of the transverse beam distributions. Two independent multiparticle simulations, one limited to the weak-strong approximation and one that models strong-strong effects in a self-consistent manner, are found in excellent agreement; both predict a percent-level shift of the absolute pp-luminosity values with respect to those assumed until recently in the physics publications of the LHC experiments. These results also provide guidance regarding further studies aimed at reducing the beam–beam-related systematic uncertainty on beam–beam corrections to absolute luminosity calibrations by the van der Meer method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Demonstration of Grid Services Using Mixed Grid-Forming and Grid-Following Technologies at the Wheatridge Renewable Energy Facility (Final Technical Report)

This Final Technical Report summarizes the analytical, modeling, and engineering work performed to evaluate grid-forming (GFM) inverter capabilities within a large hybrid renewable energy facility. The project focused on assessing the ability of advanced inverter-based resources to provide essential reliability services through coordinated operation of grid-forming (GFM) and grid-following (GFL) technologies. During Budget Period 1, the project developed a comprehensive framework of GFM performance metrics, including angle support, voltage regulation, frequency response, damping behavior, and current-limiting performance. Extensive electromagnetic transient (EMT) studies and hardware-in-the-loop (HIL) testing were conducted to validate GFM behavior under a range of operating conditions, including weak-grid scenarios, voltage disturbances, and multi-resource interactions. The project also established validated modeling approaches, plant-level control integration strategies, commissioning frameworks, and high-speed measurement infrastructure to support future field demonstration. Although the project concluded prior to field demonstration, the results provide a utility-scale foundation for evaluating, modeling, and deploying grid-forming technologies. The methodologies and tools developed contribute to industry understanding of inverter-based resource behavior and support future power system reliability under increasing renewable penetration.

14 SOLAR ENERGY↗
DOE OSTI2022

Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma

The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1990

Rheology of solid methane and nitrogen - Applications of Triton

The existing information about the rheological properties of solid methane and nitrogen is reviewed. The main limitation is the absence of creep data for solid nitrogen, but estimated properties are derived on the basis of simple scaling arguments for van der Waals bonded materials. On Triton solid methane and most likely nitrogen are too soft to support any topography over billions of years. Topographic support by water ice-rich materials seems almost certain. Nitrogen and methane do not intermix by solid-state diffusion. The weakness (absence?) of the opposition effect on Triton can be explained by annealing of nitrogen grains, provided the latter have sizes characteristic of freshly deposited frost.

Eluszkiewicz, Janusz↗
NASA NTRS1993

Ion mixing in the plasma sheet boundary layer by drift instabilities

The linear stability properties of collisionless drift instabilities are analyzed in a Harris equilibrium model of the plasma sheet boundary layer (PSBL). The strearmng ions with drift-type instabilities driven in the PSBL are considered. The fluid approximation leads to growth but predicts that the mode width approaches the gyroradius of the energetic ions. Thus an integral equation theory for the modes is developed taking into account that in the PSBL the curvature drift is weak compared with the grad-B drift. The exact wave particle resonance is kept in the nonlocal response functions. Plasma density, temperature, and magnetic gradient drift motions are taken into account. The drift modes produce an anomalous cross-field momentum transport mixing the PSBL ions on the time scale of tens of seconds. A nonlinear simulation is performed which shows the coalescence of the small scale, fast growing modes into large-scale vortices. The relation between these collective modes and plasma sheet transport phenomena is discussed including the comparison with the competing plasma mixing from single-particle stochasticity.

Horton, W.↗
NASA NTRS1993

The multiplicity of T Tauri stars in the star forming regions Taurus-Auriga and Ophiuchus-Scorpius: A 2.2 micron speckle imaging survey

We present the results of a magnitude limited (K less than = 8.5 mag) speckle imaging survey of 69 T Tauri stars in the star forming regions Taurus-Auriga and Ophiuchus-Scorpius. Thirty-three companion stars were found with separations ranging from 0.07 sec to 2.5 sec, nine are new detections. This survey reveals a distinction between the classical T Tauri stars (CTTS) and the weak-lined T Tauri stars (WTTS) based on the binary star frequency as a function of separation: the WTTS binary star distribution is enhanced at the closer separations (less than 50 AU) relative to the CTTS binary star distribution. We suggest that the nearby companion stars shorten the accretion time scale in multiple star systems, thereby accounting for the presence of WTTS that are coeval with many CTTS. The binary star frequency in the projected linear separation range 16 to 252 AU for T Tauri stars (60 (+/- 17)%) is a factor of 4 greater than that of the solar-type main-sequence stars (16(+/- 3)%). Given the limited separation range of this survey, the rate at which binaries are detected suggests that most, if not all, T Tauri stars have companions. We propose that the observed overabundance of companions of T Tauri stars is an evolutionary effect, in which triple and higher order T Tauri stars are disrupted by close encounters with another star or system of stars.

Ghez, A. M.↗
NASA NTRS2002

Sensitivities of Modeled Tropical Cyclones to Surface Friction and the Coriolis Parameter

In this investigation the sensitivities of a 2-D tropical cyclone (TC) model to surface frictional coefficient and the Coriolis parameter are studied and their implication is discussed. The model used is an axisymmetric version of the latest version of the Goddard cloud ensemble model. The model has stretched vertical grids with 33 levels varying from 30 m near the bottom to 1140 m near the top. The vertical domain is about 21 km. The horizontal domain covers a radius of 962 km (770 grids) with a grid size of 1.25 km. The time step is 10 seconds. An open lateral boundary condition is used. The sea surface temperature is specified at 29C. Unless specified otherwise, the Coriolis parameter is set at its value at 15 deg N. The Newtonian cooling is used with a time scale of 12 hours. The reference vertical temperature profile used in the Newtonian cooling is that of Jordan. The Newtonian cooling models not only the effect of radiative processes but also the effect of processes with scale larger than that of TC. Our experiments showed that if the Newtonian cooling is replaced by a radiation package, the simulated TC is much weaker. The initial condition has a temperature uniform in the radial direction and its vertical profile is that of Jordan. The initial winds are a weak Rankin vortex in the tangential winds superimposed on a resting atmosphere. The initial sea level pressure is set at 1015 hPa everywhere. Since there is no surface pressure perturbation, the initial condition is not in gradient balance. This initial condition is enough to lead to cyclogenesis, but the initial stage (say, the first 24 hrs) is not considered to resemble anything observed. The control experiment reaches quasi-equilibration after about 10 days with an eye wall extending from 15 to 25 km radius, reasonable comparing with the observations. The maximum surface wind of more than 70 m/s is located at about 18 km radius. The minimum sea level pressure on day 10 is about 886 hPa. Thus the overall simulation is considered successful and the model is considered adequate for our investigation.

Chao, Winston C.↗
NASA NTRS2010

Global Patterns of Lightning Properties Derived by LIS

The Lightning Imaging Sensor LIS aboard the TRMM satellite provides unmatched empirical data of the global lightning distribution (up to approx.35deg S/N) since end of 1997. Climatological flash rate densities derived from LIS are standard references, e.g. for flash rate parameterizations used in GCMs. It is known that flash characteristics are quite variable, and that various quantities (like the flash energy or the NOx production per flash) vary considerably, statistically as well as systematically on regional and seasonal scales. LIS provides information beyond flash counts, in particular radiance and flash footprint. Here we present an analysis of global patterns of various lightning properties derived from LIS, in relation to the number of flashes. These normalized flash characteristics show consistent spatial patterns of regions with "strong" versus regions with "weak" lightning. Most striking is a clear land-ocean contrast, with oceanic flashes being "stronger" than continental flashes. But also over continents, flash strength shows systematic variations. Highest continental values are found over the US, while values over South America and India are quite low. These regional variations cannot be simply parameterized as function of latitude. Information on spatial patterns of mean flash "strength", though rather qualitative up to now, is potentially a valuable input for improving empirical parameterizations based on flash counts (like precipitation or lightning NOx). Further investigation is in progress to come to a more physical and quantitative understanding of the spatial patterns of the different LIS properties. In particular, it has to be checked how far they could be related to established lightning properties (like energy or the fraction of intra-cloud to cloud-to-ground flashes) or to meteorological quantities (like CAPE).

Beirle, S.↗
NASA NTRS2012

Global CO2 Distributions over Land from the Greenhouse Gases Observing Satellite (GOSAT)

January 2009 saw the successful launch of the first space-based mission specifically designed for measuring greenhouse gases, the Japanese Greenhouse gases Observing SATellite (GOSAT). We present global land maps (Level 3 data) of column-averaged CO2 concentrations (X(sub CO2)) derived using observations from the GOSAT ACOS retrieval algorithm, for July through December 2009. The applied geostatistical mapping approach makes it possible to generate maps at high spatial and temporal resolutions that include uncertainty measures and that are derived directly from the Level 2 observations, without invoking an atmospheric transport model or estimates of CO2 uptake and emissions. As such, they are particularly well suited for comparison studies. Results show that the Level 3 maps for July to December 2009 on a lO x 1.250 grid, at six-day resolution capture much of the synoptic scale and regional variability of X(sub CO2), in addition to its overall seasonality. The uncertainty estimates, which reflect local data coverage, X(sub CO2) variability, and retrieval errors, indicate that the Southern latitudes are relatively well-constrained, while the Sahara Desert and the high Northern latitudes are weakly-constrained. A probabilistic comparison to the PCTM/GEOS-5/CASA-GFED model reveals that the most statistically significant discrepancies occur in South America in July and August, and central Asia in September to December. While still preliminary, these results illustrate the usefulness of a high spatiotemporal resolution, data-driven Level 3 data product for direct interpretation and comparison of satellite observations of highly dynamic parameters such as atmospheric CO2.

Hammerling, Dorit M.↗
NASA NTRS2015

Quantitative Measurements of Electronically Excited CH Concentration in Normal Gravity and Microgravity Coflow Laminar Diffusion Flames

With the conclusion of the SLICE campaign aboard the ISS in 2012, a large amount of data was made available for the analysis of the effect of microgravity on laminar coflow diffusion flames. Previous work focused on the study of sooty flames in microgravity as well as the ability of numerical models to predict its formation in a simplified buoyancy-free environment. The current work shifts the investigation to soot-free flames, putting an emphasis on the chemiluminescence emission from electronically excited CH (CH*). This radical species is of significant interest in combustion studies: it has been shown that the electronically excited CH spatial distribution is indicative of the flame front position and, given the relatively simple diagnostic involved with its measurement, several works have been done trying to understand the ability of electronically excited CH chemiluminescence to predict the total and local flame heat release rate. In this work, a subset of the SLICE nitrogen-diluted methane flames has been considered, and the effect of fuel and coflow velocity on electronically excited CH concentration is discussed and compared with both normal gravity results and numerical simulations. Experimentally, the spectral characterization of the DSLR color camera used to acquire the flame images allowed the signal collected by the blue channel to be considered representative of the electronically excited CH emission centered around 431 nm. Due to the axisymmetric flame structure, an Abel deconvolution of the line-of-sight chemiluminescence was used to obtain the radial intensity profile and, thanks to an absolute light intensity calibration, a quantification of the electronically excited CH concentration was possible. Results show that, in microgravity, the maximum flame electronically excited CH concentration increases with the coflow velocity, but it is weakly dependent on the fuel velocity; normal gravity flames, if not lifted, tend to follow the same trend, albeit with different peak concentrations. Comparisons with numerical simulations display reasonably good agreement between measured and computed flame lengths and radii, and it is shown that the integrated electronically excited CH emission scales proportionally to the computed total heat release rate; the two-dimensional electronically excited CH spatial distribution, however, does not appear to be a good marker for the local heat release rate.

chemiluminescence↗
NASA NTRS2016

Possible Mechanisms for Turbofan Engine Ice Crystal Icing at High Altitude

A thermodynamic model is presented to describe possible mechanisms of ice formation on unheated surfaces inside a turbofan engine compression system from fully glaciated ice crystal clouds often formed at high altitude near deep convective weather systems. It is shown from the analysis that generally there could be two distinct types of ice formation: (1) when the "surface freezing fraction" is in the range of 0 to 1, dominated by the freezing of water melt from fully or partially melted ice crystals, the ice structure is formed from accretion with strong adhesion to the surface, and (2) when the "surface melting fraction" is the range of 0 to 1, dominated by the further melting of ice crystals, the ice structure is formed from accumulation of un-melted ice crystals with relatively weak bonding to the surface. The model captures important qualitative trends of the fundamental ice-crystal icing phenomenon reported earlier (Refs. 1 and 2) from the research collaboration work by NASA and the National Research Council (NRC) of Canada. Further, preliminary analysis of test data from the 2013 full scale turbofan engine ice crystal icing test (Ref. 3) conducted in the NASA Glenn Propulsion Systems Laboratory (PSL) has also suggested that (1) both types of ice formation occurred during the test, and (2) the model has captured some important qualitative trend of turning on (or off) the ice crystal ice formation process in the tested engine low pressure compressor (LPC) targeted area under different icing conditions that ultimately would lead to (or suppress) an engine core roll back (RB) event.

Crystal Icing↗
DOE OSTI2023

Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo

Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative—particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods.

Chemistry↗
DOE OSTI2025

Magnon spectroscopy in the electron microscope

Abstract The miniaturization of transistors is approaching its limits owing to challenges in heat management and information transfer speed 1 . To overcome these obstacles, emerging technologies such as spintronics 2 are being developed, which make use of the electron’s spin as well as its charge. Local phenomena at interfaces or structural defects will greatly influence the efficiency of spin-based devices, making the ability to study spin-wave propagation at the nanoscale and atomic scale a key challenge 3,4 . The development of high-spatial-resolution tools to investigate spin waves, also called magnons, at relevant length scales is thus essential to understand how their properties are affected by local features. Here we detect bulk THz magnons at the nanoscale using scanning transmission electron microscopy (STEM). By using high-resolution electron energy-loss spectroscopy with hybrid-pixel electron detectors, we overcome the challenges posed by weak signals to map THz magnon excitations in a thin NiO nanocrystal. Advanced inelastic electron scattering simulations corroborate our findings. These results open new avenues for detecting magnons and exploring their dispersions and their modifications arising from nanoscale structural or chemical defects. This marks a milestone in magnonics and presents exciting opportunities for the development of spintronic devices.

Science & Technology - Other Topics↗
DOE OSTI2023

Uncovering dark matter density profiles in dwarf galaxies with graph neural networks

Dwarf galaxies are small, dark-matter-dominated galaxies, some of which are embedded within the Milky Way. Their lack of baryonic matter (e.g., stars and gas) makes them perfect test beds for probing the properties of dark matter—understanding the spatial dark matter distribution in these systems can be used to constrain microphysical dark matter interactions that influence the formation and evolution of structures in our Universe. Here, we introduce a new method that leverages simulation-based inference and graph-based machine learning in order to infer the dark matter density profiles of dwarf galaxies from observable kinematics of stars gravitationally bound to these systems. Our approach aims to address some of the limitations of established methods based on dynamical Jeans modeling. We show that this novel method can place stronger constraints on dark matter profiles and, consequently, has the potential to weigh in on some of the ongoing puzzles associated with the small-scale structure of dark matter halos, such as the core-cusp discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Combined Coherent Manipulation and Single-Shot Measurement of an Electron Spin in a Quantum Dot (Final Technical Report)

Semiconductor quantum dots (QDs) are promising candidates to act as single-photon sources and/or quantum bits in future optical quantum information applications. Their excellent optical properties such as high brightness, single-photon purity, and narrow linewidth have potential utility in many areas. One challenge is to control the energy levels of the QD without using a magnetic field. The AC Stark effect offers the opportunity to do this. We have demonstrated record-large AC Stark shifts of the energy states of a single charged QD. We showed that the shifts can be applied in a spin-selective manner by controlling the polarization of the laser producing the AC Stark effect. In order to characterize the effect, we developed a novel spectral filtering scheme to discriminate the high-power AC Stark laser from the QD fluorescence. We also developed a compact, low-cost, homemade polarimeter to help control the polarization of the AC Stark laser, which enabled the spin selectivity of the AC Stark effect. With the practical capabilities thus developed, we learned that the spin-selective AC Stark effect causes electron spin pumping, which in turn causes nuclear spin pumping via the hyperfine interaction. The nuclear spin pumping causes a mean field Zeeman interaction between the nuclear spin ensemble and the electron trapped in the QD, resulting in the so-called Overhauser shift. The magnitude of the Overhauser shift and the measured linewidth of the QD’s optical transitions enabled characterization of the mean nuclear spin polarization and fluctuations. The potentially rapid (ns-scale) control of the QD energy levels and nuclear spin polarization will enable measurements of electron and nuclear spin polarization in the absence of any real magnetic field, which is a completely unique capability. The ability to rapidly apply a spin-selective AC Stark effect in the presence of a weak real magnetic field will allow control of the polarization selection rules of the transitions, enabling all single-qubit operations on the electron spin, i.e., initialization, coherent manipulation, and quantum non-demolition measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Tensile Modeling PVC Gels for Electrohydraulic Actuators

Polyvinyl chloride (PVC)-dibutyl adipate (DBA) gels are a fascinating dielectric elastomer actuator showing promise in soft robotics. When actuated with high voltages, the gel deforms towards the anode. A recent application of PVC gels in electrohydraulic actuators motivates elastic and hyperelastic constitutive relationships for tensile loading modes. PVC gels with plasticizer-to-polymer weight ratios of 2:1, 4:1, 6:1, and 8:1 w/w were evaluated. PVC gels exhibit a linear elastic region up to 25% strain. The elastic modulus decreased with increasing plasticizer content from 288.8 kPa, 56.1 kPa, 24.7 kPa, to 11 kPa. Poisson’s ratio also decreased with increasing plasticizer content from 0.42, 0.43, 0.39, to 0.35. We suggest that the decrease in polymer concentration facilitates a weakly interconnected polymer network susceptible to chain slippage that hinders the network response, thus lowering Poisson’s ratio. Our work suggests that PVC gels can be treated as isotropic and incompressible for large strains and hyperelastic modeling; however, highly plasticized gels tend to act less incompressible at small strains. The power scaling law between the elastic modulus and plasticizer weight ratio showed high agreement, making the elastic modulus deterministic for any plasticizer content. The Neo–Hookean, Mooney–Rivlin, Yeoh, Gent, Ogden, and extended tube hyperelastic constitutive models are investigated. The Yeoh model shows the highest feasibility when evaluated up to 3.5 stretch, showing a maximum normalized root-mean-square-error of 6.85%. Together, these findings establish a constitutive basis for PVC-DBA gels, incorporating small strain elasticity, large strain non-linear behavior, and network analysis while providing suggestive insight into the network structure required for accurately modeling the EPIC.

Polymer Science↗