Search NASA⌕ Search

SEARCH · Search NASA

Results for “PDF equation”

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.

141 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Predictive drug dosage control through a Fokker–Planck observer

During the years, various techniques have been presented for optimal control of tumor growth for which its stochastic behavior has rarely been considered. This article uses the well-known Gompertz stochastic model, for describing tumor growth, and by receding horizon model predictive control (RH-MPC) scheme computes the drug dosage for minimizing the tumor-cell population. To do this, a predictive control problem is defined based on the difference between the probability density function of tumor-cell population and the desired probability density function. In the model, both the drug dosage limitation and the Fokker–Planck equation have been considered as constraints. By solving this problem, the drug dosage which is considered as an external input to tumor model, is computed. In this article, the Fokker–Planck equation is used (1) as a nonlinear observer of probability density function of tumor-cell population and (2) a mapping vehicle from stochastic to deterministic domain. The Fokker–Planck equation is written for the Gompertz stochastic model of tumor growth. The equation is then solved through the path integral method. In this way, the probability density function of tumor-cell population, which contains the entire stochastic characteristics of the tumor growth, is obtained for any instance of time. The simulation results have also been presented for the evaluation of the suggested approach. The results show that the tumor-cell population can be controlled within a number of time windows (15 time windows in our case study) if an appropriate desired PDF is selected.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Additive Molar Volumes in Amorphous Ca/Sr Carbonate Solid Solutions

The development of predictive models of minor element incorporation in crystalline carbonate end products requires an understanding of the fundamental controls on metastable intermediate phase composition. In this study, we used small-angle X-ray scattering (SAXS), X-ray pair distribution function (PDF), thermogravimetric analysis (TGA), inductively coupled plasma mass spectrometry (ICP-MS), and transmission electron microscopy (TEM) to determine the composition-dependent density of an amorphous calcium–strontium carbonate (ACSC) solid solution. The amorphous calcium carbonate (ACC) and strontium carbonate (ASC) measured densities were ρ ACC = 2.19 ± 0.04 g/cm 3 and ρ ASC = 2.97 ± 0.05 g/cm 3 . Throughout our work, the experiments showed a dependence of the water content of the amorphous solid solution on the Sr mole fraction. An equation that relates the molar volume to the average cation radius, the carbonate ion radius, and the water volume was parameterized for hydrated crystalline carbonates and predicted well the molar volume of the ACSC solid solution. This finding indicates that, as for hydrated crystalline carbonates, the molar volumes of amorphous carbonates are additive and that water is a structural component of ACSC. Ab initio molecular dynamics (AIMD) simulations of the ACSC solid solutions showed strong linear correlations between calculated molar volumes and Sr and H 2 O contents, thus supporting the experimental results. In conclusion, our findings highlight the need to consider the full CaCO 3 –MeCO 3 –H 2 O ternary when quantifying metal cation incorporation in ACC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

SU-F-BRD-07: Empirical Derivation of Site-Specific Margin Formulas

Purpose: To empirically derive margin formulas from existing clinical radiotherapy plans accounting for respiratory motion and setup uncertainties. Methods: We simulated realistic treatment scenarios, including respiratory motion and setup errors. Individual probability density functions (PDF) from respiratory data were used to simulate respiratory motion. Random (σ) and systematic (Σ) setup errors were modeled as Gaussian distributions. One-dimensional dose profiles were extracted from existing radiotherapy plans and convolved with respiratory PDFs and random error distributions to produce blurred dose profiles. Each blurred dose profile was then shifted 1000 times by randomly sampling the simulated systematic error distribution. Margins were determined from the distance between the simulated treatment and the original 95% isodose level. An equation was fit for each (σ, Σ) combination to derive margin formulas for 90% of the population receiving 95% dose. This methodology can be applied to different tumor sites. Here, dose profiles were extracted from partial breast 3DCRT plans in the anterior-posterior (AP) and superior-inferior (SI) directions. Respiratory motion data was from healthy volunteers, and a clinically relevant range of random and systematic setup errors (standard deviations 1 – 4 mm) was determined from the literature. Results: The PBI margin formulas in the AP and SI directions for 95% dose coverage for 90% of the population were very similar: M= 0.68σ + 1.54Σ and M= 0.72σ + 1.50Σ, respectively. Systematic setup errors had the largest influence on required margin size, whereas realistic respiratory amplitude had minimal impact. The derived formulas resulted in a smaller systematic component than commonly-used theoretical margin recipes. Conclusion: We have demonstrated a method to derive empirical margin formulas from existing patient radiotherapy plans, incorporating realistic respiratory motion and appropriate ranges of random and systematic error. Site-specific margin formulas may be better suited to patient populations than theoretically derived, ideal geometry general margin formulas.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A Novel Feature of Valence Quark Distributions in Hadrons

Examining the evolution of the maximum of valence quark distribution, qV, weighted by Bjorken x, h(x,t)≡xqV(x,t), it is observed that h(x,t) at the peak becomes a one-parameter function; h(xp,t)=Φ(xp(t)), where xp is the position of the peak, t=logQ2, and Q2 is the resolution scale. This observation is used to derive a new model-independent relation which connects the partial derivative of the valence parton distribution functions (PDFs) in xp to the quantum chromodynamics (QCD) evolution equation through the xp derivative of the logarithm of the function Φ(xp(t)). A numerical analysis of this relation using empirical PDFs results in an observation of the exponential form of the Φ(xp(t))=h(xp,t)=CeDxp(t) for leading to next-to-next leading order approximations of PDFs for the range of Q2, covering four orders in magnitude. The exponent, D, of the observed “height-position” correlation function converges with the increase in the order of approximation. This result holds for all the PDF sets considered. A similar relation is observed also for the pion valence quark distribution, indicating that the obtained relation may be universal for any non-singlet partonic distribution. The observed “height-position” correlation is used also to indicate that no finite number of exchanges can describe the analytic behavior of the valence quark distribution at the position of the peak at fixed Q2.

Leon, Christopher↗
DOE OSTI2024

Interpolating the ‘t Hooft model between Instant and Light-Front dynamics in the Coulomb Gauge

The 1+1D model of quantum chromodynamics (QCD) in the infinite number of colors, or 't Hooft model, can be interpolated between the instant form dynamics (IFD) and the light-front dynamics (LFD) using an interpolation parameter 0 (IFD) ≤ δ ≤ π/4 (LFD). This was realized in the interpolating axial gauge which links the axial gauge (A 1 = 0) in IFD and the light-front gauge (A + = 0) [1]. In this presentation, we discuss the corresponding realization in the interpolating Coulomb gauge which links the temporal gauge (A 0 = 0) in IFD and the light-front gauge (A + = 0) and its benefit of resolving the issue associated with the absence of the conjugate field to the gauge field A 0 in the axial gauge. In both gauges, all degrees of freedom are physical making these gauge choices ideal for finding the bound-state equations and for renormalizability. Although the gauge independence of the physical observables such as the meson mass spectra following Regge trajectories may be guaranteed due to the gauge symmetry of QCD, the realization and interpretation of the identical physical results may depend on the gauge choices. Here, we discuss such difference in the realization of the confinement phenomena ala linear potential in the two different gauges, Coulomb vs. Axial, and highlight the gauge independent physical results expected. We also comment on the utility of the interpolation which leads to an alternative quasi-PDF that can be implemented in the lattice QCD without suffering from the large momentum boost.

Duggin, Hunter↗
DOE OSTI2024

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Fluid age-based analysis of a lifted turbulent DME jet flame DNS [Residence Time-Based Analysis of a Lifted Turbulent DME Jet Flame]

The link between the distribution of fluid residence time and the distribution of reactive scalars is analysed using Direct Numerical Simulation data. Information about the reactive scalar distribution is needed in order to model the reaction terms that appear in Large Eddy and Reynolds-Averaged simulations of turbulent reacting flows. The lifted flame is simulated taking account of multi-step chemistry for dimethyl-ether fuel and differential diffusion. Due to autoignition and flame propagation, the reaction progress increases with residence time. The variation of fluid residence time is evaluated by solving an Eulerian transport equation for the fluid age. The fluid age is a passive scalar with a spatially-uniform source term, meaning that its moments and dissipation rates in turbulent flows can be modelled using closures already established for conserved scalars such as mixture fraction. In combination with the mixture fraction, the fluid age serves as a useful mapping variable to distinguish younger less-reacted fluid near the inlet from older more-reacted fluid downstream. The local fluctuations of mixture fraction and fluid age have strong negative correlation and, building upon established presumed-pdf models for mixture fraction, this feature can be used to construct an accurate presumed-pdf model for the joint mixture fraction/fluid age pdf. It is demonstrated that the double-conditional first-order moment closure combined with the proposed presumed model for the joint pdf of mixture fraction and fluid age gives accurate predictions for unconditional reaction rates – both for pre-ignition radical species produced by low-temperature processes upstream of the flame base, and for major species that are produced at the flame front.

42 ENGINEERING↗
DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

A mapping closure for turbulent scalar mixing using a time-evolving reference field

A general mapping-closure approach for modeling scalar mixing in homogeneous turbulence is developed. This approach is different from the previous methods in that the reference field also evolves according to the same equations as the physical scalar field. The use of a time-evolving Gaussian reference field results in a model that is similar to the mapping closure model of Pope (1991), which is based on the methodology of Chen et al. (1989). Both models yield identical relationships between the scalar variance and higher-order moments, which are in good agreement with heat conduction simulation data and can be consistent with any type of epsilon(phi) evolution. The present methodology can be extended to any reference field whose behavior is known. The possibility of a beta-pdf reference field is explored. The shortcomings of the mapping closure methods are discussed, and the limit at which the mapping becomes invalid is identified.

Girimaji, Sharath S.↗
DOE OSTI2020

Three-dimensional imaging of proton in basis light-front quantization

Here, we employ an effective Hamiltonian that includes the transverse and longitudinal confinement and the one-gluon exchange interaction with fixed coupling constant. By solving the eigenvalue equation in basis light-front quantization (BLFQ), we generate the light-front wavefunctions (LFWFs) for the nucleon in the valence quark Fock space. Fitting the model parameters, we obtain high quality descriptions of electromagnetic form factors and radius for proton while the results deviate somewhat from experimental data for neutron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2013

The impedance problem of wave diffraction by a strip with higher order boundary conditions

This work is devoted to analyse an impedance boundary-transmission problem for the Helmholtz equation originated by a problem of wave diffraction by an infinite strip with higher order imperfect boundary conditions. A constructive approach of operator relations is built, which allows a transparent interpretation of the problem in an operator theory framework. In particular, different types of operator relations are exhibited for different types of operators acting between Lebesgue and Sobolev spaces on a finite interval and the positive half-line. All this has consequences in the understanding of the structure of this type of problems. In particular, a Fredholm characterization of the problem is obtained in terms of the initial space order parameters. At the request of the author and the Proceedings Editor the above article has been replaced with a corrected version. The original PDF file supplied to AIP Publishing contained an error in the title of the article. The original title appeared as: 'The Impedance Problem of Wave Diffraction by a trip with Higher Order Boundary Conditions.' This article has been replaced and the title now appears correctly online. The corrected article was published on 8 November 2013.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Gluon splitting at small x: a unified derivation for the JIMWLK, DGLAP and CSS equations

We revisit the calculation of the next-to-leading order (NLO) corrections to dijet production in electron-ion collisions at small x. We focus on the back-to-back configuration where the relative transverse momentum P ⊥ of the measured jets is much larger than both their momentum imbalance K ⊥ and the target saturation momentum Q s (x, A). In this regime, we present for the first time a complete calculation of the real NLO corrections, associated with gluon emissions outside the jet cones, at leading power in 1/P ⊥ . Our result exhibits TMD factorisation, with the same hard factor as at tree-level and a NLO correction to the Weiszäcker-Williams (WW) gluon transverse momentum dependent (TMD) distribution which involves four Wilson-line operators. By studying different kinematical regimes for K ⊥ and for the radiated gluon, we recover all the quantum evolutions that were previously identified for this process at NLO: the B-JIMWLK high-energy evolution and the CSS evolution of the gluon WW TMD, and the DGLAP evolution of the gluon PDF. When both K ⊥ and the transverse momentum transferred by the target are large compared to Q s , all the Wilson-line operators boil down to the unintegrated gluon distribution and our NLO result for the gluon TMD can be used to isolate the transverse-momentum dependent gluon splitting function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Computation of probabilistic hazard maps and source parameter estimation for volcanic ash transport and dispersion

Volcanic ash advisory centers are charged with forecasting the movement of volcanic ash plumes, for aviation, health and safety preparation. Deterministic mathematical equations model the advection and dispersion of these plumes. However initial plume conditions – height, profile of particle location, volcanic vent parameters – are known only approximately at best, and other features of the governing system such as the windfield are stochastic. These uncertainties make forecasting plume motion difficult. As a result of these uncertainties, ash advisories based on a deterministic approach tend to be conservative, and many times over/under estimate the extent of a plume. This paper presents an end-to-end framework for generating a probabilistic approach to ash plume forecasting. This framework uses an ensemble of solutions, guided by Conjugate Unscented Transform (CUT) method for evaluating expectation integrals. This ensemble is used to construct a polynomial chaos expansion that can be sampled cheaply, to provide a probabilistic model forecast. The CUT method is then combined with a minimum variance condition, to provide a full posterior pdf of the uncertain source parameters, based on observed satellite imagery. The April 2010 eruption of the Eyjafjallajökull volcano in Iceland is employed as a test example. The puff advection/dispersion model is used to hindcast the motion of the ash plume through time, concentrating on the period 14–16 April 2010. Variability in the height and particle loading of that eruption is introduced through a volcano column model called bent. Output uncertainty due to the assumed uncertain input parameter probability distributions, and a probabilistic spatial-temporal estimate of ash presence are computed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Bayesian quantification of observability and equation of state of twin stars

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equations of state (EOSs) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a metamodel with nine parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12 and 18% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. Here, we show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared 𝑐$^2_{qm}$ has two peaks, one below and another above the conformal limit 𝑐$^2_{qm}$ = 1/3 predicted by perturbative QCD.

QCD in nuclear reactions↗
DOE OSTI2025

SPRUCE Root Production Assessed with Manual Minirhizotrons Resolved to Plant Functional Type, 2015-2021

This dataset contains raw root length and diameter for individual roots and estimated root population production measurements from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experimental site within the Marcell Experimental Forest in northern Minnesota, USA. Measurements started at the beginning of whole ecosystem warming manipulations in 2015 through 2021 (2015-05-26 to 2021-09-01). Root morphology and estimated production were quantified throughout the peat profile with manual minirhizotrons deployed within SPRUCE plots. Images were processed using commercial software to quantify the length and diameter of individual roots. Roots were visually assigned to a plant functional type (PFT) of either (ericaceous) shrub, herb (sedges and Maianthemum trifolium), or tree (Larix laricina, Picea mariana) based on expert opinion. The biomass of individual roots was estimated using PFT-specific allometric equations (Iversen et al., 2018). Production per day was estimated as the length of new roots produced between imaging sessions, divided by the number of days between imaging sessions. These values were placed on a m2 aboveground area basis and scaled to a standard depth of 1m (roots are not evenly distributed with depth, do not interpret value as being on a m3 basis). Maximum and average (weighted by production length) depth of each PFT were also estimated within each minirhizotron tube. Annual production was interpolated as the average of four methods to scale these data (see Weber et al, 2026). Standing crop of roots was estimated for each tube as the maximum visible amount (both length and mass) of roots of that PFT for that year. These data expand the ability of researchers to accurately estimate the belowground dynamics of peatland vegetation, as well as the role that fine roots may play in impacting the fluxes of carbon within peatlands. This dataset contains three data files in comma-separate values (*.csv) format. This dataset contains one data file in comma-separate values (.csv) format. Additional metadata are provided: three data dictionaries and a file-level metadata file in comma-separate values (.csv) format and a user guide in PDF (*.pdf) format.

54 ENVIRONMENTAL SCIENCES↗