Search NASA⌕ Search

SEARCH · Search NASA

Results for “process 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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Axial Tensile Failure Analysis of SiC{sub f}/Ti Composite Based on Continuum Cohesive Zone Model

In this paper, the failure behavior of SiC fiber-reinforced Ti-6Al-4V matrix composite (SiC{sub f}/Ti) under longitudinal loads was studied experimentally and theoretically. Experimental results show that the interface of SiC{sub f}/Ti composite is continuous and integral. Through SEM observations, it was found that transverse matrix cracks initiated from the debonding region of the interface. A three-dimensional representative volume element was developed to simulate the failure process of SiC{sub f}/Ti-6Al-4V composite. A continuum cohesive zone model is employed to describe the debonding behavior at the interface, compared with a discrete one. More accurate simulated results were obtained by a continuum cohesive zone model, corresponding to experimental results. The transverse cracking behavior of the matrix was simulated. The stress–strain curves predicted by the continuum cohesive zone model show a good agreement with the experimental results. The fiber volume fraction, interface strength and interface toughness are investigated to exhibit their influences upon the longitudinal mechanical characteristics of SiC{sub f}/Ti-6Al-4V composite.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Radionuclide gas transport through nuclear explosion-generated fracture networks

Underground nuclear weapon testing produces radionuclide gases which may seep to the surface. Barometric pumping of gas through explosion-fractured rock is investigated using a new sequentially-coupled hydrodynamic rock damage/gas transport model. Fracture networks are produced for two rock types (granite and tuff) and three depths of burial. The fracture networks are integrated into a flow and transport numerical model driven by surface pressure signals of differing amplitude and variability. There are major differences between predictions using a realistic fracture network and prior results that used a simplified geometry. Matrix porosity and maximum fracture aperture have the greatest impact on gas breakthrough time and window of opportunity for detection, with different effects between granite and tuff simulations highlighting the importance of accurately simulating the fracture network. In particular, maximum fracture aperture has an opposite effect on tuff and granite, due to different damage patterns and their effect on the barometric pumping process. From stochastic simulations using randomly generated hydrogeologic parameters, normalized detection curves are presented to show differences in optimal sampling time for granite and tuff simulations. In conclusion, seasonal and location-based effects on breakthrough, which occur due to differences in barometric forcing, are stronger where the barometric signal is highly variable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Estimation of heat transfer coefficients for biomass particles by direct numerical simulation using microstructured particle models in the Laminar regime

Here, direct numerical simulation of convective heat transfer from hot gas to isolated biomass particle models with realistic morphology and explicit microstructure was performed over a range of conditions with laminar flow of hot gas (500 degrees C). Steady-state results demonstrated that convective interfacial heat transfer is dependent on the wood species. The computed heat transfer coefficients were shown to vary between the pine and aspen models by nearly 20%. These differences are attributed to the species-specific variations in the exterior surface morphology of the biomass particles. We also quantify variations in heat transfer experienced by the particle when positioned in different orientations with respect to the direction of fluid flow. These results are compared to previously reported heat transfer coefficient correlations in the range of 0.1 < Pr < 1.5 and 10 < Re < 500. Comparison of these simulation results to correlations commonly used in the literature (Gunn, Ranz-Marshall, and Bird-Stewart-Lightfoot) shows that the Ranz-Marshall (sphere) correlation gave the closest h values to our steady-state simulations for both wood species, though no existing correlation was within 20% of both species at all conditions studied. In general, this work exemplifies the fact that all biomass feedstocks are not created equal, and that their species-specific characteristics must be appreciated in order to facilitate accurate simulations of conversion processes.

09 BIOMASS FUELS↗
DOE OSTI2016

New insights into the nonadiabatic state population dynamics of model proton-coupled electron transfer reactions from the mixed quantum-classical Liouville approach

In a previous study [F. A. Shakib and G. Hanna, J. Chem. Phys. 141, 044122 (2014)], we investigated a model proton-coupled electron transfer (PCET) reaction via the mixed quantum-classical Liouville (MQCL) approach and found that the trajectories spend the majority of their time on the mean of two coherently coupled adiabatic potential energy surfaces. This suggested a need for mean surface evolution to accurately simulate observables related to ultrafast PCET processes. In this study, we simulate the time-dependent populations of the three lowest adiabatic states in the ET-PT (i.e., electron transfer preceding proton transfer) version of the same PCET model via the MQCL approach and compare them to the exact quantum results and those obtained via the fewest switches surface hopping (FSSH) approach. We find that the MQCL population profiles are in good agreement with the exact quantum results and show a significant improvement over the FSSH results. All of the mean surfaces are shown to play a direct role in the dynamics of the state populations. Interestingly, our results indicate that the population transfer to the second-excited state can be mediated by dynamics on the mean of the ground and second-excited state surfaces, as part of a sequence of nonadiabatic transitions that bypasses the first-excited state surface altogether. This is made possible through nonadiabatic transitions between different mean surfaces, which is the manifestation of coherence transfer in MQCL dynamics. We also investigate the effect of the strength of the coupling between the proton/electron and the solvent coordinate on the state population dynamics. Drastic changes in the population dynamics are observed, which can be understood in terms of the changes in the potential energy surfaces and the nonadiabatic couplings. Finally, we investigate the state population dynamics in the PT-ET (i.e., proton transfer preceding electron transfer) and concerted versions of the model. The PT-ET results confirm the participation of all of the mean surfaces, albeit in different proportions compared to the ET-PT case, while the concerted results indicate that the mean of the ground- and first-excited state surfaces only plays a role, due to the large energy gaps between the ground- and second-excited state surfaces.

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

Shock-Driven Hydrodynamic Instability Growth Near Phase Boundaries and Material Property Transitions: Final Report

Activities for this grant included: 1) Development of dynamic impact experiments to probe strength and phase transition influence on dynamic deformation, 2) development of modern strength and phase aware simulation capabilities, 3) and post-processing of experimental data with simulation and closed form analytical techniques. Two different dynamic experiments were developed to probe material strengths in solid metals (largely copper and iron in this effort). In the first experiment a flyer plate impacts a flat target with an opposite rippled surface that is partially supported by a weaker window material. Post mortem analysis of the target sample showed a strong and repeatable residual plastic deformation dependence on grain orientation. Yield strengths for strain rates near 10 5 s -1 and plastic strains near ~50% were estimated to be around 180 to 240 MPa, varying in this range with grain orientation. Unfortunately dynamic real-time measurements were difficult with this setup due to diagnostic laser scattering; hence, an additional experimental setup was developed to complement these results. In the second set of experiments a rippled surface was ablated by a controlled laser pulsed, which launched a rippled shock front to an opposite initially flat diagnostic surface that was monitored in real-time with spatially resolved velocimetry techniques, e.g., line VISAR in addition to Transient Imaging Displacement Interferometry (TIDI) displacement measurements. This setup limited the displacements at the diagnostic surface to a reasonable level for TIDI measurements (~ less than one micrometer). These experiments coupled with analytical and numerical solutions provided evidence that viscous and elastic deviatoric strength affect shock front perturbation evolution in clearly different ways. Particularly, normalized shock front perturbation amplitudes evolve with viscosity (η) and perturbation wavelength (λ) as η/λ, such that increasing viscosity (or decreasing the initial wavelength) delays the perturbation decay. Conversely our experimental data, analysis and simulations show that for materials with elastic yield strength Y the normalized shock perturbation amplitude evolves with Yλ/A 0 , which shows wavelength increases have the opposite effect as in viscous materials and perturbation decay is also dependent on initial amplitude A 0 (viscous materials are independent of this parameter). Materials where strength had clear strain rate dependence, e.g., such as a PTW material law, behaved similarly to materials with only an effective yield stress (elastic-perfectly plastic) in the shock front perturbation studies obeying a Y eff λA 0 relationship where Y eff was a constant (near ~400 MPa for Cu for strain rates around 10 6 s -1 ). Magnitude changes in strain rate would increase Y eff as would be expected from the PTW behavior, but small perturbations (typical of regions behind the shock front) near a mean had little effect. Additional work based on simulations showed that phase transformation kinetics can affect the behavior of the perturbed shock front as well as the evolution of the RM-like instability that develops due to the imprint of the perturbed shock front on the initially flat surface as the shock breaks out.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Tracing the Rain Formation Pathways in Numerical Simulations of Deep Convection

Quantifying the microphysical process contributions to surface precipitation in numerical simulations can be challenging. This is due to the fact that many microphysical processes contribute to the formation and depletion of rain drops and there is almost always a spatial/temporal mismatch between where/when rain is formed and where/when it strikes the surface. In this work, we develop a tracing method that tracks the sources and sinks of raindrop mass and number as they are advected by the Weather Research and Forecasting model. Applying the method to an idealized squall line confirms that convective precipitation is dominated by warm rain processes (autoconversion and accretion) while stratiform precipitation is dominated by the melting of rimed and unrimed ice crystals. Sensitivity experiments in which the prescribed cloud drop number concentration is increased confirm the conventional wisdom that weakened autoconversion increases the fraction of raindrops originating from cold rain processes. The method also reveals that when applied to deep convection the Khairoutdinov and Kogan autoconversion scheme produces an excessive number of raindrops which are subsequently clipped in P3 microphysics to keep the rain size distribution within prescribed limits. This problem can mostly be mitigated by increasing the assumed radius for raindrops created by autoconversion.

58 GEOSCIENCES↗
DOE OSTI2021

Data from Technical and Economic Feasibility of an Integrated Ethanol and Anthocyanin Coproduction Process Using Purple Corn Stover

The coproduction of high-value anthocyanin extract in the cellulosic ethanol process would diversify the co-product market, increase revenue, and potentially improve the economics of the process. The high anthocyanin concentration in the cob and structural carbohydrates in residual stover make purple corn stover an attractive source for anthocyanin and ethanol coproduction. This study aimed to develop simulation models for processes integrating ethanol production and anthocyanin extraction using purple corn stover, to evaluate their techno-economic feasibility, and to compare their performance with the conventional ethanol production process using corn stover. The annual ethanol production for plants processing 2000 MT dry feedstock / day was 148.6 million L/year for the integrated processes compared with 222.6 million L/year for the conventional process. Anthocyanin production in the modified processes using dilute acid-based and water-based anthocyanin extraction processes was 1779 and 1099 MT/year, respectively. Capital investments for the integrated processes ($448.1 to $443.8 million) were higher than the conventional process ($371.9 million). Due to high revenue from anthocyanin extract, the ethanol production cost for the integrated process using acid-based anthocyanin extraction ($0.36/L) was 34.5% lower than conventional ethanol production ($0.55/L). The ethanol production cost for the integrated process using water-based anthocyanin extraction ($0.68/L) was higher than conventional ethanol production due to low ethanol and anthocyanin yields. The minimum ethanol selling price for the integrated process using acid-based anthocyanin extraction ($0.65/L) was also lower than the conventional process ($0.72/L), indicating an improvement in economic performance.

Conversion↗
DOE OSTI2016

Predicting Correlation Coefficients for Monte Carlo Eigenvalue Simulations

Monte Carlo methods are most often considered as a reference for neutron transport simulations since very limited approximations are made about nuclear data and system geometry. To report uncertainty of any tally evaluated as generation averages, the sample variance is divided by the number of active generations, which is based on the assumption that the neutron generations are independent. Correlation effects between neutrons in multiplying systems, particularly when performing power iteration to evaluate eigenvalues have been observed in previous work. Neglecting the correlation effect results in an underestimate of uncertainty reported by Monte Carlo calculations. Previous work has also proposed methods to predict the underestimation ratio. Yamamoto et al expanded the fission source distribution with diffusion equation modes, performed numerical simulation of the AR (autoregressive) process of the expansion coefficients and used the correlation of the AR process to predict that of the Monte Carlo eigenvalue simulation. Sutton applied the discretized phase space (DPS) approach to predict the underestimation ratio but the method cannot predict the ratio when one neutron generates offspring in different phase space regions or generates a random number of offspring. This paper presents a method to predict the correlation effect with the model of multi-type branching processes (MBP). The method requires simulations for one generation of neutrons without knowing the source distribution and can predict the underestimation ratio for the cases where the traditional DPS approach does not work. The generation-to-generation correlation determines the convergence rate of active generations, the bias of variance estimator for each generation and the underestimation ratio of variance estimator for tallies averaged over active generations. The generation-to-generation correlation is characterized by the Auto-Correlation Coefficients (ACC) between tallies from different generations. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Numerical simulation of the beam self-cleaning process in a multimode graded-index fibre during propagation of a pump wave and a Stokes component

A model of coupled modes is proposed, which makes it possible to describe the propagation of a pump wave and a Stokes component in a multimode graded-index fibre. A decrease in the depression in the centre of the profile of the pump wave intensity distribution as a result of the influence of Kerr effects and random linear coupling of modes is demonstrated. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Laser-induced ultrafast demagnetization time and spin moment in ferromagnets: First-principles calculation

When a laser pulse excites a ferromagnet, its spin undergoes a dramatic change. The initial demagnetization process is very fast. Experimentally, it is found that the demagnetization time is related to the spin moment in the sample. In this study, we employ the first-principles method to directly simulate such a process. We use the fixed spin moment method to change the spin moment in ferromagnetic nickel, and then we employ the Liouville equation to couple the laser pulse to the system. We find that in general the dependence of demagnetization time on the spin moment is nonlinear: It decreases with the spin moment up to a point, after which an increase with the spin moment is observed, followed by a second decrease. To understand this, we employ an extended Heisenberg model, which includes both the exchange interaction and spin-orbit coupling. The model directly links the demagnetization rate to the spin moment itself and demonstrates analytically that the spin relaxes more slowly with a small spin moment. A future experimental test of our predictions is needed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Data-driven prediction of geometry- and toolpath sequence-dependent intra-layer process conditions variations in laser powder bed fusion

Geometrical features and toolpath sequence are two important factors that cause process condition variations, such as variations in the meltpool temperature or meltpool size, that might lead to undesired material properties in the laser powder bed fusion (LPBF) process. Due to the high dynamics and complex physics of the LPBF process, it is difficult to predict variations in process conditions with simulations alone. Advances in measurement technology and computational technologies open up new possibilities for smart manufacturing. In this paper, a data-driven method to predict intra-layer variations in the processing conditions that source from the toolpath sequence and part geometry is presented. The approach is demonstrated using two-color on-axis pyrometer measurements. Three demonstration cases are presented in which it is demonstrated (1) how the trained predictive model can be used as a filter to ease the interpretation of process variations and discover patterns related to toolpath and part geometry, and (2) how to generate predictions that can be used for feedforward control, i.e., for adjusting laser power or scanning speed along the toolpath using a meltpool temperature prediction model generated based on on-axis measurements. Results show that the developed prediction model is able to meaningfully predict process variations resulted from toolpath sequence and geometry. Predictions are aligned with the results from the related work of others and for the case of 180° laser path turnarounds in our high-speed X-ray imaging experiments. In conclusion, the potential issues related to the current maturity status of the process and measuring equipment that could in practice affect the performance of the proposed solutions are also discussed.

42 ENGINEERING↗
DOE OSTI2022

Formation Process of Particles and Cloud Condensation Nuclei Over the Amazon Rainforest: The Role of Local and Remote New-Particle Formation

Understanding the formation processes of particles and cloud condensation nuclei (CCN) in pristine environments is a major challenge in accurately assessing anthropogenic radiative forcing. Here, using a state-of-the-art model that systematically simulates the chemical processes controlling new-particle formation (NPF) and multi-scale chemical transport, we find that NPF contributes ~90% of the particle number and ~80% of the CCN at 0.5% supersaturation (CCN0.5%) in the pristine Amazon boundary layer during the wet season. The corresponding contributions are only ~30% and ~20% during the dry season because of prevalent biomass burning. In both seasons, about half of the NPF-induced particles and ~85% of the NPF-induced CCN0.5% in the boundary layer originate from the long-range transport of new particles formed hundreds to thousands of kilometers away. Moreover, about 50–65% of the NPF-induced particles and 35–50% of the NPF induced CCN0.5% originate from the downward transport of new particles formed in the free troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Global variability in atmospheric new particle formation mechanisms

A key challenge in aerosol pollution studies and climate change assessment is to understand how atmospheric aerosol particles are initially formed. Although new particle formation (NPF) mechanisms have been described at specific sites, in most regions, such mechanisms remain uncertain to a large extent because of the limited ability of atmospheric models to simulate critical NPF processes. Here we synthesize molecular-level experiments to develop comprehensive representations of 11 NPF mechanisms and the complex chemical transformation of precursor gases in a fully coupled global climate model. Combined simulations and observations show that the dominant NPF mechanisms are distinct worldwide and vary with region and altitude. Previously neglected or underrepresented mechanisms involving organics, amines, iodine oxoacids and HNO 3 probably dominate NPF in most regions with high concentrations of aerosols or large aerosol radiative forcing; such regions include oceanic and human-polluted continental boundary layers, as well as the upper troposphere over rainforests and Asian monsoon regions. These underrepresented mechanisms also play notable roles in other areas, such as the upper troposphere of the Pacific and Atlantic oceans. Accordingly, NPF accounts for different fractions (10–80%) of the nuclei on which cloud forms at 0.5% supersaturation over various regions in the lower troposphere. The comprehensive simulation of global NPF mechanisms can help improve estimation and source attribution of the climate effects of aerosols.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of the dynamics of gas mixtures during plasma processing in the C75 Cavity

Plasma processing using a mixture of noble gas and oxygen is a technique that is currently being used to reduce field emission and multipacting in accelerating cavities. Plasma is created inside the cavity when the gas mixture is exposed to an electromagnetic field that is generated by applying RF power through the fundamental power or higher-order mode couplers. Oxygen ions and atomic oxygen are created in the plasma which breaks down the hydrocarbons on the surface of the cavity and the residuals from this process are removed as part of the process gas flow. Removal of hydrocarbons from the surface increases the work function and reduces the secondary emission coefficient. This work describes the initial results of plasma simulation, which provides insight into the ignition process, distribution of different species, and interactions of free oxygen and oxygen ions with the cavity surfaces. The simulations have been done with an Ar/?2 plasma using COMSOL® multiphysics. These simulations help in understanding the dynamics and control of plasma inside the cavity and the exploration of different gas mixtures.

Dhakal, Pashupati↗
DOE OSTI2025

A Novel Manufacturing Process of Lightweight Automotive Seats (Integration of Additive Manufacturing and Reinforced Polymer Composite)

Lightweight automotive seats offer multiple benefits to original equipment manufacturers in terms of cost savings from various aspects, including less material usage, more integrated processes, and compliance with Corporate Average Fuel Economy Standards. Original equipment manufacturers have been focusing on innovative ways to produce light weight automotive seats. The commercially available automotive seats are currently made of multiple metal components combined through welding and fasteners. The use of additive manufacturing and composite structures is particularly useful for light weighting the automotive components. Additive manufacturing (AM) offers multiple advantages over traditional manufacturing processes such as freedom of design thereby enabling complex structural geometries, mass customization and waste minimization, and control over the fiber alignment through deposition in a predetermined pattern. Combining metal inserts with polymer composites through a novel manufacturing process allows design of lightweight and high-performance materials for automotive components. However, fabricating these metal polymer composite structures through traditional manufacturing processes limits their mechanical properties due to limited design freedom, lack of control over fiber orientation in composite parts, and poor interfacial bonding between the constituent materials. It is essential to develop a novel manufacturing process to enable high throughput production of lightweight automotive seats using metal and polymer composites. As such it is important to design the automotive seat suitable for manufacturing via this process and perform mechanical characterization on various subcomponents of the seat to ensure that the design and performance requirements provided by the auto manufacturer are met. The aim of this project is to develop a novel manufacturing technique to produce lightweight automotive seat by combining AM with conventional manufacturing processes. The car seat back panel will be designed via topology optimization and numerical simulations to minimize the overall weight while ensuring it meets all the performance requirements. The optimization of the seat back structure will be based on computational stress analysis to maximize the stiffness and minimize the weight. Materials currently used by Ford Motor Company will be adopted for a few subcomponents while the in-house composite materials will be used for the rest of the seat back. The composite and metallic materials will be tested to determine their mechanical properties as these are necessary for simulations. A novel manufacturing process will be developed to integrate AM metal inserts with discontinuous reinforced composite through large scale additive manufacturing and compression overmolding processes. The developed manufacturing technique will be used to fabricated various subcomponents suitable for the seat back design and mechanically tested to determine their properties. The manufacturing of the lightweight seat back design through this process involves integrated AM metal inserts with the composite structure for recliner connection. The manufacturing of the entire seat back which is lightweight through the novel manufacturing process will be discussed. The performance of the designed seat back will be investigated through numerical simulations and shown to meet all the requirements provided by the auto manufacturer. The final goal of developing a novel manufacturing process for lightweight automotive seats is met through design optimization of seat back, manufacturing of subcomponents, mechanical characterization, and validation through numerical simulations. The routes to achieve the final goal of the project and the depth in which they were investigated changed throughout the project due to personnel changes and the COVID-19 pandemic. The project resulted in the development of a novel manufacturing process to integrate metal inserts with tailored polymer composite preforms through overmolding. Leveraging this proven manufacturing process, a lightweight seat back was designed through topology optimization and numerical simulations. The designed seat back uses AM metal inserts and compression overmolding of tailored polymer composite preforms obtained via large scale additive manufacturing. The metal polymer composite structures fabricated through this process exhibited enhancement in stiffness and improved ductility upon testing. Overall, the project provided an alternative design and manufacturing technique for automotive seat back that enables weight saving while meeting the safety and performance requirements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Technical and economic feasibility of an integrated ethanol and anthocyanin coproduction process using purple corn stover

The coproduction of high-value anthocyanin extract in the cellulosic ethanol process would diversify the co-product market, increase revenue, and potentially improve the economics of the process. The high anthocyanin concentration in the cob and structural carbohydrates in residual stover make purple corn stover an attractive source for anthocyanin and ethanol coproduction. In this study, we aimed to develop simulation models for processes integrating ethanol production and anthocyanin extraction using purple corn stover, to evaluate their techno-economic feasibility, and to compare their performance with the conventional ethanol production process using corn stover. The annual ethanol production for plants processing 2000 MT dry feedstock / day was 148.6 million L/year for the integrated processes compared with 222.6 million L/year for the conventional process. Anthocyanin production in the modified processes using dilute acid-based and water-based anthocyanin extraction processes was 1779 and 1099 MT/year, respectively. Capital investments for the integrated processes ($\$448.1$ to $\$443.8$ million) were higher than the conventional process ($\$371.9$ million). Due to high revenue from anthocyanin extract, the ethanol production cost for the integrated process using acid-based anthocyanin extraction ($\$0.36$/L) was 34.5% lower than conventional ethanol production ($\$0.55$/L). The ethanol production cost for the integrated process using water-based anthocyanin extraction ($\$0.68$/L) was higher than conventional ethanol production due to low ethanol and anthocyanin yields. The minimum ethanol selling price for the integrated process using acid-based anthocyanin extraction ($\$0.65$/L) was also lower than the conventional process ($\$0.72$/L), indicating an improvement in economic performance.

09 BIOMASS FUELS↗
DOE OSTI2021

Pretreatment with Sodium Methyl Mercaptide Increases Carbohydrate Yield during Kraft Pulping

Kraft pulping is the predominant technology in the pulp and paper industry for removing lignin from wood carbohydrates to produce paper, board, packaging, tissue, and specialty cellulose. However, the kraft process is energy intensive and expensive, and its yield is limited by the degradation of carbohydrates. Pretreatment can increase carbohydrate yield by limiting degradation via primary peeling of reducing end groups. However, protection of galactoglucomannan (GGM), the primary hemicellulose component of softwood, is minimal when conventional pretreatments are used. Here, we investigated the effectiveness of sodium methyl mercaptide pretreatment on southern pine wood chips under a range of experimental conditions. We found that pretreatment of biomass with 4.38% sodium methyl mercaptide at pH 12 and 105 °C for 60 min provided small but significant increases in xylan and cellulose yields relative to control conditions, but preservation of GGM was minimal. Furthermore, to provide insight into molecular-scale details of primary peeling, pretreatment, and alkaline hydrolysis, we performed classical molecular dynamics (MD) simulations under selected process conditions and quantum mechanical (QM) calculations of selected reactions. MD simulations showed that C1 of the GGM reducing end is more readily accessible by HO – and CH 3 S – ions than in cellulose. The free energy barrier for peeling calculated with QM is lower for GGM than for cellulose, indicating increased susceptibility to peeling. In addition, we found that GGM may be more susceptible to internal chain cleavage than cellulose. Thus, even though reducing end groups may be protected initially through pretreatment, new unprotected reducing end groups may be generated through alkaline hydrolysis. Taken together, these findings show the promise of methyl mercaptide as a pretreatment technology for cellulose retention and also provide molecular insight for improving its effectiveness toward GGM.

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

Effects of strain-dependent surface stress on the adhesive contact of a rigid sphere to a compliant substrate

Recent experiments have reported that the surface stress of soft elastic solids can increase rapidly with surface strain. For example, when a small hard sphere in adhesive contact with a soft silicone gel is slowly retracted from its rest position, it was found that the retraction force versus displacement relation cannot be explained either by the Johnson–Kendall–Roberts (JKR) theory or a recent indentation theory based on an isotropic surface stress that is independent of surface strain. In this paper, we address this problem using a finite element method to simulate the retraction process. Our numerical model does not have the restrictions of the aforementioned theories; that is, it can handle large nonlinear elastic deformation as well as a surface-strain-dependent surface stress. Our simulation is in good agreement with experimental force versus displacement data with no fitting parameters. Therefore, our results lend further support to the claim that significant strain-dependent surface stresses can occur in simple soft elastic gels. However, significant challenges remain in the reconciliation of theory and experiments, particularly regarding the geometry of the contact and substrate deformation.

36 MATERIALS SCIENCE↗