Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “edge effects”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Effects of Revegetation on Soil Physical and Chemical Properties in Solar Photovoltaic Infrastructure

Solar photovoltaic (PV) technology is being deployed at an unprecedented rate. However, utility-scale solar energy development is land intensive and its large-scale installation can have negative impacts on the environment. In particular, solar energy infrastructure can require extensive landscape modification that transforms soil ecological functions, thereby impacting hydrologic, vegetative, and carbon dynamics. However, reintroducing native vegetation to solar PV sites may be a means of restoring their soils. To this end, we investigated critical soil physical and chemical parameters at a revegetated photovoltaic array and an adjacent reference grassland in Colorado, United States. Seven years after revegetation, we found that carbon and nitrogen remained lower in the PV soil than in the reference soil and contained a greater fraction of coarse particles. We also found that the PV modules introduced heterogeneity in the soil moisture distribution, with precipitation accumulating along the lower edges of panels. The redistribution of soil moisture by panel arrays could potentially be used in concert with planting strategies to maximize plant growth or minimize soil erosion, and should be considered when evaluating the potential to co-locate vegetation with solar infrastructure.

14 SOLAR ENERGY↗
DOE OSTI2024

From anti-Arrhenius to Arrhenius behavior in a dislocation-obstacle bypass: Atomistic simulations and theoretical investigation

Dislocations are the primary carriers of plasticity in metallic materials. Understanding the basic mechanisms for dislocation movement is paramount to predicting the material mechanical response. Relying on atomistic simulations, we observe a transition from non-Arrhenius to Arrhenius behavior in the rate of an edge dislocation overcoming the long-range elastic interaction with a prismatic loop in tungsten. Close to the critical resolved shear stress, the process shows a non-Arrhenius behavior at low temperatures. However, as the temperature increases, the activation entropy starts to dominate, leading to a traditional Arrhenius-like behavior. Here, we have computed the activation entropy analytically along the minimum energy path following Schoeck’s method, which captures the cross-over between anti-Arrhenius and Arrhenius domains. Also, the Projected Average Force Integrator (PAFI), another simulation method to compute free energies along an initial transition path, exhibits considerable concurrence with Schoeck’s formalism. We conclude that entropic effects need to be considered to understand processes involving dislocations bypassing elastic barriers close to the critical resolved shear stress. More work needs to be performed to fully understand the discrepancies between Schoeck’s and PAFI compared to molecular dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

X-ray and gamma-ray tomographic imaging of fuel relocation inside sodium fast reactor test assemblies during severe accidents

The present work reports on x-ray and gamma-ray high-spatial resolution computerized tomography measurements of the Pin Bundle Metallic Fuel Relocation (PBR) assemblies tested in the Metallic Uranium Safety Experiment (MUSE) facility at Argonne National Lab (ANL). The aim of the study was to characterize fuel relocation structures that develop during severe core accidents pertaining to SFR assemblies; these include but are not limited to advanced core disruption recreated in the PBR-1 assembly, and cladding breach recreated in the PBR-2 assembly. We report the x-ray tomography measurements were able to resolve small quantities of relocation fuel; with increased presence of relocation fuel, the x-ray 1measurements spatially mapped the material but could not resolve the inner regions of these. The gamma-tomography measurements showed improved results, resolving the relocation structures in great detail. The upper plenum of the PBR-1 assembly where the molten uranium was initially inserted presented high structural damage, reflected by the partial and complete disintegration of the central rods. Relocation fuel filled the subchannels, adhering to surviving cladding walls and the assembly casing. In the lower portion of the measured section, the tomogram degrades due to photon starvation effects hinting at the increased amount of relocation fuel potentially plugging the assembly; flow blockage in this section was difficult to determine due to the tomogram’s degradation from photon starvation. Small fragments were observed further down the assembly, dislodged from the initial insertion of the molten material. This section was used as an unperturbed assembly reference, with a calculated blockage of less than 1% from the present fragments. The PBR-2 assembly was characterized by columnar relocation structures propagating through the subchannels. Three of the relocation structures were captured in the measured section, with evidence of cross migration on to adjacent subchannels. The measured section in this assembly captures the leading edge of two structures. The calculated flow blockage was 16% in the planes where the three relocation structures are present, but this quickly decreases to approximately 5% past the leading edge of two of the structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Node-degree aware edge sampling mitigates inflated classification performance in biomedical random walk-based graph representation learning

Motivation: Graph representation learning is a family of related approaches that learn low-dimensional vector representations of nodes and other graph elements called embeddings. Embeddings approximate characteristics of the graph and can be used for a variety of machine-learning tasks such as novel edge prediction. For many biomedical applications, partial knowledge exists about positive edges that represent relationships between pairs of entities, but little to no knowledge is available about negative edges that represent the explicit lack of a relationship between two nodes. For this reason, classification procedures are forced to assume that the vast majority of unlabeled edges are negative. Existing approaches to sampling negative edges for training and evaluating classifiers do so by uniformly sampling pairs of nodes. Results: We show here that this sampling strategy typically leads to sets of positive and negative examples with imbalanced node degree distributions. Using representative heterogeneous biomedical knowledge graph and random walk-based graph machine learning, we show that this strategy substantially impacts classification performance. If users of graph machine-learning models apply the models to prioritize examples that are drawn from approximately the same distribution as the positive examples are, then performance of models as estimated in the validation phase may be artificially inflated. We present a degree-aware node sampling approach that mitigates this effect and is simple to implement. Availability and implementation: Our code and data are publicly available at https://github.com/monarch-initiative/negativeExampleSelection.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Hardware-in-the-Loop Evaluation of Grid-Edge DER Chip Integration Into Next-Generation Smart Meters

To facilitate the implementation of distributed energy resource management systems (DERMS), we propose to insert a grid-edge distributed energy resource (DER) chip hosting a DERMS algorithm into the next generation of smart meters. This will create a pathway for the wide adoption of DERMS technology because many utilities plan to invest in advanced metering infrastructure in the near future. This will also bridge the gap between an electrical power utility and DERs behind the meter. The DER chip is designed to follow power direction signals from the DERMS coordinator while balancing its local objectives. We tested the chip using a controller- and power-hardware-in-the-loop evaluation under three scenarios that a DERMS could face in the real world. The DER chip was capable of and effective at directing four heterogeneous DERs to respond to a DERMS coordinator for grid services (e.g., voltage regulation and a virtual power plant).

distributed energy resource management system↗
DOE OSTI2023

Hardware-in-the-Loop Evaluation of Grid-Edge DER Chip Integration Into Next-Generation Smart Meters: Preprint

To facilitate the implementation of distributed energy resource management systems (DERMS), we propose to insert a grid-edge distributed energy resource (DER) chip hosting a DERMS algorithm into the next generation of smart meters. This will create a pathway for the wide adoption of DERMS technology because many utilities plan to invest in advanced metering infrastructure in the near future. This will also bridge the gap between an electrical power utility and DERs behind the meter. The DER chip is designed to follow power direction signals from the DERMS coordinator while balancing its local objectives. We tested the chip using a controller- and power-hardware-in-the-loop evaluation under three scenarios that a DERMS could face in the real world. The DER chip was capable of and effective at directing four heterogeneous DERs to respond to a DERMS coordinator for grid services (e.g., voltage regulation and a virtual power plant).

distributed energy resource management system↗
DOE OSTI2025

Tungsten-dioxo single-site heterogeneous catalyst on carbon: synthesis, structure, and catalysis

This study investigates the application of a novel third-row metal, tungsten, to carbon-supported single-site metal-oxo heterogeneous catalysis. Tungsten is a green and earth-abundant metal, but an unexplored candidate in this role. The carbon (AC = activated carbon)-supported tungsten dioxo complex, AC/WO 2 was prepared via grafting of (DME)WO 2 Cl 2 (DME = 1,2-dimethoxyethane) onto high-surface-area activated carbon. AC/WO 2 was fully characterized by ICP-OES, XPS, EXAFS, XANES, SMART-EM, and DFT. W 4d 7/2 XPS and W L III -Edge XANES assign the oxidation state as W(VI), while EXAFS reveals two W=O double and two W–O single bonds at distances of 1.73 and 1.92 Å, respectively. These data align well with DFT computational results, supporting the structure as Carbon(–μ-O–) 2 M(=O) 2 . SMART-EM verifies that single W(VI) catalytic sites are bonded in an out-of-plane manner. The catalytic performance of air- and water-stable AC/WO 2 is compared to that of AC/MoO 2 . AC/WO 2 is more active and selective than the molybdenum analog in mediating alcohol dehydration of various substrates, and is recyclable. Notably, AC/WO 2 is an effective and recyclable catalyst for primary aliphatic alcohol dehydration and forms no dehydrogenation side products in contrast to AC/MoO 2 . However, AC/WO 2 is less effective in epoxidation and PET depolymerization. Overall, this work demonstrates the potential of carbon-supported third row metals for future studies.

02 PETROLEUM↗
DOE OSTI2020

Short-term load demand forecasting through rich features based on recurrent neural networks

With the emerging penetration of renewables and dynamic loads, the understanding of grid edge loading conditions becomes increasingly substantial. Load modelling researches commonly consist of explicitly expressed load models and non-explicitly expressed techniques, of which artificial intelligence approaches turn out to be the major path. This paper reveals the artificial intelligence-based load modelling technique to enhance the knowledge of current and future load information considering geographical and weather dependencies. This paper presents a recurrent neural network based sequence to sequence (Seq2Seq) model to forecast the short-term power loads. Also, a feature attention mechanism, which is along channel and time directions, is developed to improve the efficiency of feature learning. The experiments over three publicly available datasets demonstrate the accuracy and effectiveness of the proposed model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Alpha particle channeling in ITER

Alpha particle channeling, the attempt to extract energy from high-energy fusion alpha particles using a combination of high- and low-frequency waves introduced in the plasma with antennae, is studied using a guiding center code. We examine the effect of the combination of cooling caused by a cyclotron mode, particle diffusion due to Alfvén modes introduced for ash removal, the inclusion of microturbulent diffusion, and the slowing down of the alpha particles on electrons. Induced cooling is achieved with a cyclotron resonance, producing additional losses of particles above 3.5 MeV in an amount about equal to the prompt loss. The removal of the low-energy cooled particles by additional modes or by microturbulent diffusion is the hardest challenge for alpha channeling. The cooled alpha particles remain deep in the plasma center, so for this removal, the method used must be capable of providing transport from the plasma center to the edge without also causing losses of high-energy particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Surface facet dependence of competing alloying mechanisms

Metal alloys are ubiquitous in many branches of heterogeneous catalysis, and it is now fairly well established that the local atomic structure of an alloy can have a profound influence on its chemical reactivity. While these effects can be difficult to probe in nanoparticle catalysts, model studies using well defined single crystal surfaces alloyed with dopants enable these structure–function correlations to be drawn. The first step in this approach involves understanding the alloying mechanism and the type of ensembles formed. In this study, we examined the atomic structure of RhCu single-atom alloys formed on Cu(111), Cu(100), and Cu(110) surfaces. Our results show a striking difference between Rh atoms alloying in Cu(111) vs the more open Cu(100) and Cu(110) surface facets. Unlike Cu(111) on which Rh atoms preferentially place-exchange with Cu atoms in the local regions above step edges leaving the majority of the Cu surface free of Rh, highly dispersed, homogeneous alloys are formed on the Cu(100) and (110) surfaces. These dramatically different alloying mechanisms are understood by quantifying the energetic barriers for atomic hopping, exchange, swapping, and vacancy filling events for Rh atoms on different Cu surfaces through theoretical calculations. Density functional theory results indicate that the observed differences in the alloying mechanism can be attributed to a faster hopping rate, relatively high atomic exchange barriers, and stronger binding of Rh atoms in the vicinity of step edges on Cu(111) compared to Cu(110) and Cu(100). Furthermore, these model systems will serve as useful platforms for examining structure sensitive chemistry on single-atom alloys.

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

Unraveling the Orbital Physics in a Canonical Orbital System KCuF 3

We explore the existence of the collective orbital excitations, orbitons, in the canonical orbital system KCuF 3 using the Cu L 3 -edge resonant inelastic x-ray scattering. We show that the nondispersive high-energy peaks result from the Cu 2+ dd orbital excitations. These high-energy modes display good agreement with the ab initio quantum chemistry calculation, indicating that the d d excitations are highly localized. At the same time, the low-energy excitations present clear dispersion. They match extremely well with the two-spinon continuum following the comparison with Müller ansatz calculations. The localized d d excitations and the observation of the strongly dispersive magnetic excitations suggest that the orbiton dispersion is below the resolution detection limit. Our results can reconcile with the strong local Jahn-Teller effect in KCuF 3 , which predominantly drives orbital ordering.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The role of B T -dependent flows on W accumulation at the edge of the confined plasma

Abstract Near-separatrix impurity accumulation between the crown and the outer midplane of tokamaks is a common feature in results from codes such as SOLPS-ITER and DIVIMP; however, experimental evidence of accumulation has only recently been obtained and is reported here. The codes find that the poloidal distribution of impurity ions in the scrape-off layer (SOL) depends primarily on toroidal field ( B T )-dependent parallel flow patterns of the background plasma and the parallel ion temperature gradient (∇ ‖ T ion ) force. Experimentally, Mach probes used in L-mode plasmas with favorable (for H-mode access) B T measure fast ( M ∼ 0.3–0.5) inner-target-directed (ITD) background plasma flows at the crown of single-null discharges. This study reports a set of DIVIMP simulations for two similar H-mode discharges from the DIII-D W metal rings campaign differing primarily in B T -direction to assess the effect that fast ITD flows have on the distribution of W ions in the SOL. It is found that for imposed ITD flows of M = 0.3, W ions that otherwise accumulate due to the ∇ ‖ T ion -force are largely flushed out. It is also found that doubling the radial diffusion coefficient from 0.3 to 0.6 m 2 s −1 prevents accumulation due to rapid cross-field transport into the far-SOL, where background plasma flows drain W ions to the divertors. Far-SOL W distributions from DIVIMP are then used to specify input to the impurity transport code 3DLIM, which is used to interpretively model collector probe (CP) deposition patterns measured in the ‘wall-SOL’. It is demonstrated that the deposition patterns are consistent with the DIVIMP predictions of near-SOL accumulation for the unfavorable- B T direction, and little/no accumulation for the favorable- B T direction. The wall-SOL CPs have thus provided the first experimental evidence, albeit indirect, of near-SOL W accumulation—finding it occurs for the unfavorable- B T direction only. For the favorable- B T direction, fast flows can largely prevent accumulation from occurring.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Crystal Chemistry, Band-Gap Red Shift, and Electrocatalytic Activity of Iron-Doped Gallium Oxide Ceramics

This work for the first time unfurls the fundamental mechanisms and set the stage for an approach to derive electrocatalytic activity, which is otherwise not possible, in a traditionally known wide band gap oxide material. Specifically, we report on the tunable optical properties, in terms of wide spectral selectivity and red-shifted band gap, and electrocatalytic behavior of iron (Fe) doped gallium oxide (ß-Ga 2 O 3 ) model system. X-ray diffraction (XRD) studies of sintered Ga 2-x Fe x O 3 (GFO) (0.0 ≤ x ≤ 0.3) compounds provide evidence for the Fe 3+ substitution at Ga 3+ site without any secondary phase formation. Rietveld refinement of XRD patterns reveal that the GFO compounds crystallize in monoclinic crystal symmetry with C2/m space group. The electronic structure of the GFO compounds probed using X-ray photoelectron spectroscopy data reveals that, at lower concentrations of Fe exhibits mixed chemical valence states (Fe 3+ , Fe 2+ ), whereas single chemical valence state (Fe 3+ ) is evident for higher Fe content (x = 0.20-0.30). The optical absorption spectra reveal a significant red shift in the optical band gap with Fe doping. The origin of the significant red shift even at low concentrations of Fe (x = 0.05) is attributed to the strong sp-d exchange interaction originated from the 3 d 5 electrons of Fe 3+ . The optical absorption edge observed at ˜450 nm with lower intensity is the characteristic of Fe doped compounds associated with Fe 3+ - Fe 3+ double excitation process. Coupled with optical band gap red shift, electrocatalytic studies of GFO compounds reveals that, interestingly, Fe doped Ga 2 O 3 compound exhibits electrocatalytic activity in contrast to intrinsic Ga 2 O 3 . Fe doped samples (GFO) demonstrated appreciable electrocatalytic activity towards the generation of H 2 through electrocatalytic water splitting. An onset potential and tafel slope of GFO compounds includes: ~900 mV, ~210 mVdec -1 (x = 0.15) and ~1036 mV, ~290 mVdec-1 (x = 0.30). The electrocatalytic activity of Fe doped Ga-oxide compounds is attributed to cumulative effect of different mechanisms such as; doping resulted new catalytic centers, enhanced conductivity and electron mobility. Hence, in this report, for the first time, we explored a new pathway; electrocatalytic behavior of Fe doped Ga 2 O 3 resulted due to Fe chemical states, red shift in optical band gap. The implications derived from this work may be applicable to a large class of compounds and further options may be available to design functional materials for electrocatalytic energy production.

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

Valley phonons and exciton complexes in a monolayer semiconductor

The coupling between spin, charge, and lattice degrees of freedom plays an important role in a wide range of fundamental phenomena. Monolayer semiconducting transitional metal dichalcogenides have emerged as an outstanding platform for studying these coupling effects. Here, we report the observation of multiple valley phonons – phonons with momentum vectors pointing to the corners of the hexagonal Brillouin zone – and the resulting exciton complexes in the monolayer semiconductor WSe2. We find that these valley phonons lead to efficient intervalley scattering of quasi particles in both exciton formation and relaxation. This leads to a series of photoluminescence peaks as valley phonon replicas of dark trions. Using identified valley phonons, we also uncover an intervalley exciton near charge neutrality. Our work not only identifies a number of previously unknown 2D excitonic species, but also shows that monolayer WSe 2 is a prime candidate for studying interactions between spin, pseudospin, and zone-edge phonons.

36 MATERIALS SCIENCE↗
DOE OSTI2024

High-magnification Faraday rotation imaging and analysis of X-pinch implosion dynamics

An X-pinch load driven by an intense current pulse (>100 kA in ~100 ns) can result in the formation of a small radius, runaway compressional micro-pinch. A micro-pinch is characterized by a hot (>1 keV), current-driven (>100 kA), high-density plasma column (near solid density) with a small neck diameter (1–10 µm), a short axial extent (<1 mm), and a short duration (≲1 ns). With material pressures often well into the multi-Mbar regime, a micro-pinch plasma often radiates an intense, sub-ns burst of sub-keV to multi-keV x rays. A low-density coronal plasma immediately surrounding the dense plasma neck could potentially shunt current away from the neck and thus reduce the magnetic drive pressure applied to the neck. Here, to study the current distribution in the coronal plasma, a Faraday rotation imaging diagnostic (1064 nm) capable of producing simultaneous high-magnification polarimetric and interferometric images has been developed for the MAIZE facility at the University of Michigan. Designed with a variable magnification (1–10×), this diagnostic achieves a spatial resolution of ~35 µm, which is useful for resolving the ~100-μm-scale coronal plasma immediately surrounding the dense core. This system has now been used on a reduced-output MAIZE (100–200 kA, 150 ns) to assess the radial distribution of drive current immediately surrounding the dense micro-pinch neck. The total current enclosed was found to increase as a function of radius, r, from a value of ≈ 50 ± 25 kA at r ≈ 140 µm (at the edge of the dense neck) to a maximal value of ≈150 ± 75 kA for r ≥ 225 µm. This corresponds to a peak magnetic drive pressure of ≈75 ± 50 kbar at r ≈ 225 µm. The limitations of these measurements are discussed in the paper.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Comparison of juvenile white sturgeon blade strike survival through a conventional and a novel hydropower turbine

Hydropower is a globally important renewable energy source with profound impacts on aquatic organisms, including direct injury and mortality to fish that move through hydropower facilities. While it is expected that the characteristics of turbine design and operation influence the type, frequency, and magnitude of injuries suffered by entrained fish, few direct comparisons exist to illustrate differences in fish survival outcomes for conventional and novel turbine designs. This study evaluated blade strike-associated injury and mortality rates for juvenile white sturgeon (Acipenser transmontanus) passed through a model-scale turbine equipped with a runner having conventional blade profiles (thin, straight leading edges) and the same turbine equipped with a runner designed for improved fish survival (thick, slanted leading edges). In both trials, head and turbine runner rotational speeds were matched to produce five blade peripheral speeds between 15.0 and 27.6 m/s. High-speed video was captured for all turbine passage events. Fish were assessed for injuries following passage, and mortalities were assessed immediately after passage and after 48 h. Conventional runner passage resulted in 42% to 78% survival after 48 h, with approximately one-third of all tested fish killed by severing. Under the same test conditions, immediate and 48 h survival rates through the novel runner were 100% except at the highest speed condition (95.6% survival at 48 h). These results clearly indicate that turbine design has a profound effect on fish survival, and suggest that thoughtful redesign of hydropower equipment could significantly ameliorate the risks of blade strike mortality.

13 HYDRO ENERGY↗
DOE OSTI2020

Time-resolved site-selective imaging of predissociation and charge transfer dynamics: the CH 3 I B-band

The predissociation dynamics of the 6s (B 2 E) Rydberg state of gas-phase CH 3 I were investigated by time-resolved Coulomb-explosion imaging using extreme ultraviolet (XUV) free-electron laser pulses. Inner-shell ionization at the iodine 4d edge was utilized to provide a site-specific probe of the ensuing dynamics. The combination of a velocity-map imaging (VMI) spectrometer coupled with the pixel imaging mass spectrometry (PImMS) camera permitted three-dimensional ionic fragment momenta to be recorded simultaneously for a wide range of iodine charge states. In accord with previous studies, initial excitation at 201.2 nm results in internal conversion and subsequent dissociation on the lower-lying A-state surface on a picosecond time scale. Examination of the time-dependent yield of low kinetic energy iodine fragments yields mechanistic insights into the predissociation and subsequent charge transfer following multiple ionization of the iodine products. The effect of charge transfer was observed through differing delay-dependencies of the various iodine charge states, from which critical internuclear distances for charge transfer could be inferred and compared to a classical over-the-barrier model. Time-dependent photofragment angular anisotropy parameters were extracted from the central slice of the Newton sphere, without Abel inversion, and highlight the effect of rotation of the parent molecule before dissociation, as observed in previous works. Our results demonstrate the ability to perform three-dimensional ion imaging at high event rates and showcase the potential benefits of this approach, particularly in relation to further time-resolved studies at free-electron laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Data-driven Modeling for Grid Edge IBRs: A Digital Twin Perspective of User-Defined Models

Recent events in Odessa have brought attention to the challenges associated with the interaction between Inverter- Based Resources (IBRs) and the transmission and distribution system. The NERC event diagnosis report has highlighted sev- eral issues, emphasizing the need for continuous performance monitoring of these IBRs by system operators. Key areas of concern include the mismatch of control and protection perfor- mance of IBRs between the original equipment manufacturer (OEM)-provided models and field measurements. The inability to replicate the realistic response can result in incorrect reliability and resilience studies. In this paper, we developed an approach on how to emulate the behavior of an IBR using measurement data obtained for system operators to utilize in real-time and long- term planning. Two experiments are conducted in the phasor domain and electromagnetic transients (EMT) domain to emulate the behavior for grid forming and grid following inverters under various operating conditions and the effectiveness of the proposed model is demonstrated in terms of accuracy and ease of utilizing user-defined models (UDMs)

Mahapatra, Kaveri [BATTELLE (PACIFIC NW LAB)]↗