Search NASA⌕ Search

SEARCH · Search NASA

Results for “initial conditions”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Investigating the use of field solvers for simulating classical systems

We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the nonlinear regime. We also investigate extensions of the Schrödinger-Poisson system that employ multiple stacked cold streams, and the von Neumann–Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. Finally, we present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Anti-Ultralocality and Plateau Models of Inflation

Anti-ultralocality refers to the growth of spatial gradient terms relative to velocity terms in the coupled Einstein--scalar field equations. It is a characteristic feature of decelerated expansion before the onset of inflation. Previous numerical relativity studies have shown that anti-ultralocality prevents the onset of inflation in models with power-law inflaton potentials. In this paper, we show that models with plateau-shaped inflaton potentials, which are considered to be the simplest way to generate a tensor-to-scalar ratio below current observational upper limits, are especially vulnerable to anti-ultralocality effects. The reasons are the flatness of the plateau and the energy density gap of $\sim 10$ orders of magnitude between the Planck density and the plateau potential energy. To study the problem, we develop a protocol for assessing the viability of inflationary models in general, and we apply it to a plateau potential using a previously validated numerical relativity code. We find that, starting from generic initial conditions, the growth of gradient terms in the Einstein equations relative to non-gradient terms either prevents inflation from lasting for enough $e$-folds or triggers a phase of quantum runaway. We show that the fine-tuning of initial conditions necessary to avoid these issues becomes more severe as the energy scale of inflation is made smaller, disfavoring common approaches for reducing the tensor-to-scalar ratio.

FOS: Physical sciences↗
DOE OSTI2022

Freeze-in, glaciation, and UV sensitivity from light mediators

Dark matter (DM) freeze-in through a light mediator is an appealing model with excellent detection prospects at current and future experiments. Light mediator freeze-in is UV-insensitive insofar as most DM is produced at late times, and thus the DM abundance does not depend on the unknown early evolution of our universe. However the final DM yield retains a dependence on the initial DM population, which is usually assumed to be exactly zero. We point out that in models with light mediators, the final DM yield will also depend on the initial conditions assumed for the light mediator population. We describe a class of scenarios we call “glaciation” where DM freezing in from the SM encounters a pre-existing thermal bath of mediators, and study the dependence of the final DM yield on the initial temperature of this dark radiation bath. To compute DM scattering rates in this cosmology, we derive for the first time an exact integral expression for the Boltzmann collision term describing interactions between two species at different temperatures. We quantify the dependence of the DM yield on the initial dark temperature and find that it can be sizeable in regions near the traditional (zero initial abundance) freeze-in curve. We generalize the freeze-in curve to a glaciation band, which can extend as much as an order of magnitude below the traditional freeze-in direct detection target, and point out that the DM phase space distribution as well as the yield can be strongly dependent on initial conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Study of The Compression of an Imploding Plasma and the Magnetic- Field Flux Using Advanced, Non-Intrusive, Spectroscopic Techniques

The main motivation of the project was to measure the magnetic field in pulsed power driven plasma compression using advanced non-intrusive spectroscopic techniques. The strength of the project is that it encompasses a combination of laboratory experiments with diagnostics for both the magnetic field and plasma properties and computational modeling. The following issues were addressed: (1) Plasma and magnetic field dynamics, (2) Field compression and diffusion, (3) Energy balance in the system, (4) Onset of instabilities, and (5) Correlation of the plasma and field evolutions with the initial conditions. The unique feature of the research was the use of advanced spectroscopy to measure both the magnetic field and various plasma properties in multiple locations using both side-on and end-on views, simultaneously. The experimental system allowed for imploding different gas mixtures (which facilitated local measurements) and for flexibility in modifying the initial conditions. Two current generators were used with significantly different current rise times (200 ns and 1.6 µs) but at the same current level of 400 kA. The slower pulsed power generator (at rise time 1.6 µs) was based at the Weizmann Institute of Science (WIS) and imploded plasma into which an axial magnetic field (up to 1 T) was initially embedded. The faster current rise time generator (at rise time 200 ns) based at University of California San Diego (UCSD) was used for plasma parameter studies initially, with collaborative efforts we successfully implemented the advanced spectroscopy techniques at UC San Diego. One graduate student (Nick Aybar) and a postdoctoral scholar (Maylis Dozeries) were trained. They both are based at the Lawrence Livermore National Laboratory now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Extremum seeking for optimal control problems with unknown time-varying systems and unknown objective functions

We consider the problem of optimal feedback control of an unknown, noisy, time-varying, dynamic system that is initialized repeatedly. Examples include a robotic manipulator which must perform the same motion, such as assisting a human, repeatedly and accelerating cavities in particle accelerators which are turned on for a fraction of a second with given initial conditions and vary slowly due to temperature fluctuations. In this paper, we present an approach that applies to systems of practical interest. The method presented here is model independent; does not require knowledge of the objective function; is robust to measurement noise; is applicable for any set of initial conditions; is applicable to simultaneously controlling an arbitrary number of parameters; and may be implemented with a broad range of continuous or discontinuous functions such as sine or square waves. For systems with convex cost functions we prove that our algorithm will produce controllers that approach the minimal cost. For linear systems we reproduce the cost minimizing linear quadratic regulator optimal controller that could have been designed analytically had the system and cost function been known. We demonstrate the effectiveness of the algorithm with simulation studies of noisy and time-varying systems.

42 ENGINEERING↗
DOE OSTI2022

Preliminary Thermal-hydraulics and Selected Safety Analysis for Holos-Quad Reactor Design

This report presents thermal hydraulics and safety analyses of the 10 MWe Holos-Quad micro-reactor design developed by HolosGen LLC. These analyses were executed under steady-state and station black-out (SBO) transient conditions using the System Analysis Module (SAM), with a focus on ensuring peak fuel temperatures do not exceed the safety thresholds of 1,250°C for steady state normal operating conditions and 1,600°C for transient conditions. For steady-state temperature predictions, a coupled 1D fluid-to-3D solid heat conduction methodology was used. The steady-state simulations included a model of the central core region where power peaking occurs to predict the maximum steady-state fuel temperature, and a coarse-mesh full-core model to provide the initial condition for the consequent transient SBO decay heat removal simulation. The maximum fuel temperature predicted from the central core region simulations under the normal operating condition was 1,222°C, which is below the design limit 1,250°C. To simulate a SBO transient, a full core model was necessary because of the non-axial symmetric core design and heat loss from the outer core structures to the environment. To reduce the required computational costs, a novel scaled subassembly approach was used to model the full core domain. With the full core steady-state temperature distribution provided as the initial condition, the passive decay heat removal of the core was simulated with the scaled subassembly full core model under SBO transient conditions. This analysis showed that the peak fuel temperatures of the Holos-Quad core remained below their steady-state values during the transient, therefore much lower than the safety limit 1,600°C. This is mainly due to the large thermal inertia of the graphite matrix, relatively low power density, and the large surface-to-volume-ratio of the core. The preliminary analyses presented in this report confirm the inherent safety of the Holos-Quad micro-reactor thermal-hydraulics design, as the peak fuel temperatures under both normal and off-normal operational conditions remain within design and safety limits.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Feasibility of DEIM for retrieving the initial field via dimensionality reduction

When parameter estimation is solved in a high-dimensional space, the dimensionality reduction strategy becomes the primary consideration for alleviating the tremendous computational cost. Here, the discrete empirical interpolation method (DEIM) is explored to retrieve the initial condition (IC) by combining the polynomial chaos (PC) based ensemble Kalman filter (i.e. PC-EnKF), where a non-intrusive PC expansion is considered as a surrogate model in place of the forward model in the prediction step of the ensemble Kalman filter, resulting in fewer forward model integrations but with a comparable accuracy as Monte Carlo-based approaches. The DEIM acts as a hyper-reduction tool to provide the low-dimensional input for the high-dimensional initial field, which can be reconstructed using the information on the sparse interpolation grid points that is adaptively obtained through PC-EnKF data assimilation method. Thus an innovative framework to reconstruct the IC is developed. The detailed procedure at each assimilation iteration includes: the determination of the spatial interpolation points, the estimation of the initial values on the interpolation locations using the optimal observations, and the reconstruction of IC in the full space. The current study uses the reconstruction field of initial conditions of the Navier-Stokes equations as an example to illustrate the efficacy of our method. The experimental results demonstrate the proposed algorithm achieves a satisfactory reconstruction for the initial field. The proposed method helps to extend the applicable area of DEIM in solving inverse problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Beware of fake ν ’s: The effect of massive neutrinos on the nonlinear evolution of cosmic structure

Massive neutrinos suppress the growth of cosmic structure on small, nonlinear scales. It is thus often proposed that using statistics beyond the power spectrum can tighten constraints on the neutrino mass by extracting additional information from these nonlinear scales. We study the information content regarding neutrino mass at the field level, quantifying how much of this information arises from the difference in nonlinear evolution between a cosmology with one fluid [cold dark matter (CDM)] and two fluids ( CDM + neutrinos ). We do so by running two N -body simulations, one with and one without massive neutrinos, both with the same phases, and matching their linear power spectrum at a given low redshift. This effectively isolates the information encoded in the linear initial conditions from the nonlinear cosmic evolution. We demonstrate that, for k ≲ 1 h / Mpc , and for a single redshift, there is negligible difference in the real-space CDM field between the two simulations. This suggests that all the information regarding neutrino mass is in the linear power spectrum set by the initial conditions. Thus, any probe based on the CDM field alone will have negligible constraining power beyond that which exists at the linear level over the same range of scales. Consequently, any probe based on the halo field will contain little information beyond the linear power. We find similar results for the matter field responsible for weak lensing. We also demonstrate that there may be much information beyond the power spectrum in the 3D matter field; however, this is not observable in modern surveys via dark matter halos or weak lensing. Finally, we show that there is additional information to be found in redshift space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Non-linear flow modes of identified particles in Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 5.02 TeV

The pT-differential non-linear flow modes, v 4,22 , v 5,32 , v 6,33 and v 6,222 for π ± , K ± , K$^0_S$ , $p + \bar{p}, Λ + \bar{Λ}$ and ϕ-meson have been measured for the first time at $\sqrt{s_{NN}}$ = 5.02 TeV in Pb-Pb collisions with the ALICE detector at the Large Hadron Collider. The results were obtained with a multi-particle technique, correlating the identified hadrons with reference charged particles from a different pseudorapidity region. These non-linear observables probe the contribution from the second and third order initial spatial anisotropy coefficients to higher flow harmonics. All the characteristic features observed in previous p T -differential anisotropic flow measurements for various particle species are also present in the non-linear flow modes, i.e. increase of magnitude with increasing centrality percentile, mass ordering at low p T and particle type grouping in the intermediate p T range. Hydrodynamical calculations (iEBE-VISHNU) that use different initial conditions and values of shear and bulk viscosity to entropy density ratios are confronted with the data at low transverse momenta. These calculations exhibit a better agreement with the anisotropic flow coefficients than the non-linear flow modes. These observations indicate that non-linear flow modes can provide additional discriminatory power in the study of initial conditions as well as new stringent constraints to hydrodynamical calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Impacts of Degradation on Water, Energy, and Carbon Cycling of the Amazon Tropical Forests

Selective logging, fragmentation, and understory fires directly degrade forest structure and composition. However, studies addressing the effects of forest degradation on carbon, water, and energy cycles are scarce. Here, we integrate field observations and high-resolution remote sensing from airborne lidar to provide realistic initial conditions to the Ecosystem Demography Model (ED-2.2) and investigate how disturbances from forest degradation affect gross primary production (GPP), evapotranspiration (ET), and sensible heat flux (H). We used forest structural information retrieved from airborne lidar samples (13,500 ha) and calibrated with 817 inventory plots (0.25 ha) across precipitation and degradation gradients in the eastern Amazon as initial conditions to ED-2.2 model. Our results show that the magnitude and seasonality of fluxes were modulated by changes in forest structure caused by degradation. During the dry season and under typical conditions, severely degraded forests (biomass loss ≥66%) experienced water stress with declines in ET (up to 34%) and GPP (up to 35%) and increases of H (up to 43%) and daily mean ground temperatures (up to 6.5°C) relative to intact forests. In contrast, the relative impact of forest degradation on energy, water, and carbon cycles markedly diminishes under extreme, multiyear droughts, as a consequence of severe stress experienced by intact forests. Our results highlight that the water and energy cycles in the Amazon are driven by not only climate and deforestation but also the past disturbance and changes of forest structure from degradation, suggesting a much broader influence of human land use activities on the tropical ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Empirically verified analysis of dual pre-cooling system for hydrogen refuelling station

The dual cooling system at the California State University Los Angeles Hydrogen Research and Fueling Facility is analysed to predict the output hydrogen temperature in different ambient conditions. Initially, the facility was built with a coil chiller T20 cooling, sufficient for occasional fuelings. A flat plate evaporator system was added in series to resolve any issues for multiple fuelings even in the hottest weather. This study was commenced to verify the set points and determine the control requirements for the dual setup in order to inform future designs. A numerical-empirical model of the cooling system was developed by reverse engineering, using the station experimental measurements and database records. The subsystems are characterized separately but validated as a whole system imposing specific and standard operating conditions. The cooling model is tested with different inputs to validate the pre-set parameter variation of the cooling components. Ten scenarios were investigated to evaluate the fueling output parameters, testing two hydrogen temperature input trends (according to the facility initial conditions) and five ambient temperatures (from 15 °C to 35 °C). In closing, the simulation results have confirmed that the hydrogen output temperature remains within the range imposed by the SAE J2601 for each scenario studied, around -22 °C and -25 °C.

08 HYDROGEN↗
DOE OSTI2023

A greedy Galerkin method to efficiently select sensors for linear dynamical systems

A key challenge in inverse problems is the selection of sensors to gather the most effective data. In this paper, we consider the problem of inferring the initial condition to a linear dynamical system and develop an efficient control-theoretical approach for greedily selecting sensors. Our method employs a Galerkin projection to reduce the size of the inverse problem, resulting in a computationally efficient algorithm for sensor selection. As a byproduct of our algorithm, we obtain a preconditioner for the inverse problem that enables the rapid recovery of the initial condition. Here, we analyze the theoretical performance of our greedy sensor selection algorithm as well as the performance of the associated preconditioner. Finally, we verify our theoretical results on various inverse problems involving partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Fate of supernova progenitors in massive binary systems

ABSTRACT How massive stars end their lives depends on the core mass, core angular momentum, and hydrogen envelopes at death. However, these key physical facets of stellar evolution can be severely affected by binary interactions. In turn, the effectiveness of binary interactions itself varies greatly depending on the initial conditions of the binaries, making the situation much more complex. We investigate systematically how binary interactions influence core–collapse progenitors and their fates. Binary evolution simulations are performed to survey the parameter space of supernova progenitors in solar metallicity binary systems and to delineate major evolutionary paths. We first study fixed binary mass ratios ($q=M_2/M_1$ = 0.5, 0.7, and 0.9) to elucidate the impacts of initial mass and initial separation on the outcomes, treating separately Type Ibc supernova, Type II supernova, accretion-induced collapse (AIC), rapidly rotating supernova (Ibc-R), black hole formation, and long gamma ray burst (long GRB). We then conduct 12 binary population synthesis model calculations, varying the initial condition distributions and binary evolution parameters, to estimate various supernova fractions. We obtain a Milky Way supernova rate $R_{\rm SN} = (1.78$–$2.47) \times 10^{-2} \, {\rm yr}^{-1}$ which is consistent with observations. We find the rates of AIC, Ibc-R, and long GRB to be $\sim 1/100$ the rate of regular supernovae. Our estimated long GRB rates are higher than the observed long GRB rate and close to the low luminosity GRB rate, although care must be taken considering our models are computed with solar metallicity. Furthering binary modelling and improving the inputs one by one will enable more detailed studies of these and other transients associated with massive stars.

Astronomy & Astrophysics↗
DOE OSTI2024

New horizons in the holographic conformal phase transition

Abstract We describe 5D dynamical cosmological solutions of the stabilized holographic dilaton and their role in completion of the conformal phase transition. This analysis corresponds, via the AdS/CFT dictionary, to a study of out-of-equilibrium dynamics where trajectories of the dilaton do not depend solely on thermodynamic quantities in the early universe, but have sensitivity also to initial conditions. Unlike the well-studied thermal transition, which requires quantum tunneling of an infrared brane through the surface of an AdS-Schwarzschild horizon, our approach instead invokes an early epoch in which the cosmology is fully 5-dimensional, with highly relativistic brane motion and with Rindler horizons obscuring the infrared brane at early times. In this context, we demonstrate the existence of a large class of natural initial conditions that seed trajectories where the brane simply passes through the Rindler horizon and into the basin of attraction of the stabilized dilaton potential. This corresponds to successful completion of the phase transition without sacrificing perturbativity of the 5D theory.

Physics↗
DOE OSTI2021

Probabilistic error estimation for non-intrusive reduced models learned from data of systems governed by linear parabolic partial differential equations

This work derives a residual-based a posteriori error estimator for reduced models learned with non-intrusive model reduction from data of high-dimensional systems governed by linear parabolic partial differential equations with control inputs. It is shown that quantities that are necessary for the error estimator can be either obtained exactly as the solutions of least-squares problems in a non-intrusive way from data such as initial conditions, control inputs, and high-dimensional solution trajectories or bounded in a probabilistic sense. Here, the computational procedure follows an offline/online decomposition. In the offline (training) phase, the high-dimensional system is judiciously solved in a black-box fashion to generate data and to set up the error estimator. In the online phase, the estimator is used to bound the error of the reduced-model predictions for new initial conditions and new control inputs without recourse to the high-dimensional system. Numerical results demonstrate the workflow of the proposed approach from data to reduced models to certified predictions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Skeletal Kinetics Reduction for Astrophysical Reaction Networks

A novel methodology is developed to extract accurate skeletal reaction models for nuclear combustion. Local sensitivities of isotope mass fractions with respect to reaction rates are modeled based on the forced optimally time-dependent (f-OTD) scheme. These sensitivities are then analyzed temporally to generate skeletal models. The methodology is demonstrated by conducting skeletal reduction of constant density and temperature burning of carbon and oxygen relevant to Type Ia supernovae (SNe Ia). The 495-isotopes Torch model is chosen as the detailed reaction network. A map of maximum production of 56 Ni in SNe Ia is produced for different temperatures, densities, and proton-to-neutron ratios. The f-OTD simulations and the sensitivity analyses are then performed with initial conditions from this map. A series of skeletal models are derived and their performances are assessed by comparison against currently existing skeletal models. Previous models have been constructed intuitively by assuming the dominance of α-chain reactions. The comparison of the newly generated skeletal models against previous models is based on the predicted energy release and 44 Ti and 56 Ni abundances by each model. The consequences of ye ≠ 0.5 in the initial composition are also explored where ye is the electron fraction. The simulated results show that 56 Ni production decreases by decreasing ye as expected, and that the 43 Sc is a key isotope in proton and neutron channels toward 56 Ni production. It is shown that an f-OTD skeletal model with 150 isotopes can accurately predict the 56Ni abundance in SNe Ia for ye ≲ 0.5 initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗