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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2003

Combining Static Analysis and Model Checking for Software Analysis

We present an iterative technique in which model checking and static analysis are combined to verify large software systems. The role of the static analysis is to compute partial order information which the model checker uses to reduce the state space. During exploration, the model checker also computes aliasing information that it gives to the static analyzer which can then refine its analysis. The result of this refined analysis is then fed back to the model checker which updates its partial order reduction. At each step of this iterative process, the static analysis computes optimistic information which results in an unsafe reduction of the state space. However we show that the process converges to a fired point at which time the partial order information is safe and the whole state space is explored.

Brat, Guillaume↗
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↗
NASA NTRS1989

Central cores of post-collapse clusters

By means of conventional static models, the modelling of post-collapse-core star clusters is shown. The rationale adopted is that core collapse will cause the formation of binary stars, which stabilize the cluster and make it static again, but with a very high central concentration. It is shown that the models require massive remnant stars at the cluster centers, and that the fitting of static models to post-collapse-core clusters determines how many such stars are needed.

King, Ivan R.↗
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↗
NASA NTRS1986

The structure of the static corona and transition region

Static models of coronal loops are investigated. For loops that are low-lying with heights above the chromosphere below about 5000 km, it is shown that a new type of solution appears to the static equations, in addition to the well-known coronal loop solution. The new solution is characterized by a maximum plasma temperature less than about 100,000 K. The structure and properties of these cool solutions are discussed. The differential emission measure Q(T) expected for a magnetic arcade, which must naturally contain both hot and cool loops, is calculated. It is shown that the cool loops have a dramatic effect on the form of Q(T) in the lower transition region. In particular, they can account for the observed rise in Q at low T, which has long been thought to be incompatible with the static-loop model. Finally, the implications of the cool loops on other observations of both the solar and stellar coronae and transition regions are discussed.

Antiochos, S. K.↗
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↗
NASA NTRS2018

A Dual Spring Modeling Approach for Static and Fatigue Failure Assessments of Carbon/Epoxy Composite Sub-Elements

A dual spring model is developed for the static and fatigue damage predictions of laminates interface in composite structures. Stress concentrations can be induced by the defects formed in the fabrication or service process. A conventional S-N based fatigue damage model may not be accurate to predict the fatigue life of a structure with high stress concentration. With the dual spring model, static delamination failure can be simulated using springs of cohesive type material model while fatigue delamination development can be predicted using linear springs, where the crack driving force is computed based on virtual crack closure technique (VCCT). A Paris law type fatigue growth law with its mode mixity is applied for fatigue crack growth prediction. After verified using benchmark examples, including Double Cantilever Beam (DCB), End-Notched Flexure (ENF) and Mix-Mode Bending (MMB), the proposed dual spring model is applied in the static and fatigue damage prediction of NASA/Boeing sub-elements and UTC sub-elements.

Cui, Xiaodong↗
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↗
NASA NTRS1979

Numerical modeling of quasi-static coronal loops. I - Uniform energy input

A quasi-static numerical model for coronal loops is considered for the case of a uniform energy input per unit volume into the loops. A line dipole model is used to represent the loop magnetic field, and the variations in loop cross section observed in X-ray photographs are parameterized by the ratio between the cross-sectional areas at the loop apex and base. The results of numerical modeling indicate that for an area ratio greater than unity, increases in the area ratio of a loop with a given length and apex area cause a general rise in electron density and a fall in the temperature gradient, leading to large increases in the differential emission factor at high temperatures. The differential function obtained is significantly different from that predicted by analytical models; however, analytical predictions for the temperature-electron density relations are comparable to numerical results. It is also concluded that even a symmetrical loop may have a maximum temperature away from the apex.

Vesecky, J. F.↗
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↗
NASA NTRS1975

Instability of a gravity gradient satellite due to thermal distortion

A nonlinear analytical model and a corresponding computer program were developed to study the influence of solar heating on the anomalous low frequency, orbital instability of the Naval Research Laboratory's gravity gradient satellite 164. The model's formulation was based on a quasi-static approach in which deflections of the satellite's booms were determined in terms of thermally induced bending without consideration of boom vibration. Calculations, which were made for variations in absorptivity, sun angle, thermal lag, and hinge stiffness, demonstrated that, within the confines of a relatively narrow stability criteria, the quasi-static model of NRL 164 not only becomes unstable, but, in a number of cases, responses were computed that closely resembled flight data.

Goldman, R. L.↗
NASA NTRS2005

Influence of Coronal Abundance Variations

During the final year of this program we concentrated on understanding the how to constrain the models with the best available observations. Work on developing accurate temperature and density diagnostics fkom TRACE and CDS together with constrained fits of non-potential force free fields will be extremely useful in the guiding the next generation of coronal models. The program has produced three fully operation numerical codes that model multi-species of ions in coronal loops: Static models and constant flow models. The time dependent numerical models have not been completed. We have extended the steady flow investigations to study the effect these flows have on coronal structure as observed with TRACE. Coronal observations derive from heavy-ion emission; thus, we focus on the extent to which flow may modify coronal abundances by examining the heavy-ion abundance stratification within long-lived loops. We discuss the magnitudes of the physical effects modeled and compare simulated results with TRACE observations. These results can have a profound effect on the interpretation of TRACE observations.

Gurman, Joseph↗
NASA NTRS1997

A Model of the THUNDER Actuator

A THUNDER actuator is a composite of three thin layers, a metal base, a piezoelectric wafer and a metal top cover, bonded together under pressure and at high temperature with the LaRC SI polyimid adhesive. When a voltage is applied between the metal layers across the PZT the actuator will bend and can generate a force. This document develops and describes an analytical model the transduction properties of THUNDER actuators. The model development is divided into three sections. First, a static model is described that relates internal stresses and strains and external displacements to the thermal pre-stress and applied voltage. Second, a dynamic energy based model is described that allows calculation of the resonance frequencies, developed force and electrical input impedance. Finally, a fully coupled electro-mechanical transducer model is described. The model development proceeds by assuming that both the thermal pre-stress and the piezoelectric actuation cause the actuator to deform in a pure bend in a single plane. It is useful to think of this as a two step process, the actuator is held flat, differential stresses induce a bending moment, the actuator is released and it bends. The thermal pre-stress is caused by the different amounts that the constituent layers shrink due to their different coefficients of thermal expansion. The adhesive between layers sets at a high temperature and as the actuator cools, the metal layers shrink more than the PZT. The PZT layer is put into compression while the metal layers are in tension. The piezoelectric actuation has a similar effect. An applied voltage causes the PZT layer to strain, which in turn strains the two metal layers. If the PZT layer expands it will put the metal layers into tension and PZT layer into compression. In both cases, if shear force effects are neglected, the actuator assembly will experience a uniform in-plane strain. As the materials each have a different elastic modulus, different stresses will develop in each layer and these stresses will induce a bending moment. When the actuator is released from its flat configuration, the differential stresses are relieved as the actuator bends.

Curtis, Alan R. D.↗
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↗
NASA NTRS1990

A model for the space shuttle main engine high pressure oxidizer turbopump shaft seal system

A simple static model is presented which solves for the flow properties of pressure, temperature, and mass flow in the Space Shuttle Main Engine pressure Oxidizer Turbopump Shaft Seal Systems. This system includes the primary and secondary turbine seals, the primary and secondary turbine drains, the helium purge seals and feed line, the primary oxygen drain, and the slinger/labyrinth oxygen seal pair. The model predicts the changes in flow variables that occur during and after failures of the various seals. Such information would be particularly useful in a post flight situation where processing of sensor information using this model could identify a particular seal that had experienced excessive wear. Most of the seals in the system are modeled using simple one dimensional equations which can be applied to almost any seal provided that the fluid is gaseous. A failure is modeled as an increase in the clearance between the shaft and the seal. Thus, the model does not attempt to predict how the failure process actually occurs (e.g., wear, seal crack initiation). The results presented were obtained using a FORTRAN implementation of the model running on a VAX computer. Solution for the seal system properties is obtained iteratively; however, a further simplified implementation (which does not include the slinger/labyrinth combination) was also developed which provides fast and reasonable results for most engine operating conditions. Results from the model compare favorably with the limited redline data available.

Paxson, Daniel E.↗