Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “higher-order calculations”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

A critical study of higher-order numerical methods for solving the boundary-layer equations

A fourth-order box method is presented for calculating numerical solutions to parabolic, partial differential equations in two variables or ordinary differential equations. The method is the natural extension of the second-order Keller Box Scheme to fourth order and is demonstrated with application to the incompressible, laminar and turbulent boundary-layer equations. The efficiency of the present method is compared with other two-point and three-point higher-order methods; namely, the Keller Box Scheme with Richardson extrapolation, the method of deferred corrections, and the three-point spline methods. For equivalent accuracy, numerical results show the present method to be more efficient than the other higher-order methods for both laminar and turbulent flows.

Wornom, S. F.↗
DOE OSTI2026

Gradient-informed Hamiltonian Monte Carlo for multicomponent CALPHAD model optimization and uncertainty quantification

CALPHAD model parameter optimization is inherently challenging due to non-smooth objective functions, high-dimensional parameter spaces, and the need for uncertainty quantification (UQ). Traditional weighted nonlinear least squares approaches are computationally efficient but local, whereas black-box global optimizers and ensemble Markov Chain Monte Carlo (MCMC) methods provide broader exploration at substantial computational cost. The objective of this work is to combine the global exploration capability of gradient-informed Hamiltonian Monte Carlo – specifically the No-U-Turn Sampler (NUTS) – with local deterministic refinement using BFGS to efficiently optimize multicomponent CALPHAD models with minimal manual intervention. Analytic gradients are computed via the Jansson derivative framework. The methodology is demonstrated on the Cr—Fe binary system and extended to the Cr—Fe—Ni ternary system with 32 degrees of freedom. For Cr—Fe, NUTS achieves comparable or superior optimality relative to ensemble MCMC while requiring over an order-of-magnitude fewer likelihood evaluations. Parameter uncertainties are quantified through NUTS sampling and propagated to thermodynamic observables using local expansion, demonstrating a novel modular approach that combines binary and ternary parameter subsets without requiring global relaxation. These results establish gradient-informed exploration as a scalable strategy for multicomponent CALPHAD optimization and provide a practical route towards efficient higher-order database development with quantified uncertainty.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct observation of key aluminum hydroxide prenucleation oligomers for gibbsite nucleation and crystallization in sodium aluminate solution by liquid ToF-SIMS

The mechanism of gibbsite (aluminum hydroxide) crystallization from highly alkaline solutions such as Bayer liquors remains poorly understood, where aluminum (Al) transforms from largely tetrahedrally coordinated aluminate monomers in sodium aluminate solutions into a network of octahedra in gibbsite crystals. A variety of traditional analytical approaches applied to this system do not readily reveal the presence of higher-order oligomeric intermediates. To overcome this limitation, we employed in-situ liquid Time-of-Flight Secondary Ion Mass Spectrometry (ToF-SIMS) to examine the Al species present in concentrated sodium aluminate solutions favoring crystallization of gibbsite or sodium aluminate. A complex mixture of Al oligomers was found. By comparing the change in the relative concentration of Al oligomers with +1 and -1 charge, we were able to identify three major Al oligomer candidates, iso-tetramers, iso-pentamers, and cyclic-hexamers, for the nucleation and crystallization of gibbsite. The concentrations of iso-tetramers and iso-pentamers significantly surpass those of cyclic-hexamers. Time-dependent in-situ Raman spectroscopy analysis indicated that the appearance of gibbsite coincided with the peak concentration of these oligomers. The Density-functional theory (DFT) calculation suggests that the formation of iso-oligomers is more favorable than that of cyclic-hexamers. The combined results suggest that iso-tetramers and iso-pentamers play the most substantial role in the nucleation and growth of gibbsite in the sodium aluminate solutions. Our findings also suggest that the oligomers that promote gibbsite crystallization are more stable in dilute sodium aluminate solutions, making these solutions particularly suitable for efficient gibbsite crystallization. In conclusion, our study fills a major knowledge gap in understanding Al speciation that leads to the nucleation and crystallization of gibbsite in concentrated sodium aluminate solutions.

Bayer liquor↗
NASA NTRS1977

Use of finite element techniques in the determination of the acoustic properties of turbofan inlets

This paper describes the application of the finite-element method in combination with Galerkin's method in the determination of the acoustic properties of turbofan inlets containing high-subsonic-Mach-number flows. An approximate solution to the steady inviscid flow field is obtained using an integral method for calculating the potential flow field in the inlet with a correction to account for compressibility effects. The accuracy of the finite-element technique in predicting the acoustic properties of annular ducts has been checked by comparison with available analytical solutions for the problems of planeand spinning-wave propagation through a hard-walled annular duct with a constant mean flow. Results are presented comparing low-frequency plane-wave propagation through a hard-walled turbofan inlet containing a one-dimensional flow with the same inlet containing a fully two-dimensional axisymmetric mean flow. It is shown that when one-dimensional mean flow is assumed to exist in the duct, the plane wave propagates with relatively little distortion. However, propagation of a plane wave through the fully two-dimensional flow field in the inlet produces severe distortion due to the excitation of higher-order modes.

Sigman, R. K.↗
DOE OSTI2024

Three-dimensional higher-order saddle-point-induced flatbands in Co-based kagome metals

The saddle point (Van Hove singularity) exhibits a divergent density of states in two-dimensional systems, leading to fascinating phenomena such as strong correlations and unconventional superconductivity, yet it is seldom observed in three-dimensional (3D) systems. In this work we find two types of 3D higher-order saddle points (HOSPs) in emerging 3D kagome metals YbCo 6⁢ Ge 6 and MgCo 6 ⁢Ge 6 . Both HOSPs exhibit a singularity in their density of states, which is significantly enhanced compared to the ordinary saddle point. The HOSP near the Fermi energy generates a flatband extending a large area in the Brillouin zone, potentially amplifying the correlation effect and fostering electronic instabilities. Two types of HOSPs exhibit distinct robustness upon element substitution and lattice distortions in these kagome compounds. Our work paves the way for engineering exotic band structures, such as saddle points and flatbands, and exploring interesting phenomena in Co-based kagome materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

EFT approach to the endpoint of muon decay-in-orbit

As upcoming experiments aim to probe muon conversion with unprecedented precision, equally precise theoretical predictions are crucial to maximize discovery potential. This applies not only to the new physics signal, muon-electron conversion, but also to its only irreducible background, muon decay-in-orbit (DIO) near the endpoint. Accurate computation of higher-order corrections in bound states is a long-standing challenge due to the difficulty of systematically organizing contributions. In previous work, we developed an Effective Field Theory framework to address this issue and applied it to muon conversion. Here, we extend this approach to the DIO endpoint, a more complex problem due to the presence of a neutrino-antineutrino pair in the final state. We present the most precise prediction to date of the background spectrum relevant for future muon conversion searches, achieving next-to-leading logarithmic prime accuracy for QED corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Lattice anharmonicity effects in fluorite oxide single crystals and anomalous increase in phonon lifetime in ceria at elevated temperature

Here, we investigate the temperature dependence of the frequency and linewidth of the triply degenerate T 2g zone-centered optical phonon in flux-grown ceria and hydrothermally synthesized thoria single crystals from room temperature to 1273 K using Raman spectroscopy. Both crystals exhibit an expected increase in the phonon linewidth with temperature due to enhanced phonon–phonon scattering. However, ceria displays an anomalous linewidth reduction in the temperature range of 1023–1123 K. First-principles phonon linewidth calculations considering cubic and quartic phonon interactions within temperature-independent phonon dispersion fail to describe this anomaly. A parameterization of the temperature-dependent second-order interatomic force constants based on previously reported phonon dispersion measured at room and high temperatures predicts a deviation from the monotonic linewidth increase, albeit at temperatures lower than those observed experimentally for ceria. The qualitative agreement in the trend of temperature-dependent linewidth suggests that lattice anharmonicity-induced phonon renormalization plays a role in phonon lifetime. Specifically, a change in the overlap between softened acoustic and optical branches in the dispersion curve reduces the available phonon scattering phase space of the Raman-active mode at the zone center, leading to an increased phonon lifetime within a narrow temperature interval. These findings provide insights into higher-order anharmonic interactions in ceria and thoria, motivating further investigations into the role of anharmonicity-induced phonon renormalization on phonon lifetimes at high temperatures.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Erratum to: Measurements of higher-order cumulants of multiplicity and net-electric charge distributions in inelastic proton-proton interactions by NA61/SHINE

This Erratum replaces, due to a discovery of coding mistakes, the following quantities: κ 3 /κ 1 of the h + − h − distribution presented in Fig. 6 and Table 6, κ 4 listed in Table 4, and $\hat{C}$ 4 presented in Fig. 7 and Table 5. All mentioned figures and tables were updated.

High-Energy Particle Collision Data Analysis↗
DOE OSTI2025

Higher-order symmetry plane correlations in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

The correlations between event-by-event fluctuations of symmetry planes are measured in Pb-Pb collisions at a center-of-mass energy per nucleon pair $\sqrt{s_{NN}}$ = 5.02 TeV recorded by the ALICE detector at the Large Hadron Collider. This analysis is conducted using the Gaussian estimator technique, which is insensitive to biases from correlations between different flow amplitudes. The study presents, for the first time, the centrality dependence of correlations involving up to five different symmetry planes. The correlation strength varies depending on the harmonic order of the symmetry plane and the collision centrality. Comparisons with measurements from lower energies indicate no significant differences within uncertainties. Additionally, the results are compared with hydrodynamic model calculations. Although the model predictions provide a qualitative explanation of the experimental results, they overestimate the data for some observables. This is particularly true for correlators that are sensitive to the nonlinear response of the medium to initial-state anisotropies in the collision system. As these new correlators provide unique information—independent of flow amplitudes—their usage in future model developments can further constrain the properties of the strongly interacting matter created in ultrarelativistic heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1978

Two-photon excitation of atomic oxygen

A standard perturbation expansion in the atom-radiation-field interaction is used to calculate the two-photon excitation cross section for the 1s2 2s2 2p4 3p to 1s2 2s2 3(4S)3p 3P transition in atomic oxygen. The lowest-order summation over bound and continuum intermediate states is handled by solving the equivalent inhomogeneous differential equation. Higher-order electron-correlation effects are also examined.

Pindzola, M. S.↗
NASA NTRS1978

An extension of MacCormack's method for flows with higher-order equations and in different configurations

The numerical scheme for the computation of a shock discontinuity developed by MacCormack has been extended to solve a number of differential equations, including cases explicitly containing higher-order derivatives: (1) Korteweg-de Vries equation with a term of third-order derivative, (2) a system of nonlinear equations governing nonsteady one-dimensional plasma flow in cylindrical coordinate, (3) equations of solar wind. Comparisons with previous results are made, if available, to illustrate the advantages of the present method. The question of convergence of the numerical calculation is discussed.

Ying, S. J.↗
NASA NTRS1979

Higher-order skewness and excess coefficients of some probability distributions applicable to optical propagation phenomena

Expressions are derived for higher-order skewness and excess coefficients using central moments and cumulants up to 8th order. These coefficients are then calculated for three probability distributions: (1) Log-normal, (2) Rice-Nakagami, and (3) Gamma distributions. Curves are given to shown the variation of skewness with excess coefficients for these distributions. These curves are independent of the particular distribution parameters. This method is useful for studying fluctuating phenomena, which obey non-Gaussian statistics.

Majumdar, A. K.↗
NASA NTRS1981

On various treatments of potential equations at shocks

The potential equation that governs transonic inviscid flowfields is elliptic in subsonic regions and hyperbolic in supersonic regions. The transition from supersonic to subsonic flow may occur either continuously or across a surface of discontinuity known as a shock. Because of the changing type of governing equation and the presence of the surfaces of discontinuity, special treatments of the governing equations are needed at shocks to preserve mass conservation and to ensure satisfaction of the governing equation. In the so-called shock-capturing techniques, shocks are captured by smearing the discontinuity over several mesh cells in the computational domain. Most existing conservative schemes conserve mass flux throughout the flowfield, including shock boundaries; however, they also introduce zero-order errors in the approximation to the governing equation. It is shown that the zero-order errors, which do not diminish as the mesh spacing approaches zero, can cause discrepancies in the prediction of shock strength and location. These parameters are extremely sensitive features of transonic inviscid flowfield calculations. It is also shown that the zero-order errors can be avoided by developing higher-order schemes that properly model the governing equation at shocks.

Chen, L. T.↗
NASA NTRS1983

Higher-order derivative correlations and the alignment of small-scale structures in isotropic numerical turbulence

In a three dimensional simulation higher order derivative correlations, including skewness and flatness factors, are calculated for velocity and passive scalar fields and are compared with structures in the flow. The equations are forced to maintain steady state turbulence and collect statistics. It is found that the scalar derivative flatness increases much faster with Reynolds number than the velocity derivative flatness, and the velocity and mixed derivative skewness do not increase with Reynolds number. Separate exponents are found for the various fourth order velocity derivative correlations, with the vorticity flatness exponent the largest. Three dimensional graphics show strong alignment between the vorticity, rate of strain, and scalar-gradient fields. The vorticity is concentrated in tubes with the scalar gradient and the largest principal rate of strain aligned perpendicular to the tubes. Velocity spectra, in Kolmogorov variables, collapse to a single curve and a short minus 5/3 spectral regime is observed.

Kerr, R. A.↗
NASA NTRS1984

The application of a second generation low-order panel method - program 'Vsaero' - to powerplant installation studies

Results from the application of a second generation low-order panel method to a number of powerplant installations are presented. The method, which retains the advantages of ease of use and low operating cost of earlier panel methods while retaining the aerodynamic modelling rigor of higher-order methods was used to predict the flow in and around typical nacelles and aircraft for a wide range of conditions. The ability of the second generation low-order analysis to calculate internal flows without the leakage problems associated with earlier programs is demonstrated and correlation between calculated and measured surface pressures and flow behavior, including predicting streamlines and locating regions of separated flow, is presented.

Clark, D. R.↗
NASA NTRS1986

Group-kinetic theory of turbulence

The two phases are governed by two coupled systems of Navier-Stokes equations. The couplings are nonlinear. These equations describe the microdynamical state of turbulence, and are transformed into a master equation. By scaling, a kinetic hierarchy is generated in the form of groups, representing the spectral evolution, the diffusivity and the relaxation. The loss of memory in formulating the relaxation yields the closure. The network of sub-distributions that participates in the relaxation is simulated by a self-consistent porous medium, so that the average effect on the diffusivity is to make it approach equilibrium. The kinetic equation of turbulence is derived. The method of moments reverts it to the continuum. The equation of spectral evolution is obtained and the transport properties are calculated. In inertia turbulence, the Kolmogoroff law for weak coupling and the spectrum for the strong coupling are found. As the fluid analog, the nonlinear Schrodinger equation has a driving force in the form of emission of solitons by velocity fluctuations, and is used to describe the microdynamical state of turbulence. In order for the emission together with the modulation to participate in the transport processes, the non-homogeneous Schrodinger equation is transformed into a homogeneous master equation. By group-scaling, the master equation is decomposed into a system of transport equations, replacing the Bogoliubov system of equations of many-particle distributions. It is in the relaxation that the memory is lost when the ensemble of higher-order distributions is simulated by an effective porous medium. The closure is thus found. The kinetic equation is derived and transformed into the equation of spectral flow.

Tchen, C. M.↗
NASA NTRS1942

On the use of residue theory for treating the subsonic flow of a compressible fluid

A new mathematical technique, due to Milne-Thomson, is used to obtain an improved form of the method of Poggi for calculating the effect of compressibility on the subsonic flow past an obstacle. By means of this new method, the difficult surface integrals of the original Poggi method can be replaced by line integrals. These line integrals are then solved by the use of residue theory. In this way an equation is obtained giving the second-order effect of compressibility on the velocity of the fluid. The method is practicable for obtaining the higher-order effects of compressibility on the velocity field. As an illustration of the general result, the flow past an elliptic cylinder is discussed.

Kaplan, Carl↗