Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “grid simulator”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Numerical simulation of Shuttle ascent transonic flow using an unstructured-grid approach

An unstructured-grid, finite-volume method has been developed for simulating the inviscid flow over spacecrafts of realistic configuration. The grid generation is accomplished by a new technique on the basis of the advancing-front concept. This simple technique is shown to be equally as powerful for a complex multibody as for a single vehicle. Second- or third-order accuracy is obtained via an innovative interpolation procedure similar to the conventional MUSCL approach. This method has been applied to the Shuttle orbiter and a representative Shuttle launch vehicle consisting of the orbiter, the external tank, and the solid rocket boosters. A comparison is discussed between the present results and other results obtained from structured- and unstructured-grid methods.

Wey, T. C.↗
DOE OSTI2010

MONTE CARLO SIMULATIONS OF GLOBULAR CLUSTER EVOLUTION. V. BINARY STELLAR EVOLUTION

We study the dynamical evolution of globular clusters containing primordial binaries, including full single and binary stellar evolution using our Monte Carlo cluster evolution code updated with an adaptation of the single and binary stellar evolution codes SSE and BSE from Hurley et al. We describe the modifications that we have made to the code. We present several test calculations and comparisons with existing studies to illustrate the validity of the code. We show that our code finds very good agreement with direct N-body simulations including primordial binaries and stellar evolution. We find significant differences in the evolution of the global properties of the simulated clusters using stellar evolution compared with simulations without any stellar evolution. In particular, we find that the mass loss from the stellar evolution acts as a significant energy production channel simply by reducing the total gravitational binding energy and can significantly prolong the initial core contraction phase before reaching the binary-burning quasi-steady state of the cluster evolution. We simulate a large grid of models varying the initial cluster mass, binary fraction, and concentration parameter, and we compare properties of the simulated clusters with those of the observed Galactic globular clusters (GGCs). We find that simply including stellar evolution in our simulations and assuming the typical initial cluster half-mass radius is approximately a few pc independent of mass, our simulated cluster properties agree well with the observed GGC properties such as the core radius and the ratio of the core radius to the half-mass radius. We explore in some detail qualitatively different clusters in different phases of their evolution and construct synthetic Hertzsprung-Russell diagrams for these clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

THE ANGULAR MOMENTUM OF MAGNETIZED MOLECULAR CLOUD CORES: A TWO-DIMENSIONAL-THREE-DIMENSIONAL COMPARISON

In this work, we present a detailed study of the rotational properties of magnetized and self-gravitating dense molecular cloud (MC) cores formed in a set of two very high resolution three-dimensional (3D) MC simulations with decaying turbulence. The simulations have been performed using the adaptative mesh refinement code RAMSES with an effective resolution of 4096{sup 3} grid cells. One simulation represents a mildly magnetically supercritical cloud and the other a strongly magnetically supercritical cloud. We identify dense cores at a number of selected epochs in the simulations at two density thresholds which roughly mimic the excitation densities of the NH{sub 3} (J - K) = (1,1) transition and the N{sub 2}H{sup +} (1-0) emission line. A noticeable global difference between the two simulations is the core formation efficiency (CFE) of the high-density cores. In the strongly supercritical simulations, the CFE is 33% per unit free-fall time of the cloud (t{sub ff,cl}), whereas in the mildly supercritical simulations this value goes down to {approx}6 per unit t{sub ff,cl}. A comparison of the intrinsic specific angular momentum (j{sub 3D}) distributions of the cores with the specific angular momentum derived using synthetic two-dimensional (2D) velocity maps of the cores (j{sub 2D}) shows that the synthetic observations tend to overestimate the true value of the specific angular momentum by a factor of {approx}8-10. We find that the distribution of the ratio j{sub 3D}/j{sub 2D} of the cores peaks at around {approx}0.1. The origin of this discrepancy lies in the fact that contrary to the intrinsic determination of j which sums up the individual gas parcels' contributions to the angular momentum, the determination of the specific angular momentum using the standard observational procedure which is based on a measurement on the global velocity gradient under the hypothesis of uniform rotation smoothes out the complex fluctuations present in the 3D velocity field. Our results may well provide a natural explanation for the discrepancy by a factor of {approx}10 observed between the intrinsic 3D distributions of the specific angular momentum and the corresponding distributions derived in real observations. We suggest that previous and future measurements of the specific angular momentum of dense cores which are based on the measurement of the observed global velocity gradients may need to be reduced by a factor of {approx}10 in order to derive a more accurate estimate of the true specific angular momentum in the cores. We also show that the exponent of the size-specific angular momentum relation is smaller ({approx}1.4) in the synthetic observations than their values derived in the 3D space ({approx}1.8).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

CALIBRATION OF THE MIXING-LENGTH THEORY FOR CONVECTIVE WHITE DWARF ENVELOPES

A calibration of the mixing-length parameter in the local mixing-length theory (MLT) is presented for the lower part of the convection zone in pure-hydrogen-atmosphere white dwarfs. The parameterization is performed from a comparison of three-dimensional (3D) CO5BOLD simulations with a grid of one-dimensional (1D) envelopes with a varying mixing-length parameter. In many instances, the 3D simulations are restricted to the upper part of the convection zone. The hydrodynamical calculations suggest, in those cases, that the entropy of the upflows does not change significantly from the bottom of the convection zone to regions immediately below the photosphere. We rely on this asymptotic entropy value, characteristic of the deep and adiabatically stratified layers, to calibrate 1D envelopes. The calibration encompasses the convective hydrogen-line (DA) white dwarfs in the effective temperature range 6000 ≤ T {sub eff} (K) ≤15, 000 and the surface gravity range 7.0 ≤ log g ≤ 9.0. It is established that the local MLT is unable to reproduce simultaneously the thermodynamical, flux, and dynamical properties of the 3D simulations. We therefore propose three different parameterizations for these quantities. The resulting calibration can be applied to structure and envelope calculations, in particular for pulsation, chemical diffusion, and convective mixing studies. On the other hand, convection has no effect on the white dwarf cooling rates until there is a convective coupling with the degenerate core below T {sub eff} ∼ 5000 K. In this regime, the 1D structures are insensitive to the MLT parameterization and converge to the mean 3D results, hence they remain fully appropriate for age determinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Triple-Moment Representation of Ice in the Predicted Particle Properties (P3) Microphysics Scheme

In the original Predicted Particle Properties (P3) bulk microphysics scheme, all ice-phase hydrometeors are represented by one or more “free” ice categories, where the physical properties evolve smoothly through changes to the four prognostic variables (per category), and with a two-moment representation of the particle size distribution. As such, the spectral dispersion cannot evolve independently, which thus results in limitations in representation of ice—in particular, hail—due to necessary constraints in the scheme to prevent excessive gravitational size sorting. To overcome this, P3 has been modified to include a three-moment representation of the size distribution of each ice category through the addition of a fifth prognostic variable, the sixth moment of the size distribution. The details of the three-moment ice parameterization in P3 are provided. The behavior of the modified scheme, with the single-ice-category configuration, is illustrated through simulations in a simple 1D kinematic model framework as well as with near large-eddy-resolving (250-m grid spacing) 3D simulations of a hail-producing supercell. It is shown that the three-moment ice configuration controls size sorting in a physically based way and leads to an improved capacity to simulate large, heavily rimed ice (hail), including mean and maximum sizes and reflectivity, and thus an overall improvement in the representation of ice-phase particles in the P3 scheme.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Resilience of the Electric Grid Through Trustable IoT-Coordinated Assets

The electricity grid has evolved from a physical system to a cyberphysical system with digital devices that perform measurement, control, communication, computation, and actuation. The increased penetration of distributed energy resources (DERs) including renewable generation, flexible loads, and storage provides extraordinary opportunities for improvements in efficiency and sustainability. However, they can introduce new vulnerabilities in the form of cyberattacks, which can cause significant challenges in ensuring grid resilience. We propose a framework in this paper for achieving grid resilience through suitably coordinated assets including a network of Internet of Things devices. A local electricity market is proposed to identify trustable assets and carry out this coordination. Situational Awareness (SA) of locally available DERs with the ability to inject power or reduce consumption is enabled by the market, together with a monitoring procedure for their trustability and commitment. With this SA, we show that a variety of cyberattacks can be mitigated using local trustable resources without stressing the bulk grid. Multiple demonstrations are carried out using a high-fidelity cosimulation platform, real-time hardware-in-the-loop validation, and a utility-friendly simulator.

distributed energy resources↗
DOE OSTI2026

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Comparison and analysis of aircraft measurements and mesoscale atmospheric chemistry model simulations of tropospheric ozone

The Regional Acid Deposition Model (RADM) has been applied to several of the field experiments which were part of the Acid Models Operational and Diagnostic Evaluation Study (Acid MODES). The experiment which was of particular interest with regards to ozone photochemistry involved horizontal zig-zag flight patterns (ZIPPER) over an area from the eastern Ohio River valley to the Adirondacks of New York. Model simulations by both the standard resolution RADM (delta x = 80 km) and the nested grid RADM (delta x = 26.7 km) compare well to measurements in the low emission regions in central Pennsylvania and upstate New York, but underestimate in the high emission upper Ohio River valley. The nested simulation does considerably better, however, than the coarse grid simulation in terms of horizontal pattern and concentration magnitudes. Analysis of NO(x) and HO(x) concentrations and photochemical products rates of ozone show that the model's response to large point source emissions is very unsystematic both spatially and temporally. This is due to the models instability to realistically simulate the small scale (subgrid) gradients in precursor concentrations in and around large point source plumes.

Pleim, Jonathan E.↗
DOE OSTI2010

NON-LOCALITY OF HYDRODYNAMIC AND MAGNETOHYDRODYNAMIC TURBULENCE

We compare non-locality of interactions between different scales in hydrodynamic (HD) turbulence and magnetohydrodynamic (MHD) turbulence in a strongly magnetized medium. We use three-dimensional incompressible direct numerical simulations to evaluate non-locality of interactions. Our results show that non-locality in MHD turbulence is much more pronounced than that in HD turbulence. Roughly speaking, non-local interactions count for more than 10% of total interactions in our MHD simulation on a grid of 512{sup 3} points. However, there is no evidence that non-local interactions are important in our HD simulation with the same numerical resolution. We briefly discuss how non-locality affects the energy spectrum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Software Requirements for a Consequence Management Sample Data Simulator for Training and Drills

This document describes the requirements for a software tool that will enable FRMAC to simulate large sets of sample result data that is based realistically on simulated radionuclide deposition grids from NARAC. The user of this tool would be scientists involved in exercise and drill planning or part of the simulation cell of an exercise controller team. A key requirement is that this tool must be able to be run with a reasonable amount of training and job aids by any person within the Assessment, Laboratory Analysis, or Monitoring and Sampling divisions of the FRMAC to support any level of exercise from the small IPX to the national level full scale exercise. This tool should be relatively lean and stand-alone so that the user can run it in the field with limited IT resources. This document will describe the desired architecture, design characteristics, order of operations, and algorithms that can be given to a software development team to assist them in project scoping, costing, and eventually, development.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1988

Three-dimensional elliptic grid generation for an F-16

A case history depicting the effort to generate a computational grid for the simulation of transonic flow about an F-16 aircraft at realistic flight conditions is presented. The flow solver for which this grid is designed is a zonal one, using the Reynolds averaged Navier-Stokes equations near the surface of the aircraft, and the Euler equations in regions removed from the aircraft. A body conforming global grid, suitable for the Euler equation, is first generated using 3-D Poisson equations having inhomogeneous terms modeled after the 2-D GRAPE code. Regions of the global grid are then designated for zonal refinement as appropriate to accurately model the flow physics. Grid spacing suitable for solution of the Navier-Stokes equations is generated in the refinement zones by simple subdivision of the given coarse grid intervals. That grid generation project is described, with particular emphasis on the global coarse grid.

Reese L. Sorenson↗
DOE OSTI2019

Seminar 13 - Advanced Methods for Grid Integration of High-Performance Residential Communities

High-performance residential communities such as net-zero energy (NZE) and NZE-ready communities are emerging across the United States as many cities and states have committed to achieve ambitious energy efficiency and renewable energy goals between now and 2050. The achievement of these goals will transform the nation's communities into the holistic paradigms consisting of pervasive solar photovoltaics (PVs) and other renewable distributed energy resources, energy storage and flexible loads. However, integrating residential communities with high PV penetration may bring many challenges to the grid. For example, during many hours of the year, electric power is flowing out of PV-enabled homes and back into the grid. When this happens simultaneously in many nearby homes, distribution feeder voltage can be raised significantly. Intermittent solar power can also cause voltage flickers, leading to potential damages on both distribution feeder and home appliances. Existing solutions, such as PV curtailment, battery storage, and direct load control, are not able to address all the challenges in a cost-effective way. Better understanding is needed on how to deploy rooftop PV, how to predict building loads, and how to coordinate the homes in a community to enhance grid reliability. This seminar consists of four integral presentations that cover different aspects of a solution to these challenges. First, we will present an automated method to identify the best surfaces to deploy PV and size PV according to the PV penetration goals in a community. Second, we will present a streamlined way to create gray-box building models that are suitable for large-scale building-grid co-simulation. Next, we will present a data-driven method to estimate the end-use energy consumption to better coordinate the behind-the-meter resources to reduce utility cost, provide grid services, and improve renewable energy integration. Lastly, we will present a hierarchical control system that integrates the building models and the estimated non-dispatchable loads into the decision-making process to enhance the reliability and resilience of the residential community. Simulation results from an NZE-ready community in Fort Collins, Colorado will be presented in all four presentations and the guidelines on better grid integration of residential communities will be provided.

14 SOLAR ENERGY↗
NASA NTRS1991

Monte Carlo simulation of ion-neutral charge exchange collisions and grid erosion in an ion thruster

A combined particle-in-cell (PIC)/Monte Carlo simulation model has been developed in which the PIC method is used to simulate the charge exchange collisions. It is noted that a number of features were reproduced correctly by this code, but that its assumption of two-dimensional axisymmetry for a single set of grid apertures precluded the reproduction of the most characteristic feature of actual test data; namely, the concentrated grid erosion at the geometric center of the hexagonal aperture array. The first results of a three-dimensional code, which takes into account the hexagonal symmetry of the grid, are presented. It is shown that, with this code, the experimentally observed erosion patterns are reproduced correctly, demonstrating explicitly the concentration of sputtering between apertures.

Peng, Xiaohang↗
DOE OSTI2019

Assessing Adaptive Mesh Refinement (AMR) in a Forced Shallow-Water Model with Moisture

Two forced shallow-water flow scenarios are explored in a 2D fourth-order finite-volume dynamical core with adaptive mesh refinement (AMR) to investigate AMR’s ability to track and resolve complex evolving features. Traditional shallow-water test cases are mainly characterized by large-scale smooth flows that do not effectively test the multiscale abilities of variable-resolution and AMR models to resolve sharp gradients and small-scale flow filaments. Therefore, adding forcing mechanisms to the shallow-water system to model key atmospheric processes adds complexity and creates small-scale phenomena. These can serve as foci for dynamic grid refinement while remaining simple enough to study the numerical design of a model’s dynamical core. The first shallow-water flow scenario represents a strengthening, tropical cyclone–like, vortex that is driven by a Betts–Miller-like convection scheme. The second shallow-water test is built upon a barotropically unstable jet with an added Kessler-like warm rain scheme that leads to precipitating frontal zones. The key feature of both tests is that there is significant sensitivity to the model grid while converging (structurally) at high resolution. Both test cases are investigated for a series of uniform resolutions and a variety of AMR tagging criteria. The AMR simulations demonstrate that grid refinement can resolve local features without requiring global high-resolution meshes. However, the results are sensitive to the refinement criteria. Criteria that trigger refinement early in a simulation reproduce the uniform-resolution reference solutions most reliably. Finally, in contrast, AMR criteria that delay refinement for several days require careful tuning of the AMR thresholds to improve results compared with uniform-resolution simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Probabilistic Risk Assessment of CO2 Trapping Mechanisms in a Sandstone CO2-EOR Field in Northern Texas, USA

This study employs a polynomial chaos expansion (PCE) technique to develop the reduced order models (ROMs) integrated with the Monte Carlo simulations to quantify uncertainties of CO 2 storage associated with enhanced oil recovery (EOR) in a commercial-scale CO 2 -EOR field. The Morrow B sandstone at the Farnsworth Unit (FWU) EOR field in northern Texas was selected as a case study. A 3-D reservoir model, including simulation of a five-spot injection/production well pattern, was constructed based on the data from well logs and cores, and an updated seismic interpretation in the FWU. Porosity and permeability were considered as sources of uncertainty. The 1000-realization heterogeneous random fields of porosity and permeability were generated based on the variogram analysis on the core measurements of the FWU. Three model outputs of interest include the amount of CO 2 trapped by three trapping mechanisms, i.e. hydrodynamic trapping, oil dissolution trapping, and aqueous dissolution trapping. The total of 25 runs were conducted with the simulations of a 20-year CO 2 -EOR period, and a 20-year monitoring period to develop the ROMs between uncertain model input parameters (i.e., porosity and permeability in this study) and model outputs at the end of simulation for each grid cell using PCE approach. The high coefficient of determination (R 2 ) and small normalized root mean square error (NRMSE) between the simulated (using reservoir modeling) and predicted (using ROMs) CO 2 storage indicate that the ROMs are acceptable and reliable to be used for the predictions. Given the 1000 Monte Carlo samplings of the model input parameters, cumulative distribution functions (CDFs) and uncertainty bounds (10 th and 90 th percentiles) of model outputs were estimated based on ROMs. At the end of simulation, the most injected CO 2 at the FWU was stored by hydrodynamic trapping (between 121,400 and 166,860 tonnes), followed by oil dissolution trapping (between 55,303 and 76,816 tonnes), and aqueous dissolution trapping (between 1979 and 2751 tonnes).

58 GEOSCIENCES↗
DOE OSTI2022

The impact of forest-controlled snow variability on late-season streamflow varies by climatic region and forest structure

Previous studies have documented how forests influence snow at fine spatial scales, but none have documented the influence that existing forest-snow variability has on streamflow. Here to test how much forest-controlled snow variability influences streamflow, a tiling parameterization based on classifications from high-resolution (1–3 m) vegetation maps was incorporated into the Distributed Hydrology and Soil Vegetation Model (DHSVM). Within each grid cell (90–150 m), the tiling parameterization simulated forest-snow variability with four independently evolving snowpacks. Each tile had unique radiation conditions to represent conditions underneath the canopy, in exposed areas, and along north- and south-facing forest edges. This tiled parameterization was used to test where and when detailed forest-snow modelling should be considered further and where and when the impacts are too small to be worth the effort. To test this, tiled model simulations of streamflow were compared to non-tiled model simulations in the Sierra Nevada, CA, the Jemez Mountains, NM, and the Eastern Cascades, WA. In Tuolumne, CA, the tiled model simulated little difference in grid cell average SWE, and late-season streamflow decreased by only 3%–4% compared to the non-tiled model. In Jemez, NM, the tiled model decreased late-season streamflow by 18% due to increased sublimation. In Chiwawa, WA, the tiled model increased late-season streamflow by 15% due to high shortwave radiation attenuation and less longwave radiation enhancement from the forest. Furthermore, within the Chiwawa, a substantial silvicultural practice was synthetically implemented to increase the north-facing edge's fractional area. This silvicultural experiment, which used the same fractional forest area in all simulations increased late-season streamflow by 35% compared to tiled model simulations that did not represent forest edges. In conclusion, representing forest-SWE variability had an effect on late-season streamflow in some watersheds but not in others based on the fractional area of the forest edges, forest characteristics, and climate conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

New Grid and Discretization Technologies for Ocean and Ice Simulations (Final Report)

The work reported is in pursuit of these goals: high-quality unstructured, non-uniform Voronoi and Delaunay grids; improved finite element and finite volume discretization schemes; and improved finite element and finite volume discretization schemes. These are sought for application to spherical and three-dimensional applications suitable for ocean, atmosphere, ice-sheet, and other climate modeling applications.

58 GEOSCIENCES↗