Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “HEIGHTS”

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 343 records · Page 19

Ion Beam Figuring System for Synchrotron X-Ray Mirrors Achieving Sub-0.2-µrad and Sub-0.5-nm Root Mean Square

Abstract Optics with high-precision height and slope are increasingly desired in numerous industrial fields. For instance, Kirkpatrick–Baez (KB) mirrors play an important role in synchrotron X-ray applications. A KB system is composed of two aspherical, grazing-incidence mirrors used to focus an X-ray beam. The fabrication of KB mirrors is challenging due to the aspherical departure of the mirror surfaces from base geometries and the high-quality requirements for slope and height residuals. In this paper, we present the process of manufacturing an elliptical cylinder KB mirror using our in-house-developed ion beam figuring (IBF) and metrology technologies. First, the key aspects of figuring and finishing processes with IBF are illustrated in detail. The effect of positioning error on the convergence of the residual slope error is highlighted and compensated. Finally, inspection and cross-validation using different metrology instruments are performed and used as the final validation of the mirror. Results confirm that relative to the requested off-axis ellipse, the mirror has achieved 0.15-µrad root mean square (RMS) and 0.36-nm RMS residual slope and height errors, respectively, while maintaining the initial 0.3-nm RMS microroughness.

36 MATERIALS SCIENCE↗

Insights into factors that affect non-Arrhenius migration of a simulated incoherent Σ3 grain boundary

Non-Arrhenius grain boundary migration, sometimes referred to as antithermal migration where temperature and GB velocity values are inversely related to each other, is examined in an incoherent twin Σ3 [111] 60° (11 8 5) nickel grain boundary. Molecular dynamics is used to simulate migration and examine the effect of various factors on the migration, including interatomic potential, system size, driving force, and variation of atomic grain boundary structure. A classical model for grain boundary migration, in its unsimplified form, is used to analyze the results. The grain boundaries exhibit migration mechanisms with very low apparent barrier heights to migration. As a result, the boundaries migrate quickly but exhibit a velocity saturation similar to that of dislocations. Furthermore, the various interatomic potentials exhibit different migration velocities, but their similarities suggest they all predict similar overall behaviors of migration. The variation of atomic structure in the same incoherent twin grain boundary leads to diverse behaviors with barrier heights that vary from non-Arrhenius to Arrhenius migration. Facet nucleation is confirmed not to be a requirement for this boundary based on an examination of simulation cell size; however, the presence and/or interaction between numerous facets does suggest a slowing and increased barrier height to migration for larger boundaries.

36 MATERIALS SCIENCE↗

Controls and process planning strategies for 5-axis laser directed energy deposition of Ti-6Al-4V using an 8-axis industrial robot and rotary motion

Interest in multi-axis additive manufacturing has grown significantly in recent years due to associated benefits that range from improved quality to support for greater geometric complexity. A method of 5-axis Laser Hot-Wire Directed Energy Deposition is developed and demonstrated in this study using an 8-axis industrial robot, comprised of a 6-axis serial link robot and a 2 degree-of-freedom positioner. Coordinated 5-axis motion is used to print Ti-6Al-4V objects predominantly through build plate rotation. Methods of slicing and toolpath generation are presented, including the use of an alternating toolpath direction philosophy and the manipulation of the CAD model surface normal vectors for determining the required build plate tilt angles for maintaining consistent deposition orientations relative to gravity. Another objective of the study was to determine the limits of conventional slicing and layer height control techniques for printing geometries with significant unsupported overhangs. A comparison is made between the behavior of normal and conformal bead alignment in this regard, and the performance of the layer height controller is reported as it relates to process stability and the limits of printing overhang geometries. Unsupported overhangs were successfully printed at 35° without build plate tilt and 45° with build plate tilt. Toolpath velocity and net-shape printing performance are reported as well, with 3D scans of the printed surfaces used for net-shape characterization. In all, the printing of three geometries is presented, including the overhang geometry, a 5-axis toolpath geometry, and a large-scale, multi-bead demonstrator part, which exceeded 500 mm in height and weighed 44 kg (97 lbs).

36 MATERIALS SCIENCE↗

Turbulence regimes in the nocturnal roughness sublayer: Interaction with deep convection and tree mortality in the Amazon

We investigated the influence of seasonality and proximity to the forest canopy on nocturnal turbulence regimes in the roughness sublayer of a Central Amazon forest. Since convective systems of different scales are common in this region, we also analyzed the effect of extreme wind gusts (propagated from convective downdrafts) on the organization of the turbulence regimes, and their potential to cause the mortality of canopy trees. Our data include high-frequency winds, temperature and ozone concentration at different heights during the dry and wet seasons of 2014. In addition, we used critical wind-speed data derived from a tree-winching experiment and a modeling study conducted in the same study site. Two different turbulence regimes were identified at three heights above the canopy: a weakly stable (WS) and a very stable regime (VS). The threshold wind speeds that mark the transition between turbulence regimes were larger during the dry season and increased as a function of the height above the canopy. The turbulent fluxes of sensible heat and momentum during the WS accounted for 88% of the entire nighttime flux. Downdrafts occurred only in the WS and favored a fully coupled state of wind flow along the canopy profile. The destructive potential of winds was four times higher than on nights without downdrafts.

54 ENVIRONMENTAL SCIENCES↗

Scale effects on core design, fuel costs, and spent fuel volume of pressurized water reactors

The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.

42 ENGINEERING↗

Evaluation of riparian enhancement actions in the Columbia River Basin

Riparian enhancement is a common restoration technique in the Columbia River Basin (CRB) and the Pacific Northwest. However, relatively few studies have evaluated its success and even fewer studies include long-term monitoring. Forty-one riparian planting projects located in the CRB, each with a paired treatment and control reach, were evaluated in the summer of 2018 and 2019 using an extensive post-treatment (EPT) design. At each reach, we quantified woody plant abundance, richness, diversity, and vegetation and canopy cover as riparian response variables, and measured terrace height and documented individuals with bud browse, deceased, and with predator protection as potential explanatory variables. Species richness, woody plant abundance (all height classes combined and within shrub height class), and the proportion of woody plants with bud browse or deceased were all higher in treatment than control reaches. However, no significant improvements were observed for any other riparian response or explanatory variables. Implementation techniques, including invasive species removal efforts, follow-up planting, watering, and the type of restoration project (floodplain restoration and planting, in-stream restoration and planting, or solely planting) were also considered as potential explanatory variables and were found to have significant impacts on restoration response for abundance and cover variables. Further, our results suggest that planting implementation methods, site level physical factors, and time since planting all influence the success of riparian planting projects and if not addressed in restoration design and implementation can contribute to a lack of detectable response.

54 ENVIRONMENTAL SCIENCES↗

A new phenomenological model to describe root-soil interactions based on percolation theory

In his paper on net primary productivity of terrestrial communities predicted from climatological data, Rosenzweig (1968) argued that variability in productivity is well accounted for by (evapo)-transpiration, and that water from transpiration is, on global scales, the most variable component in the photosynthesis reaction. The goal of this paper is to investigate whether variability in plant growth on local scales and within species is primarily related to transpiration under several scenarios including different terrain curvature, slope aspect, soil characteristics, and climate ranges. Here, we test the hypothesis that this relationship exists because root growth into the surface soil layers (0–2 m) tends to follow paths with minima in resistance, which in turn maximizes water flow and nutrient delivery rates that regulate growth. The set of all connected paths with individual pore-to-pore flow resistances less than a critical, percolating, value forms a cluster with mass fractal dimensionality, d f . We propose that roots follow paths through the 2D percolation cluster, defining the set of all optimal flow paths, making the 2D value of d f from percolation relevant to root fractal dimensionality. The tortuosity of such optimal paths as defined in percolation theory should then relate root length to root radial extent, linking the parameters of root tortuosity and plant productivity. Our analysis of large data sets across species implies that root radial extent and tree height are both proportional to cumulative transpiration until trees approached maximum height, and their growth rates are proportional to the transpiration rate, not to the moisture content. Local variations in tree height as functions of the variables investigated appear generally consistent with deduced variations in transpiration. Here this correlation is investigated more closely in the context of studies addressing individual tree species.

54 ENVIRONMENTAL SCIENCES↗

Process design of the piperazine advanced stripper for a 460 MW NGCC

The Texas Carbon Management Program is conducting a Front-End Engineering Design (FEED) for Piperazine with the Advanced Stripper™ (PZAS™) at a natural gas combined cycle power plant in west Texas. Here, this paper presents the design of the advanced stripper system with an emphasis on the selection of packing type, packing height, column diameter, and heat exchanger size, demonstrates the effect of ambient temperature swing on the regeneration energy requirement, and investigates two alternate process configurations that provide better energy performance than PZAS™. The design uses 2 stripper trains, and each train uses two cold cross exchangers, one hot cross exchanger, a stripper column, a CO 2 exchanger, and a condenser. Each column consists of two sections of random packing with RSR 2 packing at the top and RSR 3 packing at the bottom, with a total packed height of 7.5 m. The cold exchanger was designed with an NTU of 6, the hot exchanger with an NTU of 2, and the CO 2 exchanger with an NTU of 5. At 0.4 rich loading and 0.2 lean loading, this design resulted in a heat duty of 2.85 GJ/tonne CO 2 . The heat duty can vary up to 3.1 GJ/tonne CO 2 due to lower rich loading at higher summer temperature, but the absolute steam requirement is 8–15% lower in the summer. An alternate design, PZAS with hot rich bypass, is proposed, which can reduce the equivalent work of stripping by 1.5–3.5% compared to PZAS at low heat exchanger size, due to a reduction in stripper driving force. Using a hot rich bypass increases the vapor fraction of the total rich bypass stream from 0.3 to 1 mol% compared to using the warm rich bypass, which also contributes to a reduction in lost work in the column. PZAS with hot rich bypass outperforms PZAS in energy requirement below a packing height of 15 m, but fails to do so above 15 m of packing due to a column pinch that results in little change in LMTD of the column. PZAS with hot rich bypass has the potential to minimize capital cost in addition to reducing energy requirement as it requires less packing to strip CO 2 than PZAS and requires only one cross exchanger with 8 heat transfer units.

03 NATURAL GAS↗

Direct metal contacts printing on 4H-SiC for alpha detectors and inhomogeneous Schottky barriers

Electrical characterization of several 4H-SiC Schottky diodes with Aerosol-Jet printed gold (Au), silver (Ag), and platinum (Pt) contacts was performed using forward and reverse current–voltage (IV) measurements. From these measurements, device parameters such as Schottky barrier height and ideality factor were determined; however, many of the devices exhibited nonideal behavior and inferior performance characterized by ideality factors significantly greater than one, disproportionate low voltage leakage current and low barrier heights. Forward current–voltage (FIV) characteristics were fitted to an inhomogeneous barrier height theory to explain the abnormal behavior exhibited by the printed devices. Here, transmission electron microscopy (TEM) of the device cross-sections was performed to investigate the printed metal and semiconductor epitaxial layer interface, which revealed the imperfections in the metal–semiconductor contact. Despite these irregularities, alpha radiation detection capability of these devices was still achieved with an energy resolution of 1.89% at 5.486 MeV, and the best achievable resolution at 0.51% with no energy degradation of 5.486 MeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Convolution-based frequency domain multiplexing of SiPM readouts using the $\mathrm{DRS4}$ digitizer

Here we present 4:1 multiplexing of organic scintillators, each coupled to a silicon photomultiplier (SiPM), to reduce the need for a large number of digitizer input channels to readout highly pixelated radiation detection systems. Frequency domain multiplexing (FDM) encodes a detector pulse by assigning it a unique frequency via convolution before combining the encoded signal into a single channel. The combined signal is then read through a digitizer input channel. We have designed an FDM system to multiplex four SiPMs using DRS4 digitizer evaluation board from Paul Scherrer Institute (PSI). We demonstrate 4:1 multiplexing of the SiPM fast output signals and pulse recovery from the digitized multiplexed signal using deconvolution. The noise in the recovered pulse introduces a bias and uncertainty in the estimate of energy and timing that changes with pulse height. The relative uncertainty in the estimated energy from the recovered pulse decreases with pulse height with a maximum uncertainty of 3.1% for the low energy pulses (corresponding to 100 keV); the uncertainty in the estimated time pick-off also decreases with pulse height with a maximum uncertainty of 110 ps for the low energy pulses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparative study with two reduced test facilities for air-cooled RCCS scaling law

Korea Atomic Energy Research Institute (KAERI) and Argonne National Laboratory (ANL) conducted steady-state tests at constant heat flux conditions to investigate the scale effect on heat removal behavior of the air-cooled Reactor Cavity Cooling System (RCCS). Two differently scaled-down test facilities were used in this study, KAERI’s ¼ scale with 4.0-m cavity height, and ANL’s ½ scale with 6.8-m cavity height. Two scaling laws were proposed to simulate the radiation across the cavity and the buoyancy-driven duct flow in the riser, respectively. The test matrix for KAERI and ANL test facilities focused on the scaling effect of the heated riser length with constant heat flux at the PMR200 RCCS design. The test results showed that mixed convection in the riser duct is an important factor for accurately extrapolating the thermo-fluid behavior in the prototype from the test results in the scale-down facilities. The system analysis code, GAMMA+, with improved heat transfer models predicted fairly well the air-cooled RCCS test data from two facilities. GAMMA+ analysis results showed that the predicted radiation fractions on the heated plate in the scale-down test conditions were larger than those on the reactor vessel in the prototype. The scaling law for air-cooled RCCS was improved by considering the mixed convection in the riser duct and the same radiation fraction on the heated plate. The mixed convection effect was calculated by the height ratio. The same radiation fraction provided a conservative extrapolation from the scale-down test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance Evaluation of an Offshore Wave Measurement Buoy in Monochromatic Waves

The accurate measurement of waves underpins marine energy resource characterization, device design, and project development. Datawell wave buoys are widely deployed and have long served as a trusted standard for wave measurements. We quantify the measurement performance, including wave elevation and energy flux estimation, of a Datawell DWR-MkIII buoy using prescribed monochromatic heave motions on a large-amplitude six-degree-of-freedom motion platform at the National Laboratory of the Rockies, assuming the buoy behaves as an ideal wave follower. Commanded motions were validated with an optical motion tracking system while buoy elevation and raw acceleration were recorded. Wave elevations were propagated to wave energy flux estimation using four methods, including one frequency-domain method and three time-domain methods. The Bayesian optimization was applied for design of experiments, and records from three test sites were also applied and evaluated in the present study. Results show two error regions within the nominal period range of 1.6 s to 30 s. For wave periods between 5 s and 25 s, the buoy provides accurate wave height measurements. For short periods less than 5 s, the 1.28 Hz sampling frequency induces sub-Nyquist artifacts that bias elevation and can drive maximum energy flux estimation errors above 100%. For long periods exceeding 25 s, the buoy reported elevation is underpredicted with error depending on period but relatively independent of wave height, with maximum wave height and wave energy flux errors reaching 64% and 87%, respectively. Furthermore, analysis of three field-derived cases shows that frequency-domain estimates at 1.28 Hz agree within 2% of the corresponding 100 Hz estimates, while larger method-dependent differences are observed for the Hilbert method.

16 TIDAL AND WAVE POWER↗

Collimated phase measuring deflectometry

Phase measuring deflectometry has been applied for free-form specular surface metrology, but its measured slope results are sensitive to the depth of sample positioning, which is also called the height-slope ambiguity. The objective of this work is to tackle this height-slope ambiguity problem. The main idea is to introduce collimated camera rays using a telecentric imaging lens and collimated structured-light illumination with a Fourier lens. This setup makes the fringe phases become only sensitive to the surface slopes and insensitive to the depth of the sample positioning. In this way, the slope calculation is theoretically independent of the sample depth. We call this new deflectometry technique Collimated Phase Measuring Deflectometry (CPMD). With our developed CPMD experimental setup, the measurement is insensitive to the depth of sample positioning, e.g., the measured height dispersion is less than 30 nm RMS within a 10 mm depth range when measuring a 50-mm-diameter spherical mirror with a 200 mm radius of curvature. In conclusion, the merits and limitations of the proposed CPMD technique are discussed, revealing its prospects in practical metrology applications and potential future investigations.

36 MATERIALS SCIENCE↗

Real-Time Defect Correction in Large-Scale Polymer Additive Manufacturing via Thermal Imaging and Laser Profilometer

Defects can result in a failed part and are costly in terms of time and material. This cost is even greater in the context of large-scale additive manufacturing where the objects can be very large. As a result, a great deal of research has focused on defect identification and mitigation. To address defects during object construction, researchers at Oak Ridge National Laboratory’s Manufacturing Demonstration Facility investigated an in-situ control system comprised of two sensors: a thermal camera and laser profilometer. This control system adjusted material flow and build speed to mitigate three types of defects: low layer times, underfill, and overfill. Several test objects were constructed. The control system was found to adjust build parameters to handle low layer times of approximately 15 seconds and height deviations from -100% underfill (the absence of a layer) to 50% overfill. Within two layers, height deviations could be returned to within 10% of the expected layer height. Further, preliminary results suggest the system can compensate for uneven build surfaces.

Borish, Michael↗

Aboveground biomass density models for NASA’s Global Ecosystem Dynamics Investigation (GEDI) lidar mission

NASA's Global Ecosystem Dynamics Investigation (GEDI) is collecting spaceborne full waveform lidar data with a primary science goal of producing accurate estimates of forest aboveground biomass density (AGBD). This paper presents the development of the models used to create GEDI's footprint-level (~25 m) AGBD (GEDI04_A) product, including a description of the datasets used and the procedure for final model selection. The data used to fit our models are from a compilation of globally distributed spatially and temporally coincident field and airborne lidar datasets, whereby we simulated GEDI-like waveforms from airborne lidar to build a calibration database. We used this database to expand the geographic extent of past waveform lidar studies, and divided the globe into four broad strata by Plant Functional Type (PFT) and six geographic regions. GEDI's waveform-to-biomass models take the form of parametric Ordinary Least Squares (OLS) models with simulated Relative Height (RH) metrics as predictor variables. From an exhaustive set of candidate models, we selected the best input predictor variables, and data transformations for each geographic stratum in the GEDI domain to produce a set of comprehensive predictive footprint-level models. We found that model selection frequently favored combinations of RH metrics at the 98th, 90th, 50th, and 10th height above ground-level percentiles (RH98, RH90, RH50, and RH10, respectively), but that inclusion of lower RH metrics (e.g. RH10) did not markedly improve model performance. Second, forced inclusion of RH98 in all models was important and did not degrade model performance, and the best performing models were parsimonious, typically having only 1-3 predictors. Third, stratification by geographic domain (PFT, geographic region) improved model performance in comparison to global models without stratification. Fourth, for the vast majority of strata, the best performing models were fit using square root transformation of field AGBD and/or height metrics. There was considerable variability in model performance across geographic strata, and areas with sparse training data and/or high AGBD values had the poorest performance. These models are used to produce global predictions of AGBD, but will be improved in the future as more and better training data become available.

54 ENVIRONMENTAL SCIENCES↗

A scale-wise analysis of intermittent momentum transport in dense canopy flows

We investigate the intermittent dynamics of momentum transport and its underlying time scales in the near-wall region of the neutrally stratified atmospheric boundary layer in the presence of a vegetation canopy. This is achieved through an empirical analysis of the persistence time scales (periods between successive zero-crossings) of momentum flux events, and their connection to the ejection–sweep cycle. Using high-frequency measurements from the GoAmazon campaign, spanning multiple heights within and above a dense canopy, the analysis suggests that, when the persistence time scales ( $t_p$ ) of momentum flux events from four different quadrants are separately normalized by $\varGamma _{w}$ (integral time scale of the vertical velocity), their distributions $P(t_p/\varGamma _{w})$ remain height-invariant. This result points to a persistent memory imposed by canopy-induced coherent structures, and to their role as an efficient momentum-transporting mechanism between the canopy airspace and the region immediately above. Moreover, $P(t_p/\varGamma _{w})$ exhibits a power-law scaling at times $t_{p}<\varGamma _{w}$ , with an exponential tail appearing for $t_{p} \geq \varGamma _{w}$ . By separating the flux events based on $t_p$ , we discover that around 80 % of the momentum is transported through the long-lived events ( $t_{p} \geq \varGamma _{w}$ ) at heights immediately above the canopy, while the short-lived ones ( $t_{p} < \varGamma _{w}$ ) only contribute marginally ( $\approx 20\,\%$ ). To explain the role of instantaneous flux amplitudes in momentum transport, we compare the measurements with newly developed surrogate data and establish that the range of time scales involved with amplitude variations in the fluxes tends to increase as one transitions from within to above the canopy.

Mechanics↗

MgO Nanostructures on Cu(111): Understanding Size- and Morphology-Dependent CO 2 Binding and Hydrogenation

To design and optimize cost-effective technologies for the capture, utilization, and storage of carbon dioxide (CO 2 ), we need a fundamental knowledge and control of chemical interactions associated with the capture and conversion of the molecule into high-value chemicals, minerals, and all kinds of materials. Bulk magnesium oxide (MgO) is frequently used for the trapping and storage of CO 2 by generation of magnesium carbonates. In this study, the growth and reactivity of MgO nanostructures on a Cu 2 O/Cu(111) substrate were investigated using scanning tunneling microscopy (STM) and synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS). For extremely small concentrations of Mg (~ 0.01 monolayer (ML)), a well-ordered film of copper oxide with small clusters (0.2-0.5 nm in width, 0.4-0.6 Å in height) of embedded MgO was seen. At a coverage of 0.1 ML, MgO nanoparticles with a width of 0.4 to 1 nm and a height of ~ 1.5 Å were randomly distributed on the copper oxide. Further, random distribution was also observed when the MgO coverage was raised to 0.25 ML, with the width of the MgO particles increasing to 2-2.5 nm and the height reaching 2 Å. These oxide nanostructures displayed a high reactivity towards CO 2 and H 2 that is not seen for bulk MgO. Dissociation of H 2 was observed at room temperature with reaction of the H adatoms with CuO x and C-containing groups. On the small MgO nanostructures (< 1 nm in width), instead of plain carbonate formation, there was dissociation of CO 2 into CO and C species, opening reaction channels for the conversion of this harmful molecule into oxygenates and light alkanes.

03 NATURAL GAS↗

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗