Search NASA⌕ Search

SEARCH · Search NASA

Results for “vacuum simulation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical analysis of gas exhaust in Wendelstein 7-X using the direct simulation Monte Carlo method

The present work is focused on a 3D numerical assessment of the Wendelstein 7-X (W7-X) particle exhaust. For all the numerical simulations the direct simulation Monte Carlo solver of the DIVGAS workflow, has been employed. The complex 3D geometry of the sub-divertor region includes the pumping gap panel, supporting structures, cooling pipes as well as the cryo-vacuum pump. All the considered flow simulations correspond to the Standard magnetic configuration of W7-X. The main conclusions, which can be extracted from the present numerical analysis could be summarized as follows; The coupling between EMC3-EIRENE and DIVGAS, which considers the fact that the incoming neutral particle flux at the sub-divertor is based on realistic plasma background, has been demonstrated. Three plasma scenarios have been considered, for which is clearly seen that by increasing the heating power, the neutral pressure as well as the resulting pumping efficiency is increased. The obtained numerical results of the neutral pressure in the sub-divertor lie within a more general scan matrix, which assumes a wider range of incoming particle flux, namely 10 19 –10 24 (s −1 ). It has been observed that, the sub-divertor neutral pressure is proportional to the incoming neutral particle flux, with the effective pumping speed to be a constant of proportionality. The influence of switching off the cryo-vacuum pump on the sub-divertor pressure is rather modest and a weak increase of the neutral pressure in the sub-divertor is expected. Correlations of the sub-divertor pressure with the total incoming particle flux as well as the individual pumped flux at each of the AEH and AEP sections have been deduced. Moreover, it has been demonstrated that the influence of the incoming neutral particle flux on the albedo coefficient at the AEH and AEP pumping gaps is rather weak. All the above numerical findings will actively support the optimization of the W7-X particle exhaust, in view of future experimental campaigns.

DSMC method↗
DOE OSTI2025

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Fast ion confinement in the presence of core magnetic islands in Wendelstein 7-X

The effect of magnetic islands in the core region of Wendelstein 7-X (W7-X) on fast ion confinement is explored through simulations with the BEAMS3D code. A magnetic configuration where the n/m = 5/5 island chain is shifted to r/a ~ 0.7 allows the exploration of core island physics in W7-X. The control coil system on W7-X allows the tuning of the island size either increasing the island width or decreasing it. A coupling of the BEAMS3D code to the FIELDLINES code provides a versatile mechanism for incorporating magnetic islands and stochastic regions into the BEAMS3D code. Collisionless simulations suggest that the presence of core islands degrade the confinement of passing particles in the region of the island chain. Full neutral beam simulations of W7-X show a similar behavior with confinement decreasing as the island width is increased. Comparisons between a vacuum magnetic field and low beta HINT2 simulation are made showing similar fast ion behavior. Measurements of lost fast ions in W7-X confirm this trend with the control coil suppressed island configuration showing lower losses than that with no control coils applied. Simulations of fast ion wall loads are performed suggesting no drastic change in loss pattern and a slight reduction in losses with minimized islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Organic Evaporation, Oxidation, and Hydrolysis Testing in Support of Hanford Sample-and-Send

The Hanford site has approximately 54 to 56 million gallons of radioactive mixed waste stored in 156 unretrieved underground storage tanks. The Hanford Waste Treatment and Immobilization Plant (WTP) is being built to treat and immobilize the tank waste. The baseline method for immobilization of Low Activity Waste (LAW) through the WTP is vitrification, but additional immobilization capacity is needed to supplement the initial LAW melters. An alternative cementitious waste form is being investigated for that future immobilization method. However, one impediment to a cementitious waste form is the presence of Land Disposal Restricted (LDR) organic chemicals in tank waste, which are regulated on a concentration based standard in the final waste form. Hence, if the quantity of organics in LAW is high enough, they must be destroyed or removed to make a waste form compatible with disposal in a mixed low level waste landfill. This work evaluates potential avenues for treatment of LDR organics to eliminate the impediment and permit possible use of a cementitious waste form. Vacuum evaporation testing to remove LDR organics consisted of preparing a non-radioactive LAW simulant, spiking that simulant with organic chemicals, and evaporating the mixture via differential distillation. The apparatus was a laboratory-scale vacuum evaporator operated at 60 ±5 torr absolute (vacuum evaporation). The LAW simulant represented the liquid expected to be retrieved from the Hanford tank farms at approximately 4.0 M [Na+] total sodium ion concentration. The concentration of the organic chemicals added was significantly higher than typically found in the tank waste samples since the higher levels were necessary to assist in analytical measurement and tracking of the spiked species.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗
DOE OSTI2022

Three-Dimensional Permeability of Thick-Section Glass Fabric Reinforced Polymer Composite by Vacuum-Assisted Resin Infusion Molding

Abstract Determination of permeability of thick-section glass fabric preforms with fabric layers of different architectures is critical for manufacturing large, thick composite structures with complex geometry, such as wind turbine blades. The thick-section reinforcement permeability is inherently three-dimensional and needs to be obtained for accurate composite processing modeling and analysis. Numerical simulation of the liquid stage of vacuum-assisted resin infusion molding (VARIM) is important to advance the composite manufacturing process and reduce processing-induced defects. In this research, the 3D permeability of thick-section E-glass fabric reinforcement preforms is determined, and the results are validated by a comparison between flow front progressions from experiments and from numerical simulations using ansys fluent software. The orientation of the principal permeability axes were unknown prior to experiments. The approach used in this research differs from those in literature in that the through-thickness permeability is determined as a function of flow front positions along the principal axes and the in-plane permeabilities and is not dependent on the inlet radius. The approach was tested on reinforcements with fabric architectures which vary through-the-thickness direction, such as those in a spar cap of a wind turbine blade. The computational simulations of the flow-front progression through-the-thickness were consistent with experimental observations.

Engineering↗
DOE OSTI2021

Coalescence and splashing threshold for head-on collisions of liquid metal nanodroplets

Head-on collisions of liquid metal nanodroplets in a vacuum are investigated through molecular dynamics simulations in order to determine the transition threshold between the coalescing and splashing regimes for six different materials (aluminum, calcium, cerium, gold, platinum, and tin). Droplets of various sizes and initial speeds are simulated, and it is found that the Reynolds and Ohnesorge numbers are able to predict the transition between the coalescing and splashing regimes. An energy balance for coalescing droplets shows that the initial energy is mainly converted to thermal energy increasing the temperature of the combined droplets by several hundred to several thousand kelvin depending on the material, and this result is confirmed in the simulations. Furthermore, when splashing occurs, the number of smaller droplets formed and the spreading rate are found to be dependent on the initial size and initial speed of the original droplets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Organic Evaporation and Oxidation Testing in Support of Hanford Sample-and-Send

The Hanford site has approximately 56 million gallons of radioactive mixed waste stored in 177 underground storage tanks. The Hanford Waste Treatment and Immobilization Plant (WTP) is being built to treat and immobilize the tank waste. The baseline method for immobilization of Low Activity Waste (LAW) through the WTP is vitrification, but additional immobilization capacity is needed to supplement the initial LAW melters. An alternative cementitious waste form is being investigated for that future immobilization method to supplement vitrification. However, one impediment to a cementitious waste form is the presence of Land Disposal Restricted (LDR) organic chemicals in tank waste. This work evaluates potential avenues to eliminate that impediment to permit possible use of a cementitious waste form and work towards a decision whether additional LDR organic pretreatment would be required. Savannah River National Laboratory (SRNL) performed testing using simulants to examine evaporation as a method to remove some prevalent organics from LAW. Spiking the caustic LAW simulant with selected regulated organic chemicals found one that clearly decomposes because of caustic instability. Oxidation testing of other organic chemicals found some LDR organics degrade as desired and others are stable in the presence of peroxide and permanganate. In addition to studies with simulants, a literature review was performed to evaluate radiological stability of LDR organics. Descriptions of the experimental details, equipment, and results are included in this report. Evaporation testing consisted of preparing the LAW simulant, spiking that simulant with organic chemicals, and evaporating the mixture via differential distillation. The apparatus was a laboratory-scale vacuum evaporator operated at 60 ±5 torr absolute (vacuum evaporation) and also at atmospheric pressure. The LAW simulant represented the liquid expected to be retrieved from the Hanford tank farms at approximately 4.0 M [Na + ] total sodium ion concentration. The concentration of the organic chemicals added was significantly higher than typically found in the tank waste samples since the higher levels were necessary to assist in analytical measurement and tracking of the spiked species. Organic chemicals were chosen for the work with a consideration of how their volatility compares with that of methanol. This was done by comparing the ratio of the pure water Henry’s law coefficient (K h ) of methanol to that of the compound in question (hereafter termed the K h ratio), where ratios above unity indicated less volatility than methanol. Methanol was chosen because it is a common regulated chemical with relatively low volatility but which has been removed by evaporation in previous laboratory work. While organic separation results depend on evaporator design, laboratory experiments verified that organic partitioning to the overhead condensate stream by evaporation is a practical process. The work reported here found difficulties in quantitative analysis of the organic chemicals in aqueous samples. Most of the time there was insufficient analysis to close a mass balance for evaporator runs, but qualitative evidence of carryover was obtained. The methods were also able to show whether organic chemicals were susceptible or resistant to solution oxidation in permanganate or hydrogen peroxide tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Study of vacuum ultraviolet emission in helium and helium/nitrogen mixtures

Helium is frequently used as a working medium for the generation of plasmas and is capable of energetic photon emissions. These energetic photon emissions are often attributed to the formation of helium excimer and subsequent photon emission. When the plasma device is exposed to another gas, such as nitrogen, this energetic photon emission can cause photoionization and further ionization wave penetration into the additional gas. Often ignored are the helium resonance emissions that are assumed to be radiation trapped and therefore not pertinent to photoionization. Here, experimental evidence for the presence of helium atomic emission in a pulsed discharge at ten’s of Torr is shown. Simulations of a discharge in similar conditions agree with the experimental measurements. In this context, the role of atomic and molecular helium light emission on photoionization of molecular nitrogen in an ionization wave is studied using a kinetic modeling approach that accounts for radiation dynamics in a developing low-temperature plasma. Three different mixtures of helium at a total pressure of 250 Torr are studied in simulation. Photoionization of the nitrogen molecule by vacuum ultraviolet helium emission is used as the only seed source ahead of the ionization front. It is found that even though radiation trapped, the atomic helium emission lines are the significant source of photoionization of nitrogen. The significant effect of radiation trapped photon emission on ionization wave dynamics demonstrates the need to consider these radiation dynamics in plasma reactors where self-absorbed radiation is ignored.

42 ENGINEERING↗
DOE OSTI2023

Outgassing behavior and heat treatment optimization of JSC-1A lunar regolith simulant

As NASA strives towards a long duration presence on the Moon, it has become increasingly important to learn how to better utilize resources from the lunar surface for everything from habitats, vehicle infrastructure, and chemical extraction. To that end, a variety of lunar simulants have been sourced from terrestrially available volcanic minerals and glass as Apollo regolith is unavailable for experimentation needing large masses. However, while mineralogy and chemical composition can approach that of lunar material in these simulants, there are still distinct non-lunar phases such as hydrates, carbonates, sulfates, and clays that can cause simulants to behave distinctly non-lunar in a variety of processing conditions that maybe applied in-situ to lunar material. Notably, severe glassy bubbling has been documented in a variety of vacuum sintering experiments on JSC-1A lunar mare simulant heated via microwaves. The origins of this outgassing have not been well understood but are normally attributed to the decomposition of non-lunar contaminates intrinsic to virtually all terrestrially sourced simulants. As such, a series of controlled environmental tests were performed to ascertain the origins of the high temperature outgassing and to develop heat treatments that can drive JSC-1A closer to lunar composition and behavior. It was found that in JSC-1A at elevated temperatures distinct gas evolutions of water, carbon dioxide, and sulfur dioxide occur in both inert gas and vacuum. Additionally, the presence of hydrogen during heat treatments was shown to dramatically change gas evolutions, leading to distinctly more lunar-like composition and behavior from JSC-1A simulant.

58 GEOSCIENCES↗
DOE OSTI2023

Design of a Nanophotonic Hybrid Coupler

This project aims to design and simulate a nanophotonic hybrid coupler for detecting vacuum ultraviolet (VUV) light with higher efficiency. VUV light is of significant interest because of its applications ranging from high energy physics to space science and electronic industry. However, current VUV photodetectors have very low efficiencies because most materials absorb in the VUV region (100-200 nm). To address this challenge, I have designed a plasmonic-dielectric coupler with wavelength shifting properties, enabling high Purcell enhancement and far-field confinement. The coupler consists of a cylindrical cSi resonator on top of a thin ZnO layer, and Ag and SiO2 are used as substrates. Simulations are done using the finite-difference time domain method to solve Maxwell’s equations, and specifically with Meep and Lumerical. The effectiveness of the coupler is characterized by the Purcell factor and far-field radiation. The Purcell factor reaches a maximum near 390 nm and is in the o rder of 10^4. Furthermore, the coupler also allows partial control of the directionality of the far-field radiation. With the presence of the resonator, the radiation is directed towards the coupled metalens. Overall, our design combines the advantages of plasmonic and dielectric couplers, along with the intrinsic wavelength shifting properties of materials to enhance VUV detection significantly, and such a device can be used in, for example, the DUNE project.

Zhou, Jiayang↗
DOE OSTI2022

SimELIT: A Novel GUI-Based Comprehensive Ion Trajectory Simulation Software for Mass Spectrometry

Ion trajectory simulation in mass spectrometry systems from injection to detection is technically challenging but very important for better understanding the ion dynamics in instrument development. Here, in this work, we present SimELIT (Simulator of Eulerian and Lagrangian Ion Trajectories), a novel ion trajectory simulation platform. SimELIT is built upon a suite of multiphysics solvers compiled into OpenFOAM (an open-source numerical solver library particularly used for computational mechanics), with a simple web-based graphical user interface (GUI) allowing users to define the details of OpenFOAM cases and run simulations. SimELIT is a modular program and can provide extensions of physics (e.g., gas flows, electrodynamic fields) and thus enable ion trajectory simulations from the ion source to detector. The current version (SimELIT) provides two numerical solvers for ion trajectory simulations–(1) a Lagrangian particle tracker in vacuum and (2) a Eulerian ion density solver in background gas in the presence of electric fields. Here, we describe the architecture of SimELIT, including its use of Docker and the React Framework, and demonstrate the computation of ion trajectories of multiple m/z values in a static/linear voltage drop in vacuum (across a 1 m long flight tube). Further, the drift motion of ions under 1 Torr pressure conditions in a static background (N 2 ) gas through a 20 V/cm static electric field is shown. The results produced from SimELIT were compared with SIMION and theoretical estimates. In addition, we report the computation of ion trajectories in electrodynamic fields within a planar FAIMS device operating at atmospheric pressure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Minimum Free-Energy Shapes of Ag Nanocrystals: Vacuum vs Solution

Here, we use two variants of replica-exchange molecular dynamics (MD) simulations, parallel tempering MD and partial replica exchange MD, to probe the minimum free-energy shapes of Ag nanocrystals containing 100–200 atoms in a vacuum, ethylene glycol (EG) solvent, and EG solvent with a PVP polymer containing 100 repeat units. Our simulations reveal a shape intermediate between a Dh and an Ih, a Dh-Ih, that has distinct structural signatures and magic sizes. We find several prominent features associated with entropy: pure FCC nanocrystals are less common than FCC crystals containing stacking faults, and crystals with the minimum potential energy are not always preferred over the range of relevant temperatures. The shapes of the nanocrystals in solution are influenced by the chemical identities of the solution-phase molecules. Comparing Ag nanocrystal shapes in EG to those in an EG+PVP solution, we find more icosahedra in EG and more decahedra in EG+PVP across all of the nanocrystal sizes probed in this study. At certain critical sizes, nanocrystal shapes can change dramatically with the addition and removal of a single atom or with a change in temperature at a fixed size. The information in our study could be useful in efforts to devise processing routes to achieve selective nanocrystal shapes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Self-driven ion deflectometry measurements using MeV fusion-driven protons and accelerated deuterons in the deuterated hybrid x-pinch on the MAIZE LTD generator

Abstract We report on the results of point-projection ion deflectometry measurements from a mid-size university z-pinch experiment. A 1 MA 8 kJ LTD generator at the University of Michigan (called MAIZE) drove a hybrid x-pinch (HXP) with a deuterated polyethylene fiber load to produce a point-like source of MeV ions for backlighting. In these experiments, 2.7 MeV protons were generated by DD beam-target fusion reactions. Due to the kinematics of beam-target fusion, the proton energies were down-shifted from the more standard 3.02 MeV proton energy that is released from the center-of-mass rest frame of a DD reaction. In addition to the 2.7 MeV protons, strongly anisotropic beams of 3 MeV accelerated deuterons were detected by ion diagnostics placed at a radial distance of 90 mm from the x-pinch. Numerical reconstruction of experimental data generated by deflected hydrogen ion trajectories evaluated the total current in the vacuum load region. Numerical ion-tracking simulations show that accelerated deuteron beams exited the ion source region at large angles with respect to the pinch current direction.

Physics↗
DOE OSTI2025

Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations

Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.

BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗