Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity Simulation”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mean Climate and Tropical Rainfall Variability in Aquaplanet Simulations Using the Model for Prediction Across Scales‐Atmosphere

Abstract Aquaplanet experiments are important tools for understanding and improving physical processes simulated by global models; yet, previous aquaplanet experiments largely differ in their representation of subseasonal tropical rainfall variability. This study presents results from aquaplanet experiments produced with the Model for Prediction Across Scales‐Atmosphere (MPAS‐A)—a community model specifically designed to study weather and climate in a common framework. The mean climate and tropical rainfall variability simulated by MPAS‐A with varying horizontal resolution were compared against results from a recent suite of aquaplanet experiments. This comparison shows that, regardless of horizontal resolution, MPAS‐A produces the expected mean climate of an aquaplanet framework with zonally symmetric but meridionally varying sea‐surface temperature. MPAS‐A, however, has a stronger signal of tropical rainfall variability driven by convectively coupled equatorial waves. Sensitivity experiments with different cumulus parameterizations, physics packages, and vertical grids consistently show the presence of those waves, especially equatorial Kelvin waves, in phase with lower‐tropospheric convergence. Other models do not capture such rainfall‐kinematics phasing. These results suggest that simulated tropical rainfall variability depends not only on the cumulus parameterization (as suggested by previous studies) but also on the coupling between physics and dynamics of climate and weather prediction models.

Rios‐Berrios, R.↗
DOE OSTI2024

Modeling a divertor with mid-leg pumping for high-power H-mode scenarios in DIII-D considering E x B drift flows

Edge-plasma simulations of a baffled, long-legged divertor in DIII-D, performed using the multi-fluid code UEDGE, indicate that the position of the detachment front is constrained to the location of the pump duct along the low-field side (LFS) baffle. Simulations including magnetic and E x B drifts were performed for 12.5 MW deuterium plasmas including intrinsic carbon and seeded neon to assess the optimal location of the LFS divertor pump to create a stable detachment front between the target and the X-point. The radiation front position in the simulations, taken to be indicative of the detachment front, can be controlled between the pump and X-point in the favorable magnetic field direction for H-mode access by moving the pump duct location upstream of the target along the LFS baffle. In the unfavorable magnetic field direction, the radial E θ x B drift flows are directed towards the pumping surface, efficiently removing the injected deuterium gas and limiting the sensitivity of the radiation front location to the gas injection rate. The role of pumping rate and drift direction on the pumping efficiency are also found to affect the divertor plasma conditions and detachment front location in UEDGE simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Numerical Analysis of Fuel Effects on Advanced Compression Ignition Using a Cooperative Fuel Research Engine Computational Fluid Dynamics Model

Growing environmental concerns and demand for a better fuel economy are driving forces that motivate the research for more advanced engines. Multi-mode combustion strategies have gained attention for their potential to provide high thermal efficiency and low emissions for light-duty applications. These strategies target optimizing the engine performance by correlating different combustion modes to load operating conditions. The extension from boosted spark ignition (SI) mode at high loads to advanced compression ignition (ACI) mode at low loads can be achieved by increasing the compression ratio and utilizing intake air heating. Further, in order to enable an accurate control of intake charge condition for ACI mode and rapid mode-switches, it is essential to gain fundamental insights into the autoignition process. Within the scope of ACI, homogeneous charge compression ignition (HCCI) mode is of significant interest. It is known for its potential benefits, operation at low fuel consumption, low NOx, and particulate matter (PM) emissions. In the present work, a virtual Cooperative Fuel Research (CFR) engine model is used to analyze fuel effects on ACI combustion. In particular, the effect of fuel octane sensitivity (S) (at constant Research Octane Number (RON)) on autoignition propensity is assessed under beyond-RON (BRON) and beyond-MON (BMON) ACI conditions. The three-dimensional CFR engine computational fluid dynamics (CFD) model employs a finite-rate chemistry approach with a multi-zone binning strategy to capture autoignition. Two binary blends with Research Octane Number (RON) of 90 are chosen for this study: primary reference fuel (PRF) with S = 0 and toluene–heptane (TH) blend with S = 10.8, representing paraffinic and aromatic gasoline surrogates. Reduced mechanisms for these blends are generated from a detailed gasoline surrogate kinetic mechanism. Simulation results with the reduced mechanisms are validated against experimental data from an in-house CFR engine, with respect to in-cylinder pressure, heat release rate, and combustion phasing. Thereafter, the sensitivity of combustion behavior to ACI operating condition (BRON versus BMON), air-fuel ratio (λ = 2 and 3), and engine speed (600 and 900 rpm) is analyzed for both fuels. It is shown that the sensitivity of a fuel’s autoignition characteristics to λ and engine speed significantly differs at BRON and BMON conditions. Moreover, this sensitivity is found to vary among fuels, despite the same RON. It is also observed that the presence of low-temperature heat release (LTHR) under BRON condition leads to more sequential autoignition and longer combustion duration than BMON condition. Finally, the study indicates that the octane index (OI) fails to capture the trend in the variation of autoignition propensity with S under the BMON condition.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Development of a U-19Pu-10Zr fuel performance benchmark case based on the IFR-1 experiment

Metallic nuclear fuels are subject to research and development for use in advanced reactors. Robust, accurate metallic fuel performance models are important for the design, safety analysis, and licensing of these reactors. However, metallic fuel performance models require additional development; they are not as mature as UO 2 fuel performance models. A benchmark case based on the IFR-1 experiment was developed to better gauge the accuracy of existing models, identify models for high-priority development, and potentially quantify any future improvements made by further model development.This work collected publicly available information on the IFR-1 experiment and used it to develop the benchmark case. Fuel behavior during the IFR-1 irradiation was simulated by using the fuel performance code BISON, and the predicted results were compared with postirradiation examination data from the IFR-1 experiment. Furthermore, a sensitivity study and tuning studies were performed as a preliminary investigation into the causes of inaccurate temperature and dimensional change predictions.The benchmark predicted reasonably accurate values for the burnup and fission gas release. There was error in the predicted temperatures, which could be explained by uncertainty in the input parameters and legacy temperatures. BISON over-predicted dimensional changes in the fuel and cladding. The sensitivity study showed that the dimensional changes were most sensitive to the fuel swelling anisotropy and the cladding void swelling model. Future benchmark and model development should focus on cladding swelling behaviors to improve dimensional change predictions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Search for CME in U+U and Au+Au collisions in STAR with different approaches of handling backgrounds

We report the chiral magnetic effect (CME) refers to charge separation along a strong magnetic field between left- and right-handed quarks, caused by interactions with topological gluon fields from QCD vacuum fluctuations. We present two approaches to handle the dominant elliptic flow (v 2 ) background in the three-particle correlator (Δγ 112 ), sensitive to CME. In the first approach, we present the Δγ 112 and Δγ 123 measurements in U+U and Au+Au collisions. While hydrodynamic simulations including resonance decays and local charge conservation predict that Δγ 112 scaled by N part /ν 2 will be similar in U+U and Au+Au collisions, the projected B-field exhibits a distinct difference between the two systems and with varying N part . Therefore, U+U and Au+Au collisions provide configurations with different expectations for both CME signal and background. Moreover, the three-particle observable Δγ 123 scaled by N part /ν 3 provide baseline measurement for only the background. In the second approach, we handle the ν 2 background by measuring Δγ 112 with respect to the planes of spectators measured by Zero Degree Calorimeters and participants measured by Time Projection Chamber. These measurements contain different amounts of contributions from CME signal (along B-field, due to spectators) and ν 2 background (determined by the participant geometry). With the two Δγ 112 measurements, the possible CME signal and the background contribution can be determined. We report such a measurement in Au+Au collisions at $\sqrt{s_{NN}}$ = 27GeV with the newly installed event plane detector, and report the new findings in U+U system where the spectator-participant plane correlations are expected to differ from those in Au+Au collisions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Temperature Trends in the Northwestern Tibetan Plateau Constrained by Ice Core Water Isotopes Over the Past 7,000 Years

Abstract The reasons for the Holocene temperature conundrum, known as the inconsistency between the reconstructed cooling and the inferred warming simulations during the Holocene, remain unclear. Temperature reconstructions from the Tibetan Plateau (TP) provide important insights for understanding the Holocene temperature conundrum due to enhanced sensitivity to climate at high altitudes. Given the significant positive correlation between air temperature and δ 18 O in precipitation over the northern TP, the stable isotopic records of ice cores recovered from this area are widely used for paleotemperature reconstruction. Here we present a new high‐resolution δ 18 O record from the Chongce ice cores to bedrock, dated back to 7 ka BP by the accelerator mass spectrometry (AMS) 14 C dating technique. Our reconstructed temperature record shows a long‐term warming trend until ~2 ka BP, followed by an abrupt change to a relatively cool period until the start of the industrial‐era warming. This record challenges the widely recognized Holocene reconstruction from the neighboring Guliya ice core. It is also different from many previous temperature reconstructions, most of which have summer biases and show a long‐term cooling trend over the past two millennia. In addition, our record shows that temperatures during the recent decades are almost the highest during the past 7 ka BP, highlighting the unusual warming forced by anthropogenic greenhouse gases.

Pang, Hongxi↗
DOE OSTI2022

Importance of the Antarctic Slope Current in the Southern Ocean Response to Ice Sheet Melt and Wind Stress Change

Abstract We use two coupled climate models, GFDL‐CM4 and GFDL‐ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre‐industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near‐surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large‐scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean‐state representation and meltwater‐driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.

Beadling, R. L.↗
DOE OSTI2022

A Landau–Devonshire analysis of strain effects on ferroelectric Al 1-x Sc x N

We present a thermodynamic analysis of the recently discovered nitride ferroelectric materials using the classic Landau–Devonshire approach. Electrostrictive and dielectric stiffness coefficients of Al 1-x Sc x N with a wurtzite structure (6 mm) are determined using a free energy density function assuming a hexagonal parent phase (6/mmm), with the first-order phase transition based on the dielectric stiffness relationships. The results of this analysis show that the strain sensitivity of the energy barrier is one order of magnitude larger than that of the spontaneous polarization in these wurtzite ferroelectrics, yet both are less sensitive to strain compared to classic perovskite ferroelectrics. These analysis results reported here explain experimentally reported sensitivity of the coercive field to elastic strain/stress in Al 1-x Sc x N films and would enable further thermodynamic analysis via phase field simulation and related methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Real-space simulations of the hyperbolic plasmon polaritons in 1T′ tungsten ditelluride

Recent discoveries of hyperbolic surface plasmon polaritons (SPPs) in 1T′ WTe 2 have garnered significant attention in 2D materials and nanophotonics. In this study, we employ finite-element simulations to investigate the real-space characteristics of hyperbolic SPPs in thin WTe 2 flakes. Our results show that the SPPs exhibit a pronounced sensitivity to excitation energy and sample thickness. By analyzing the plasmonic field patterns, we extract key plasmonic parameters including plasmon wavelengths, hyperbolic angles, and plasmonic figures of merit. In addition, we examine SPP modes in stacked WTe 2 flakes with varying twist angles, demonstrating that the plasmonic field patterns can be effectively tuned by adjusting the twist angle. Notably, as the twist angle increases, the SPPs undergo a topological transition from open hyperbolic modes to closed elliptic modes. This twist-engineering capability offers promising potential for the development of tunable plasmonic devices based on WTe 2 .

2D materials↗
DOE OSTI2024

Study of synchronous RF pulsing in dual frequency capacitively coupled plasma

Low-pressure multi-frequency capacitively coupled plasmas (CCPs) are used for numerous etch and deposition applications in the semiconductor industry. Pulsing of the radio-frequency (RF) sources enables control of neutral and charged species in the plasma on a millisecond timescale. The synchronous (i.e. simultaneous, in-phase) pulsing of both power sources in a dual frequency CCP is examined in this article. Due to the low gas pressure, modeling has been done using the electrostatic particle-in-cell/Monte Carlo collision method. Here, the objective of this work is to investigate the sensitivity of the plasma properties to small changes in timing during synchronous pulsing of the two RF sources. It is demonstrated that small deviations in the on and off times of the two RF sources can lead to major changes in the plasma characteristics. This high sensitivity is of concern for process repeatability but can be utilized to enable better control of the dynamics of plasma-surface interaction. In the simulations, the pulsing parameters (on and off times and ramp rates) are varied and the temporal evolution of plasma characteristics such as electron density (n e ), species current at the electrode, and electron temperature are examined. It is demonstrated that if the low-frequency (LF) source is turned off a few μs before (or after) the high-frequency source, n e during the off-state is significantly higher (or lower) due to the frequency coupling effect. Similarly, turning on the LF source with a small delay results in a sharp increase in the plasma density when the RF sources are turned on.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Optimizing detachment control using the magnetic configuration of divertors

Abstract As tokamak research moves to reactor conditions, the control of a stable, optimally-detached divertor plasma has become increasingly relevant. Simple predictions of such detachment control have been performed previously using the detachment location sensitivity (DLS) model (Lipschultz et al 2016 Nucl. Fusion 56 056007). In this study the DLS model is extended and combined with SOLPS-ITER simulations of isolated divertor grids to study the effects of alternate divertor magnetic field properties on detachment control. The DLS model predicts that divertors can achieve easier access to detachment through a long connection length, a high total flux expansion, and a high average magnetic field in the divertor compared with that at the x -point. SOLPS-ITER simulations suggest an even stronger impact of total flux expansion and connection length on detachment access than the DLS model predicts. In terms of detachment evolution, both simulation and modelling show a high gradient in the total magnetic field and low B p o l B are able to more easily keep a detachment front at a desired poloidal location. Regions with no stable detachment front locations can arise for high magnetic field gradients directed towards the target. Significant differences have been found between impurity scan simulations and the DLS model. These differences may be attributed to sources and sinks of power and electron pressure, and should be explored further in future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Redshift-space streaming velocity effects on the Lyman- α forest baryon acoustic oscillation scale

The baryon acoustic oscillation (BAO) scale acts as a standard ruler for measuring cosmological distances and has therefore emerged as a leading probe of cosmic expansion history. However, any physical effect that alters the length of the ruler can lead to a bias in our determination of distance and expansion rate. One of these physical effects is the streaming velocity, the relative velocity between baryons and dark matter in the early Universe, which couples to the BAO scale due to their common origin in acoustic waves at recombination. In this work, we investigate the impact of streaming velocity on the BAO feature of the Lyman- α forest auto-power spectrum, one of the main tracers being used by the recently commissioned Dark Energy Spectroscopic Instrument (DESI). To do this, we develop a perturbative model for Lyman- α flux fluctuations which is complete to second order for a certain set of fields, and applicable to any redshift-space tracer of structure since it is based only on symmetry considerations. We find that there are 8 biasing coefficients through second order. We find streaming velocity-induced shifts in the BAO scale of 0.081%–0.149% (transverse direction) and 0.053%–0.058% (radial direction), depending on the model for the biasing coefficients used. These are smaller than, but not negligible compared to, the DESI Lyman- α BAO error budget, which is 0.46% on the overall scale. Overall, the sensitivity of these results to our choice of bias parameters underscores the need for future work to measure the higher-order biasing coefficients from simulations, especially for future experiments beyond DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermal effects mediating the flow induced by laser-induced optical breakdown

Here, a short focused laser pulse can generate a hot plasma, which expands then contracts and can eject a hot jet, the speed and direction of which is sensitive to the details of the plasma kernel. The coupling of thermal and chemical nonequilibrium plasma mechanisms with this hydrodynamic development is assessed with detailed simulations of a two-temperature, three-species plasma model. Time scales for electron recombination, thermal relaxation, and diffusion are compared to that of the plasma expansion to anticipate conditions in which these mechanisms might affect the vorticity generation that leads to the ultimate flow pattern. The effect of these mechanisms are analyzed through comparison with corresponding inert-gas and equilibrium models. Thermal-nonequilibrium effects are found to be weak due to rapid relaxation of the heavy-particle and electronic temperatures. In contrast, chemical equilibration occurs at a rate comparable to the expansion and thereby enhance both it and subsequent hydrodynamic mechanisms as the energy stored in ion formation during the preceding breakdown is released by electron recombination. Thermal conduction, enhanced by high-temperature free electrons, weakens the ejection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probing ferroelectricity by x-ray absorption spectroscopy in molecular crystals

We carry out x-ray absorption spectroscopy experiment at the oxygen $K$ edge in croconic acid (C 5 H 2 O 5 ) crystal as a prototype of ferroelectric organic molecular solid, whose electric polarization is generated by proton transfer. The experimental spectrum is well reproduced by the electron-hole excitation theory simulations from configuration generated by ab initio molecular dynamics simulation. When inversion symmetry is broken in the ferroelectric state, the hydrogen bonding environment on the two bonded molecules become inequivalent. Such a difference is sensitively probed by the bound excitation in the pre-edge, which is strongly localized on the excited molecules. Our analysis shows that a satellite peak in the pre-edge will emerge at higher excitation energy, which serves as a clear signature of ferroelectricity in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimal Voltage Control in Low-Observable Unbalanced Distribution Systems

The increased integration of distributed energy resources (DERs) in distribution systems brings both advantages and technical challenges. High levels of DERs can often cause over/under voltage problems. Classical voltage control algorithms are based on full knowledge of voltage states across all nodes in the system. However, this may not be a practical assumption since many locations in distribution systems are unobservable. Therefore, this paper proposes a new model predictive control (MPC) based control algorithm that accounts for system unobservability to efficiently eliminate voltage violations with as low as 50% fraction of observable nodes. Additionally, an analytical voltage sensitivity framework is employed to quickly determine the change in voltage states due to PV injections. The effectiveness of the proposed method is validated via simulations on the unbalanced IEEE 37 node test system.

Abujubbeh, Mohammad↗