Search NASA⌕ Search

SEARCH · Search NASA

Results for “FORMULAS - BUCKLING”

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 OSTI2016

Numerically stable, scalable formulas for parallel and online computation of higher-order multivariate central moments with arbitrary weights

Formulas for incremental or parallel computation of second order central moments have long been known, and recent extensions of these formulas to univariate and multivariate moments of arbitrary order have been developed. Such formulas are of key importance in scenarios where incremental results are required and in parallel and distributed systems where communication costs are high. We survey these recent results, and improve them with arbitrary-order, numerically stable one-pass formulas which we further extend with weighted and compound variants. We also develop a generalized correction factor for standard two-pass algorithms that enables the maintenance of accuracy over nearly the full representable range of the input, avoiding the need for extended-precision arithmetic. We then empirically examine algorithm correctness for pairwise update formulas up to order four as well as condition number and relative error bounds for eight different central moment formulas, each up to degree six, to address the trade-offs between numerical accuracy and speed of the various algorithms. Finally, we demonstrate the use of the most elaborate among the above mentioned formulas, with the utilization of the compound moments for a practical large-scale scientific application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Annotation of DOM metabolomes with an ultrahigh resolution mass spectrometry molecular formula library

Current approaches to analyzing metabolomic data often rely on matching MS/MS fragmentation data to sparse libraries or databases. This approach results in limited identification of features, often with less than 10% of the dataset being annotated. A complementary approach is to assign molecular formula to features based on accurate mass measurements, but the platforms commonly used for metabolomics do not have the needed accuracy or resolving power to do this robustly, particularly for larger molecules. Using our newly modified analysis tool, CoreMS, we generated a library of molecular formula from pooled samples analyzed with LC-21T FT-ICR MS. This library successfully annotated approximately 53.2% of features identified from the exometabolome of marine diatom Phaeodactylum tricornutum – a nearly ten-fold increase over the 5.9% annotation rate achieved using a conventional MS/MS library matching approach. Using this FT-ICR MS library approach, we were able to differentiate differences in the exometabolome of P. tricornutum in iron replete and iron limited conditions, with 668 metabolites being differentially expressed (p < 0.05, 2 x intensity difference) under these conditions. The traditional MS/MS fragmentation-based annotation approach only annotated 61 of these metabolites, while our novel pipeline annotated 450 metabolites and revealed 12 metabolites that were significantly more abundant under low iron conditions. Our results demonstrate the utility of ultrahigh resolution mass spectrometry for generating more comprehensive and confident molecular annotations.

21T-FTICR-MS, CoreMS↗
DOE OSTI2014

Compact formulas for bounce/transit averaging in axisymmetric tokamak geometry

Compact formulas for bounce and transit orbit averaging of the fluctuation-amplitude eikonal factor in axisymmetric tokamak geometry, which is frequently encountered in bounce-gyrokinetic description of microturbulence, are given in terms of the Jacobi elliptic functions and elliptic integrals. These formulas are readily applicable to the calculation of the neoclassical susceptibility in the framework of modern bounce-gyrokinetic theory. In the long-wavelength limit for axisymmetric electrostatic perturbations, we recover the expression for the Rosenbluth-Hinton residual zonal flow [M. N. Rosenbluth and F. L. Hinton, Phys. Rev. Lett. 80, 724 (1998)] accurately.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Generalized Lapse Rate Formulas for Use in Entraining CAPE Calculations

Abstract Common assumptions in temperature lapse rate formulas for lifted air parcels include neglecting mixing, hydrostatic balance, the removal of all condensate once it forms (pseudoadiabatic), and/or the retention of all condensate within the parcel (adiabatic). These formulas are commonly derived from the conservation of entropy, which leads to errors when nonequilibrium mixed-phase condensate is present. To evaluate these assumptions, a new general lapse rate formula is derived from an expression for energy conservation, rather than entropy conservation. This new formula incorporates mixing of the parcel with its surroundings, relaxes the hydrostatic assumption, allows for nonequilibrium mixed-phase condensate, and can be formulated for pseudoadiabatic or adiabatic ascent. The new formula is shown to exactly conserve entropy for reversible ascent. Predictions by the new formula are compared to that of older and less general formulas. The errors in previous formulas arise from the assumption of hydrostatic balance, which results in considerable warm biases due to the neglect of the energy sink from buoyancy. Predictions of ascent with entrainment using the new formula are then compared to parcel properties along trajectories in large eddy simulations. Simulated parcel properties are better predicted by the formula using a diluted analogy to adiabatic ascent, wherein condensate is diluted at the same rate as other parcel properties, than by the diluted analogy to pseudoadiabatic ascent, wherein all condensate is removed. These results suggest that CAPE should be computed with adiabatic, rather than pseudoadiabatic, parcel ascent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

New simplified formulae of estimation the total number of fissions based on the Nordheim-Fuchs model - 14297

This article presents simplified formulae for the estimation of the total number of fissions for criticality accidents. One of these formulae is based on the Nordheim-Fuchs model. An additional formula takes into account the boiling for solution systems. These formulae are adapted to homogeneous media and can be used for solution and metal accidents. The principle of these formulae, associated with their assumptions, is described. Then, these formulae are compared with practical applications for SILENE and CRAC reactors (solution systems) and past criticality accidents with metal assemblies as well as a previous simplified formula. The proposed new formulae give bounding estimates of the total number of fissions, provided that their area of applicability are respected and that the simulations performed to calculate the necessary input data are representative of the criticality accident configuration. These new formulae require the determination of more parameters than previous simplified formulae, which make them more complicated to use, but may point out and quantify the impact of some physical phenomenon (for instance for the solution boiling). In particular, the comparison with past criticality accidents with metal assemblies shows that more detailed simulations (in terms of geometry and reflection conditions) might be performed in future work to improve the estimate of the number of fissions based on the Nordheim-Fuchs model. (authors)

42 ENGINEERING↗
DOE OSTI2021

Dispersion formulas in QFTs, CFTs and holography

We study momentum space dispersion formulas in general QFTs and their applications for CFT correlation functions. We show, using two independent methods, that QFT dispersion formulas can be written in terms of causal commutators. The first derivation uses analyticity properties of retarded correlators in momentum space. The second derivation uses the largest time equation and the defining properties of the time-ordered product. At four points we show that the momentum space QFT dispersion formula depends on the same causal double-commutators as the CFT dispersion formula. At n -points, the QFT dispersion formula depends on a sum of nested advanced commutators. For CFT four-point functions, we show that the momentum space dispersion formula is equivalent to the CFT dispersion formula, up to possible semi-local terms. We also show that the Polyakov-Regge expansions associated to the momentum space and CFT dispersion formulas are related by a Fourier transform. In the process, we prove that the momentum space conformal blocks of the causal double-commutator are equal to cut Witten diagrams. Finally, by combining the momentum space dispersion formulas with the AdS Cutkosky rules, we find a complete, bulk unitarity method for AdS/CFT correlators in momentum space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Chapter Six - Derivation, cross-validation, and comparison of analytic formulas for electrostatic deflector aberrations

Here, we derived first and second order analytic aberration formulas in the horizontal transverse plane for an electrostatic deflector specified by the reference orbit radius, the central angle spanning the deflector, and its inhomogeneity coefficients. The derivation was performed using an iterative order-by-order perturbation method in a Frenet–Serret beamline coordinate system. We produced a C program edabrt for calculation of the first and second order aberrations using the formulas that we derived. The electrostatic deflector aberration formulas from Wollnik (1965) disagree with COSY INFINITY and the aberration formulas that we derived here, and they also deviate from the first and second order symplecticity conditions. The reason for this discrepancy is that the aberration formulas from Wollnik (1965) consider only motion in the main field of the deflector and do not properly account for fringe field effects, in particular the necessarily occurring change of kinetic energy due to the change in potential. The code GIOS uses electrostatic aberration formulas from Wollnik (1965) and exhibits the same discrepancy. We performed a comparison of first and second order electrostatic deflector aberrations for (1) the analytic formulas derived in this work, (2) differential-algebraic (DA) numerical integration of the equations of motion using the code COSY INFINITY, (3) the aberration formulas from Wollnik (1965), (4) aberrations computed using the code GIOS, and (5) the aberration formulas from Wollnik (1965) adjusted to account for the occurring change in potential. An electrostatic spherical deflector and an electrostatic cylindrical deflector were used as test cases for this comparison. There is excellent agreement between methods (1), (2), and (5), and these three methods satisfy the first and second order symplecticity conditions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

SU-E-P-45: An Analytical Formula for Deriving Mechanical Iso-Center of Rotational Gantry Treatment Unit Rotational Gantry Treatment Unit

Purpose: To present an analytical formula for deriving mechanical isocenter (MIC) of a rotational gantry treatment unit. The input data to the formula is obtained by a custom-made device. The formula has been implemented and used in an operational proton therapy facility since 2005. Methods: The custom made device consisted of 3 mutually perpendicular dial indicators and 5 clinometers, to obtain displacement data and gantry angle data simultaneously. During measurement, a steel sphere was affixed to the patient couch, and the device was attached to the snout rotating with the gantry. The displacement data and angle data were obtained simultaneously at angular increments of less than 1 degree. The analytical formula took the displacement and angle as input and derived the positions of dial indicator tips (DIT) position in room-fixed coordinate system. The formula derivation presupposes trigonometry and 3-dimentional coordinate transformations. Due to the symmetry properties of the defining equations, the DIT position can be solved for analytically without using mathematical approximations. We define the mean of all points in the DIT trajectory as the MIC. The formula was implemented in computer code, which has been employed during acceptance test, commissioning, as well as routine QA practice in an operational proton facility since 2005. Results: It took one minute for the custom-made device to acquire the measurement data for a full gantry rotation. The DIT trajectory and MIS are instantaneously available after the measurement. The MIC Result agrees well with vendor’s Result, which came from a different measurement setup, as well as different data analysis algorithm. Conclusion: An analytical formula for deriving mechanical isocenter was developed and validated. The formula is considered to be absolutely accurate mathematically. Be analyzing measured data of radial displacements as function of gantry angle, the formula calculates the MI position in room coordinate.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Feynman formulae and phase space Feynman path integrals for tau-quantization of some Lévy-Khintchine type Hamilton functions

Evolution semigroups generated by pseudo-differential operators are considered. These operators are obtained by different (parameterized by a number τ) procedures of quantization from a certain class of functions (or symbols) defined on the phase space. This class contains Hamilton functions of particles with variable mass in magnetic and potential fields and more general symbols given by the Lévy-Khintchine formula. The considered semigroups are represented as limits of n-fold iterated integrals when n tends to infinity. Such representations are called Feynman formulae. Some of these representations are constructed with the help of another pseudo-differential operator, obtained by the same procedure of quantization; such representations are called Hamiltonian Feynman formulae. Some representations are based on integral operators with elementary kernels; these are called Lagrangian Feynman formulae. Langrangian Feynman formulae provide approximations of evolution semigroups, suitable for direct computations and numerical modeling of the corresponding dynamics. Hamiltonian Feynman formulae allow to represent the considered semigroups by means of Feynman path integrals. In the article, a family of phase space Feynman pseudomeasures corresponding to different procedures of quantization is introduced. The considered evolution semigroups are represented as phase space Feynman path integrals with respect to these Feynman pseudomeasures, i.e., different quantizations correspond to Feynman path integrals with the same integrand but with respect to different pseudomeasures. This answers Berezin’s problem of distinguishing a procedure of quantization on the language of Feynman path integrals. Moreover, the obtained Lagrangian Feynman formulae allow also to calculate these phase space Feynman path integrals and to connect them with some functional integrals with respect to probability measures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

3D small-gain formula allowing strong focusing and harmonic lasing for a ring-based x-ray free electron laser oscillator

We present a detailed derivation of a formula for the small-gain calculation for an x-ray free electron laser oscillator (XFELO) based on a medium-energy (3–4 GeV) storage ring. We found harmonic lasing and strong focusing are essential for this beam energy range. Taking the small-signal low-gain formula developed by Kim and his colleagues, we modified it in such a way that the gain can be calculated without the “no focusing approximation,” and a strong focusing can be applied, as well as harmonic lasing. In this formula, the gain is represented as a product of two factors with one of them depending only on the harmonic number, undulator period, and gap. Using this factor, we show that it is favorable to use harmonic lasing to achieve hard x-ray FEL working in the small-signal low-gain regime with the medium-energy electron beam. Our formula also allows FEL optimization by varying the vertical gradient of the undulator, the vertical dispersion, and the horizontal and vertical focusing, independently. As an example, we applied this formula to study the feasibility of an XFELO option for the National Synchrotron Light Source II (NSLS-II) upgrade. Since a quite high peak current is required for the FEL, collective effects of beam dynamics in medium-energy synchrotrons significantly affect the electron beam parameters. We carried out a multiparameter optimization taking collective effects into account. Note, even though our example is for a ring-based XFELO at 3 to 4 GeV, the formula and, in particular, the approach developed here may be applied to other types of FELs. Published by the American Physical Society 2024

36 MATERIALS SCIENCE↗
DOE OSTI2023

Suzuki type estimates for exponentiated sums and generalized Lie-Trotter formulas in JB-algebras

Lie-Trotter-Suzuki product formulas are ubiquitous in quantum mechanics, computing, and simulations. Approximating exponentiated sums with such formulas are investigated in the JB-algebraic setting. We show that the Suzuki type approximation for exponentiated sums holds in JB-algebras, we give explicit estimation formulas, and we deduce three generalizations of Lie-Trotter formulas for arbitrary number elements in such algebras. In conclusion, we also extended the Lie-Trotter formulas in a Jordan Banach algebra from three elements to an arbitrary number of elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Analysis and fringe field scaling of a legacy set of electrostatic deflector aberration formulas [Chapter Seven]

In this work we performed an analysis of the derivation of first and second order analytic aberration formulas by Wollnik (1965) for the case of electrostatic deflectors in the horizontal plane. We found that these aberration formulas are valid for the main field of the deflector; however, they do not account for any fringe field effects. To address this issue, we modified the aberration formulas from Wollnik (1965) for a hard edge fringe field by scaling the particle inclinations at the entrance and exit edges of the deflector. These modified aberration formulas fully agree with electrostatic aberration formulas derived by Valetov and Berz (2020), as well as with differential-algebraic (DA) calculations of electrostatic deflector aberrations using the code COSY INFINITY.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

A modified formula for dose calculations of stereotactic ablative body radiotherapy for non–small cell lung cancer

To provide a modified formula consistent with the Monte Carlo (MC) algorithm for dose calculations during stereotactic ablative body radiotherapy for non–small cell lung cancer. Seventy CyberKnife treatment plans were calculated and analyzed by MC and ray-tracing (RT) algorithms, separately. Parameters of treatment plans were compared, and those associated with differences of dose distributions were analyzed to establish a modified formula. Gross tumor volume and tumor tracking volume (TTV) were defined as the evident disease on the sequences of the window width and level of the lung and the mediastinum. Additionally, the formula was validated by another 20 plans. The prescription dose of the 90 patients was 60 Gy/5f. The RT algorithm overestimated the planning target volume (PTV) D{sub 95} by an average of 8.59 Gy and the gross tumor volume D{sub 99} by an average of 5.84 Gy. The homogeneity index of PTV was underestimated by 0.11 on average, whereas the conformity index and new conformity index was underestimated by 0.05. The RT algorithm overestimated the dose distribution to the spinal cord by 2.23 Gy, the esophagus by 1.96 Gy, the trachea by 1.89 Gy, the left-sided bronchus by 1.77 Gy, the right-sided bronchus by 1.64 Gy, and the heart by 2.16 Gy. The average whole-lung dose volumes of lung tissues and dose volumes of V{sub 5} were overestimated by 2.69 Gy and 7.52%, respectively. A power function distribution (R{sup 2} = 0.8626) was confirmed between PTV D{sub 95} and TTV volumes. PTV D{sub 95} calculated by the MC algorithm could be computed easily with TTV and PTV D{sub 95} calculated by the RT algorithm based on the formula. The modified equation was more consistent with MC algorithm than with other formula, which could be a reference to those not accessible to the MC algorithm.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Nearly Optimal Lattice Simulation by Product Formulas

We consider simulating an n-qubit Hamiltonian with nearest-neighbor interactions evolving for time t on a quantum computer. We show that this simulation has gate complexity (nt)1+o(1) using product formulas, a straightforward approach that has been demonstrated by several experimental groups. While it is reasonable to expect this complexity—in particular, this was claimed without rigorous justification by Jordan, Lee, and Preskill—we are not aware of a straightforward proof. Furthermore, our approach is based on an analysis of the local error structure of product formulas, as introduced by Descombes and Thalhammer and significantly simplified here. We prove error bounds for canonical product formulas, which include well-known constructions such as the Lie-Trotter-Suzuki formulas. We also develop a local error representation for time-dependent Hamiltonian simulation, and we discuss generalizations to periodic boundary conditions, constant-range interactions, and higher dimensions. Combined with a previous lower bound, our result implies that product formulas can simulate lattice Hamiltonians with nearly optimal gate complexity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Precision of approximate Kramers formulas for the nuclear-fission rate in the presence of shell corrections

The precision of approximate Kramers formulas for the nuclear-fission rate in the presence of shell corrections in a quasistationary state is studied. The rates in question obtained analytically are compared with their numerical counterparts resulting from a dynamical simulation based on the Smoluchowski equation. It is shown that the integral Kramers formula is applicable even in the presence of shell corrections. As the excitation energy of the nucleus being considered decreases, the inaccuracy of this formula decreases monotonically, reaching 2% (inaccuracy of the dynamical simulation) in which case the dimensionless controlling parameter (entropy barrier) becomes greater than four. At small and intermediate values of this parameter, the accuracy of the integral Kramers formula is higher than the accuracy of the approximate Kramers formula obtained from it.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗