Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Taylor expansion”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Uncertainty quantification for misspecified machine learned interatomic potentials

The use of high-dimensional regression techniques from machine learning has significantly improved the quantitative accuracy of interatomic potentials. Atomic simulations can now plausibly target quantitative predictions in a variety of settings, which has brought renewed interest in robust means to quantify uncertainties. In many practical settings where model complexity is constrained (e.g., due to performance considerations), misspecification — the inability of any one choice of model parameters to exactly match all training data — is a key contributor to errors that is often disregarded. Here, we employ a recent misspecification-aware regression technique to quantify parameter uncertainties, which is then propagated to a broad range of phase and defect properties in tungsten. The propagation is performed through both brute-force resampling and implicit Taylor expansion. The propagated misspecification uncertainties robustly quantify and bound errors on a broad range of material properties. We demonstrate application to recent foundational machine learning interatomic potentials, accurately predicting and bounding errors in MACE-MPA-0 energy predictions across the diverse materials project database.

36 MATERIALS SCIENCE↗
DOE OSTI2024

QCD material parameters at zero and non-zero chemical potential from the lattice

Using an eighth-order Taylor expansion in baryon-chemical potential, we recently obtained the (2+1) flavor QCD equation of state (EoS) at non zero conserved charge chemical potentials from the lattice. We focused on strangeness-neutral, isospin-symmetric QCD matter, which closely resembles the situation encountered in heavy-ion collision experiments. Using this EoS, we present here results on various QCD material parameters; in particular we compute the specific heat, speed of sound, and compressibility along appropriate lines of constant physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗
DOE OSTI2024

Compiled properties of nucleonic matter and nuclear and neutron star models from nonrelativistic and relativistic interactions

Here, this paper compiles the model parameters and zero-temperature properties of an extensive collection of published theoretical nuclear interactions, including 255 nonrelativistic (Skyrme-like) forces, 270 relativistic mean field (RMF) and point-coupling (RMF-PC) forces, and 13 Gogny-like forces. This forms the most exhaustive tabulation of model parameters to date. The properties of uniform symmetric matter and pure neutron matter at the saturation density are determined. Symmetry properties found from the second-order term of a Taylor expansion in neutron excess are compared with the energy difference of pure neutron and symmetric nuclear matter at the saturation density. Selected liquid-droplet model parameters, including the surface tension and surface symmetry energy, are determined for semi-infinite surfaces. Theoretical liquid droplet model neutron skin thicknesses of the neutron-rich closed-shell nuclei 48 Ca and 208 Pb are compared with published theoretical Hartree-Fock and experimental results. A similar comparison is made between theoretical liquid-droplet and experimental values of dipole polarizabilities. In addition, radii, binding energies, moments of inertia and tidal deformabilities of 1.2⁢M ⊙ , 1.4⁢M ⊙ and 1.6⁢M ⊙ neutron stars are computed. An extensive correlation analysis of bulk matter, nuclear structure, and low-mass neutron star properties is performed and compared with nuclear experiments and astrophysical observations.

nuclear astrophysics↗
DOE OSTI2026

Generalized parton distributions from symbolic regression

AI/ML informed Symbolic Regression is the next stage of scientific modeling. We utilize a highly customizable symbolic regression package “PySR” to model the x and t dependence of the flavor isovector combination Hu-d(x, t, ) at = 0. These PySR models were trained on GPD results provided by both Lattice QCD and phenomenological sources GGL, GK, and VGG. We demonstrate, for the first time, the consistency and systematic convergence of Symbolic Regression by quantifying the disparate models through their Taylor expansion coefficients. In addition to PySR penalizing models with higher complexity and mean-squared error, we implement schemes that test specific physics hypotheses, including force-factorized x and t dependence and Regge behavior in PySR GPDs. We show that PySR can identify factorizing GPD sources based on their response to the Force-Factorized model. Here, knowing the precise behavior of the GPDs, and their uncertainties in a wide range in x and t, crucially impacts our ability to concretely and quantitatively predict hadronic spatial distributions and their derived quantities.

Dotson, Andrew [New Mexico State University, Las C↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
DOE OSTI2024

cppTPSA/pyTPSA: a C++/Python package for truncated power series algebra

The truncated power series algebra (TPSA), also referred to as differential algebra (DA), is a well-established and widely used method in particle accelerator physics and astronomy. The most straightforward usage of TPSA/DA is to calculate the Taylor expansion of a given function at a specific point up to order ?. In recent years, as the application of TPSA/TA has been extended to other fields, a reusable implementation of TPSA/DA as a modern C++ library or other high level programming language like Python has become desirable. The cppTPSA package implements TPSA/DA in C++11 and provides developers a convenient library with which to build advanced TPSA/DA-based methods. A Python 3 library, pyTPSA, has also been developed based on the C++ lib.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

Stability of nonuniform rotor blades in hover using a mixed formulation

A mixed formulation for calculating static equilibrium and stability eigenvalues of nonuniform rotor blades in hover is presented. The static equilibrium equations are nonlinear and are solved by an accurate and efficient collocation method. The linearized perturbation equations are solved by a one step, second order integration scheme. The numerical results correlate very well with published results from a nearly identical stability analysis based on a displacement formulation. Slight differences in the results are traced to terms in the equations that relate moments to derivatives of rotations. With the present ordering scheme, in which terms of the order of squares of rotations are neglected with respect to unity, it is not possible to achieve completely equivalent models based on mixed and displacement formulations. The one step methods reveal that a second order Taylor expansion is necessary to achieve good convergence for nonuniform rotating blades. Numerical results for a hypothetical nonuniform blade, including the nonlinear static equilibrium solution, were obtained with no more effort or computer time than that required for a uniform blade.

Stephens, W. B.↗
NASA NTRS1982

A second-order Budkyo-type parameterization of landsurface hydrology

A simple, second order parameterization of the water fluxes at a land surface for use as the appropriate boundary condition in general circulation models of the global atmosphere was developed. The derived parameterization incorporates the high nonlinearities in the relationship between the near surface soil moisture and the evaporation, runoff and percolation fluxes. Based on the one dimensional statistical dynamic derivation of the annual water balance, it makes the transition to short term prediction of the moisture fluxes, through a Taylor expansion around the average annual soil moisture. A comparison of the suggested parameterization is made with other existing techniques and available measurements. A thermodynamic coupling is applied in order to obtain estimations of the surface ground temperature.

Andreou, S. A.↗
NASA NTRS1984

The solution of non-linear hyperbolic equation systems by the finite element method

A finite-element method for the solution of nonlinear hyperbolic systems of equations, such as those encountered in non-self-adjoint problems of transient phenomena in convection-diffusion or in the mixed representation of wave problems, is developed and demonstrated. The problem is rewritten in moving coordinates and reinterpolated to the original mesh by a Taylor expansion prior to a standard Galerkin spatial discretization, and it is shown that this procedure is equivalent to the time-discretization approach of Donea (1984). Numerical results for sample problems are presented graphically, including such shallow-water problems as the breaking of a dam, the shoaling of a wave, and the outflow of a river; compressible flows such as the isothermal flow in a nozzle and the Riemann shock-tube problem; and the two-dimensional scalar-advection, nonlinear-shallow-water, and Euler equations.

Loehner, R.↗
NASA NTRS1986

Embedded mesh solution of the 2-D Euler equations - Evaluation of interface formulations

Solution of the steady 2-D Euler equations using mesh embedding, or local grid refinement, with a cell-centered finite volume scheme is investigated. Embedded regions which are topologically similar to the global grid are considered. An isoenergetic model for the governing equations is used in Jameson's finite volume multistage scheme with modifications to the boundary conditions and smoothing. A detailed study of the embedding interface flux and smoothing formulations is conducted. Taylor expansion analysis reveals that local second order spatial accuracy is not possible if a conservative interface flux formulation is used. The analysis also gives constraints for local first order accuracy. An energy stability analysis indicates that downwind weighting of interface fluxes causes local instabilities. Analysis shows that conservative interface smoothing formulations must have a locally convective component, but that correct interface formulations allow globally dissipative smoothing. Embedded mesh solutions obtained with this scheme are presented for a transonic airfoil. They show that if embedding interfaces are close to the shocks, then small modifications in the interface location can have large effects on converge and solution accuracy.

Allmaras, S. R.↗
NASA NTRS1988

Application of the Thurston bifurcation solution strategy to problems with modal interaction

The solution of bifurcation problems with closely-spaced critical points is achieved by first separating the singular part of the equation system encountered during a Newton iteration and carrying the Taylor expansion of the reduced system out to higher order. This separation is accomplished by transforming the equation system into an equivalent system in which some of the original unknowns are replaced with an equal number of modal amplitude coefficients. This method was used to continue the analysis of two significant example problems well past multiple fiburcation points, allowing a detailed examination of postbuckling behavior in the presence of modal interaction.

Rankin, C. C.↗
NASA NTRS1988

Tunneling near the peaks of potential barriers - Consequences of higher-order Wentzel-Kramers-Brillouin corrections

The quantum-mechanical tunneling transmission coefficient is derived for energies near the peak of a general potential barrier, to fourth order in the WKB approximation, using a method based on an eighth-order Taylor expansion of the potential near the peak. The result agrees with contour-integral formulations of WKB tunneling. For the Poeschl-Teller potential the WKB series converges rapidly to the known exact transmission coefficient. Higher-order WKB corrections may be important in attempts to determine the internuclear potential by inversion using low-energy fusion cross-section data.

Will, Clifford M.↗
NASA NTRS1992

The free and forced vibrations of structures using the finite dynamic element method

In addition to an extensive review of the literature on exact and corrective displacement based methods of vibration analysis, a few theorems are proven concerning the various structural matrices involved in such analyses. In particular, the consistent mass matrix and the quasi-static mass matrix are shown to be equivalent, in the sense that the terms in their respective Taylor expansions are proportional to one another, and that they both lead to the same dynamic stiffness matrix when used with the appropriate stiffness matrix.

Fergusson, Neil J.↗
NASA NTRS1993

On the assumption of vanishing temperature fluctuations at the wall for heat transfer modeling

Boundary conditions for fluctuating wall temperature are required for near-wall heat transfer modeling. However, their correct specifications for arbitrary thermal boundary conditions are not clear. The conventional approach is to assume zero fluctuating wall temperature or zero gradient for the temperature variance at the wall. These are idealized specifications and the latter condition could lead to an ill posed problem for fully-developed pipe and channel flows. In this paper, the validity and extent of the zero fluctuating wall temperature condition for heat transfer calculations is examined. The approach taken is to assume a Taylor expansion in the wall normal coordinate for the fluctuating temperature that is general enough to account for both zero and non-zero value at the wall. Turbulent conductivity is calculated from the temperature variance and its dissipation rate. Heat transfer calculations assuming both zero and non-zero fluctuating wall temperature reveal that the zero fluctuating wall temperature assumption is in general valid. The effects of non-zero fluctuating wall temperature are limited only to a very small region near the wall.

Sommer, T. P.↗
NASA NTRS1993

Convection effects on radial segregation and crystal melt interface in vertical Bridgman growth

We analytically study the influence of convection caused by horizontal heat transfer through the sides of a vertical Bridgman apparatus. We consider the case when the heat transfer across the side walls is small so that the resulting interfacial deformation and fluid velocities are also small. This allows us to linearize the Navier-Stokes equations and express the interfacial conditions about a planar interface through a Taylor expansion. Using a no tangential stress conditions on the side walls, asymptotic expressions for both the interfacial slope, and radial segregation at the crystal-melt interface are obtained in closed form in the limit of large thermal Rayleigh number. It is suggested that these can be reduced by appropriately controlling a specific heat transfer property at the edge of the insulation zone in the solid side.

Tanveer, S.↗