Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “static models”

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

Static Thermochemical Model of COREX Melter Gasifier

COREX is one of the commercial smelting reduction processes. It uses the finer size ore and semi-soft coal instead of metallurgical coke to produce hot metal from iron ore. The use of top gas with high calorific value as a by-product export gas makes the process economical and green. The predictive thermochemical model of the COREX process presented here enables rapid computation of process parameters such as (1) required amount of ore, coal, and flux; (2) amount of slag and gas generated; and (3) gas compositions (based on the raw material and desired hot metal quality). The model helps in predicting the variations in process parameters with respect to the (1) degree of metallization and (2) post-combustion ratio for given raw material conditions. In general reduction in coal, flux, and oxygen, the requirement is concomitant with an increase in the degree of metallization and post-combustion ratio. The model reported here has been benchmarked using industrial data obtained from the JSW Steel Plant, India.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗
DOE OSTI2012

Sub-Hour Solar Data for Power System Modeling From Static Spatial Variability Analysis: Preprint

High penetration renewable integration studies need high quality solar power data with spatial-temporal correlations that are representative of a real system. This paper will summarize the research relating sequential point-source sub-hour global horizontal irradiance (GHI) values to static, spatially distributed GHI values. This research led to the development of an algorithm for generatingcoherent sub-hour datasets that span distances ranging from 10 km to 4,000 km. The algorithm, in brief, generates synthetic GHI values at an interval of one-minute, for a specific location, using SUNY/Clean Power Research, satellite-derived, hourly irradiance values for the nearest grid cell to that location and grid cells within 40 km.

14 SOLAR ENERGY↗
DOE OSTI2017

Convolution Correction Factor Adjustments on Static PRA Models for Event Assessment

Current probabilistic risk assessment (PRA) models will contain cut sets that have time-constrained basic events such as switched components, recovery of failed components, and failures in time. For these cut sets that have multiple time-constrained events that interact (e.g., offsite power recover at the same time as multiple diesel generators fail to run), there is the possibility that the minimal cut set quantified result is too large. Since current PRA models typically assume that the accident scenario starts at “time zero,” the fails-to-run basic events use this starting time as the basis for their respective reliability models. Over time, the standard practice was to recognize this issue and accept it as part of the conservatism built into the event tree/fault tree PRA models. However, now that these same PRA models are being used in a regulatory forum, this conservatism needs to be addressed. Further, this issue exists not only for the baseline PRA results, but for deficiencies that have been observed or those that have the potential of occurrence, whether the deficiency is complete failure or just a degraded condition. We have found that this conservatism becomes very important when evaluating scenarios such as station blackout sequences and other loss of offsite power sequences. This paper will address the theory and implemented process used by the Systems Analysis Programs for Hands-on Integrated Reliability Evaluations (SAPHIRE) PRA software. Starting with SAPHIRE version 8, the process that was developed is to calculate a “convolution correction factor” based on three distinct calculations. The first calculation is the cut set based frequency calculated via the standard event tree/fault tree PRA model (assuming all run failures start at time equal to zero). The second calculation is to perform the exact time-constrained calculation by integrating the component failure models over the time period of interest (e.g., time to recover, a mission time). For example, this exact calculation for the case of the diesel generators is to convolve their failure models with the offsite power recovery model to obtain the exact frequency. The third calculation is to produce the convolution correction factor by dividing the exact frequency calculation by the standard cut set frequency. This convolution correction factor is then applied to the PRA cut sets in order to remove the time-constrained conservatism. The process used to obtain the convolution correction factor for the baseline PRA is important; however, in this paper we also describe what happens when there is degradation of a component or a potential deficiency. These issues become very important in order to ensure that the convolution correction factor is re-calculated correctly. If the correction factor does not get re-calculated and the same baseline convolution correction factor is used, the final analysis result could be over conservative or under conservative. This paper will discuss how these issues are addressed and provide examples on how SAPHIRE automatically makes the necessary adjustments. MathCad verification results can be computed using the provided equations and input data.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2019

Online Static Load Model Estimation in Distribution Systems

This paper focuses on parameter identification for time-varying voltage-dependent polynomial load models in distribution systems—namely, the constant impedance (Z), constant current (I), and constant power (P). The paper assumes the presence of measurements at the distribution nodes from various sources, each at varying temporal resolutions. The least- squares estimation method with bounded variables is used to convert the measurements into ZIP parameter estimations for loads at various locations on the feeder. The results and discussions focus on some of the estimation issues faced in a distribution system caused by varying voltage sensitivities of loads along a feeder, measurement resolution, impact of window size and sampling rate of measurement data. Finally, the paper uses dynamic simulations data (synthetic high-resolution measurements) to further improve the estimation results.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Doppler-free, Stark broadened profiles at low plasma densities in helium

This paper utilizes Doppler-free saturation spectroscopy to measure Doppler-free, Stark broadened spectral profi les for the π-polarization of the HeI 2¹P → 6¹D transition in a low density (10¹⁰ ₋ 10¹² cm⁻³) helium plasma. The measurements were performed in an electron-cyclotron resonance discharge at 23 mTorr with electron density being diagnosed using a combination of Langmuir probe and microwave interferometry techniques. The Doppler-free pro files were observed to be nearly symmetric at n e < 10¹⁰ cm⁻³ but markedly asymmetric above this transition point. Electron density is extracted from the spectral data via fi tting to a spectral model based on quasi-static Stark broadening and compared to the measured n e . The quasi-static model is shown to be valid in the low-density regime for 2¹P → 6¹D helium transition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Generalized Quasi-Static Mooring System Modeling with Analytic Jacobians

This paper presents a generalized and efficient method for quasi-static analysis of mooring systems, including complex scenarios such as when shared mooring lines interconnect multiple floating wind or wave energy devices. While quasi-static mooring models are well established, most published formulations are focused on specific applications, and no publicly available implementations provide efficient handling of large mooring system networks. The present formulation addresses these gaps by: (1) formulating solutions for edge cases not typically supported by quasi-static models; (2) creating a fully generalized model structure such that any combination of mooring lines, point masses, and floating bodies can be assembled; and (3) deriving analytic expressions for the system Jacobians (stiffness matrices) so that systems with many degrees of freedom can be solved efficiently. These techniques form the theory basis of MoorPy, an open-source mooring analysis library. The model is demonstrated on nine scenarios of increasing complexity with features of interest for offshore renewable energy applications. When compared with steady-state results from a lumped-mass dynamic model, the results show that the quasi-static formulation accurately calculates profiles and tensions and that its analytic approach provides more efficient and reliable computation of system stiffness matrices than finite-differencing methods. These results verify the accuracy of the MoorPy model.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Economic and operational investigation of CO 2 sequestration through enhanced oil recovery in unconventional reservoirs in Colorado, USA

The ongoing CCUS commercial projects are highly relied on the support of government incentives due to massive capital investment. Here, this study analyzes the economics of carbon capture utilization and sequestration (CCUS) projects, shows a state-wide CCUS deployment exercise, followed by simulation results of enhanced oil recovery (EOR) based CO 2 storage in unconventional reservoirs. The comprehensive economic analysis of capture, transportation, sequestration costs, enhanced 45Q tax credits, and EOR revenue implies the practicality of CO 2 -EOR to offset the high CCUS costs. With the economics understanding, we study the top CO 2 sources, existing CO 2 pipelines, and sequestration sinks in the state of Colorado, USA. This paper next presents results from EOR simulation in one section of the unconventional Denver-Julesburg (DJ) Basin Niobrara and Codell reservoirs. The simulation model is based on a geological static model, incorporated with hydraulic fracture stimulation, history matched to production, and calibrated to the microseismic and time-lapse surface seismic data. The CO 2 -EOR simulation results show that oil production can be increased and more CO 2 stored with: a longer primary production period; the presence of a shut-in period; higher injection rates; and multi-well injectors. The modeling results show that about 7–10 Mscf of CO 2 will be stored when recovering 1 stb of EOR oil. By adding the enhanced oil revenue and the carbon credits together, it is estimated that the most economic case can generate $\$13$ MM when oil price is assumed to be $\$80$/stb, and the EOR oil revenue is 3.4 times greater than that generated from 45Q incentives. It corresponds to the scenario that a five-year primary production is followed by CO 2 injection into four wells with the sequence of injection (4 MMscf/day for 6 months), shut-in (6 months) and production (12 months). The best practices in this study will provide valuable insights for similar CCUS projects in other unconventional fields. Furthermore, this study defines a term named “Carbon Neutrality Index (CNI)” by comparing the amount of CO 2 stored with that burned by EOR oil. The CNI value of 0 indicates the enhanced oil is carbon neutral; a negative CNI value implies there is a net reduction in carbon emission. The 4-year huff-n-puff (HnP) simulation leads to a positive CNI value, indicating that the EOR oil generated in this process is not carbon neutral yet.

03 NATURAL GAS↗
DOE OSTI2017

Modeling of Hydraulic Fracture Propagation at the kISMET Site Using a Fully Coupled 3D Network-Flow and Quasi- Static Discrete Element Model

Aimed at supporting the design of hydraulic fracturing experiments at the kISMET site, ~1500 m below ground in a deep mine, we performed pre-experimental hydraulic fracturing simulations in order to estimate the breakdown pressure, propagation pressure, fracture geometry, and the magnitude of induced seismicity using a newly developed fully coupled three-dimensional (3D) network flow and quasi-static discrete element model (DEM). The quasi-static DEM model, which is constructed by Delaunay tessellation of the rock volume, considers rock fabric heterogeneities by using the “disordered” DEM mesh and adding random perturbations to the stiffness and tensile/shear strengths of individual DEM elements and the elastic beams between them. A conjugate 3D flow network based on the DEM lattice is constructed to calculate the fluid flow in both the fracture and porous matrix. One distinctive advantage of the model is that fracturing is naturally described by the breakage of elastic beams between DEM elements. It is also extremely convenient to introduce mechanical anisotropy into the model by simply assigning orientation-dependent tensile/shear strengths to the elastic beams. In this paper, the 3D hydraulic fracturing model was verified against the analytic solution for a penny-shaped crack model. We applied the model to simulate fracture propagation from a vertical open borehole based on initial estimates of rock mechanical properties and in-situ stress conditions. The breakdown pressure and propagation pressure are directly obtained from the simulation. In addition, the released elastic strain energies of individual fracturing events were calculated and used as a conservative estimate for the magnitudes of the potential induced seismic activities associated with fracturing. The comparisons between model predictions and experimental results are still ongoing.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

An automated and portable method for selecting an optimal GPU frequency

Power consumption poses a significant challenge in current and emerging graphics processing unit (GPU) enabled high-performance computing systems. In modern GPUs, dynamic voltage frequency scaling (DVFS) appears to be a reliable control to regulate power consumption and performance. However, the DVFS design space is large - hence, brute-force approaches are infeasible to select the optimal frequency. Furthermore, no single frequency can be universally optimal for applications with varying computational intensities. Thus, the application's complexity and the availability of a wide range of frequency settings are a challenge in selecting the optimal frequency configuration for a given GPU workload. To that end, this paper proposes a systematic approach that consists of three steps. The feature characterization study identifies the fine-grain GPU utilization metrics that influence the power consumption and execution time of a given workload. To understand the performance, power, and energy consumption behaviors of a workload across GPU's DVFS design space, we derived analytical power and performance models using the identified fine-grain features. Here, it is shown that the same set of GPU utilization metrics can estimate both the power consumption and execution time while being agnostic of changes to frequency and input sizes. Applying a power control with the single objective of reducing power may cause performance degradation, leading to more energy consumption. A multi-objective approach is proposed to select the optimal GPU DVFS configuration for a workload that reduces power consumption with negligible degradation in performance. The evaluation was conducted using SPEC ACCEL benchmarks and three real applications - NAMD LAMMPS, and LSTM on NVIDIA GV100, GA100, and AMD MI210 GPUs. On average, real applications showed 29.6% energy savings with a performance loss of 5.2% on GA100 and 22.6% energy savings with a performance loss of 4.7% on GV100. Moreover, the proposed models are portable to real applications, GPU architectures, and vendors, and require metric collection at only the default frequency rather than all supported DVFS configurations. Additionally, we conducted a comparison between our models and the GPU assembly instructions (PTX)-based static models. The results revealed a significant reduction in the average error rates, with a decrease from 19.7% to 3.1% for power models and from 29.4% to 5.2% for performance models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Data-Driven Dynamic Modeling in Power Systems: A Fresh Look on Inverter-Based Resource Modeling

Here, this article introduces ways to identify dynamic system models using measurement data. In power system analysis, a static model represents the time-invariant input-output relationship of a system, while a dynamic model describes the behavior of the system over time. For example, how will a system transit from one steady-state operation point to another?

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Kimberlina - Regional Static 3D Geologic Model of the Southern San Joaquin Basin

Regional static 3D geologic model of the southern San Joaquin Basin for the Westcarb Kimberlina demonstration project created in Earthvision (DGI, Inc.). Also contains a smaller-scale model from the larger basin domain for the geometry used in Birkholzer et al. (2011). Data sets derived from: Wagoner, J., 2009, 3D geologic modeling of the Southern San Joaquin Basin for the Westcarb Kimberlina demonstration project—A status report: Lawrence Livermore National Laboratory LLNL-TR-410813, 27 p. Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011, A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, pp.4371-4378. LLNL contact and file inventory creator: Kayyum Mansoor (mansoor1@llnl.gov)

3D Geological Model↗
DOE OSTI2014

The HNC/HCN ratio in star-forming regions

HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce. We have undertaken a detailed parameter study on which environmental characteristics and chemical reactions affect the HNC/HCN ratio and can thus contribute to the observed dependence. Using existing gas and gas-grain models updated with new reactions and reaction barriers, we find that in static models the H + HNC gas-phase reaction regulates the HNC/HCN ratio under all conditions, except for very early times. We quantitatively constrain the combinations of H abundance and H + HNC reaction barrier that can explain the observed HNC/HCN temperature dependence and discuss the implications in light of new quantum chemical calculations. In warm-up models, gas-grain chemistry contributes significantly to the predicted HNC/HCN ratio and understanding the dynamics of star formation is therefore key to model the HNC/HCN system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Validation of the Modelling of Corrosion Driven Expansion in Grouted Waste-forms - 16125

One of the current preferred solutions for the disposal of Intermediate Level radioactive waste in the UK is grout encapsulation inside steel boxes for long term storage. These boxes are then to be entombed in a Geological Disposal Facility once one is available. However, both prior to and during geological disposal, waste contents can react with the water present in the encapsulant grout and begin to corrode. This corrosion product occupies a larger volume (by several hundred percent) than the original material and the corrosion occurs over the course of decades and centuries, long after the encapsulant grout has set and hardened. The expanding waste then creates stresses in the encapsulant causing both it and the surrounding box to deform and, in extreme circumstances, and after long times, for the structural integrity of the box to fail. This behaviour needs to be quantified before the box can be approved for use as any excessive deformation could lead to a loss of containment or, since the boxes are stacked on top of each other, to a significant accident in the stores. Due to the long timescales involved, physical testing of the phenomenon in full scale is impossible. Instead, a combination of chemical modelling and finite element analysis is used to predict the extent of the corrosion in the waste and its effect of the integrity of the box. There are a number of difficulties in modelling the corrosion-driven expansion and subsequent deformation of the boxes. Although the corrosion occurs over a very long time scale, the event is highly non-linear and dynamic. All of the material models, including the steel and grout are non-linear. The deformation is in the order of 10s of centimetres, well beyond the static elastic range. Finally, the components must be allowed to move relative to each other which involve a lot of part to part contacts. This means that simple static modelling, which would solve quite quickly, cannot be used. Instead a dynamic analysis code, LS-DYNA, has been used to expand the product over time to simulate the corrosion-driven expansion. This software is generally used for modelling car crashes or scenarios where waste containers have been dropped. The analyses therefore include a significant inertial component. This paper will detail the NNL's approach to modelling this long-term, quasi-static phenomenon in the dynamic, impact oriented software code LS-DYNA. It will cover the problems inherent in this approach, the errors that it creates as well as the tests and solutions that NNL worked on in order to overcome them and be confident that the solution could be validated. This work will ultimately be used to determine the maximum inventory which can be stored in each box based on the waste content and the potential extent of expansion as well as any modification that could be made to the box design or the processing streams that may improve the storage solution. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗