Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Results”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Selecting durable building envelope systems with machine learning assisted hygrothermal simulations database

Hygrothermal simulations provide insight into the energy performance and moisture durability of building envelope components under dynamic conditions. The inputs required for hygrothermal simulations are extensive, and carrying out simulations and analyses requires expert knowledge. An expert system, the Building Science Advisor (BSA), has been developed to predict the performance and select the energy-efficient and durable building envelope systems for different climates. The BSA consists of decision rules based on expert opinions and thousands of parametric simulation results for selected wall systems. The number of potential wall systems results in millions, too many to simulate all of them. We present how machine learning can help predict durability data, such as mold growth, while minimizing the number of simulations needed to run. The simulation results are used for training and validation of machine learning tools for predicting wall durability. We tested Artificial Neural Network (ANN) and Gradient Boosted Decision Trees (GBDT) for their applicability and model accuracy. Models developed with both methods showed adequate prediction performance (root mean square error of 0.195 and 0.209, respectively). Finally, we introduce how the information supports guidance for envelope design via an easy-to-use web-based tool that does not require the end-user to run hygrothermal simulations.

Salonvaara, Mikael↗
DOE OSTI2020

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3, documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation is to describe the application of the hydraulic property fields and recharge parameters documented in ECF-HANFORD-20-0027, Null Space Monte Carlo Evaluation of the Plateau to River Model, to the CA flow and fate and transport simulation results to quantify the uncertainty in the simulated results due to input parameter selection. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, models lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. The result of NSCM analysis is a set of F&T simulations that provide an estimate of the range of possible outcomes that are used to quantify the uncertainty in simulated concentrations produced using the base case simulations. The simulated concentrations from all simulations will support calculation of the uncertainty of the total dose calculated in a separate calculation.

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

Simulation of Plasma Emission in Magnetized Plasmas

The recent Parker Solar Probe observations of type III radio bursts show that the effects of the finite background magnetic field can be an important factor in the interpretation of data. In the present paper, the effects of the background magnetic field on the plasma-emission process, which is believed to be the main emission mechanism for solar coronal and interplanetary type III radio bursts, are investigated by means of the particle-in-cell simulation method. The effects of the ambient magnetic field are systematically surveyed by varying the ratio of plasma frequency to electron gyrofrequency. The present study shows that for a sufficiently strong ambient magnetic field, the wave–particle interaction processes lead to a highly field-aligned longitudinal mode excitation and anisotropic electron velocity distribution function, accompanied by a significantly enhanced plasma emission at the second-harmonic plasma frequency. For such a case, the polarization of the harmonic emission is almost entirely in the sense of extraordinary mode. On the other hand, for moderate strengths of the ambient magnetic field, the interpretation of the simulation result is less clear. The underlying nonlinear-mode coupling processes indicate that to properly understand and interpret the simulation results requires sophisticated analyses involving interactions among magnetized plasma normal modes, including the two transverse modes of the magneto-active plasma, namely, the extraordinary and ordinary modes, as well as electron-cyclotron-whistler, plasma oscillation, and upper-hybrid modes. At present, a nonlinear theory suitable for quantitatively analyzing such complex-mode coupling processes in magnetized plasmas is incomplete, which calls for further theoretical research, but the present simulation results could provide a guide for future theoretical efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Numerical modeling of impurity powder injection in W7-X

In this work, we present numerical simulation results of powder injection in W7-X using the EMC3-EIRENE and DIS codes. First, we model powder injection experiments performed in W7-X with the Probe Mounted Powder Injector. The simulation results qualitatively agree with visible imaging measurements. Secondly, we perform predictive simulations to guide the installation of an Impurity Powder Dropper in W7-X, allowing to choose in between several available non-vertical ports to maximize the amount of powder penetrating into the plasma, as well as the verticality of the port, to minimize sticking of the powders in the in-vessel stainless steel guiding tube. Port AEM41 is selected as the best candidate for IPD installation. The robustness of the simulation results has been verified for different plasma densities, powder materials and sizes, powder friction coefficient and changes in the plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The Biomechanics of Indirect Traumatic Optic Neuropathy Using a Computational Head Model With a Biofidelic Orbit

Indirect traumatic optic neuropathy (ITON) is an injury to the optic nerve due to head trauma and usually results in partial or complete loss of vision. In order to advance a mechanistic understanding of the injury to the optic nerve, we developed a head model with a biofidelic orbit. Head impacts were simulated under controlled conditions of impactor velocity. The locations of impact were varied to include frontal, lateral, and posterior parts of the head. Impact studies were conducted using two types of impactors that differed in their rigidity relative to the skull. The simulated results from both the impactors suggest that forehead impacts are those to which the optic nerve is most vulnerable. The mode and location of optic nerve injury is significantly different between the impacting conditions. Simulated results using a relatively rigid impactor (metal cylinder) suggest optic nerve injury initiates at the location of the intracranial end of the optic canal and spreads to the regions of the optic nerve in the vicinity of the optic canal. In this case, the deformation of the skull at the optic canal, resulting in deformation of the optic nerve, was the primary mode of injury. On the other hand, simulated results using a relatively compliant impactor (soccer ball) suggest that primary mode of injury comes from the brain tugging upon the optic nerve (from where it is affixed to the intracranial end of the optic canal) during coup countercoup motion of the brain. This study represents the first published effort to employ a biofidelic simulation of the full length of the optic nerve in which the orbit is integrated within the whole head.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Comparison of Simulation of Dual-Active Bridge in Different Simulators and Using Different Simulation Methods

The accuracy of states measured from simulation of a dual-active bridge (DAB) converter with an open-loop controller has been used for comparing simulation software and algorithms. The accuracy of the measured states is sensitive to simulation settings and parameters chosen. In this paper, the simulation of the DAB converter with open-loop control is run in the MATLAB/Simscape Electrical, PLECS, and PSCAD software environments, with various simulation settings. The simulation results are also compared with a new simulation algorithm proposed for DAB simulations. Different simulation settings are evaluated, including different solver types, operating conditions, and the on-state resistance values of the MOSFET switches. The simulation results are compared with the theoretical calculations. It is shown that unexpected divergence and DC offset may occur in some cases. This paper provides guidelines to simulate the DAB converter with open-loop control.

Yin, Shiyuan↗
DOE OSTI2022

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗
DOE OSTI2020

Towards the design of a wideband reflective long period grating distributed sensor

Abstract In this paper, we computationally investigate the effects of metal coating length and coating coverage on the reflected spectrum of a long period grating (LPG) over a broad bandwidth. Simulation results indicate that coating the tail end of the fiber between the LPG and the end facet of the fiber provides a reflected spectrum that mimics the LPG transmission spectrum shape over a 400 nm bandwidth. Based on single LPG simulation results, we present the design of a distributed LPG structure containing a multiple number ( n ) of LPGs in reflection mode for the first time. Simulation results for n = 1, 2, and 3 are presented here to demonstrate the concept of a distributed reflective LPG design. It is expected that such a sensor will open a new window for distributed sensing using reflective LPGs.

Rana, Sohel (ORCID:0000000262813555)↗
DOE OSTI2020

Coupled confined phase behavior and transport of methane in slit nanopores

The proximity of the order of magnitude molecular mean free path and pore sizes in nanopores leads to remarkable interactions between molecules and walls. In such systems, the thermodynamic property and the transport behavior of fluids deviate from those at bulk conditions. Molecular dynamics (MD) simulation may be used to investigate the effects of confinement on fluid physics in nanosize pores. However, translating these subscale observations into a larger scale of interest remains a challenging task. In this work, we propose a modified extension of Peng-Robinson equation of state (PR-EOS) motivated and guided by MD simulation results. The shift of critical properties, i.e., pressure and temperature, are evaluated independently. A temperature dependent parameter is introduced to account for the existence of capillary pressure in the two-phase region. Here, this formulation is capable of capturing the shift in critical properties as well as density phase diagrams under various confinement scenarios. We incorporate the proposed EOS in the lattice Boltzmann method (LBM) to study coupled confined phase behavior and transport of methane in nano-size slit pores. Adsorptive strength is determined such that the resulting density ratio matches that from MD simulation. Results indicate that the use of the proposed EOS leads to smaller fluid velocities compared to PR-EOS. This effect is due to stronger interactions between fluid particles under confinement. Also, results show that transport characteristics are impacted by the pressure of the system: in systems with a relatively low pressure, transport seems to be dominated by Knudsen diffusion; however, at a higher pressure, the contribution from viscous flow increases as the pore widens while the influence of Knudsen and surface diffusion diminishes.

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

Investigating the Impact of Thickness, Calendering and Channel Structures of Printed Electrodes on the Energy Density of LIBs - 3D Simulation and Validation

Current lithium ion batteries (LIBs) are expensive and bulky, limited by relatively low charging rates. To increase the rate of charging and reduce weight, thin electrodes with high energy density are required. The increase in energy density can be achieved by several techniques including boosting electrolyte transport, high loading/utilization of active material, employing high conductive electrolytes and electrodes with advanced architectures, and increasing cell temperature. In this paper, a 3D physics-based electrochemical model of LIBs is developed in COMSOL simulation software for different thickness, calendering steps as well as channel structures (conical, cylindrical) to optimize the electrode design and in turn maximize volumetric energy density. The simulation results demonstrated that calendering the electrodes with high initial porosity increases the volumetric energy density of the cell. In addition, cylindrical channel structures with relatively lower edge-to-edge distance also results in increased volumetric energy density. The simulation results of the 3D model was validated by comparing it with experimental results.

improving volumetric energy density↗
DOE OSTI2023

Model Predictive Voltage Control of Large-Scale PV or Hybrid PV-BESS Plants

Increased penetration level of inverter-based resources (IBRs) and renewable energy in the power grid has called for more requirements from the control of IBRs. In this paper, the voltage-reactive power control of a photovoltaic (PV) power plant or hybrid PV-battery energy storage system (BESS) plant connected to a bulk power grid, whilst meeting the grid requirements, is studied. In this paper, a continuous-set model predictive control (MPC) formulation is proposed for voltage-reactive power control. For the same, an aggregated dynamic PV plant model is developed based on the recursive system equivalencing method. The formulation is then implemented in PSCAD simulation on a 125 MW PV plant or hybrid PV-BESS plant. The MPC implementation achieves 10.13% voltage improvement based on PSCAD simulation results in a balanced fault case study. Simulation results further demonstrate that MPC provides improved voltage support over the conventional proportional-integral (PI) controller post-fault occurrence.

Abu rub, Omar↗
DOE OSTI2025

Microstructure-based modeling of inner oxygen pressure in solid oxide electrolysis cells: Analysis of electrode delamination and mitigation

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here a microstructure-based electrochemical model, which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

25 ENERGY STORAGE↗
DOE OSTI2025

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. Furthermore, the simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

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

Design Considerations of a Coordinative Demand Charge Mitigation Strategy

This paper presents a coordinative demand charge mitigation (DCM) strategy for reducing electricity consumption during system peak periods. Available DCM resources include batteries, diesel generators, controllable loads, and conservation voltage reduction. All resources are directly controlled by load serving entities. A mixed integer linear programming based energy management algorithm is developed to optimally coordinate of DCM resources considering the load payback effect. To better capture system peak periods, two different kinds of load forecast are used: the day-ahead load forecast and the peak-hour probability forecast. Five DCM strategies are compared for reconciling the discrepancy between the two forecasting results. The DCM strategies are tested using actual utility data. Simulation results show that the proposed algorithm can effectively mitigate the demand charge while preventing the system peak from being shifted to the payback hours. We also identify the diminishing return effect, which can help load serving entities optimize the size of their DCM resources.

Hu, Rongxing↗
DOE OSTI2022

A Laplace-Domain Circuit Model for Fault and Stability Analysis Considering Unbalanced Topology

For systems subject to unbalanced faults, analytical model building for stability assessment is a challenging task. This letter presents a straightforward modeling approach. A generalized dynamic circuit representation is achieved by use of the Laplacian transform variable s . Here, we translate the voltage and current relationship at the fault location into the relationship of three subsystems. The final circuit model is an interconnected sequence network with impedances in the Laplace domain. This circuit can be directly converted from a steady-state sequence network. This modeling procedure is illustrated by an example case of an induction motor served by a grid through a series compensated line. Electromagnetic transient simulation results demonstrate that sub-synchronous oscillations can be mitigated when a single-line to ground fault is applied at the motor terminal. Stability analysis results based on the dynamic circuit corroborate the simulation results. What's more, the derived circuit effortlessly reveals why unbalance can enhance stability.

42 ENGINEERING↗
DOE OSTI2022

Isogeometric large-eddy simulations of turbulent particle-laden flows

In recent years, isogeometric analysis (IGA) has attracted significant attention from the computational mechanics community due to its ability to integrate design and analysis. Besides, IGA is also a higher-order discretization technique for solving partial differential equations, showing high approximation capability per degree of freedom. In this paper, we extend the application realm of IGA to particle-laden flows based on Eulerian–Eulerian description that couples Navier–Stokes equations with a density transport equation through a Boussinesq approximation. The coupled systems are solved by using quadratic non-uniform rational B-spline (NURBS) functions and a recently developed residual-based variational multiscale (VMS) formulation, which introduces coupling between the fine velocity scales and density equation residuals. We deploy the proposed approach to perform large-eddy simulations (LES) of dilute particle-laden flows over a flat surface at Reynolds number = 10,000. We compare the simulation results against direct numerical simulation (DNS) results from the literature. We find that combining VMS and IGA, the proposed approach enables accurate prediction of a wide range of flow/particle statistics with a relatively lower mesh resolution.

Mathematics↗
DOE OSTI2022

DSO+T: Valuation Methodology and Economic Metrics (DSO+T Study: Volume 4)

This report summarizes a rigorous valuation analysis methodology used by the Distribution System Operator with Transactive (DSO+T) study to estimate the financial benefits and costs of adopting Transactive Energy coordination of distributed energy resources for key stakeholders (for example distribution system operators and customers). This was achieved by modeling the value exchanges between stakeholders and determining the annualized costs and revenues experienced by stakeholders, enabling the evaluation of overall impact on stakeholder’s annualized cash flow. Extensive work was conducted developing methods to estimate the operating costs of distribution system operators at a level of granularity that would allow the financial impact of implementing a transactive energy approach to be estimated. This work included developing parametric models for labor and software costs, distribution system capital and maintenance costs, growth rates, and factors to determine annualized costs of capital investments. Simulation results were used to calculate wholesale energy costs and revenues from retail sales. Valuation analysis methods were also developed for other stakeholders including customers, the transmission system operator, independent system operator, and generators. The resulting capability allows a complete mapping of the flow of financial value between stakeholders that can be directly integrated with the results of demand flexibility simulations. Example results are provided for a business-as-usual case and compared to a transactive energy case as well as to actual cost data.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Microstructure-Based Modeling of Inner Oxygen Pressure in Solid Oxide Electrolysis Cells

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here, a microstructure-based electrochemical model which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

Xue, Fei↗