Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid modeling”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

ChatGPT and Other Large Language Models for Cybersecurity of Smart Grid Applications

Cybersecurity breaches targeting electrical substations constitute a significant threat to the integrity of the power grid, necessitating comprehensive defense and mitigation strategies. Any anomaly in information and communication technology (ICT) should be detected for secure communications between devices in digital substations. This paper proposes large language models (LLMs), e.g., ChatGPT, for the cybersecurity of IEC 61850-based communications. Multi-cast messages such as generic object oriented system events (GOOSE) and sampled values (SV) are used for case studies. The proposed LLM-based cybersecurity framework includes, for the first time, data pre-processing of communication systems and human-in-the-loop (HITL) training (considering the cybersecurity guidelines recommended by humans). The results show a comparative analysis of detected anomaly data carried out based on the performance evaluation metrics for different LLMs. A hardware-in-the-loop (HIL) testbed is used to generate and extract a dataset of IEC 61850 communications.

ChatGPT↗
NASA NTRS1980

Investigation of transonic flow in a cascade using an adaptive mesh

The solution of two-dimensional full potential equation for the analysis of steady transonic flow through cascades is investigated. Finite element method is employed in the analysis. Accuracy and efficiency of the obtained numerical solutions are discussed in terms of the employed computational grid. Accurate modeling of subsonic and supersonic flow regions together with the shock is discussed. The choice of artificial viscosity and relaxation factors are examined and related to the design of a computational grid. Shock capturing and shock fitting procedures are compared for improved accuracy and efficiency. Numerical results include cascades of Gostelow and NACA 0012 airfoils.

Ecer, A.↗
DOE OSTI2021

Stochastic parametric skeletal dosimetry model for humans: General approach and application to active marrow exposure from bone-seeking beta-particle emitters

The objective of this study is to develop a skeleton model for assessing active marrow dose from bone-seeking beta-emitting radionuclides. This article explains the modeling methodology which accounts for individual variability of the macro- and microstructure of bone tissue. Bone sites with active hematopoiesis are assessed by dividing them into small segments described by simple geometric shapes. Spongiosa, which fills the segments, is modeled as an isotropic three-dimensional grid (framework) of rod-like trabeculae that “run through” the bone marrow. Randomized multiple framework deformations are simulated by changing the positions of the grid nodes and the thickness of the rods. Model grid parameters are selected in accordance with the parameters of spongiosa microstructures taken from the published papers. Stochastic modeling of radiation transport in heterogeneous media simulating the distribution of bone tissue and marrow in each of the segments is performed by Monte Carlo methods. Model output for the human femur at different ages is provided as an example. The uncertainty of dosimetric characteristics associated with individual variability of bone structure was evaluated. An advantage of this methodology for the calculation of doses absorbed in the marrow from bone-seeking radionuclides is that it does not require additional studies of autopsy material. The biokinetic model results will be used in the future to calculate individual doses to members of a cohort exposed to 89,90 Sr from liquid radioactive waste discharged to the Techa River by the Mayak Production Association in 1949–1956. Further study of these unique cohorts provides an opportunity to gain more in-depth knowledge about the effects of chronic radiation on the hematopoietic system. In addition, the proposed model can be used to assess the doses to active marrow under any other scenarios of 90 Sr and 89 Sr intake to humans.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1974

The surface geometry of inherited joint and fracture trace patterns resulting from active and passive deformation

Hypothetical models are considered for detecting subsurface structure from the fracture or joint pattern, which may be influenced by the structure and propagated to the surface. Various patterns of an initially orthogonal fracture grid are modeled according to active and passive deformation mechanisms. In the active periclinal structure with a vertical axis, fracture frequency increased both over the dome and basin, and remained constant with decreasing depth to the structure. For passive periclinal features such as a reef or sand body, fracture frequency is determined by the arc of curvature and showed a reduction over the reefmound and increased over the basin.

Podwysocki, M. H.↗
NASA NTRS1980

A theoretical model of phase changes of a klystron due to variation of operating parameters

A mathematical model for phase changes of the VA-876 CW klystron amplifier output is presented and variations of several operating parameters are considered. The theoretical approach to the problem is based upon a gridded gap modeling with inclusion of a second order correction term so that actual gap geometry is reflected in the formulation. Physical measurements are contrasted to theoretical calculations.

Kupiszewski, A.↗
NASA NTRS2010

Methodology for the Assessment of 3D Conduction Effects in an Aerothermal Wind Tunnel Test

This slide presentation reviews a method for the assessment of three-dimensional conduction effects during test in a Aerothermal Wind Tunnel. The test objectives were to duplicate and extend tests that were performed during the 1960's on thermal conduction on proturberance on a flat plate. Slides review the 1D versus 3D conduction data reduction error, the analysis process, CFD-based analysis, loose coupling method that simulates a wind tunnel test run, verification of the CFD solution, Grid convergence, Mach number trend, size trends, and a Sumary of the CFD conduction analysis. Other slides show comparisons to pretest CFD at Mach 1.5 and 2.16 and the geometries of the models and grids.

Oliver, Anthony Brandon↗
NASA NTRSDate not supplied

Recent Progress on Rans-Based Transition Model Verification

The current efforts to assess and improve the Reynolds-averaged Navier-Stokes (RANS)-coupled transition models in the NASA FUN3D and OVERFLOW codes are summarized in this study. The first AIAA Transition Modeling Workshop and the NATO AVT-313 Transition Workshop both emphasized the need for code verification for transport equations based transition models as a top priority. We discuss the methods used for the model verification, the resulting grid families, the flow solutions, and other supporting information collected with at least two established NASA flow solvers, namely, FUN3D and OVERFLOW. These results, which will be uploaded onto the NASA Turbulence Modeling Resource, should assist other members of the computational fluid dynamics (CFD) community in verifying their own implementations of various transition models, such as the Langtry-Menter (LM2009) model, the one-equation γ model, and Coder’s amplification factor transport (AFT) model. Grid convergence is assessed using both global and local flow metrics of interest such as lift and drag as well as local skin-friction coefficients. We also explore the anisotropic unstructured metric-based adaptive mesh refinement library known as refine with the NASA FUN3D solver to determine if this capability can achieve the same accuracy as handcrafted structured grids with a significantly smaller node count and to learn the characteristics of the resulting grid distribution, especially in the vicinity of the transition zone.

Transition↗
NASA NTRSDate not supplied

Recent Progress on RANS-Based Transition Model Verification

The current efforts to assess and improve the Reynolds-averaged Navier-Stokes (RANS)-coupled transition models in the NASA FUN3D and OVERFLOW codes are summarized in this study. The first AIAA Transition Modeling Workshop and the NATO AVT-313 Transition Workshop both emphasized the need for code verification for transport equations based transition models as a top priority. We discuss the methods used for the model verification, the resulting grid families, the flow solutions, and other supporting information collected with at least two established NASA flow solvers, namely, FUN3D and OVERFLOW. These results, which will be uploaded onto the NASA Turbulence Modeling Resource, should assist other members of the computational fluid dynamics (CFD) community in verifying their own implementations of various transition models, such as the Langtry-Menter (LM2009) model, the one-equation γ model, and Coder’s amplification factor transport (AFT) model. Grid convergence is assessed using both global and local flow metrics of interest such as lift and drag as well as local skin-friction coefficients. We also explore the anisotropic unstructured metric-based adaptive mesh refinement library known as refine with the NASA FUN3D solver to determine if this capability can achieve the same accuracy as handcrafted structured grids with a significantly smaller node count and to learn the characteristics of the resulting grid distribution, especially in the vicinity of the transition zone.

RANS↗
NASA NTRS1989

A numerical study of a class of TVD schemes for compressible mixing layers

At high Mach numbers the two-dimensional time-developing mixing layer develops shock waves, positioned around large-scale vortical structures. A suitable numerical method has to be able to capture the inherent instability of the flow, leading to the roll-up of vortices, and also must be able to capture shock waves when they develop. Standard schemes for low speed turbulent flows, for example spectral methods, rely on resolution of all flow-features and cannot handle shock waves, which become too thin at any realistic Reynolds number. The performance of a class of second-order explicit total variation diminishing (TVD) schemes on a compressible mixing layer problem was studied. The basic idea is to capture the physics of the flow correctly, by resolving down to the smallest turbulent length scales, without resorting to turbulence or sub-grid scale modeling, and at the same time capture shock waves without spurious oscillations. The present study indicates that TVD schemes can capture the shocks accurately when they form, but (without resorting to a finer grid) have poor accuracy in computing the vortex growth. The solution accuracy depends on the choice of limiter. However a larger number of grid points are in general required to resolve the correct vortex growth. The low accuracy in computing time-dependent problems containing shock waves as well as vortical structures is partly due to the inherent shock-capturing property of all TVD schemes. In order to capture shock waves without spurious oscillations these schemes reduce to first-order near extrema and indirectly produce clipping phenomena, leading to inaccuracy in the computation of vortex growth. Accurate simulation of unsteady turbulent fluid flows with shock waves will require further development of efficient, uniformly higher than second-order accurate, shock-capturing methods.

Sandham, N. D.↗
NASA NTRSDate not supplied

Wall-modeled LES of Turbulent Flow Over a Two-dimensional Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a nominally two-dimensional Gaussian-shaped bump geometry to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on high-fidelity numerical simulation of Uzun & Malik (AIAA Journal 2022). The oncoming flow Mach number is 0.2, with the bump length-based Reynolds number (ReL) of 2 million. In our previous study, while WMLES with the constant coefficient Vreman subgrid scale model and equilibrium wall model performed satisfactorily at lower ReL = 1 million, it failed to predict flow separation at the higher Re, contrary to the experimental observations. We investigate the sensitivity of WMLES to different subgrid scale models, wall models, and grid resolution and topology on flow separation by comparing with available data.

Computational Fluid Dynamics↗
DOE OSTI2018

Implementation of a Grid Connected Battery-Inverter Fleet Model

Batteries are designed to store electrical energy. The increasing variation in time value of energy has driven the use of batteries as controllable distributed energy resources (DER). This is enabled though low-cost power electronic inverters that are able to precisely control charge and discharge. This paper describes the software implementation of an open-source battery inverter fleet models in python. The Sandia BatterylnverterFleet class model can be used by scientists, researchers, and engineers to perform simulations of one or more fleets of similar battery-inverter systems connected to the grid. The program tracks the state- of-charge of the simulated batteries and ensures that they stay within their limits while responding to separately generated service requests to charge or discharge. This can be used to analyze control and coordination, placement and sizing, and many other problems associated with the integration of batteries on the power grid. The development of these models along with their python implementation was funded by the Grid Modernization Laboratory Consortium (GMLC) project 1.4.2. Definitions, Standards and Test Procedures for Grid Services from Devices.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2015

Low-cloud characteristics over the tropical western Pacific from ARM observations and CAM5 simulations

Here, this study evaluates the ability of the Community Atmospheric Model version 5 (CAM5) to reproduce low clouds observed by the Atmospheric Radiation Measurement (ARM) cloud radar at Manus Island of the tropical western Pacific during the Years of Tropical Convection. Here low clouds are defined as clouds with their tops below the freezing level and bases within the boundary layer. Low-cloud statistics in CAM5 simulations and ARM observations are compared in terms of their general occurrence, mean vertical profiles, fraction of precipitating versus nonprecipitating events, diurnal cycle, and monthly time series. Other types of clouds are included to put the comparison in a broader context. The comparison shows that the model overproduces total clouds and their precipitation fraction but underestimates low clouds in general. The model, however, produces excessive low clouds in a thin layer between 954 and 930 hPa, which coincides with excessive humidity near the top of the mixed layer. This suggests that the erroneously excessive low clouds stem from parameterization of both cloud and turbulence mixing. The model also fails to produce the observed diurnal cycle in low clouds, not exclusively due to the model coarse grid spacing that does not resolve Manus Island. Lastly, this study demonstrates the utility of ARM long-term cloud observations in the tropical western Pacific in verifying low clouds simulated by global climate models, illustrates issues of using ARM observations in model validation, and provides an example of severe model biases in producing observed low clouds in the tropical western Pacific.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

The relationship between large-scale and convective states in the tropics - Towards an improved representation of convection in large-scale models

This report summarises an investigation into the relationship of tropical thunderstorms to the atmospheric conditions they are embedded in. The study is based on the use of radar observations at the Atmospheric Radiation Measurement site in Darwin run under the auspices of the DOE Atmospheric Systems Research program. Linking the larger scales of the atmosphere with the smaller scales of thunderstorms is crucial for the development of the representation of thunderstorms in weather and climate models, which is carried out by a process termed parametrisation. Through the analysis of radar and wind profiler observations the project made several fundamental discoveries about tropical storms and quantified the relationship of the occurrence and intensity of these storms to the large-scale atmosphere. We were able to show that the rainfall averaged over an area the size of a typical climate model grid-box is largely controlled by the number of storms in the area, and less so by the storm intensity. This allows us to completely rethink the way we represent such storms in climate models. We also found that storms occur in three distinct categories based on their depth and that the transition between these categories is strongly related to the larger scale dynamical features of the atmosphere more so than its thermodynamic state. Finally, we used our observational findings to test and refine a new approach to cumulus parametrisation which relies on the stochastic modelling of the area covered by different convective cloud types.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Model dynamics for quantum computing

A model master equation suitable for quantum computing dynamics is presented. In an ideal quantum computer (QC), a system of qubits evolves in time unitarily and, by virtue of their entanglement, interfere quantum mechanically to solve otherwise intractable problems. In the real situation, a QC is subject to decoherence and attenuation effects due to interaction with an environment and with possible short-term random disturbances and gate deficiencies. The stability of a QC under such attacks is a key issue for the development of realistic devices. We assume that the influence of the environment can be incorporated by a master equation that includes unitary evolution with gates, supplemented by a Lindblad term. Lindblad operators of various types are explored; namely, steady, pulsed, gate friction, and measurement operators. In the master equation, we use the Lindblad term to describe short time intrusions by random Lindblad pulses. The phenomenological master equation is then extended to include a nonlinear Beretta term that describes the evolution of a closed system with increasing entropy. An external Bath environment is stipulated by a fixed temperature in two different ways. Here we explore the case of a simple one-qubit system in preparation for generalization to multi-qubit, qutrit and hybrid qubit–qutrit systems. This model master equation can be used to test the stability of memory and the efficacy of quantum gates. The properties of such hybrid master equations are explored, with emphasis on the role of thermal equilibrium and entropy constraints. Several significant properties of time-dependent qubit evolution are revealed by this simple study. - Highlights: • General realistic density matrix evolution model shown for one qubit. • Unitary pulses & bias yield quantum gates without harmful phase accumulations. • Gate friction, dissipation included by random Lindblad noise pulses. • Bath & entropy terms yield equilibrium; test algorithms, gates act on Larmor grid. • Measurement modeled by strong, fast Lindblad pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Bilevel Nodal Behind-the-meter Solar Disaggregation Under Unexpected Extreme Weather Conditions

As the power grid undergoes significant paradigm shift due to the increasing penetration of renewable generation, the ever-growing installation of behind-the-meter (BTM) solar generation in the power grid also has a significant impact on nodal loads, posing challenges on transmission operators. Furthermore, increasing frequent and severe extreme weather events intertwine with ubiquitous BTM solar generations and have amplified the challenges of accurately model nodal load profiles, especially under the lack of ground-truth information for verification. To tackle these challenges, this paper introduces a bilevel model that utilizes year-long data (e.g., proxy solar, zonal load, and individual node load profiles) to disaggregate metered profiles into actual demand and BTM solar generation at each transmission node. The proxy solar not only scales the BTM solar generation of individual nodes but also create a compensation term for enhancing performance on days with unexpected extreme weather events. The proposed algorithm is validated with real-world PJM Interconnection data during unexpected events like the recent Winter Storm Elliott. For quantitative evaluations, a novel Score error is introduced, which is based on mean percentages and load scales and offers a universal assessment method suitable for all nodes and different data formats (e.g., normalized or raw values).

behind-the-meter solar, load disaggregations, load↗
NASA NTRS2003

Static Aeroelastic Predictions for a Transonic Transport Model Using an Unstructured-Grid Flow Solver Coupled With a Structural Plate Technique

An equivalent-plate structural deformation technique was coupled with a steady-state unstructured-grid three-dimensional Euler flow solver and a two-dimensional strip interactive boundary-layer technique. The objective of the research was to assess the extent to which a simple accounting for static model deformations could improve correlations with measured wing pressure distributions and lift coefficients at transonic speeds. Results were computed and compared to test data for a wing-fuselage model of a generic low-wing transonic transport at a transonic cruise condition over a range of Reynolds numbers and dynamic pressures. The deformations significantly improved correlations with measured wing pressure distributions and lift coefficients. This method provided a means of quantifying the role of dynamic pressure in wind-tunnel studies of Reynolds number effects for transonic transport models.

Allison, Dennis O.↗
NASA NTRS1989

Theoretical studies of the infrared emission from circumstellar dust shells: the infrared characteristics of circumstellar silicates and the mass-loss rate of oxygen-rich late-type giants

We have modeled the infrared emission of spherically symmetric, circumstellar dust shells with the aim of deriving the infrared absorption properties of circumstellar silicate grains and the mass-loss rates of the central stars. As a basis for our numerical studies, a simple semianalytical formula has been derived that illustrates the essential characteristics of the infrared emission of such dust shells. A numerical radiative transfer program has been developed and applied to dust shells around oxygen-rich late-type giants. Free parameters in such models include the absorption properties and density distribution of the dust. An approximate, analytical expression is derived for the density distribution of circumstellar dust driven outward by radiation pressure from a central source. A large grid of models has been calculated to study the influence of the free parameters on the emergent spectrum. These results form the basis for a comparison with near-infrared observations. Observational studies have revealed a correlation between the near-infrared color temperature, Tc, and the strength of the 10 micrometers emission or absorption feature, A10. This relationship, which essentially measures the near-infrared optical depth in terms of the 10 micrometers optical depth, is discussed. Theoretical A10-Tc relations have been calculated and compared to the observations. The results show that this relation is a sensitive way to determine the ratio of the near-infrared to 10 micrometers absorption efficiency of circumstellar silicates. These results as well as previous studies show that the near-infrared absorption efficiency of circumstellar silicate grains is much higher than expected from terrestrial minerals. We suggest that this enhanced absorption is due to the presence of ferrous iron (Fe2+) color centers dissolved in the circumstellar silicates. By using the derived value for the ratio of the near-infrared to 10 micrometers absorption efficiency, the observed A10-Tc relation can be calibrated in terms of the total dust column density of the circumstellar shell and thus the mass-loss rate of late-type giants can easily be derived. Detailed models have been made of the infrared emission of three well-studied Miras: R Cas, IRC 10011, and OH 26.5+0.6, with the emphasis on the shape of the 10 micrometers emission or absorption feature. The results show that the intrinsic shape of the 10 micrometers resonance varies from a very broad feature in R Cas to a relatively narrower feature in OH 26.5+0.6, with IRC 10011 somewhere in between. Possible origins of this variation are discussed. The mass-loss rates from these objects are calculated to be 3 x 10(-7), 2 x 10(-5), and 2 x 10(-4) M Sun yr-1 for R Cas, IRC 10011, and OH 26.5+0.6, respectively. These results are compared to other determinations in the literature.

NASA Discipline Exobiology↗
DOE OSTI2021

Self-Security for Grid-Interactive Smart Inverters Using Steady-State Reference Model

Smart inverters exchange information with other devices through a shared communication link, making the inverters more prone to receive harmful commands from external parties. This erroneous data can be received due to an anomaly in the system, such as a device fault, unintentional utility operator action, or a cyber-attack. In this paper, a device-level self-security strategy is implemented using reference models for a grid-interactive inverter to examine the incoming power setpoints, detect the anomalies, and protect the system accordingly. The PQ setpoints received from the utility supervisory controller are autonomously examined using the inverter’s normal and stable operating regions before engaging the setpoints to the inverter’s local controller. Grid parameters are estimated in real-time during the examination process. The efficacy of the self-security algorithm is tested using a three-phase 3-kVA SiC-MOSFET inverter and a 12-kW NHR 9410 regenerative grid emulator. The results verify that the proposed method can detect harmful PQ setpoints that can cause abnormal or unstable inverter operation.

Gursoy, Mehmetcan↗