SEARCH · Engineering Papers

Results for “conservative interpolation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Optimal interpolation of global dissolved oxygen: 1965–2015

Oxygen inventory of the global ocean has declined in recent decades potentially due to the warming-induced reduction in solubility as well as the circulation and biogeochemical changes associated with ocean warming and increasing stratification. Earth System Models predict continued oxygen decline for this century with profound impacts on marine ecosystem and fisheries. Observational constraint on the rate of oxygen loss is crucial for assessing the ability of models to accurately simulate these changes. There are only a few observational assessments of the global oceanic oxygen inventory reporting a range of oxygen loss. This study develops a gridded data set of dissolved oxygen for the global oceans using optimal interpolation method. The resulting gridded product includes full-depth map of dissolved oxygen as 5-year moving average from 1965 to 2015 with uncertainty estimates. The uncertainty can come from unresolved small-scale and high-frequency variability and mapping errors. The multi-decadal trend of global dissolved oxygen is in the range of -281 to -373 Tmol/decade. This estimate is more conservative than previous works. In this study, the grid points far from the observations are essentially set equal to zero anomaly from the climatology. Calculating global inventory with this approach produces a relatively conservative estimate; thus, the results from this study likely provide a useful lower bound estimate of the global oxygen loss.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2021

An Efficient and Accurate Algorithm for Computing Grid-Averaged Solar Fluxes for Horizontally Inhomogeneous Clouds

A computationally efficient method is presented to account for the horizontal cloud inhomogeneity by using a radiatively equivalent plane parallel homogeneous (PPH) cloud. The algorithm can accurately match the calculations of the reference (rPPH) independent column approximation (ICA) results, but use only the same computational time required for a single plane parallel computation. The effective optical depth of this synthetic sPPH cloud is derived by exactly matching the direct transmission to that of the inhomogeneous ICA cloud. The ffective9 scattering asymmetry factor is found from a pre-calculated albedo inverse look-up-table that is allowed to vary over the range from -1.0 to 1.0. In the special cases of conservative scattering and total absorption, the synthetic method is exactly equivalent to the ICA, with only a small bias (about 0.2% in flux) relative to ICA due to imperfect interpolation in using the look-up tables. In principle, the ICA albedo can be approximated accurately regardless of cloud inhomogeneity. For a more complete comparison, the broadband shortwave albedo and transmission calculated from the synthetic sPPH cloud and averaged over all incident directions, have the RMS biases of 0.26% and 0.76%, respectively, for inhomogeneous clouds over a wide variation of particle size. The advantages of the synthetic PPH method are that (1) it is not required that all the cloud subcolumns have uniform microphysical characteristic, (2) it is applicable to any 1D radiative transfer scheme, and (3) it can handle arbitrary cloud optical depth distributions and an arbitrary number of cloud subcolumns with uniform computational efficiency.

cloud inhomogeneity↗

NASA NTRS2000

Hydrological Balance in Tropical Cyclones with Scatterometer and TRMM Data

Precipitation over oceans can be estimated from the radar and the microwave radiometer of the Tropical Rain Measuring Mission (TRMM). It can also be estimated from the divergence of the vertically integrated water vapor transport, through the conservation principle, assuming evaporation is relatively small. In tropical cyclones, the divergence of vertically integrated water vapor is highly dependent on the vertical transport and, therefore, on the wind divergence. Spaceborne scatterometers provide surface wind velocity and, therefore, surface wind divergence at spatial resolutions that are much higher than products of numerical weather prediction (NWP). In this study, ocean surface winds derived from the observations of space-based scatterometers and surface precipitation measured by TRMM were objectively interpolated to the same time and location during the passage of a tropical cyclone. Surface precipitation distribution was derived from wind and humidity profiles provided by NWP. When the surface level winds of NWP were replaced by the scatterometer winds, the surface precipitation patterns computed with the conservation method were found to be significantly changed and the new patterns are much closer in agreement with the patterns observed by TRMM.

Liu, W. Timothy↗

DOE OSTI2025

First Look at Quartic-in-Spin Binary Dynamics at Third Post-Minkowskian Order

We compute the conservative and radiation-reaction contributions to classical observables in the gravitational scattering between a spinning and a spinless black hole to the fourth order in spin and third order in the gravitational constant. The conservative results are obtained from two-loop amplitudes for the scattering process of a massive scalar with a massive spin-𝑠 field (𝑠=0, 1, 2) minimally coupled to gravity, employing the recently introduced spin interpolation method to resolve all spin-Casimir terms. The two-loop amplitude exhibits a spin-shift symmetry in both probe limits, which we conjecture to be a sign of yet unknown integrability of Kerr orbits through the quartic order in spin and to all orders in the gravitational constant. We obtain the radial action from the finite part of the amplitude and use it to compute classical observables, including the impulse and spin kick. This is done using the recently introduced covariant Dirac brackets, which allow for the computation of classical scattering observables for general (nonaligned) spin configurations. Finally, employing the radiation-reaction amplitude proposed by Alessio and Di Vecchia, together with the Dirac brackets, we obtain radiation-reaction contributions to observables at all orders in spin and beyond the aligned-spin limit. We find agreement with known results up to the quadratic order in spin for both conservative and radiation-reaction contributions. Our results advance the state of the art in the understanding of spinning binary dynamics in general relativity and demonstrate the power and simplicity of the Dirac bracket formalism for relating scattering amplitudes to classical observables.

general relativity↗

DOE OSTI2025

Spinning black hole scattering at $ \mathcal{O} $(G 3 S 2 )(G 3 S 2 ): Casimir terms, radial action and hidden symmetry

We resolve subtleties in calculating the post-Minksowskian dynamics of binary systems, as a spin expansion, from massive scattering amplitudes of fixed finite spin. In particular, the apparently ambiguous spin Casimir terms can be fully determined from the gradient of the spin-diagonal part of the amplitudes with respect to S 2 = −s(s+1)ħ 2 , using an interpolation between massive amplitudes with different spin representations. From two-loop amplitudes of spin-0 and spin-1 particles minimally coupled to gravity, we extract the spin Casimir terms in the conservative scattering angle between a spinless and a spinning black hole at $ \mathcal{O} $(G 3 S 2 ), finding agreement with known results in the literature. This completes an earlier study [Phys. Rev. Lett. 130 (2023), 021601] that calculated the non-Casimir terms from amplitudes. We also illustrate our methods using a model of spinning bodies in electrodynamics, finding agreement between scattering amplitude predictions and classical predictions in a root-Kerr electromagnetic background up to $ \mathcal{O} $(α 3 S 2 ). For both gravity and electrodynamics, the finite part of the amplitude coincides with the two-body radial action in the aligned spin limit, generalizing the amplitude-action relation beyond the spinless case. Surprisingly, the two-loop amplitude displays a hidden spin-shift symmetry in the probe limit, which was previously observed at one loop. We conjecture that the symmetry holds to all orders in the coupling constant and is a consequence of integrability of Kerr orbits in the probe limit at the first few orders in spin.

Classical Theories of Gravity↗

NASA NTRS1985

Numerical solution of the two-dimensional Euler equations by second-order upwind difference schemes

Two time-level, five-point explicit and implicit upwind difference schemes based on the characteristic flux difference splitting concept have been developed for the two-dimensional Euler equations. The method is conservative, second-order accurate in time and space, and general coordinate systems are used to treat complex geometries. Nonlinear flux limiters are employed to yield oscillation free sharp shock profiles. Upstream interpolation is used to yield a class of higher-order upwind schemes which closely mimic the locally one-dimensional method of characteristics (with fixed time intervals) through operator splitting. Numerical results have been obtained for a plane shock reflection and for flow over a circular arc in a channel. Characteristics of upwind TVD schemes, as applied to two-dimensional flows with embedded shocks are discussed.

Yang, J. Y.↗

NASA NTRS1985

Finite element algorithms for compressible flow computation on a supercomputer

Methods of applying computational fluid dynamics techniques to engineering problems are discussed. The linking factors between theoretical and applied research are the shape and weight functions which produce the Galerkin finite element schemes. The flow is unsteady, viscous, heat conducting and compressible, and steady-state flows are the asymptotic limit of unsteady flows. All flows are described by the time-averaged Navier-Stokes equations (NSE) with mass and energy conservation. Local curvilinear intrinsic coordinates are applied to discretize the NSE in arbitrary geometric domains. A modified weighted residuals approach defines the discrete analogs of the physical systems modeled and a general interpolants method is used to derive families of numerical models, both implicit and explicit, finite difference and finite elemnt, from a single point of departure. Time-dependent element approximations are achieved with a progressive assembly of generalized elements method. Hyperbolic steady-state Euler algorithms and quasi-parabolic are spatial marching algorithms for solving the discretized equations. Application of the techniques is illustrated through calculation of the compressible flow around a hypersonic flight vehicle at 100,000 ft at zero angle of attack at speeds of Mach 4-10.

Spradley, L. W.↗

DOE OSTI2022

Field validation of data-driven BSDF and peak extraction models for light-scattering fabric shades

Shading and daylighting systems affect cooling, heating, and lighting energy use by modulating solar radiation through the building façade. Characterizing shading systems holistically and accurately helps designers and engineers evaluate shading systems to achieve energy and non-energy performance goals. These complex fenestration systems can be modeled using Bidirectional Scattering Distribution Functions (BSDF), which map incident radiation to hemispherical distributions of outgoing radiation. Data-driven, tabulated BSDFs are derived from interpolated goniophotometer measured data, then sampled during the raytracing calculation. A peak extraction (PE) algorithm was developed to circumvent limits in BSDF angular resolution, where the specular peak is extracted during simulation by evaluating the BSDF in the through direction and surrounding region. The objective of this study was to validate this measurement and modeling workflow using field monitored data from a full scale testbed with eleven installed fabrics of different weaves, openness factors, and colors and assess the accuracy of the workflow under different adaptation and contrast conditions. Test conditions were limited to clear sky conditions with the sun in the field of view. Results showed that, for tensor tree datasets, vertical illuminance, solar luminance (2.5° apex), and daylight glare probability (DGP) were predicted to within a mean bias error (MBE) error of -456 lx (-12.3%), -3.46e5 (-38.4%), and -0.042 (-7.8%) when full PE occurred. With a binary classification of glare/ no glare, DGP was predicted accurately with a true positive rate of 0.98 and true negative rate of 1.0 using tensor tree data and less accurately with Klems BSDF data, particularly for cases of no glare. The workflow may be of insufficient accuracy to distinguish borderline performance between fabrics using the four-point glare scale, particularly under low adaptation, high contrast daylit conditions. Errors were due to reductions in peak shape and intensity across the BSDF interpolation and data reduction workflow. Future work is needed to better preserve measurement fidelity during interpolation and sampling, which in turn will improve PE performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1990

Triangle based TVD schemes for hyperbolic conservation laws

A triangle based total variation diminishing (TVD) scheme for the numerical approximation of hyperbolic conservation laws in two space dimensions is constructed. The novelty of the scheme lies in the nature of the preprocessing of the cell averaged data, which is accomplished via a nearest neighbor linear interpolation followed by a slope limiting procedures. Two such limiting procedures are suggested. The resulting method is considerably more simple than other triangle based non-oscillatory approximations which, like this scheme, approximate the flux up to second order accuracy. Numerical results for linear advection and Burgers' equation are presented.

Durlofsky, Louis J.↗

NASA NTRS1992

Triangle based adaptive stencils for the solution of hyperbolic conservation laws

Durlofsky, Louis J.↗

DOE OSTI2020

A coupled polygonal DEM-LBM technique based on an immersed boundary method and energy-conserving contact algorithm

Here we present a framework of coupling polygonal discrete elements and the lattice Boltzmann method using a direct forcing immersed boundary scheme. In this technique, an energy-conserving contact algorithm is utilized to handle the interactions between convex and concave polygonal particles. The surface of a polygon is represented by discrete boundary points which includes vertices of polygonal particles and/or points interpolated from vertices. The fluid-particle coupling is obtained through the interactions of the boundary points and the imaginary fluid particles using a direct-forcing immersed boundary method. Validations of the proposed technique are made by single particle and multiple arbitrarily-shaped particle sedimentation tests, and the effect of particle shape is illustrated using a drafting-kissing-tumbling benchmark.

42 ENGINEERING↗

NASA NTRS1986

A two-dimensional numerical study of the flow inside the combustion chambers of a motored rotary engine

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust.

Shih, T. I. P.↗

NASA NTRS1986

A two-dimensional numerical study of the flow inside the combustion chamber of a motored rotary engine

Shih, T. I-P.↗

DOE OSTI2021

A Constrained Data Assimilation Algorithm Based on GSI Hybrid 3D-EnVar and Its Application

Abstract Data assimilation (DA) at mesoscales is important for severe weather forecasts, yet the techniques of data assimilation at this scale remain a challenge. This study introduces dynamical constraints in the Gridpoint Statistical Interpolation (GSI) three-dimensional ensemble variational (3D-EnVar) data assimilation algorithm to enable the use of high-resolution surface observations of precipitation to improve atmospheric analysis at mesoscales. The constraints use the conservations of mass and moisture. Mass constraint suppresses the unphysical high-frequency oscillation, while moisture conservation constrains the atmospheric states to conform with the observed high-resolution precipitation. We show that the constrained data assimilation (CDA) algorithm significantly reduced the spurious residuals of the mass and moisture budgets compared to the original data assimilation (ODA). A case study is presented for a squall line over the Southern Great Plains on 20 May 2011 during Midlatitude Continental Convective Clouds Experiment (MC3E) of the Atmospheric Radiation Measurement (ARM) program by using ODA or CDA analysis as initial condition of forecasts. The state variables, and the location and intensity of the squall line are better simulated in the CDA experiment. Results show how surface observation of precipitation can be used to improve atmospheric analysis through data assimilation by using the dynamical constraints of mass and moisture conservations.

Wang, Jia↗

DOE OSTI2022

Downscaling atmospheric chemistry simulations with physically consistent deep learning

Abstract. Recent advances in deep convolutional neural network (CNN)-based super resolution can be used to downscale atmospheric chemistry simulations with substantially higher accuracy than conventional downscaling methods. This work both demonstrates the downscaling capabilities of modern CNN-based single image super resolution and video super-resolution schemes and develops modifications to these schemes to ensure they are appropriate for use with physical science data. The CNN-based video super-resolution schemes in particular incur only 39 % to 54 % of the grid-cell-level error of interpolation schemes and generate outputs with extremely realistic small-scale variability based on multiple perceptual quality metrics while performing a large (8×10) increase in resolution in the spatial dimensions. Methods are introduced to strictly enforce physical conservation laws within CNNs, perform large and asymmetric resolution changes between common model grid resolutions, account for non-uniform grid-cell areas, super-resolve lognormally distributed datasets, and leverage additional inputs such as high-resolution climatologies and model state variables. High-resolution chemistry simulations are critical for modeling regional air quality and for understanding future climate, and CNN-based downscaling has the potential to generate these high-resolution simulations and ensembles at a fraction of the computational cost.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1979

M2, S2, K1 models of the global ocean tide

Ocean tidal signals appear in many geophysical measurements. Geophysicists need realistic tidal models to aid in interpretation of their data. Because of the closeness to resonance of dissipationless ocean tides, it is difficult for numerical models to correctly represent the actual open ocean tide. As an approximate solution to this problem, test functions derived by solving Laplace's Tidal Equations with ocean loading and self gravitation are used as a basis for least squares dynamic interpolation of coastal and island tidal data for the constituents M2, S2, and Kl. The resulting representations of the global tide are stable over at least a ?5% variation in the mean depth of the model basin, and they conserve mass. Maps of the geocentric tide, the induced free space potential, the induced vertical component of the solid earth tide, and the induced vertical component of the gravitational field for each contituent are presented.

Parke, M. E.↗

DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗

DOE OSTI2024

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. Here, we use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗