Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 1,063 records · Page 59

Composite Analysis Vadose Zone Transport of Selected Hanford Site Composite Analysis: Vadose Zone Transport of Selected Radionuclides Released from DOE O 435.1 Waste Sites

The main objectives of the composite analysis (CA) compliance case vadose zone (VZ) facet are to simulate transport of radionuclide releases from the surface to the water table, and to provide radionuclide transfer rates to groundwater for inclusion into the Plateau-to-River (P2R) saturated zone (SZ) model (CP-57037, Model Package Report: Plateau to River Groundwater Model, Version 8.3). The parallel exascale version of the Subsurface Transport Over Multiple Phases (eSTOMP) simulator is used to simulate flow and transport for the CA VZ models. The CA compliance case approach to address the scale and distribution of contaminant sources in the VZ facet is to subdivide the Hanford Site Central Plateau into 26 models that contain contaminant sources and liquid discharges likely to commingle during migration through the VZ to the water table.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Wake steering of multirotor wind turbines

Abstract In this paper, wake steering is applied to multirotor turbines to determine whether it has the potential to reduce wind plant wake losses. Through application of rotor yaw to multirotor turbines, a new degree of freedom is introduced to wind farm control such that wakes can be expanded, channelled or redirected to improve inflow conditions for downstream turbines. Five different yaw configurations are investigated (including a baseline case) by employing large‐eddy simulations (LES) to generate a detailed representation of the velocity field downwind of a multirotor wind turbine. Two lower‐fidelity models from single‐rotor yaw studies (curled‐wake model and analytical Gaussian wake model) are extended to the multirotor case, and their results are compared with the LES data. For each model, the wake is analysed primarily by examining wake cross‐sections at different downwind distances. Further quantitative analysis is carried out through characterisations of wake centroids and widths over a range of streamwise locations and through a brief analysis of power production. Most significantly, it is shown that rotor yaw can have a considerable impact on both the distribution and magnitude of the wake velocity deficit, leading to power gains for downstream turbines. The lower‐fidelity models show small deviation from the LES results for specific configurations; however, both are able to reasonably capture the wake trends over a large streamwise range.

17 WIND ENERGY↗

The Simplified Aircraft-Based Paired Approach With the ALAS Alerting Algorithm

This paper presents the results of an investigation of a proposed concept for closely spaced parallel runways called the Simplified Aircraft-based Paired Approach (SAPA). This procedure depends upon a new alerting algorithm called the Adjacent Landing Alerting System (ALAS). This study used both low fidelity and high fidelity simulations to validate the SAPA procedure and test the performance of the new alerting algorithm. The low fidelity simulation enabled a determination of minimum approach distance for the worst case over millions of scenarios. The high fidelity simulation enabled an accurate determination of timings and minimum approach distance in the presence of realistic trajectories, communication latencies, and total system error for 108 test cases. The SAPA procedure and the ALAS alerting algorithm were applied to the 750-ft parallel spacing (e.g., SFO 28L/28R) approach problem. With the SAPA procedure as defined in this paper, this study concludes that a 750-ft application does not appear to be feasible, but preliminary results for 1000-ft parallel runways look promising.

Perry, Raleigh B.↗

Particle‐In‐Cell Simulations of Starfish Prime

The Starfish Prime high altitude nuclear test created a transient diamagnetic cavity in the Earth's magnetic field above Johnston Island and launched an electromagnetic pulse (EMP) that was detected around the globe. For this work, we use the ion‐kinetic particle‐in‐cell code Topanga to simulate diamagnetic cavity evolution and the E3 EMP signal for Starfish Prime out to over a minute of physical time. The simulation domain has a longitudinal and latitudinal extent of 60° X 60° and a vertical extent of 2,000 km from the surface of the Earth. We compare our simulated results to magnetometer measurements taken in space and on the ground, finding good agreement in both cases. The diamagnetic cavity in the simulation forms in about a second, while the associated debris flux tube takes approximately 30 s to decay. The debris flux tube undergoes significant motion during this timeframe, rising upward. The measured E3 EMP signal on the ground consists of several components, all of which are present in our simulations. We discuss the physical origin of these components in relation to E3a (blast) and E3b (heave).

Belyaev, Mikhail A. [Lawrence Livermore National L↗

Requirements for Large Eddy Simulation Computations of Variable-Speed Power Turbine Flows

Variable-speed power turbines (VSPTs) operate at low Reynolds numbers and with a wide range of incidence angles. Transition, separation, and the relevant physics leading to them are important to VSPT flow. Higher fidelity tools such as large eddy simulation (LES) may be needed to resolve the flow features necessary for accurate predictive capability and design of such turbines. A survey conducted for this report explores the requirements for such computations. The survey is limited to the simulation of two-dimensional flow cases and endwalls are not included. It suggests that a grid resolution necessary for this type of simulation to accurately represent the physics may be of the order of Delta(x)+=45, Delta(x)+ =2 and Delta(z)+=17. Various subgrid-scale (SGS) models have been used and except for the Smagorinsky model, all seem to perform well and in some instances the simulations worked well without SGS modeling. A method of specifying the inlet conditions such as synthetic eddy modeling (SEM) is necessary to correctly represent the inlet conditions.

Grid↗

Direct Simulation of a Self-Similar Turbulent Mixing Layer

Three direct numerical simulations of incompressible turbulent plane mixing layers have been performed. All the simulations were initialized with the same two velocity fields obtained from a direct numerical simulation of a turbulent boundary layer with a momentum thickness Reynolds number of 300 computed by Spalart (J. Fluid Mech. 187, 61, 1988). In addition to a baseline case with no additional disturbances, two simulations were begun with two-dimensional disturbances of varying strength in addition to the boundary layer turbulence. After a development stage, the baseline case and the case with weaker additional two-dimensional disturbances evolve self-similarly, reaching visual thickness Reynolds numbers of up to 20 000. This self-similar period is characterized by a lack of large-scale organized pairings, a lack of streamwise vortices in the 'braid' regions, and scalar mixing that is characterized by 'marching' Probability Density Functions (PDFs). The case begun with strong additional two-dimensional disturbances only becomes approximately self-similar, but exhibits sustained organized large-scale pairings, clearly defined braid regions with streamwise vortices that span them, and scalar PDFs that are 'nonmarching.' It is also characterized by much more intense vertical velocity fluctuations than the other two cases. The statistics and structures in several experiments involving turbulent mixing layers are in better agreement with those of the simulations that do not exhibit organized pairings.

Rogers, Michael M.↗

Kinetic modeling of hot tail runaway electron generation during plasma disruptions using the JOREK code

The generation of runaway electrons (REs) during disruptions poses a significant challenge for the operation of tokamaks. The production of these high-energy electrons can cause substantial damage, particularly when the plasma current is high, making it a critical concern for ITER. For the high-temperature plasmas anticipated in ITER, the primary generation of REs may be dominated by the hot tail mechanism, which consists of the acceleration of hot electrons from the pre-disruption population which have not yet thermalized with the bulk following the rapid cooling of the plasma. To account for the significant 3D effects on RE production, a hot tail modeling framework has been developed within the non-linear 3D extended MHD code JOREK. This paper presents the structure of this framework, which is based on test electrons evolving in MHD fields. The verification of the method shows good agreement with the reference DREAM code for 0D test cases, as well as for axisymmetric simulations of 15 MA ITER H-mode disruption scenarios. Furthermore, a proof-of-principle application to a DIII-D case demonstrates the framework’s capability to capture for the first time the hot tail generation in 3D MHD simulations in realistic geometry. Preliminary results suggest that the production of REs is significantly reduced by stochastic losses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modelling of electrostatic ion-scale turbulence in divertor tokamaks with the gyrokinetic code $\mathrm{COGENT}$

Here, continuum gyrokinetic simulations of electrostatic ion scale turbulence are presented for the case of a diverted (single-null) tokamak geometry. The simulation model, implemented in the finite-volume code COGENT, solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to a vorticity equation for electrostatic potential variations, where a fluid model is used for an electron response. The model describes the ion scale ion temperature gradient (ITG) and resistive drift modes as well as neoclassical ion physics effects. Different turbulence regimes are observed depending on the plasma profiles, and the roles of a self-consistent background electric field and an X-point geometry are explored. In particular, increasing the pedestal density gradient and the corresponding radial electric field is demonstrated to suppress the ITG turbulence, whereas the same edge plasma background can still be destabilized by the resistive modes when the plasma resistivity is increased. The effects of X-point geometry are assessed by comparing cross-separatrix simulations with counterpart calculations performed for a toroidal annulus geometry. For the simulation parameters considered, similar global behaviour is observed in both cases, whereas strong local suppression of turbulence fluctuations is demonstrated near the X-point for the case of a single-null geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of Modeled to Measured Spectra using MCNP and GADRAS to Benchmark and Contrast Modeling Limitations

The desire to improve nuclear material detection through portal monitors and other low resolution detectors has led to interest in benchmarking the performance of radiation transport codes. These codes can be used to generate a variety and quantity of spectra that may be cost prohibitive to measure directly. Particular characteristics of typical detection scenarios were isolated to compare the performance of radiation transport codes to laboratory experiments. These benchmark experiments were performed to validate simulations in a number of key areas. These experiments included high areal density configurations in both one- and three-dimensional configurations. In addition, off-angle measurements were conducted to investigate the impact of detector response assumptions on complex geometries, for example where an unknown source is not physically collocated with an apparent hot spot. Furthermore, by examining both simple shielding and backscatter configurations, isolation of the scatter emanating from the object and environmental background scatter contributions to the spectra could be elucidated. A series of experiments were performed in which the order of shielding layers and the thickness were evaluated to examine the ability of the simulations to address these variables. Experiments with increasing thickness of polyethylene around a neutron source also allows for examinations of the ability of the transport simulations to properly model n,γ reactions. Finally, depleted uranium (DU) measurements have more complicated spectra than the cobalt-60 and consequently allow for investigations into the ability of transport simulations to model self-shielding. In each case, the results of the experiments were compared to MCNP simulations to judge the performance of the code, as was done previously with GADRAS 18.7.9 simulations [7]. All of the benchmark measurements were taken with a liquid nitrogen cooled 140% high purity germanium (HPGe) detector with a bismuth collimator that had the front tin filter removed. The measurement location and detector configuration were subsequently used for all experimental configurations. The remainder of this report details the experimental measurements, Section 2 and the MCNP model, Section 3. Comparisons between the experimental results and MCNP-generated spectra are performed, Section 4. Section 5 compares the results obtained with MCNP with the previously reported results obtained with GADRAS [2]. Finally, Section 6 presents the conclusions.

61 RADIATION PROTECTION AND DOSIMETRY↗

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Universal energy-speed-accuracy trade-offs in driven nonequilibrium systems

The connection between measure theoretic optimal transport and dissipative nonequilibrium dynamics provides a language for quantifying nonequilibrium control costs, leading to a collection of thermodynamic speed limits, which rely on the assumption that the target probability distribution is perfectly realized. This is almost never the case in experiments or numerical simulations, so here we address the situation in which the external controller is imperfect. We obtain a lower bound for the dissipated work in generic nonequilibrium control problems that (1) is asymptotically tight and (2) matches the thermodynamic speed limit in the case of optimal driving. Along with analytically solvable examples, we refine this imperfect driving notion to systems in which the controlled degrees of freedom are slow relative to the nonequilibrium relaxation rate, and identify independent energy contributions from fast and slow degrees of freedom. Furthermore, we develop a strategy for optimizing minimally dissipative protocols based on optimal transport flow matching, a generative machine learning technique. Furthermore, this latter approach ensures the scalability of both the theoretical and computational framework we put forth. Crucially, we demonstrate that we can compute the terms in our bound numerically using efficient algorithms from the computational optimal transport literature and that the protocols we learn saturate the bound.

59 BASIC BIOLOGICAL SCIENCES↗

An Experimental and Computational Study of Directional Solidification in Transparent Materials

An experimental and numerical study of the horizontal Bridgman growth of pure succinonitrile (SCN) and of a succinonitrile- 1.0 mol.% acetone alloy (SCN- 1.0 mol.% ACE) has been performed. Experiments at growth rates of 0, 2 and 40 micron/s were investigated. The solid/liquid interface was stable (non-dendritic and non-cellular); however, it was not flat. Rather, it was significantly distorted by the influence of convection in the melt and, for the growth cases, by the moving temperature boundary conditions along the ampoule. For the alloy, the interface was.determined to be unstable at growth rates greater than 2.8 micron/s, but stable for the no-growth and 2 micron/s growth cases. When compared to the pure SCN interface, the alloy interface forms closer to the cold zone, indicating that the melting temperature has been suppressed by the addition of the alloying element. Two-dimensional computer simulations were performed for the no-growth case for both the pure and alloy materials. These simulations indicate that a primary longitudinal convective cell is formed in the melt. The maximum magnitude of velocity was calculated to be 1.515 mm/s for pure SCN and 1.724 mm/s for the alloy. The interface shape predicted by the computer simulation agrees well with the experimentally determined shape for the pure SCN case. In ongoing work, numerical simulations of the process during growth conditions are being performed.

Simpson, James E.↗

CFD Approaches for Simulation of Wing-Body Stage Separation

A collection of computational fluid dynamics tools and techniques are being developed and tested for application to stage separation and abort simulation for next-generation launch vehicles. In this work, an overset grid Navier-Stokes flow solver has been enhanced and demonstrated on a matrix of proximity cases and on a dynamic separation simulation of a belly-to-belly wing-body configuration. Steady cases show excellent agreement between Navier-Stokes results, Cartesian grid Euler solutions, and wind tunnel data at Mach 3. Good agreement has been obtained between Navier-Stokes, Euler, and wind tunnel results at Mach 6. An analysis of a dynamic separation at Mach 3 demonstrates that unsteady aerodynamic effects are not important for this scenario. Results provide an illustration of the relative applicability of Euler and Navier-Stokes methods to these types of problems.

Buning, Pieter G.↗

Synthetic Scenarios from CMIP5 Model Simulations for Climate Change Impact Assessments in Managed Ecosystems and Water Resources: Case Study in South Asian Countries

Increasing population, urbanization, and associated demand for food production compounded by climate change and variability have important implications for the managed ecosystems and water resources of a region. This is particularly true for south Asia, which supports one quarter of the global population, half of whom live below the poverty line. This region is largely dependent on monsoon precipitation for water. Given the limited resources of the developing countries in this region, the objective of our study was to empirically explore climate change in south Asia up to the year 2099 using monthly simulations from 35 global climate models (GCMs) participating in the fifth phase of the Climate Model Inter-comparison Project (CMIP5) for two future emission scenarios (representative concentration pathways RCP4.5 and RCP8.5) and provide a wide range of potential climate change outcomes. This was carried out using a three-step procedure: calculating the mean annual, monsoon, and non-monsoon precipitation and temperatures; estimating the percent change from historical conditions; and developing scenario funnels and synthetic scenarios. This methodology was applied for the entire south Asia region; however, the percent change information generated at 1.5deg grid scale can be used to generate scenarios at finer spatial scales. Our results showed a high variability in the future change in precipitation (-23% to 52%, maximum in the non-monsoon season) and temperature (0.8% to 2.1%) in the region. Temperatures in the region consistently increased, especially in the Himalayan region, which could have impacts including a faster retreat of glaciers and increased floods. It could also change rivers from perennial to seasonal, leading to significant challenges in water management. Increasing temperatures could further stress groundwater reservoirs, leading to withdrawal rates that become even more unsustainable. The high precipitation variability (with higher propensity for localized intense rainfall events) observed in the region can be a key factor for managed ecosystems and water management and could also lead to more incidence of severe urban flooding. The results could be used to assess both mitigation and adaptation alternatives to reduce vulnerabilities in managed ecosystems (agricultural and urban) and water resources.

agriculture↗

Multiclustering needlet ILC for CMB B -mode component separation

ABSTRACT The Cosmic Microwave Background (CMB) primordial B-mode signal is predicted to be much lower than the polarized Galactic emission (foregrounds) in any region of the sky pointing to the need for sophisticated component separation methods. Among them, the blind Needlet Internal Linear Combination (NILC) has great relevance given our current poor knowledge of the B-mode foregrounds. However, the expected level of spatial variability of the foreground spectral properties complicates the NILC subtraction of the Galactic contamination. We therefore propose a novel extension of the NILC approach, the Multiclustering NILC (MC-NILC), which performs NILC variance minimization on separate regions of the sky (clusters) properly chosen to have similar spectral properties of the B-mode Galactic emission within them. Clusters are identified thresholding either the ratio of simulated foregrounds-only B modes (ideal case) or the one of cleaned templates of Galactic emission obtained from realistic simulations. In this work we present an application of MC-NILC to the future LiteBIRD satellite, which targets the observation of both reionization and recombination peaks of the primordial B-mode angular power spectrum with a total error on the tensor-to-scalar ratio δr < 0.001. We show that MC-NILC provides a CMB solution with residual foreground and noise contamination that is significantly lower than the NILC one and the primordial signal targeted by LiteBIRD at all angular scales for the ideal case and at the reionization peak for a realistic ratio. Thus, MC-NILC will represent a powerful method to mitigate B-mode foregrounds for future CMB polarization experiments.

Astronomy & Astrophysics↗

Direct simulation Monte Carlo prediction of on-orbit contaminant deposit levels for HALOE

A three-dimensional version of the direct simulation Monte Carlo method is adapted to assess the contamination environment surrounding a highly detailed model of the Upper Atmosphere Research Satellite. Emphasis is placed on simulating a realistic, worst-case set of flow field and surface conditions and geometric orientations for the satellite in order to estimate an upper limit for the cumulative level of volatile organic molecular deposits at the aperture of the Halogen Occultation Experiment. A detailed description of the adaptation of this solution method to the study of the satellite's environment is also presented. Results pertaining to the satellite's environment are presented regarding contaminant cloud structure, cloud composition, and statistics of simulated molecules impinging on the target surface, along with data related to code performance. Using procedures developed in standard contamination analyses, along with many worst-case assumptions, the cumulative upper-limit level of volatile organic deposits on HALOE's aperture over the instrument's 35-month nominal data collection period is estimated at about 13,350 A.

Woronowicz, Michael S.↗