Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Consider Parameter”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Cross-machine comparison of runaway electron generation during tokamak start-up for extrapolation to ITER

A cross-machine comparison of global parameters that determine the runaway electron (RE) generation and loss process during tokamak start-up was carried out with the aim to extrapolate these to ITER. The study found that all considered discharges, also those that do not show signs of RE, are non-thermal at the start, i.e. have a streaming parameter larger than 0.1. During the current ramp-up the electric field, E, remains above the critical value, E c , that allows RE in the plasma. The distinction to be made is not if RE can form but, if sufficient RE can form fast enough such that they are detected or start to dominate the dynamics of the tokamak discharge. The dynamics of the value of E, density and temperature during tokamak are key to the formation of RE. It was found that larger devices operate with E closer to E c , due to their higher temperatures, hence the RE generation is relatively slower. The slower time scales for the formation of RE, estimated to be of the order of 100s of ms in ITER simplifies the development of avoidance schemes. The RE confinement time is also an important determinant of the entire process and is found to increase with the device size. The study also revealed that drift orbit losses, a mechanism often attributed as the main RE loss mechanism during the early tokamak discharge, are actually more difficult to achieve. RE losses might be more likely attributed to RE diffusion due to magnetic turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet

We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
DOE OSTI2026

Understanding the Influence of Chain Architecture on the Transport Quantities of Polymer Electrolytes with Covalently Bonded Anions

Here, we use a combination of experiments and coarse-grained molecular dynamics simulations to elucidate the structure–property relationships in polymer electrolytes obtained by the copolymerization of poly(vinyl ethylene carbonate─lithium styrene bis(trifluoromethanesulfonyl)imide) or p(VEC-LiSTFSI). Experiments show that the conductivity reduces with increasing anion (i.e., STFSI) fraction on the chain, and the cation transference number (t + ) is found to be dependent on the anion fraction. Furthermore, a significant fraction of unpolymerized VEC monomers are observed. Since it is inherently difficult to experimentally control the chain architecture and the amount of unpolymerized VEC in these systems, we perform coarse-grained molecular dynamics simulations on model polymer systems with different chain architectures to mimic the plausible experimental systems. Specifically, we look at the differences in transference numbers arising from (i) a random copolymer of VEC and STFSI monomers; (ii) a blend of VEC-STFSI copolymer with VEC monomers; and (iii) a ternary blend of the VEC homopolymer, STFSI homopolymers, and VEC monomers. The ternary blend model demonstrates the closest resemblance with the experimental transference numbers and diffusivities. The lithium diffusivity obtained from the coarse-grained models with VEC monomers (plasticizers) is about 1.5 times that of the model without VEC monomers, showing that the plasticizing effect of VEC monomers is modest. We rationalize the experimental observations based on aggregate and cluster analyses obtained from molecular simulations. This work reveals that polymer electrolyte chain architecture and plasticizers can critically influence the transport properties, and these parameters should be considered when designing single ion conducting polymeric electrolytes.

cluster distribution↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gravitational wave probes of massive gauge bosons at the cosmological collider

We extend the reach of the "cosmological collider" for massive gauge boson production during inflation from the CMB scales to the interferometer scales. Considering a Chern-Simons coupling between the gauge bosons and the pseudoscalar inflaton, one of the transverse gauge modes is efficiently produced and its inverse decay leaves an imprint in the primordial scalar and tensor perturbations. We study the correlation functions of these perturbations and derive the updated constraints on the parameter space from CMB observables. We then extrapolate the tensor power spectrum to smaller scales consistently taking into account the impact of the gauge field on inflationary dynamics. Our results show that the presence of massive gauge fields during inflation can be detected from characteristic gravitational wave signals encompassing the whole range of current and planned interferometers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

PIC Simulations of One-side Multipactor on Dielectric

Breakdown of dielectric RF windows is an important issue for particle accelerators and high-power RF sources. One of the generally considered reasons for the RF windows failure is the multipactor on dielectric surface. The multipactor may be responsible for excessive heating of dielectric and discharge of charges that accumulated in ceramic due to secondary emission. In this study the comprehensive self-consistent PIC simulations with space charge effect were performed in order to better understand the dynamic of one-side multipactor development and floating potential on dielectric induced by the emission. The important correlations between the multipactor parameters at saturation and the secondary emission properties of dielectric and the applied RF field parameters were found and are reported in the paper.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Chaotic saddles and interior crises in a dissipative nontwist system

Here, we consider a dissipative version of the standard nontwist map. Nontwist systems present a robust transport barrier, called the shearless curve, that becomes the shearless attractor when dissipation is introduced. This attractor can be regular or chaotic depending on the control parameters. Chaotic attractors can undergo sudden and qualitative changes as a parameter is varied. These changes are called crises, and at an interior crisis the attractor suddenly expands. Chaotic saddles are nonattracting chaotic sets that play a fundamental role in the dynamics of nonlinear systems; they are responsible for chaotic transients, fractal basin boundaries, and chaotic scattering, and they mediate interior crises. In this work we discuss the creation of chaotic saddles in a dissipative nontwist system and the interior crises they generate. We show how the presence of two saddles increases the transient times and we analyze the phenomenon of crisis induced intermittency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluation of Converter Performance Considering Static and Dynamic Device Part-to-Part Variability

This paper presents a methodology to incorporate and analyze the impact of semiconductor device part-to-part variation on power converter performance. By integrating extensive static and dynamic device characterization data with an automated compact model generation process that reflects manufacturing variability, device models with inherent variability features are utilized in converter simulations for a comprehensive assessment of performance impacts. The traditional converter performance evaluation process typically yields fixed efficiency values, often dismissing the inherent part-to-part variability caused by the manufacturing process of semiconductor devices. To address this limitation, a large population of devices was characterized to capture variations in static parameters-such as transfer, output, and capacitance characteristics-as well as dynamic behaviors, including switching losses. This data-driven approach enables the development of individual compact models, which were then integrated into converter simulations to evaluate efficiency ranges rather than single point estimated values. The converter simulation results show that part-to-part component variation can lead to significant efficiency deviations, exceeding several percentage points in high-power conversion applications. By offering a more accurate representation of converter behavior under real-world manufacturing conditions, this methodology enables designers to anticipate performance variability, improving the robustness of power converter designs.

device characterization↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The Study on the Shock Wave Propagation Rule of a Gas Explosion in a Gas Compartment

Combined with the k-ε turbulence model of general application, a refined finite element model of a utility tunnel’s gas compartment filled with the methane/air mixture is developed. A series of analyses are made by using the powerful industry-leading computational fluid dynamics (CFD) software flame acceleration simulator (FLACS) to study the shock wave propagation rule in the gas compartment. The longitudinal and transversal distribution laws of the explosion shock wave are gained taking into consideration the spatial characteristics of the gas compartment. The influences of a few parameters, such as initial conditions and section size of the gas compartment, on the shock wave propagation rule are further discussed. The basic procedure for predicting the peak pressure of the blast wave is provided by considering the initial conditions and the gas compartment, and the corresponding injury effect of the explosion wave on the living beings is assessed. The investigation demonstrates that the peak pressure by the coupled effect between the initial conditions is significantly influenced, especially at the upper and lower gas explosion limits. The peak pressure increases gradually as the width or height increases, and both basically meet the linear relation. The proposed method can forecast the peak pressure of the explosion shock wave in the gas compartment accurately. According to the peak pressure longitudinal and transversal distributions of the blast wave, the peak pressure is far greater than the killing pressure threshold in the underground and closed space; consequently, it is not safe for the living beings in the gas compartment.

Lin, Shu-Chao↗
DOE OSTI2023

Electron kinetics in a high- Z plasmoid

The problem of the electron dynamics on a closed magnetic field line passing through a high-Z plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the H theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mechanism-Based Energy Regularization in Computational Modeling of Quasibrittle Fracture

Abstract Quasibrittle materials are featured by a strain-softening constitutive behavior under many loading scenarios, which could eventually lead to localization instability. It has long been known that strain localization would result in spurious mesh sensitivity in finite element (FE) simulations. Previous studies have shown that, for the case of fully localized damage, the mesh sensitivity can be mitigated through energy regularization of the material constitutive law. However, depending on the loading configuration and structural geometry, quasibrittle structures could exhibit a complex damage process, which involves both localized and diffused damage patterns at different stages of loading. This study presents a generalized energy regularization method that considers the spatial and temporal evolution of damage pattern. The method introduces a localization parameter, which describes the local damage pattern. The localization parameter governs the energy regularization of the constitutive model, which captures the transition from diffused to localized damage during the failure process. The method is cast into an isotropic damage model, and is further extended to rate-dependent behavior. The energy regularization scheme is directly incorporated into the kinetics of damage growth. The model is applied to simulate static and dynamic failures of ceramic specimens. It is shown that the present model is able to effectively mitigate the spurious mesh sensitivity in FE simulations of both types of failure. The present analysis demonstrates the essential role of mechanism-based energy regularization of constitutive relation in FE simulations of quasibrittle fracture.

Mechanics↗
DOE OSTI2023

Thermal acclimation of plant photosynthesis and autotrophic respiration in a northern peatland

Peatlands contain one-third of global soil carbon (C), but the responses of peatland ecosystems to long-term warming are not well understood. Here, we pursue an emergent understanding of warming effects on ecosystem C fluxes at peatlands by constraining a process-oriented model, the terrestrial ECOsystem model, with observational data from a long-term warming experiment at the Spruce and Peatland Responses Under Changing Environments site. Model-based assessments show that ecosystem-level photosynthesis and autotrophic respiration exhibited significant thermal acclimation, with temperature sensitivities being linearly decreased with warming. Using the thermal-acclimated parameter values, simulated gross primary production, net primary production, and plant autotrophic respiration (R a ), were all lower than those simulated with non-thermal acclimated parameter values. In contrast, ecosystem respiration simulated with thermal acclimated parameter values was higher than that simulated with non-thermal acclimated parameter values. Net ecosystem CO 2 exchange was much higher after constraining model parameters with observational data from the warming treatments, releasing C at a rate of 28.3 g C m -2 yr -1 °C -1 . Our data-model integration study suggests that peatlands are likely to release more C than previously estimated. Earth system models may overestimate C uptake by peatlands under warming if physiological thermal acclimation of plants is not incorporated. Thus, it is critical to consider the long-term physiological thermal acclimation of plants in the models to better predict global C dynamics under future climate and their feedback to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ice-shelf ocean boundary layer dynamics from large-eddy simulations

Abstract. Small-scale turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small-scale processes, which regulate heat and salt transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

DynaShape: Dynamic Shaping of Grid Response with Inverters

We consider the problem of controlling the frequency of low-inertia power systems via inverter-based resources (IBRs) that are weakly connected to the grid. We propose a novel grid-forming control strategy, the so-called frequency shaping control, that aims to shape the frequency response of synchronous generators (SGs) to load perturbations so as to efficiently arrest sudden frequency drops. Our solution relaxes several existing assumptions in the literature and is able to navigate trade-offs between peak power requirements and maximum frequency deviations. Finally, we analyze the robustness to imperfect knowledge of network parameters, while particularly highlighting the importance of accurate estimation of these parameters.

frequency-shaping control↗
DOE OSTI2021

Experiment and computational fluid dynamics investigation of biochar elutriation in fluidized bed

Here, in fluidized bed biomass fast pyrolysis, the biomass is converted to biochar and elutriated. The elutriation rate is a key parameter in reactor designs and operations. This research presents a video-based continuous measurement of biochar elutriation rate in a fluidized bed with sands and biomass as bed materials. The fluidized bed is simulated with the Computational Fluid Dynamics - Coarse-Grained Discrete Element Method (CFD-CGDEM) in MFiX. The fluidization behavior of non-spherical sands can be more accurately captured when a rolling friction model is used. The predicted elutriation rate is close to the experimental measurement when the particle size distributions are considered and the filtered drag with a shape correction is used. These results validated the accuracy of the MFiX based CFD framework for the prediction of biochar elutriations in the fluidized bed biomass fast pyrolysis reactor.

09 BIOMASS FUELS↗
DOE OSTI2023

Switching of Hybrid Improper Ferroelectricity in Oxide Double Perovskites

In ABO 3 -type perovskite oxides with Pnma symmetry, rotation (Q R+ , a 0 a 0 c + ) and tilt (Q T , a – a – c 0 ) of BO 6 octahedra are the two primary order parameters. These order parameters establish an inherent trilinear coupling with anti-ferroelectric A-site displacement (Q AFE ) to form the low-symmetry phase. The symmetry is further lowered in double perovskite oxides (DPOs) due to A/A' cation ordering. It in turn makes these systems polar via hybrid improper ferroelectric mechanism, primarily driven by Q R+ and Q T . Naturally, it has been believed that functionalities such as polarization can also be switched by tuning these primary order parameters. However, mystery around finding switching mechanism still remains. Our study based on density functional theory calculations combined with finite-temperature molecular dynamics simulations shows that the polarization switching is a two-step process, driven by out-of-phase rotation (Q R– , a 0 a 0 c – when Q T = 0 or, a – a – b – when Q T ≠ 0). A series of polar DPOs such as KLnFeOsO 6 [Ln = Sm, Gd, Dy, Tm (lanthanides) and Y (rare earth)], all belonging to P2 1 symmetry, are considered in this investigation. The polarization switching P ($\overrightarrow{P}$) occurs at a very high temperature of ~1150 K through a phase transition, from a polar (P2 1 ) phase with $\overrightarrow{P}$(+) to $\overrightarrow{P}$(-) via a non-polar P4/n phase. The switching itself is metastable in nature. The switching (both polarization and spin state) is only observed for a very short period of time (~23 ps) that poses limitation on using such a mechanism in memory device realization. We demonstrate a concurrent heating–cooling procedure to overcome such shortcoming. In conclusion, simulations conducted at 600 K further imply that long lasting switching can be achieved, at least for 1.2 ns for 600 K, and ideally for an infinite time, if the material is heated just above the T c followed by rapid cooling to a temperature below T c .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗