Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation technique”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of biomass liquefaction on glucose and xylose prices predicted by National Renewable Energy Laboratory biochemical sugar model

Abstract The National Renewable Energy Laboratory (NREL) published a model in 2017 that enables the minimum selling price of lignocellulosic sugar to be calculated. The model can be modified to suit any biomass feedstock and operational design. In the present case, the model is used to understand the economics of a process configuration that incorporates liquefaction as a preprocessing step in corn‐stover‐fed biorefinery. This study demonstrates a quantitative approach utilizing an existing biorefinery setup to simulate the biomass liquefaction technique while estimating the price of the resulting sugar product. The objective is to understand whether the addition of liquefaction methodology – and the substitution of acid pretreatment – can reduce the cost of lignocellulosic sugars. The reason for setting up a liquefaction unit at the start is to reduce the yield stress of biomass slurry so that the flow remains unrestricted downstream. The liquefaction process is also unique because it uses no chemicals and saves the cost of pretreatment. The pretreatment step is bypassed because the output of liquefaction can be fed into enzyme hydrolysis just after simple cooking. The liquefaction process can be performed in three modes: enzyme, enzyme mimetic, and a combination of enzyme and enzyme mimetic. The results from the BC1707 model indicate the minimum cost for the enzyme liquefaction route. © 2022 The Authors. Biofuels, Bioproducts and Biorefining published by Society of Industrial Chemistry and John Wiley & Sons Ltd.

09 BIOMASS FUELS↗
DOE OSTI2023

Enhancing CO 2 Transport Across a PEEK‐Ionene Membrane and Water‐Lean Solvent Interface

Abstract Efficient direct air capture (DAC) of CO 2 will require strategies to deal with the relatively low concentration in the atmosphere. One such strategy is to employ the combination of a CO 2 ‐selective membrane coupled with a CO 2 capture solvent acting as a draw solution. Here, the interactions between a leading water‐lean carbon‐capture solvent, a polyether ether ketone (PEEK)‐ionene membrane, CO 2 , and combinations were probed using advanced NMR techniques coupled with advanced simulations. We identify the speciation and dynamics of the solvent, membrane, and CO 2 , presenting spectroscopic evidence of CO 2 diffusion through benzylic regions within the PEEK‐ionene membrane, not spaces in the ionic lattice as expected. Our results demonstrate that water‐lean capture solvents provide a thermodynamic and kinetic funnel to draw CO 2 from the air through the membrane and into the bulk solvent, thus enhancing the performance of the membrane. The reaction between the carbon‐capture solvent and CO 2 produces carbamic acid, disrupting interactions between the imidazolium (Im + ) cations and the bistriflimide anions within the PEEK‐ionene membrane, thereby creating structural changes through which CO 2 can diffuse more readily. Consequently, this restructuring results in CO 2 diffusion at the interface that is faster than CO 2 diffusion in the bulk carbon‐capture solvent.

Walter, Eric D.↗
DOE OSTI2023

Vision-based localization for cooperative robot-CNC hybrid manufacturing

Wire and arc additive manufacturing (WAAM) has shown promise in recent years for producing large-scale parts with higher deposition rates than other additive processes. WAAM is often combined with subtractive machining to form a hybrid manufacturing process. This hybrid process can be realized by retrofitting computer numerical control (CNC) machines with deposition heads, adding spindles and deposition heads to robots, or developing part localization methods to transfer parts from an additive cell to a CNC machine. Here, a novel, robot-CNC hybrid configuration is introduced where a maneuverable robot is placed in front of a CNC machine to deposit material within the machine envelop. Furthermore, this method removes the need for part localization and the extensive machine modifications required for retrofitting; however, the problem of robot localization is also added. In this work, the effects of error in vision-based, contactless robot localization on machining parameters in a robot-machine hybrid process were studied. Performance was characterized on an implementation of this system using classical computer vision techniques. In addition, machining simulations were conducted to evaluate the effects of image-induced error on chip thickness, material removal rate, and machining allowance. Initial tests show that computer vision could adequately locate a robot for the hybrid WAAM process within .5 mm.

Hybrid manufacturing↗
DOE OSTI2024

Efficient sensitivity analysis of the thermal profile in powder bed fusion of metals using hypercomplex automatic differentiation finite element method

Rapid cyclic temperature fluctuation occurring in powder bed fusion of metals using a laser beam (PBF-LB/M) influences the formation of flaws in printed parts. Consequently, there is a pressing need to enhance the quality of printed parts by developing innovative methodologies that can predict thermal histories and help uncover the intricate relationships between process parameters and thermal profiles. Sensitivity Analysis (SA) emerges as an essential tool for this, offering the potential for process optimization and enhanced quality control. Nonetheless, conventional SA methodologies often incur in excessive computational costs and potential numerical approximation errors. Here, to address this technical challenge, we present a novel method for SA that integrates the HYPercomplex-based Automatic Differentiation (HYPAD) technique with transient thermal simulations conducted via the finite element method (FEM). Leveraging this methodology, we efficiently and accurately perform SA for PBF-LB/M processes in a post-processing step. Compared to traditional methods like Finite Differences (FD), HYPAD-FEM required 96 % less computational time for obtaining sensitivities for 22 process parameters, under a comparative study conducted within the context of the 2018–02 AM benchmark of the National Institute of Standards and Technology. In summary, HYPAD-FEM offers superior efficiency and accuracy in SA over conventional methods, delivering the best sensitivity of a model without the need for step-size selection and problem or parameter-based implementations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A robust offering strategy for wind producers considering uncertainties of demand response and wind power

This paper proposes a risk-constrained decision-making approach for a wind power producer participating in the day-ahead market. In the developed model, a flexible demand response trading scheme between the wind power producer and different customers is employed. Through the proposed demand response mechanism, the wind power producer is able to trade demand response resource internally with different customers, and then trade energy externally with the market to increase the expected profit and the wind energy utilization. The uncertainties in the wind power and demand response are modeled by using the information gap decision theory approach from risk averse (robust) and risk-seeking (opportunistic) perspectives. The objective of the robust model is to maximize the robust level while satisfying the desired profit, whereas the opportunistic model aims to evaluate the possibility of achieving windfall profits with favorable uncertainties. The overall offering strategy problem is modeled as a bi-objective mixed integer nonlinear programming, which is linearized by proper techniques and solved efficiently by using the normal boundary intersection technique. In this work, simulation results show that utilizing demand response resource to mitigate wind power deviations can increase a wind power producer's profit and reduce potential risks. In addition, the results demonstrate that the proposed bi-objective optimization approach enables the wind power producer to select appropriate offering decisions with respect to uncertainties.

17 WIND ENERGY↗
DOE OSTI2020

Comparative investigations of multi-fidelity modeling on performance of electrostatically-actuated cracked micro-beams

Silicon is a commonly used material for the fabrication of beams for use in micro-electrical-mechanical systems (MEMS). Although silicon is a brittle material, it has been shown to accumulate fatigue damage at the micro-scale. Understanding the effect this has on the overall device performance is critical to the design of reliable devices. Analytical methods for modeling damage provide expedient results but are limited by broad modeling assumptions. Numerical models account for more detailed physical phenomena but can be computationally intensive. In this work, two different crack scenarios are modeled using both analytical techniques and 3D computational simulations. First, the effects of a single surface crack on the static deflection and natural frequency of an electrostatically actuated micro-beam are formulated and compared. Then, a new method for approximating damage associated with realistic distributed crack networks is formulated for use in an analytical model and numerical simulations. A method for utilizing experimentally derived crack statistics to inform the analytical and numerical distributed crack models is developed. Good agreement between the analytical and numerical models is obtained for both crack scenarios. Altogether, these models can be used to effectively simulate a variety of damage and fatigue behaviors in silicon-based MEMS devices.

42 ENGINEERING↗
DOE OSTI2021

The effects of ACRT on melt undercooling during the growth of CZT via the traveling heater method: Ekman versus Taylor-Görtler flows

For the first time, a comprehensive model for the traveling heater method, including detailed species transport and phase-change thermodynamics, is applied to understand the effects of accelerated crucible rotation (ACRT) on the growth of cadmium zinc telluride. Here, a rotation schedule that is designed from classical considerations does not completely mix the melt, rather undercooled regions of localized melt are swept periodically across the growth interface. Even without melt homogenization, these ACRT-driven flows disrupt the persistent undercooling that occurs during THM growth without rotation. An ACRT schedule designed to promote Ekman flows results in lower levels of melt undercooling than a schedule designed to promote Taylor-Görtler flows and greater mixing.

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

Stability-based optimization of ACRT for the growth of CZT by the traveling heater method

Numerical modeling is employed to study the effects of several accelerated crucible rotation technique (ACRT) rotation schedules on the growth of cadmium zinc telluride (CZT) by the traveling heater method (THM). In conjunction with these analyses, a measure is developed to estimate the amplitude of disturbances associated with cellular interface growth morphologies that develop over time and across the melt-solid interface. Furthermore, this disturbance amplitude measure is used to assess the ability of ACRT to stabilize the growth interface and thereby minimize inclusion formation processes. Notably, our analyses reveal that interfacial instability is not well correlated with melt undercooling, explaining why classical ACRT approaches to reduce undercooling by mixing may not stabilize growth. Computations show that interface stability is best achieved via rotation schedules that accentuate outward Ekman flows during spin-up and, during spin-down, eliminate Taylor-Görtler flows and minimize inward Ekman flows.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct sensitivity analysis on the parameterization of crystal plasticity models

Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.

Elementary Effects Test↗
DOE OSTI2022

High strain-rate compression behavior of polymeric rod and plate Kelvin lattice structures

The compressive high strain-rate behavior of polymeric Kelvin lattice structures with rod-based or plate-based unit cells was investigated through experimental techniques and finite element simulations. Polymeric lattice structures with 5x5x5 unit cell geometries were manufactured on the millimeter scale using vat polymerization additive manufacturing and tested at low (0.001/s) and high (1000/s) strain-rates. High strain-rate experiments were performed and validated for a viscoelastic split-Hopkinson (Kolsky) pressure bar system (SHPB) coupled with high-speed imaging and digital image correlation (DIC). Experimental results at both low and high strain-rates show the formation of a localized deformation band which was more prevalent in low relative density specimens and low strain-rate experiments. Strain-rate effects of lattice specimens strongly correlate with effects of the base polymer material; both bulk polymer and lattice specimen demonstrated strain-rate hardening, strain-rate stiffening, and decreased fracture strain under dynamic loading. Results show mechanical failure properties and energy absorption depended strongly on the relative density of the lattice specimen and exhibited distinct scaling between relative density and geometry type (rod, plate) and loading rate. High relative density plate-lattices demonstrated inferior mechanical properties to rod-lattices; however, there exists a critical relative density for a given mechanical property (17%- 28%) below which plate-lattices outperform rod-lattices of similar mass. As a result, high strain-rate explicit finite element simulations were performed and showed good agreement with the mechanical failure trends and deformation modes observed in the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nonlinear optics from hybrid dispersive orbits

Here in this paper we expand the technique of characterizing nonlinear optics from off-energy closed orbits (NOECO) to cover harmonic sextupoles in storage rings. The existing NOECO technique has been successfully used to correct the chromatic sextupole errors on the MAX-IV machine (Olsson et al., 2020), however, it does not account for harmonic sextupoles, which are still widely used on many other machines. By generating vertical dispersion with chromatic skew quadrupoles, a measurable chromatic function dependence on harmonic sextupoles can be observed on hybrid dispersive orbits. Proof of concept of this technique was accomplished by simulations and beam measurements on the National Synchrotron Light Source II (NSLS-II) storage ring.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Residual stress accumulation in large-scale Ti-6Al-4V wire-arc additive manufacturing

Large-scale additive manufacturing (AM) is of interest in the manufacturing industry to produce components of dimensions on the order of hundreds of millimeters to meters in scale. Wire-arc AM (WAAM) of Ti-6Al-4V (Ti64) is an attractive technique for large-scale AM in the aerospace industry, give the high strength-to-weight ratio of the material and high deposition rate of the process. However, due to the large scale, significant distortion and residual stresses are developed in the material during deposition and cooling, potentially leading to part failure. WAAM of a prototypical large-scale Ti64 machine tool component is studied in this work via finite element analysis (FEA). Element activation/deactivation technique is employed to simulate deposition and the resulting distortion and residual stress (RS) predictions are analyzed to evaluate the possibility of crack formation. Significant distortion (∼10 mm) and RS (>1300 MPa) is predicted, suggesting the formation of cracks and possible crack propagation into the build region. Incorporation of fillets significantly reduces RS concentration regions, thus reducing the possibility of part failure. Material deposition sequence also affects the RS pattern in the build.

Mathews, Ritin↗
DOE OSTI2025

Comparative Analysis of Rear Irradiance Modeling Methods for Bifacial PV Systems on Single-Axis Trackers Under Varying Albedo Conditions

This study compares three rear-side irradiance modeling methods for bifacial PV systems on single-axis trackers: (i) the 2D View Factor (VF) model in PVsyst(R), (ii) the open-source PVFactors VF model, and (iii) the Ray Tracing (RT) technique using bifacial_radiance. Simulations were benchmarked against field measurements from a pilot PV plant with rear-side sensors at the torque tube height. Results show that all models underestimated the non-uniformity of rear irradiance along the module length and overestimated the total rear irradiance incident on the module. However, since bifacial gain represents only a fraction of the system's total energy, the resulting energy yield differences among methods remained within +- 2 % of measured values, which is typical for such simulations. While the overall energy impact is limited, this study characterizes the specific limitations of each modeling approach, supporting further improvements in bifacial PV performance assessment methods.

14 SOLAR ENERGY↗
DOE OSTI2023

Characterizing Surface Ice-Philicity Using Molecular Simulations and Enhanced Sampling

The formation of ice, which plays an important role in diverse contexts ranging from cryopreservation to atmospheric science, is often mediated by solid surfaces. Although surfaces that interact favorably with ice (relative to liquid water) can facilitate ice formation by lowering nucleation barriers, the molecular characteristics that confer ice-philicity to a surface are complex and incompletely understood. To address this challenge, here we introduce a robust and computationally efficient method for characterizing surface ice-philicity that combines molecular simulations and enhanced sampling techniques to quantify the free energetic cost of increasing surface–ice contact at the expense of surface–water contact. Using this method to characterize the ice-philicity of a family of model surfaces that are lattice matched with ice but vary in their polarity, we find that the nonpolar surfaces are moderately ice-phobic, whereas the polar surfaces are highly ice-philic. In contrast, for surfaces that display no complementarity to the ice lattice, we find that ice-philicity is independent of surface polarity and that both nonpolar and polar surfaces are moderately ice-phobic. Furthermore, our work thus provides a prescription for quantitatively characterizing surface ice-philicity and sheds light on how ice-philicity is influenced by lattice matching and polarity.

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

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. Here we generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

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

Interfacial Ice Density Fluctuations Inform Surface Ice-Philicity

The propensity of a surface to nucleate ice or bind to ice is governed by its ice-philicity─its relative preference for ice over liquid water. However, the relationship between the features of a surface and its ice-philicity is not well understood, and for surfaces with chemical or topographical heterogeneity, such as proteins, their ice-philicity is not even well-defined. In the analogous problem of surface hydrophobicity, it has been shown that hydrophobic surfaces display enhanced low water-density (vapor-like) fluctuations in their vicinity. To interrogate whether enhanced ice-like fluctuations are similarly observed near ice-philic surfaces, here we use molecular simulations and enhanced sampling techniques. Using a family of model surfaces for which the wetting coefficient, k , has previously been characterized, we show that the free energy of observing rare interfacial ice-density fluctuations decreases monotonically with increasing k . By utilizing this connection, we investigate a set of fcc systems and find that the (110) surface is more ice-philic than the (111) or (100) surfaces. By additionally analyzing the structure of interfacial ice, we find that all surfaces prefer to bind to the basal plane of ice, and the topographical complementarity of the (110) surface to the basal plane explains its higher ice-philicity. Using enhanced interfacial ice-like fluctuations as a measure of surface ice-philicity, we then characterize the ice-philicity of chemically heterogeneous and topologically complex systems. In particular, we study the spruce budworm antifreeze protein (sbwAFP), which binds to ice using a known ice-binding site (IBS) and resists engulfment using nonbinding sites of the protein (NBSs). We find that the IBS displays enhanced interfacial ice-density fluctuations and is therefore more ice-philic than the two NBSs studied. We also find the two NBSs are similarly ice-phobic. By establishing a connection between interfacial ice-like fluctuations and surface ice-philicity, our findings thus provide a way to characterize the ice-philicity of heterogeneous surfaces.

crystals↗
DOE OSTI2026

Molecular Tuning of Ether Cosolvent Chemistry for High-Voltage Sodium-Ion Batteries

Ethers as electrolyte cosolvents in sodium-ion batteries (SIBs) provide favorable Na + solvation and interfacial properties, but their low oxidative stability limits their use in high-voltage SIBs. Herein, we address this limitation via molecular tuning of ether cosolvents for high-voltage (4.2 V) hard carbon || NaNi 0.33 Fe 0.33 Mn 0.33 O 2 full cells. Tetrahydropyran (THP) is functionalized with a nitrile group to form tetrahydropyran-4-carbonitrile (THPCN). To delineate the effect of nitrile functionalization and benchmark ether against a conventional carbonate, THP, THPCN, and diethyl carbonate (DEC) are evaluated as cosolvents with ethylene carbonate. Nitrile functionalization lowers the HOMO energy of the ether, extends the electrolyte stability window, and alters Na⁺ solvation. Spectroscopic techniques and molecular dynamics simulations reveal that THPCN exhibits predominantly aggregate-dominated solvation (95.1 %) with weakened Na + -solvent interactions, producing the most anion-rich environment relative to DEC and THP cosolvents. THPCN-modified solvation promotes the formation of highly conductive, fluorine-enriched interphases that suppress parasitic reactions. Pouch full cells with THPCN sustained ~ 600 cycles at 4.2 V, outperforming THP and DEC. Operando gas analysis reveals that THPCN reduces CO 2 generation by 45% and H 2 generation by 30% relative to THP. Furthermore, the findings demonstrate nitrile functionalization as a molecular design strategy to stabilize ethers and enable high-voltage SIBs.

25 ENERGY STORAGE↗