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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Influence of plastic properties on the grain size effect on twinning in Ti and Mg

Deformation twinning is an inevitable mode of plastic deformation in hexagonal close packed (HCP) crystals and has been considered a hinderance to achieving formability and ductility. It has long been recognized that reductions in grain size can lower the propensity for twinning. In prior work, EBSD statistical analyses on big data sets of twins in pure Ti and Mg were carried out for the same twin type, loading conditions, initial texture, and range of grain sizes. Here, we revisit these data sets to investigate the role of material plastic properties on grain size effects on deformation twinning. We show that while reductions in grain size dramatically lower the number of twins per grain and the twin thickness in both metals, these effects are several times stronger in Ti than that for Mg. To rationalize this difference, full-field, 3D crystal plasticity calculations are performed. The analysis indicates that the stronger grain size effect in Ti arises due to the larger critical resolved shear stress (CRSS) of its easiest slip system, making Ti less accommodating of the twin shear than Mg. Overall, this result points to plastic properties as additional variables that can be tuned towards developing HCP materials for structural applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fast dynamic aperture optimization with forward-reversal integration

A fast dynamic aperture (DA) optimization method for storage rings has been developed through the use of reversal integration. Even if dynamical systems have an exact reversal symmetry, a numerical forward integration differs from its reversal. For a chaotic trajectory, cumulative round-off errors are scaled, which results in an exponential growth on the difference. The exponential effect is a generic chaos indicator which represents the sensitivity of the chaotic motion to its initial condition. The chaos indicator of the charged particle motion can be obtained by comparing the forward integrations of particle trajectories with corresponding reversals, a.k.a. “backward integrations.” The indicator is observable even through short-term particle tracking simulations. Therefore, adopting it as an objective function could speed up optimization. Finally, the DA of the National Synchrotron Light Source II storage ring, and another test diffraction-limited light source ring, were optimized using this method for the purpose of demonstration.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Charge sharing in pixelated semiconductor sensors

The charge sharing between neighboring pixels in pixelated sensors can be used to measure particle or x-ray coordinates with accuracy better than the pixel pitch. The accurate model of the charge distribution shape is essential to achieve ultimate coordinate accuracy. The charge sharing is caused by charge carriers' diffusion on the path from the generation point to pixels. This paper is focused on the diffusion of the initially compact charge cloud in the field free region. The diffusion equation solutions are obtained using separation of variable and Fourier synthesis method for different initial conditions and resulting charge distributions are integrated over pixel areas. In conclusion, the look up table containing pre-calculated values for pixel charge fractions is proposed to speed up numerical calculations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Expanded analysis of machine learning models for nuclear transient identification using TPOT

Industries around the world are becoming more and more data driven. The nuclear field is no exception with several different applications being proposed. One popular area of research is the use of machine learning in transient detection. This paper seeks to build upon a previous study which made use of the AutoML package TPOT to train traditional machine learning models to classify transient events occurring with a reactor. Synthetic data was once again collected using a GPWR reactor simulator. Data on 12 different events was collected using 15 different initial conditions. Here, a dataset consisting of over 100,000 data points was compiled and used to train 7 different machine learning models using a pre-defined TPOT dictionary with 12 different preprocessing techniques. Three of the trained models were able to produce validation results in the 90s with the expanded dataset. Once the models were trained, it was possible to look into where during the simulation, misclassifications occurred. Using these three models, analysis was done to determine if TPOT could be used to train models that were effective if important features were missing. The results from this were positive with the newly trained models scoring close to the original models. Finally, to conclude this study, the three high performing models were retrained using different random states to see if there was any major variation when different states were used.

42 ENGINEERING↗
DOE OSTI2020

Helicity-dependent extension of the McLerran–Venugopalan model

We construct a generalization of the McLerran-Venugopalan (MV) model including helicity effects for a longitudinally polarized target (a proton or a large nucleus). The extended MV model can serve as the initial condition for the helicity generalization of the JIMWLK evolution equation constructed in our previous paper, as well as for calculation of helicity-dependent observables in the quasi-classical approximation to QCD.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A potential model study of the nucleon's charge and mass radius

We study the charge and mass distributions within a nucleon and compute the associated squared radii based on a potential model approach. Different constituent quark configurations such as Δ, Y, and quark-diquark are considered and compared, with model parameters calibrated by experimental measurements of the proton and neutron charge radius. The results suggest that while the charge radius is dictated by quark dynamics, the mass radius is strongly influenced by nonperturbative QCD contributions to a nucleon's mass that are not sensitive to the constituent quarks. As a result, the mass radius could become substantially different from the charge radius. The obtained nucleon mass distributions of different configurations are further used for simulations of the initial conditions in heavy ion collisions. Furthermore, the computed eccentricities ε 2 and ε 3 are found to demonstrate a considerable sensitivity to the input nucleon profiles, especially to the mass radius in the peripheral region as well as for systems with fewer participants.

Gallimore, Daniel [Indiana University, Bloomington↗
DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗
DOE OSTI2019

Parametric investigation of the transition to turbulence in Rayleigh-Taylor mixing

In this study, a database is developed consisting of highly resolved large-eddy simulations of the transition process in a Rayleigh-Taylor mixing layer. Over 300 simulations are performed with parametric variations in Grashof number, Schmidt number, Atwood number, and initial perturbation spectrum. The database is then analyzed to explore how the factors that influence the transition-to-turbulence process impact mixing evolution and growth rate prior to the development of late-time, steady-state growth. The competition between entrainment-driven mixing and diffusion-driven mixing is observed, and two different mechanisms are suggested by which the homogeneity of mixing might increase, rather than decrease, at early time. The impact of parametric variation on growth rate is additionally examined, and all variations in initial conditions are found to converge to approximately the same growth parameter around α ≈ 0.025.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A crisis for the verification and validation of turbulence simulations

Three algorithms have been proposed for solution of the Rayleigh–Taylor turbulent mixing problem. They are based upon three different physical principles governing the Euler equations for fluid flow. The principles serve to select the physically relevant solution from among many nonunique solutions. The admissibility principle is in dispute. The three different algorithms, expressing the three physical admissibility principles can be formulated in terms of the three energy dissipation rates or the entropy production rates, as selected by the size of the sub grid scale coefficients. These have maximal values or less than maximal values. The resulting solutions are markedly different. We find strong validation evidence that supports the maximum rate principle, based on a review of prior results and on new results presented here. We review experimental data used for validation and sufficient to discriminate among the three. We present a new analysis of this data. We show that the hypothesized long wave length perturbations in the initial conditions are not significant, so that validation can be based on this data in a straight forward manner. One of the algorithms is labeled direct numerical simulation, but is not, and as a consequence, the two algorithms with less than maximal SGS coefficients are variants of one another. Finally recommendations for the numerical modeling of the deflagration to detonation transition in type Ia supernova are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Advection versus diffusion in Richtmyer-Meshkov mixing

The Richtmyer-Meshkov (RM) instability is one of the most severe degradation mechanisms for inertial confinement fusion (ICF), and mitigating it has been a priority for the global ICF effort. In this Letter, the instability's ability to atomically mix is linked to its background decay of residual turbulent energy. Here, we show how recently derived inequalities from the mathematical theory of PDEs constrain the evolution. A model RM process at leading order may diffusively mix or retain imprints of its initial structures indefinitely, depending on initial conditions, and there exists a theoretical range of zero-mixing for certain values of parameters. The results may apply to other systems resembling scalar transport in decaying turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Two-particle azimuthal correlations in γp interactions using pPb collisions at s NN = 8.16 TeV

The first measurements of the Fourier coefficients ( V n Δ ) of the azimuthal distributions of charged hadrons emitted from photon-proton (γp) interactions are presented. The data are extracted from 68.8 nb −1 of ultra-peripheral proton-lead (pPb) collisions at s NN = 8.16 TeV using the CMS detector. The high energy lead ions produce a flux of photons that can interact with the oncoming proton. This γp system provides a set of unique initial conditions with multiplicity lower than in photon-lead collisions but comparable to recent electron-positron and electron-proton data. The V n Δ coefficients are presented in ranges of event multiplicity and transverse momentum ( p T ) and are compared to corresponding hadronic minimum bias pPb results. For a given multiplicity range, the mean p T of charged particles is smaller in γp than in pPb collisions. For both the γp and pPb samples, V 1 Δ is negative, V 2 Δ is positive, and V 3 Δ consistent with 0. For each multiplicity and p T range, V 2 Δ is larger for γp events. The γp data are consistent with model predictions that have no collective effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Pseudorapidity dependence of anisotropic flow and its decorrelations using long-range multiparticle correlations in Pb–Pb and Xe–Xe collisions

The pseudorapidity dependence of elliptic (v 2 ), triangular (v 3 ), and quadrangular (v 4 ) flow coefficients of charged particles measured in Pb–Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{S_{NN}}$ = 5.02 TeV and in Xe–Xe collisions at $\sqrt{S_{NN}}$ = 5.44 TeV with ALICE at the LHC are presented. The measurements are performed in the pseudorapidity range − 3.5 < η < 5 for various centrality intervals using two- and multi-particle cumulants with the subevent method. The flow probability density function (p.d.f.) is studied with the ratio of flow coefficient v 2 calculated with four- and two-particle cumulant, and suggests that the variance of flow p.d.f. is independent of pseudorapidity. The decorrelation of the flow vector in the longitudinal direction is probed using two-particle correlations. The results measured with respect to different reference regions in pseudorapidity exhibit differences, argued to be a result of saturating decorrelation effect above a certain pseudorapidity separation, in contrast to previous publications which assign this observation to non-flow effects. The results are compared to 3 + 1 dimensional hydrodynamic and the AMPT transport model calculations. Neither of the models is able to simultaneously describe the pseudorapidity dependence of measurements of anisotropic flow and its fluctuations. The results presented in this work highlight shortcomings in our current understanding of initial conditions and subsequent system expansion in the longitudinal direction. Therefore, they provide input for its improvement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Skewness and kurtosis of mean transverse momentum fluctuations at the LHC energies

The first measurements of skewness and kurtosis of mean transverse momentum ($\langle$$p_T\rangle$) fluctuations are reported in Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV, Xe–Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV and pp collisions at $\sqrt{s}$ = 5.02 TeV using the ALICE detector. The measurements are carried out as a function of system size $\langle$dN ch /dη$\rangle$$^{1/3}_{|η|<0.5}$, using charged particles with transverse momentum ($p_T$) and pseudorapidity (η), in the range 0.2<$p_T$<3.0 GeV/c and |η|<0.8, respectively. In Pb–Pb and Xe–Xe collisions, positive skewness is observed in the fluctuations of $\langle$$p_T\rangle$ for all centralities, which is significantly larger than what would be expected in the scenario of independent particle emission. This positive skewness is considered a crucial consequence of the hydrodynamic evolution of the hot and dense nuclear matter created in heavy-ion collisions. Furthermore, similar observations of positive skewness for minimum bias pp collisions are also reported here. Kurtosis of $\langle$$p_T\rangle$ fluctuations is found to be in good agreement with the kurtosis of Gaussian distribution, for most central Pb–Pb collisions. Hydrodynamic model calculations with MUSIC using Monte Carlo Glauber initial conditions are able to explain the measurements of both skewness and kurtosis qualitatively from semicentral to central collisions in Pb–Pb system. Color reconnection mechanism in PYTHIA8 model seems to play a pivotal role in capturing the qualitative behavior of the same measurements in pp collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Jet cone radius dependence of R AA and v 2 at PbPb 5.02 TeV from JEWEL+T R ENTo+v-USPhydro

We combine, for the first time, event-by-event T R ENTo initial conditions with the relativistic viscous hydrodynamic model v-USPhydro and the Monte Carlo event generator JEWEL to make predictions for the nuclear modification factor R AA and jet azimuthal anisotropies v n { 2 } in $\sqrt{s_{NN}}$ = 5.02 TeV PbPb collisions for multiple centralities and values of the jet cone radius R. The R-dependence of R AA and v 2 { 2 } strongly depends on the presence of recoiling scattering centers. We find a small jet v 3 { 2 } in mid-central collisions and consistent results in wide jet p T regions and centralities with ATLAS data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ab-initio nucleon-nucleon correlations and their impact on high energy 16 O+ 16 O collisions

Investigating nucleon-nucleon correlations inherent to the strong nuclear force is one of the core goals in nuclear physics research. We showcase the unique opportunities offered by collisions of 16 O nuclei at high-energy facilities to reveal detailed many-body properties of the nuclear ground state. We interface existing knowledge about the geometry of 16 O coming from ab-initio calculations of nuclear structure with transport simulations of high-energy 16 O+ 16 O collisions. Bulk observables in these processes, such as the elliptic flow or the fluctuations of the mean transverse momentum, are found to depend significantly on the input nuclear model and to be sensitive to realistic clustering and short-range repulsive correlations, effectively opening a new avenue to probe these features experimentally. This finding demonstrates collisions of oxygen nuclei as a tool to elucidate initial conditions of small collision systems while fostering connections with effective field theories of nuclei rooted in quantum chromodynamics (QCD).

Zhang, Chunjian [Fudan University, Shanghai (China↗
DOE OSTI2025

Proton spin from small-x with constraints from the valence quark model

We apply the valence quark model [1] to constrain the non-perturbative initial condition for the small-x helicity evolution. The remaining free parameters are constrained by performing a global analysis akin to [2] to the available polarized small-x deep inelastic scattering data. A good description of the world data is obtained with only 8 free parameters. The model parameters are tightly constrained by the data, allowing us to predict the proton polarized structure function $g^p_1$ to be negative at small x. Furthermore, we obtain the small-x quark and gluon spins to give a contribution $\int^{0.1}_{10^{-5}} dx (\frac{1}{2}\Delta\Sigma + \Delta{G}) = 0.63 ± 0.10$ or $1.35 ± 0.16$ (or to the proton spin, depending on the applied running coupling prescription.

Color-class condensate↗
DOE OSTI2025

Symmetric-asymmetric collision comparison: Disentangling nuclear structure and subnucleonic structure effects for small system flow

Previous flow measurements in small collision systems were mostly based on highly asymmetric collisions (p+Pb, p+Au, d+Au, 3 He+Au), where both nuclear structure and subnucleonic fluctuations are important. Comparing these asymmetric systems with the newly available symmetric 16 O+ 16 O collisions at RHIC and LHC provides a unique opportunity to disentangle these two contributions. Using Glauber models incorporating both nucleon and quark-level substructure, we analyze multiplicity distributions and initial-state estimators: eccentricities ε n for anisotropic flow ν n and inverse transverse size d ⟂ for radial flow. We find that subnucleonic fluctuations impact O+O collisions differently from asymmetric systems, creating specific patterns in flow observables that enable disentangling the competing contributions. Such experimental comparisons will reduce uncertainties in the initial conditions and improve our understanding of the properties of the QGP-like medium produced in small systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis

Uncertainty in predicting solar energy resources introduces major challenges in power system management and necessitates the development of reliable probabilistic solar forecasts. As the first part of the development of probabilistic forecasts based on the Weather Research and Forecasting model with solar extensions (WRF-Solar), this study presents a tangent linear approach to identify input variables responsible for the largest uncertainties in predicting surface solar irradiance and clouds. A tangent linear analysis is capable of efficiently investigating sensitivities of output variables with respect to various input variables of WRF-Solar because this approach avoids the computational burden of perturbing the initial conditions of individual input variables. We develop tangent linear models (TLMs) for six WRF-Solar physics packages that control the formation and dissipation of clouds and solar radiation, and we evaluate the validity of TLMs using a linearity test. The tangent linear sensitivity analysis is conducted under various scenarios based on satellite observations and model simulations to consider realistic input conditions. A simple method is used to quantify the impact of the uncertainty of input variables on the output variables from the TLMs. The results demonstrate that uncertainties in the output variables that are the focus of this study—including global horizontal irradiance, direct normal irradiance, cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes—are highly sensitive to uncertainties in 14 input variables. This study indicates that the tangent linear method can identify key variables of physics modules in WRF-Solar that can be stochastically perturbed to generate ensemble-based probabilistic forecasts.

14 SOLAR ENERGY↗