Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “optimization analysis”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Numerical Simulation of Ion Transport in a Nano-Electrospray Ion Source at Atmospheric Pressure

Understanding ion transport properties from the ion source to the mass spectrometer (MS) is essential for optimizing device performance. Numerical simulation helps in understanding of ion transport properties and, furthermore, facilitates instrument design. In contrast to previously reported numerical studies, ion transport simulations in a continuous injection mode whilst considering realistic space-charge effects have been carried out. The flow field was solved using Reynolds-averaged Navier-Stokes (RANS) equations, and a particle-in-cell (PIC) method was applied to solve a time-dependent electric field with local charge density. A series of ion transport simulations were carried out at different cone gas flow rates, ion source currents, and capillary voltages. A force evaluation analysis reveals that the electric force, the drag force, and the Brownian force are the three dominant forces acting on the ions. Both the experimental and simulation results indicate that cone gas flow rates of ≤250 slph (standard liter per hour) are important for high ion transmission efficiency, as higher cone gas flow rates reduce the ion signal significantly. The simulation results also show that the ion transmission efficiency reduces exponentially with an increased ion source current. Additionally, the ion loss due to space-charge effects has been found to be predominant at a higher ion source current, a lower capillary voltage, and a stronger cone gas counterflow. The interaction of the ion driving force, ion opposing force, and ion dispersion is discussed to illustrate ion transport mechanism in the ion source at atmospheric pressure. .

42 ENGINEERING↗
DOE OSTI2013

Optimization of beam configuration in laser fusion based on the laser beam pattern

A simple method based on the laser beam pattern is proposed and numerically demonstrated to optimize a beam configuration for direct drive laser fusion. In this method, both the geometrical factor G{sub l} and the single beam factor B{sub l} are considered. By diminishing the product of B{sub l}·G{sub l}, the irradiation nonuniformity can be decreased to the order of 10{sup −5}. This optimization method can be applied on the design of irradiation systems for an arbitrary number of beams and any axially symmetric beam patterns.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Crystallographic Insight of Reduced Lattice Volume Expansion in Mesoporous Cu 2+ -Doped TiNb 2 O 7 Microspheres during Li + Insertion

TiNb 2 O 7 represents a promising anode material for lithium-ion batteries (LIBs), but its practical applications are currently hampered by the non-negligible volumetric expansion and contraction during the charge/discharge process and the sluggish ion/electron kinetics. Here a combination technique is reported by systematically optimizing the porous and spherical morphology, crystal structure, and surface decoration of mesoporous Cu 2+ -doped TiNb 2 O 7 microspheres to enhance the electrochemical Li + storage performance and stability simultaneously. The Cu 2+ dopants preferentially replace Ti 4+ in crystal lattices, which decreases the Li + diffusion barrier and increases the electronic conductivity, as confirmed by density functional theory (DFT) calculation and demonstrated by diverse electrochemical characterizations. The successful Cu 2+ doping significantly reduces the lattice expansion coefficient from 7.26% to 4.61% after Li + insertion along the b-axis of TiNb 2 O 7 , as visualized from in situ and ex situ XRD analysis. The optimal 5% Cu 2+ -doped TiNb 2 O 7 with surface coating of N-doped carbon exhibits significantly enhanced specific capacity and rate and cyclic performances in both half- and full-cell configurations, demonstrating an excellent electrochemical behavior for fast-charging LIB applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Designing Monte Carlo Simulation and an Optimal Machine Learning to Optimize and Model Space Missions

This paper investigates applying artificial intelligence (AI) algorithms to attitude control system of satellites to optimally tune the controller using high performance computing. This methodology is applied to the Virtual Telescope for X-ray Observation mission, which is a precise formation of two separate spacecraft observing multiple objects in the space in the X-ray domain. The mission is divided into phases based on the instrumentation and the mission goal. To reach an stable precise formation robust to stochastic slew and slew rate (i.e., Euler angles and angular velocities) in a minimal constrained time T , consumed energy of the attitude control system, denoted as E, and root-mean-square state error of attitude control system, denoted as e, are minimized. Monte-Carlo simulation is used for the sensitivity analysis of optimization and designing a controller. Deep neural networks (DNN), Gaussian processes (GP), and support vector regression (SVR) learn this optimization as a surrogate model, while their hyperparameters are optimized in a novel approach. THETA supercomputer at Argonne Leadership Computing Facility (ALCF) is used for optimizing the hyperparameters of DNN. The surrogate model meets the requirements of the mission, and it shows a better performance over the optimization and Monte-Carlo. The optimal DNN can satisfy the mission requirements e and T while reducing E for 90% compared to the other given methods.

42 ENGINEERING↗
DOE OSTI2019

Phosphate glass matrix composites incorporating trisilanol phenyl polyhedral oligomeric silsesquioxane prepared by viscous flow sintering method with enhanced benefits

The effect of mixing and sintering processes to prepare tin fluorophosphate glass (Pglass) matrix composites incorporating trisilanol phenyl polyhedral oligomeric silsesquioxane (TSP-POSS) was investigated by comparing manual and suspension mixing, one-step and stepwise sintering processes to explore the structure dynamics and physical properties in the composites as a function of the different process conditions used. Energy Dispersive X-ray analysis confirmed optimal homogeneous dispersion of the TSP-POSS molecules in the composites prepared by the suspension method. The observed increase of glass transition temperature and the reduction of non-bridging bonds in the composites are believed to be the reason for the effective dispersion of TSP-POSS molecules in the composites. The chemical reaction between the TSP-POSS and Pglass was strongly influenced by the mixing/dispersion and sintering processes investigated. 13 C cross polarized magic angle spinning (CP MAS) solid state nuclear magnetic resonance (NMR) spectroscopy confirmed the chemical stability of the TSP-POSS during the sintering process at elevated temperatures. In addition, a chemical reaction between the TSP-POSS and Pglass was evidenced by 29 Si CP MAS NMR analysis. Here, this study will provide a better fundamental understanding of the effective dispersion mechanism of the TSP-POSS molecules in the Pglass matrix that will facilitate tailoring the physicochemical properties of the composites with addition of various small concentrations of TSP-POSS for a number of applications where the pure Pglass is not applicable due to the intrinsic properties of the Pglass.

36 MATERIALS SCIENCE↗
DOE OSTI2022

System Modeling Frameworks for Wind Turbines and Plants: Review and Requirements Specifications

System modeling frameworks for wind turbines and plants are used by research groups and industry to design wind energy systems that take into account key trade-offs across performance, cost, and reliability at both the turbine and plant level. The frameworks are exercised using a variety of multi-disciplinary design, analysis and optimization (MDAO) methods. To improve inter-operability and foster collaboration, this report proposes a classification system for the frameworks along dimensions of model fidelity and scope. The classification system is first motivated with reviews the state-of-the-art in the development of software frameworks for integrated wind turbine and plant simulation. Within each major wind turbine and power plant subsystem, a matrix is developed for the disciplines used and the fidelity levels with which each discipline can be modeled. The existing frameworks are then classified according to the matrix. Next, an ontology is proposed that will allow for standardizing how data is transferred between the most common discipline-fidelity combinations used in the frameworks. A common representation of data creates the ability to 1) share system descriptions and analysis results, supporting more transparent benchmarks and comparison, and 2) integrate models together into workflows within and across organizations for improving the efficiency and performance of wind turbine and power plant design processes. Ultimately, this integration leads to better overall wind energy system designs with high performance and low costs.

17 WIND ENERGY↗
DOE OSTI2019

Effect of Yagi–Uda nano-antenna element shape on the directivity and radiation efficiency

In this paper, different optimal designs of Yagi–Uda nano-antennas (NAs) are introduced and analyzed based on different element shapes to maximize the directivity and radiation efficiency. The studied element shapes include triangular, hexagonal, cube, square, rotated square, rectangular, ellipsoid, and elliptical cylinder shapes. The NAs parameters for each design are optimized using particle swarm optimization algorithm to achieve high directivity at a wavelength of 500 nm for wireless point-to-point applications. The designs are numerically analysed using 3-D finite difference time domain method (3-D FDTD) via Lumerical software package. To illustrate the performance of the proposed antennas, different radiation parameters such as radiation pattern, directivity, and radiation efficiency have been studied. Moreover, a tolerance of ± 5% is calculated to consider the imperfections that may appear in the fabrication process at a nanoscale. The optimized rectangular prism NA achieves the highest directivity of 24.21 with an efficiency of 50.74%. However, the optimized triangular prism NA has the highest efficiency of 73% with a directivity of 20.89. This enhancement is mainly attributed to the efficient coupling between array elements along with low side lobe level. However, a trade-off between directivity and radiation efficiency is obtained by the hexagonal prism NA with directivity of 21.9 and radiation efficiency of 69.1%.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Optimizing fused silica debris shield use in the National Ignition Facility

The National Ignition Facility (NIF) deliberately and routinely operates with its final optics exposed to fluences likely to induce and grow damage. This choice has enabled the NIF to operate at energies previously unachievable at the expense of limitations to delivered power and energy contingent on the ability to repair damage on the final optics. Prior to the introduction of the fused silica debris shield (FSDS), the grating debris shield (GDS) was the optic type on the NIF most prone to damage and hence a primary limiter to the delivered energy. The introduction of the FSDS has reduced the damage initiation rate on the GDS by about two orders of magnitude to date. Despite the FSDS being a disposable optic, it is important to understand how the FSDS damages and when its damage requires it to be exchanged to maximize the FSDS lifetime while still protecting the GDS from damage. We will define the FSDS exchange criteria and evaluate the damage initiation mechanisms of the FSDS. The surface damage morphologies provide insight into the impact of the various damage mechanisms that limit the FSDS lifetime. These results are utilized in our analysis to optimize the FSDS exchange criteria to further extend the GDS lifetime.

42 ENGINEERING↗
DOE OSTI2017

Topology-Aware Performance Optimization and Modeling of Adaptive Mesh Refinement Codes for Exascale

Here, we introduce a topology-aware performance optimization and modeling workflow for AMR simulation that includes two new modeling tools, ProgrAMR and Mota Mapper, which interface with the BoxLib AMR framework and the SSTmacro network simulator. ProgrAMR allows us to generate and model the execution of task dependency graphs from high-level specifications of AMR-based applications, which we demonstrate by analyzing two example AMR-based multigrid solvers with varying degrees of asynchrony. Mota Mapper generates multiobjective, network topology-aware box mappings, which we apply to optimize the data layout for the example multigrid solvers. While the sensitivity of these solvers to layout and execution strategy appears to be modest for balanced scenarios, the impact of better mapping algorithms can be significant when performance is highly constrained by network hop latency. Furthermore, we show that network latency in the multigrid bottom solve is the main contributing factor preventing good scaling on exascale-class machines.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Low-order preconditioning of the Stokes equations

A well-known strategy for building effective preconditioners for higher-order discretizations of some PDEs, such as Poisson's equation, is to leverage effective preconditioners for their low-order analogs. In this work, we show that high-quality preconditioners can also be derived for the Taylor–Hood discretization of the Stokes equations in much the same manner. In particular, we investigate the use of geometric multigrid based on the Q 1 iso Q 2 /Q 1 discretization of the Stokes operator as a preconditioner for the Q 2 /Q 1 discretization of the Stokes system. We utilize local Fourier analysis to optimize the damping parameters for Vanka and Braess–Sarazin relaxation schemes and to achieve robust convergence. Furthermore, these results are then verified and compared against the measured multigrid performance. While geometric multigrid can be applied directly to the Q 2 /Q 1 system, our ultimate motivation is to apply algebraic multigrid within solvers for Q 2 /Q 1 systems via the Q 1 iso Q 2 /Q 1 discretization, which will be considered in a companion paper.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Design Optimization Study on the Single Tank Packed-bed Thermal Energy Storage System

Idaho National Laboratory (INL) is establishing the Dynamic Energy Transport and Integration Lab (DETAIL) as part of its commitment to research on nuclear-renewable hybrid energy systems and associated advanced reactor technologies. DETAIL is designed to allow several different energy systems to work in unison. Current plans include a PWR simulator, high-temperature steam electrolysis (HSTE) unit, and a renewable energy system (e.g., photovoltaics). DETAIL will provide the real world basis for studies on the systems integration and system configurations to be completed. Encompassed in the DETAIL program is the Thermal Energy Storage (TES) system. The TES system is one of the key components in DETAIL because it allows for the delayed release of energy and can be used to simulate storage capabilities currently being considered by utility providers. Of the various TES concepts, this paper deals with the single-tank packed bed TES system. Particular attention is given to the transient thermal behavior of fluid and solid particles within the packed bed thermocline tank and heat storage efficiency influenced by various design parameters. The effects of tank geometry (height-to-diameter ratio), filler size, filler packing ratio and operating temperature differentials are investigated. Based on the parametric study and cost analysis, the optimal TES tank design for DETAIL is discussed.

25 ENERGY STORAGE↗
DOE OSTI2015

Low cost back contact heterojunction solar cells on thin c-Si wafers. Integrating laser and thin film processing for improved manufacturability

An interdigitated back contact (IBC) Si wafer solar cell with deposited a-Si heterojunction (HJ) emitter and contacts is considered the ultimate single junction Si solar cell design. This was confirmed in 2014 by both Panasonic and Sharp Solar producing IBC-HJ cells breaking the previous record Si solar cell efficiency of 25%. But manufacturability at low cost is a concern for the complex IBC-HJ device structure. In this research program, our goals were to addressed the broad industry need for a high-efficiency c-Si cell that overcomes the dominant module cost barriers by 1) developing thin Si wafers synthesized by innovative, kerfless techniques; 2) integrating laser-based processing into most aspects of solar cell fabrication, ensuring high speed and low thermal budgets ; 3) developing an all back contact cell structure compatible with thin wafers using a simplified, low-temperature fabrication process; and 4) designing the contact patterning to enable simplified module assembly. There were a number of significant achievements from this 3 year program. Regarding the front surface, we developed and applied new method to characterize critical interface recombination parameters including interface defect density Dit and hole and electron capture cross-section for use as input for 2D simulation of the IBC cell to guide design and loss analysis. We optimized the antireflection and passivation properties of the front surface texture and a-Si/a-SiN/a-SiC stack depositions to obtain a very low (< 6 mA/cm2) front surface optical losses (reflection and absorption) while maintaining excellent surface passivation (SRV<5 cm/s). We worked with kerfless wafer manufacturers to apply defect-engineering techniques to improve bulk minority-carrier lifetime of thin kerfless wafers by both reducing initial impurities during growth and developing post-growth gettering techniques. This led insights about the kinetics of nickel, chromium, and dislocations in PV-grade silicon and to achieving millisecond lifetimes in kerfless silicon materials. Laser fired contacts to n-Si were developed for the first time using a Al/Sb/Ti metal stack giving contact resistances < 5 mOhm-cm2 when fired through several different dielectric layers. A new 2 step laser+chemical etch isolation technique was developed using a sacrificial top coating which avoids laser damage to Si passivation. Regarding the heterojunction emitter, analysis of front FHJ (1D) and IBC (2D) cells with range of p-layer conditions found that a 2-stage high/low doped p-layer was optimum: the low doped region has lower defects giving higher Voc and the high doped region gave a better contact to the metal. A significant effort was spent studying the patterning process and its contribution to degradation of passivation and reproducibility. Several promising new cleaning, contact and deposition patterning and processing approaches were implemented leading to fabrication of several runs with cells having 19-20% efficiency which were stable over several months. This program resulted in the training and support of 12 graduate students, publication of 21 journal papers and 14 conference papers.

14 SOLAR ENERGY↗
DOE OSTI2015

Low cost back contact heterojunction solar cells on thin c-Si wafers. integrating laser and thin film processing for improved manufacturability

An interdigitated back contact (IBC) Si wafer solar cell with deposited a-Si heterojunction (HJ) emitter and contacts is considered the ultimate single junction Si solar cell design. This was confirmed in 2014 by both Panasonic and Sharp Solar producing IBC-HJ cells breaking the previous record Si solar cell efficiency of 25%. But manufacturability at low cost is a concern for the complex IBC-HJ device structure. In this research program, our goals were to addressed the broad industry need for a high-efficiency c-Si cell that overcomes the dominant module cost barriers by 1) developing thin Si wafers synthesized by innovative, kerfless techniques; 2) integrating laser-based processing into most aspects of solar cell fabrication, ensuring high speed and low thermal budgets ; 3) developing an all back contact cell structure compatible with thin wafers using a simplified, low-temperature fabrication process; and 4) designing the contact patterning to enable simplified module assembly. There were a number of significant achievements from this 3 year program. Regarding the front surface, we developed and applied new method to characterize critical interface recombination parameters including interface defect density Dit and hole and electron capture cross-section for use as input for 2D simulation of the IBC cell to guide design and loss analysis. We optimized the antireflection and passivation properties of the front surface texture and a-Si/a-SiN/a-SiC stack depositions to obtain a very low (< 6 mA/cm2) front surface optical losses (reflection and absorption) while maintaining excellent surface passivation (SRV<5 cm/s). We worked with kerfless wafer manufacturers to apply defect-engineering techniques to improve bulk minority-carrier lifetime of thin kerfless wafers by both reducing initial impurities during growth and developing post-growth gettering techniques. This led insights about the kinetics of nickel, chromium, and dislocations in PV-grade silicon and to achieving millisecond lifetimes in kerfless silicon materials. Laser fired contacts to n-Si were developed for the first time using a Al/Sb/Ti metal stack giving contact resistances < 5 mOhm-cm2 when fired through several different dielectric layers. A new 2 step laser+chemical etch isolation technique was developed using a sacrificial top coating which avoids laser damage to Si passivation. Regarding the heterojunction emitter, analysis of front FHJ (1D) and IBC (2D) cells with range of p-layer conditions found that a 2-stage high/low doped p-layer was optimum: the low doped region has lower defects giving higher Voc and the high doped region gave a better contact to the metal. A significant effort was spent studying the patterning process and its contribution to degradation of passivation and reproducibility. Several promising new cleaning, contact and deposition patterning and processing approaches were implemented leading to fabrication of several runs with cells having 19-20% efficiency which were stable over several months. This program resulted in the training and support of 12 graduate students, publication of 21 journal papers and 14 conference papers.

14 SOLAR ENERGY↗
DOE OSTI2014

Geometry optimization for micro-pressure sensor considering dynamic interference

Presented is the geometry optimization for piezoresistive absolute micro-pressure sensor. A figure of merit called the performance factor (PF) is defined as a quantitative index to describe the comprehensive performances of a sensor including sensitivity, resonant frequency, and acceleration interference. Three geometries are proposed through introducing islands and sensitive beams into typical flat diaphragm. The stress distributions of sensitive elements are analyzed by finite element method. Multivariate fittings based on ANSYS simulation results are performed to establish the equations about surface stress, deflection, and resonant frequency. Optimization by MATLAB is carried out to determine the dimensions of the geometries. Convex corner undercutting is evaluated. Each PF of the three geometries with the determined dimensions is calculated and compared. Silicon bulk micromachining is utilized to fabricate the prototypes of the sensors. The outputs of the sensors under both static and dynamic conditions are tested. Experimental results demonstrate the rationality of the defined performance factor and reveal that the geometry with quad islands presents the highest PF of 210.947 Hz{sup 1/4}. The favorable overall performances enable the sensor more suitable for altimetry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

FY17 Status Report on the Initial Development of a Constitutive Model for Grade 91 Steel

Grade 91 is a candidate structural material for high temperature advanced reactor applications. Existing ASME Section III, Subsection HB, Subpart B simplified design rules based on elastic analysis are setup as conservative screening tools with the intent to supplement these screening rules with full inelastic analysis when required. The Code provides general guidelines for suitable inelastic models but does not provide constitutive model implementations. This report describes the development of an inelastic constitutive model for Gr. 91 steel aimed at fulfilling the ASME Code requirements and being included into a new Section III Code appendix, HBB-Z. A large database of over 300 experiments on Gr. 91 was collected and converted to a standard XML form. Five families of Gr. 91 material models were identified in the literature. Of these five, two are potentially suitable for use in the ASME code. These two models were implemented and evaluated against the experimental database. Both models have deficiencies so the report develops a framework for developing and calibrating an improved model. This required creating a new modeling method for representing changes in material rate sensitivity across the full ASME allowable temperature range for Gr. 91 structural components: room temperature to 650° C. On top of this framework for rate sensitivity the report describes calibrating a model for work hardening and softening in the material using genetic algorithm optimization. Future work will focus on improving this trial model by including tension/compression asymmetry observed in experiments and necessary to capture material ratcheting under zero mean stress and by improving the optimization and analysis framework.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Stochastic approximation methods-Powerful tools for simulation and optimization: A survey of some recent work on multi-agent systems and cyber-physical systems

Stochastic approximation methods have found extensive and diversified applications. Recent emergence of networked systems and cyber-physical systems has generated renewed interest in advancing stochastic approximation into a general framework to support algorithm development for information processing and decisions in such systems. This paper presents a survey on some recent developments in stochastic approximation methods and their applications. Using connected vehicles in platoon formation and coordination as a platform, we highlight some traditional and new methodologies of stochastic approximation algorithms and explain how they can be used to capture essential features in networked systems. Distinct features of networked systems with randomly switching topologies, dynamically evolving parameters, and unknown delays are presented, and control strategies are provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Dynamic Event Tree advancements and control logic improvements

The RAVEN code has been under development at the Idaho National Laboratory since 2012. Its main goal is to create a multi-purpose platform for the deploying of all the capabilities needed for Probabilistic Risk Assessment, uncertainty quantification, data mining analysis and optimization studies. RAVEN is currently equipped with three different sampling categories: Forward samplers (Monte Carlo, Latin Hyper Cube, Stratified, Grid Sampler, Factorials, etc.), Adaptive Samplers (Limit Surface search, Adaptive Polynomial Chaos, etc.) and Dynamic Event Tree (DET) samplers (Deterministic and Adaptive Dynamic Event Trees). The main subject of this document is to report the activities that have been done in order to: start the migration of the RAVEN/RELAP-7 control logic system into MOOSE, and develop advanced dynamic sampling capabilities based on the Dynamic Event Tree approach. In order to provide to all MOOSE-based applications a control logic capability, in this Fiscal Year an initial migration activity has been initiated, moving the control logic system, designed for RELAP-7 by the RAVEN team, into the MOOSE framework. In this document, a brief explanation of what has been done is going to be reported. The second and most important subject of this report is about the development of a Dynamic Event Tree (DET) sampler named “Hybrid Dynamic Event Tree” (HDET) and its Adaptive variant “Adaptive Hybrid Dynamic Event Tree” (AHDET). As other authors have already reported, among the different types of uncertainties, it is possible to discern two principle types: aleatory and epistemic uncertainties. The classical Dynamic Event Tree is in charge of treating the first class (aleatory) uncertainties; the dependence of the probabilistic risk assessment and analysis on the epistemic uncertainties are treated by an initial Monte Carlo sampling (MCDET). From each Monte Carlo sample, a DET analysis is run (in total, N trees). The Monte Carlo employs a pre-sampling of the input space characterized by epistemic uncertainties. The consequent Dynamic Event Tree performs the exploration of the aleatory space. In the RAVEN code, a more general approach has been developed, not limiting the exploration of the epistemic space through a Monte Carlo method but using all the forward sampling strategies RAVEN currently employs. The user can combine a Latin Hyper Cube, Grid, Stratified and Monte Carlo sampling in order to explore the epistemic space, without any limitation. From this pre-sampling, the Dynamic Event Tree sampler starts its aleatory space exploration. As reported by the authors, the Dynamic Event Tree is a good fit to develop a goal-oriented sampling strategy. The DET is used to drive a Limit Surface search. The methodology that has been developed by the authors last year, performs a Limit Surface search in the aleatory space only. This report documents how this approach has been extended in order to consider the epistemic space interacting with the Hybrid Dynamic Event Tree methodology.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Hydrogen adsorption and desorption with 3D silicon nanotube-network and film-network structures: Monte Carlo simulations

Hydrogen is clean, sustainable, and renewable, thus is viewed as promising energy carrier. However, its industrial utilization is greatly hampered by the lack of effective hydrogen storage and release method. Carbon nanotubes (CNTs) were viewed as one of the potential hydrogen containers, but it has been proved that pure CNTs cannot attain the desired target capacity of hydrogen storage. In this paper, we present a numerical study on the material-driven and structure-driven hydrogen adsorption of 3D silicon networks and propose a deformation-driven hydrogen desorption approach based on molecular simulations. Two types of 3D nanostructures, silicon nanotube-network (Si-NN) and silicon film-network (Si-FN), are first investigated in terms of hydrogen adsorption and desorption capacity with grand canonical Monte Carlo simulations. It is revealed that the hydrogen storage capacity is determined by the lithium doping ratio and geometrical parameters, and the maximum hydrogen uptake can be achieved by a 3D nanostructure with optimal configuration and doping ratio obtained through design optimization technique. For hydrogen desorption, a mechanical-deformation-driven-hydrogen-release approach is proposed. Compared with temperature/pressure change-induced hydrogen desorption method, the proposed approach is so effective that nearly complete hydrogen desorption can be achieved by Si-FN nanostructures under sufficient compression but without structural failure observed. The approach is also reversible since the mechanical deformation in Si-FN nanostructures can be elastically recovered, which suggests a good reusability. This study may shed light on the mechanism of hydrogen adsorption and desorption and thus provide useful guidance toward engineering design of microstructural hydrogen (or other gas) adsorption materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗