Search NASA⌕ Search

SEARCH · Search NASA

Results for “analytical models”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probabilistic inference of the structure and orbit of Milky Way satellites with semi-analytic modelling

Semi-analytic modelling furnishes an efficient avenue for characterizing dark matter haloes associated with satellites of Milky Way-like systems, as it easily accounts for uncertainties arising from halo-to-halo variance, the orbital disruption of satellites, baryonic feedback, and the stellar-to-halo mass (SMHM) relation. We use the SatGen semi-analytic satellite generator, which incorporates both empirical models of the galaxy–halo connection as well as analytic prescriptions for the orbital evolution of these satellites after accretion onto a host to create large samples of Milky Way-like systems and their satellites. By selecting satellites in the sample that match observed properties of a particular dwarf galaxy, we can infer arbitrary properties of the satellite galaxy within the cold dark matter paradigm. For the Milky Way’s classical dwarfs, we provide inferred values (with associated uncertainties) for the maximum circular velocity v max and the radius r max at which it occurs, varying over two choices of baryonic feedback model and two prescriptions for the SMHM relation. While simple empirical scaling relations can recover the median inferred value for v max and r max , this approach provides realistic correlated uncertainties and aids interpretability. We also demonstrate how the internal properties of a satellite’s dark matter profile correlate with its orbit, and we show that it is difficult to reproduce observations of the Fornax dwarf without strong baryonic feedback. Furthermore, the technique developed in this work is flexible in its application of observational data and can leverage arbitrary information about the satellite galaxies to make inferences about their dark matter haloes and population statistics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Analytical Model for the Higher Order Moments of Midlatitude Atmospheric Temperature Distributions

Observed distributions of atmospheric temperature are non-Gaussian. Therefore, moments beyond variance are necessary in determining the frequency of extreme temperature events. Here we propose a simple kinematic model for atmospheric mid-latitude temperature variability based on symmetric advection from a non-symmetric background temperature profile. We then use this model to derive analytical expressions for the higher order moments of temperature distributions. Our results show that nonzero skewness and kurtosis arise due to the nonlinearity of the time-mean meridional temperature profile. The analytical model matches an idealized Held-Suarez atmospheric model, indicating nonlinearity of time-mean temperature in latitude is the dominant contribution to nonzero skewness and kurtosis in synoptic temperature variations. Model analysis further shows decrease in higher order moments due to climate change come roughly equally from changes in mixing length and changes in the background temperature profiles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Triboplasma assisted chemical conversion in granular systems: a semi-analytic model

We present a semi-analytic model for a novel plasma-assisted chemical conversion pathway using triboplasmas generated in granular flows. Triboelectric charge relaxation is a well known phenomena where the potential generated from contact charging of particles exceeds the breakdown voltage of the background gas. In this work, we extend the triboelectric charge relaxation theory to include non equilibrium plasma energy and particle balance equations to predict the formation of dissociated and excited species that act as precursors to chemical conversion, for example in plasma-assisted ammonia synthesis. Our example case study with nitrogen background gas and teflon/aluminum tribomaterial system yielded high excited nitrogen species densities per collision that are comparable to current plasma-assisted conversion pathways. We also present a regime diagram for various gases where Paschen breakdown parameters are used to determine whether triboplasmas can be formed for a given effective work-function difference between two materials. Our sensitivity studies indicate particle velocity, particle radius, solids fraction and space charge effects play a critical role in overall plasma densities and excited species production.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A predictive analytical model of electrical transport in multi-principal-element alloys

A predictive analytical model is presented for the electrical conductivity of multi-principal-element alloys (MPEAs), including those containing aluminum, transition metals, and refractory metals. Given that the lattice parameter of the Wigner-Seitz cell of an MPEA is similarly variable to a bulk metallic glass, it is postulated that electron scattering can be approximated by a series of two-level systems. Here, the resulting reduced-order model enabled an accurate determination of electrical resistivity and electron thermal conductivity based on the scattering of electrons in a two-level system across a Bloch-potential-based virtual crystal approximation. Model results are compared to experimental four-point probe electrical resistivity measurements between 300 K and 700 K for Al 0.3 CoCrCuFeNi, CoCrFeMnNi, (CoCrFeMnNi) 0.98 W 0.02 , (CoCrFeMnNi) 0.95 W 0.05 , and Nb 4 Ta 4 V 3 Ti, for model validation.

Analytical model↗
DOE OSTI2026

Analytic model for neutral penetration and plasma fueling

Neutral atoms recycled from wall interaction interact with confined plasma, thereby refuelling it, most strongly in the region closest to the wall. This occurs near the X-point in diverted configurations, or else near the wall itself in limited configurations. A progression of analytic models is developed for neutral density in the vicinity of a planar or linear source in an ionising domain. First-principles neutral transport simulations with DEGAS2 are used throughout to test the validity and limits of the model when using equivalent sources. The model is further generalised for strong plasma gradients or the inclusion of charge exchange. An important part of the problem of neutral fuelling from recycling is thereby isolated and solved with a closed-form analytic model. A key finding is that charge exchange with the confined plasma can be significantly simplified with a reasonable sacrifice of accuracy by treating it as a loss. The several assumptions inherent to the model (and the simulations with which it is compared) can be adapted according to the particular behaviour of neutrals in the divertor and the manner in which they cross the separatrix.

fusion plasma↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
DOE OSTI2025

A Unified Analytical Model for Pressure Solution With Fully Coupled Diffusion and Reaction

Abstract Geophysical models for pressure solution are typically developed for diffusion‐controlled or reaction‐controlled scenarios. We present a unified analytical model that considers fully coupled diffusion and reaction during pressure solution. The model recovers the diffusion‐controlled and reaction‐controlled models in the literature as specific limiting cases. When diffusion and reaction exhibit comparable influences, we validate the proposed model against independent numerical simulations. The proposed model is then employed in interpreting experimental measurements, demonstrating a better agreement compared to previous interpretations.

Wang, Ziyan [Energy Geosciences Division Lawrence ↗
DOE OSTI2026

Boron Coordination in Multicomponent Glasses: Analytical Models and Machine Learning With Uncertainty

Borosilicate glasses are extensively used in a variety of applications from kitchenware to nuclear waste immobilization due to the strong network formed by the Si-O-B bond that makes it resistant to chemical corrosion and gives it a low thermal expansion. Boron, however, exists in both trigonal BO3 and tetrahedral BO4 bonds in glass systems, which impacts the chemical durability and thermal resistance of the glass, amongst other properties. Boron coordination (N4), or the ratio of the amount of BO4 to BO3 within a glass, may aid in predicting these properties but is difficult to derive without experimental data due to the complexity of impacts from varied glass compositions and processing factors. For this reason, compositional models have been developed to predict boron coordination, but the models typically include a limited number of glass components. To help fill this gap in the models, in this work, a diverse multicomponent glass dataset of 809 glasses is compiled from a literature search, and then a number of analytical and machine learning (ML) models are trained on the dataset. Previously developed modified Bernstein and modified Du Stebbins analytical models were fitted to update parameters with the new dataset. Then, partially Bayesian neural networks, Gaussian process regressor, and heteroskedastic deterministic neural networks were evaluated. The ML models examined all have different strategies to overcome the potential for overfitting as a result of a limited training dataset, and return results that account for model uncertainty, which can be valuable for understanding model reliability. For the first time, cooling rate is introduced as an input parameter for ML models, showing consistent improvements in performance and solidifying the importance of including parameters outside of composition alone for N4 prediction. The machine learning models examined here show promise in accurate predictions of boron coordination in borosilicate glasses, all achieving R2 values of 0.91.

boron coordination↗
DOE OSTI2026

Comprehensive analytical model of the dynamic 𝑍 pinch

In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Analytical models of hydrogen transport in graphite

The importance of graphite-hydrogen chemical reactions to fusion, fission, and hydrogen storage applications, combined with the rapidly evolving knowledge on the underlying mechanisms, has led to the development of multiple models to describe hydrogen transport in graphite. Significant differences exist among these models, resulting from discrepancies in the modeling assumptions, intended degree of fidelity, and conditions of applicability. This paper attempts at reconciling these apparent differences by providing a comprehensive description of the constitutive equations governing hydrogen transport in graphite at high-temperature, identifying outstanding gaps in knowledge, illustrating how these different models approach them, and proposing alternative analytical formulations grounded on experimental results from hydrogen-graphite studies. Governing equations, closing relations, and simplifying assumptions are discussed for hydrogen transport at the inter-granular and intra-granular level, accompanied by compiled experimental data and illustrated energy diagrams associated to the proposed transport mechanisms. Analytical formulations are provided to reproduce competing hypotheses on the mechanisms, supporting the development of a range of computational models that can enable resolution of outstanding knowledge gaps through comparative testing against experimental data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Transport coefficient sensitivities in a semi-analytic model for magnetized liner inertial fusion

Performance of magnetized liner inertial fusion (MagLIF) experiments is highly dependent on transport processes including magnetized heat flows and magnetic flux losses. Magnetohydrodynamic simulations used to model these experiments require a choice of model for the transport coefficients, which are the constants of proportionality relating driving terms, such as temperature gradients and currents, to the associated heat and magnetic field transport. The coefficients have been the subject of repeated recalculation using various methods throughout the years. Using a semi-analytic MagLIF model, we compare models for the transport coefficients. The choice of model modifies magnetic-flux losses caused by the Nernst thermoelectric effect and thermal conduction losses. We present simulated results from parameter scans conducted in order to compare the effects of the different models on parameters of interest in MagLIF. In some regions of parameter space, discrepancies of up to 38% are found in integrated quantities like the fusion yield. These results may serve as a guide for experimental validation of the various models, particularly as laser preheat energies and initial axial field strengths are increased on MagLIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Analytical model for balling defects in laser melting using rivulet theory and solidification

In laser welding and additive manufacturing communities, the balling (humping) defect is primarily attributed to the Plateau–Rayleigh fluid instability (PRI) with a few authors suggesting fluid jetting and volume conservation as alternative mechanisms. As analytical descriptions of these mechanisms are unavailable, combining them into a single formalism is unfeasible. We present a new model of PRI with higher accuracy, accounting for competition with solidification, to compare the expected behavior with known experimental trends when fluid jetting is neglected. We adapt a rivulet instability model from fluid physics to account for the stabilizing effects of the substrate, which the traditional cylindrical-jet geometry does not account for and estimate the instability growth rate. Our model yields a continuous transition from non-balling to balling hitherto lacking in current literature and predicts instability growth at higher wavelengths with strong sensitivity to the solidification front curvature. While the fluid surface is most unstable for shallow melt pools, the absolute magnitude of balling relevant to printing defects scales with melt pool depth and has a maximum for a given melt pool geometry. Synchrotron-based x-ray radiography of thin samples indicates that PRI growth rates and solidification can be comparable in magnitude and thus compete, as we find in this work. We predict that deviations between model predictions and our experimental results demonstrate the importance of fluid flows and heat transport in the balling process. Our experiments further demonstrate at least one mechanism by which the melt pool length and the balling wavelength are not equivalent, as commonly claimed.

Taylor, Z. [Stanford Univ., CA (United States); SL↗
DOE OSTI2024

Analytical Model for Atomic Relaxation in Twisted Moiré Materials

By virtue of being atomically thin, the electronic properties of heterostructures built from two-dimensional materials are strongly influenced by atomic relaxation. The atomic layers behave as flexible membranes rather than rigid crystals. Here we develop an analytical theory of lattice relaxation in twisted moiré materials. We obtain analytical results for the lattice displacements and corresponding pseudo gauge fields, as a function of twist angle. We benchmark our results for twisted bilayer graphene and twisted WSe 2 bilayers using large-scale molecular dynamics simulations. Our single-parameter theory is valid in graphene bilayers for twist angles 𝜃 ≳ 0.7°, and in twisted WSe 2 for 𝜃 ≳ 1.6°. Furthermore, we also investigate how relaxation alters the electronic structure in twisted bilayer graphene, providing a simple extension to the continuum model to account for lattice relaxation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of errors in analytic modeling of permanent magnet electron spectrometers for laser-plasma accelerators

Here, we assess the error in three treatments of a critical electron beam diagnostic for laser-plasma acceleration (LPA) experiments: a permanent magnet electron spectrometer. Since LPA electron spectrometers are often difficult to calibrate due to the scarcity of well-characterized, tunable electron beam sources in the appropriate energy range and the mechanical complexity of electron spectrometers, the standard of calibration is Hall probe measurements. We first compare the electron spectrometer performance between SIMION calculations and the Hall probe measurements. We find up to a 60% (<8%) error in determining the absolute energy for electrons below (above) 80 MeV when using the SIMION-modeled magnetic field vs the measured magnetic field due to SIMION overpredicting fringe field strength. The difference in spectrometer energy resolution is within ∼10% between the modeled and measured fields. We then assess a simple block model, which is commonly used. We find that the simple block model presented here sufficiently captures the predicted as-built spectrometer performance for any application provided that the user can tolerate uncertainties on absolute energy determination up to 5% and errors on energy resolution up to 1%.

Electromagnetic optics↗