SEARCH · Engineering Papers

Results for “algebraic approach”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Slip velocity boundary conditions for the lattice Boltzmann modeling of microchannel flows

Abstract Slip flows in ducts are important in numerous engineering applications, most notably in microchannel flows. Compared to the standard no‐slip Dirichlet condition, the case of slip formulates as a Robin‐type condition for the fluid tangential velocity. Such an increase in mathematical complexity is accompanied by a more challenging numerical transcription. The present work concerns with this topic, addressing the modeling of the slip velocity boundary condition in the lattice Boltzmann method (LBM) applied to steady slow viscous flows inside ducts of nontrivial shapes. As novelty, we extend the newly revised local second‐order boundary (LSOB) Dirichlet fluid flow method [ Philos. Trans. R. Soc. A 378, 20190404 (2020)] to implement the slip velocity condition within the two‐relaxation‐time (TRT) framework. The LSOB follows an in‐node philosophy where its operation principle seeks to explicitly reconstruct the unknown boundary populations in the form of a third‐order accurate Chapman–Enskog expansion, where the wall slip condition is built‐in as a normal Taylor‐type condition. The key point of this approach is that the required first‐ and second‐order momentum derivatives, rather than computed through nonlocal finite difference approximations, are locally determined through a simple local linear algebra procedure, whose formulation is particularly aided by the TRT symmetry argument. To express the obtained derivatives, two approaches are considered, called and , which operate with node and wall variables, respectively. These two formulations are developed to prescribe the physical slip condition over plane and curved walls, including the corners. Their consistency and accuracy characteristics are examined against alternative linkwise strategies to impose the wall slip velocity, such as the kinetic‐based diffusive bounce‐back scheme, the central linear interpolation slip scheme, and the multireflection slip scheme. The several slip schemes are tested over different 3D microchannel configurations, with walls not conforming with the LBM uniform mesh. Numerical tests confirm the advanced accuracy characteristics of the proposed LSOB slip boundary scheme, revealing the added challenge of the wall slip modeling, and that parabolic accuracy is a necessary requirement to reach second‐order accuracy within this problem class.

Silva, Goncalo↗

NASA NTRS1974

The modified equation approach to the stability and accuracy analysis of finite-difference methods

The stability and accuracy of finite-difference approximations to simple linear partial differential equations are analyzed by studying the modified partial differential equation. Aside from round-off error, the modified equation represents the actual partial differential equation solved when a numerical solution is computed using a finite-difference equation. The modified equation is derived by first expanding each term of a difference scheme in a Taylor series and then eliminating time derivatives higher than first order by certain algebraic manipulations. The connection between 'heuristic' stability theory based on the modified equation approach and the von Neumann (Fourier) method is established. In addition to the determination of necessary and sufficient conditions for computational stability, a truncated version of the modified equation can be used to gain insight into the nature of both dissipative and dispersive errors.

Warming, R. F.↗

NASA NTRS2021

Sulfur Molecules in Space by X-rays: A Computational Study

X-ray astronomy lacks high resolution spectra of interstellar dust analogues and molecules, severely hampering interstellar medium studies based on upcoming X-ray missions. Various theoretical approaches may be used to address this problem, but they must first be shown to reproduce reliable spectra compared to the experiment. In this work, we calculate the sulfur Kedge X-ray absorption spectra of H2S, SO2, and OCS, whose spectra are already known from X-ray experiments and predict the X-ray spectrum of CS, which as far as we are aware has not been measured, thereby hampering its detection by X-ray telescopes. We chose these four molecules as the astrochemistry of sulfur is an unsolved problem and as the four molecules are already known to exist in space. We consider three types of methods for modeling the X-ray spectra: more accurate calculations with the algebraic-diagrammatic construction (ADC) and the CC2, CCSD, and CC3coupled cluster (CC) approaches as well as more affordable ones with transition potential density functional theory (TP-DFT). A comparison of our computational results to previously reported experimental spectra shows that the core−valence separation (CVS)approaches CVS-ADC(2)-x and CVS-CC3 generally yield a good qualitative level of agreement with the experiment, suggesting that they can be used for interpreting measured spectra, while the TP-DFT method is not reliable for these molecules. However, quantitative agreement with the experiment is still outside the reach of the computational methods studied in this work.

Goranka Bilalbegovic↗

NASA NTRS1987

An algebraic criterion for the onset of chaos in nonlinear dynamic systems

The correspondence between iterated integrals and a noncommutative algebra is used to recast the given dynamical system from the time domain to the Laplace-Borel transform domain. It is then shown that the following algebraic criterion has to be satisfied for the outset of chaos: the limit (as tau approaches infinity and x sub 0 approaches infinity) of ((sigma(k=0) (tau sup k) / (k* x sub 0 sup k)) G II G = 0, where G is the generating power series of the trajectories, the symbol II is the shuffle product (le melange) of the noncommutative algebra, x sub 0 is a noncommutative variable, and tau is the correlation parameter. In the given equation, symbolic forms for both G and II can be obtained by use of one of the currently available symbolic languages such as PLI, REDUCE, and MACSYMA. Hence, the criterion is a computer-algebraic one.

Unal, A.↗

DOE OSTI2022

Generalizations of the double-copy: the KLT bootstrap

We formulate a new program to generalize the double-copy of tree amplitudes. The approach exploits the link between the identity element of the “KLT algebra” and the KLT kernel, and we demonstrate how this leads to a set of KLT bootstrap equations that the double-copy kernel has to satisfy in addition to locality constraints. We solve the KLT bootstrap equations perturbatively to find the most general higher-derivative corrections to the 4- and 5-point field theory KLT kernel. The new kernel generalizes the string KLT kernel and its associated monodromy relations. It admits new color-structures in the effective theories it double-copies. It provides distinct generalized KK and BCJ relations for the left and right single-color theories and is in that sense a ‘heterotic’-type double-copy. We illustrate the generalized double-copy in detail for 4d Yang-Mills theory with higher-derivative corrections that produce dilaton-axion-gravity with local operators up order ∇ 10 R 4 . Finally, we initiate a search for new double-copy kernels.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Tensor Decompositions for Count Data that Leverage Stochastic and Deterministic Optimization

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the global maximum likelihood estimator from local minima. Simultaneously, a recent trend in theoretical computer science and numerical linear algebra leverages randomization to solve very large, hard problems. The typical approach is to use randomization for a fast approximation and determinism for refinement to yield effective algorithms with theoretical guarantees. Two popular algorithms for Poisson CPD reflect that emergent dichotomy: CP Alternating Poisson Regression is a deterministic algorithm and Generalized Canonical Polyadic decomposition makes use of stochastic algorithms in several variants. This work extends recent work to develop two new methods that leverage randomized and deterministic algorithms for improved accuracy and performance.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1972

Elimination on sparse symmetric systems of a special structure.

Consideration of the problem of finding a permutation of rows and columns and an algorithm for solving ordered systems of linear algebraic equations with sparse matrices having a certain regular structure. Two approaches to the solution of this problem, in which the sparsity is used to some extent, are outlined. One of them is a very general approach where optimal (or nearly optimal) ordering is sought and the algorithm for solving the ordered system treats the matrix element by element to perform only necessary operations. The other approach involves the use of band matrices. After comparing these two approaches, a third approach is then suggested which involves the use of pipe matrices, and a means of ordering the rows and columns to obtain this type of matrix is presented. Examples of matrices reordered by the proposed procedure are cited.

Segethova, J.↗

NASA NTRS2002

Anytime query-tuned kernel machine classifiers via Cholesky factorization

We recently demonstrated 2 to 64-fold query-time speedups of Support Vector Machine and Kernel Fisher classifiers via a new computational geometry method for anytime output bounds (DeCoste,2002). This new paper refines our approach in two key ways. First, we introduce a simple linear algebra formulation based on Cholesky factorization, yielding simpler equations and lower computational overhead. Second, this new formulation suggests new methods for achieving additional speedups, including tuning on query samples. We demonstrate effectiveness on benchmark datasets.

Kernel Cholesky MNIST↗

NASA NTRS2016

Comparison of Turbulent Thermal Diffusivity and Scalar Variance Models

In this study, several variable turbulent Prandtl number formulations are examined for boundary layers, pipe flow, and axisymmetric jets. The model formulations include simple algebraic relations between the thermal diffusivity and turbulent viscosity as well as more complex models that solve transport equations for the thermal variance and its dissipation rate. Results are compared with available data for wall heat transfer and profile measurements of mean temperature, the root-mean-square (RMS) fluctuating temperature, turbulent heat flux and turbulent Prandtl number. For wall-bounded problems, the algebraic models are found to best predict the rise in turbulent Prandtl number near the wall as well as the log-layer temperature profile, while the thermal variance models provide a good representation of the RMS temperature fluctuations. In jet flows, the algebraic models provide no benefit over a constant turbulent Prandtl number approach. Application of the thermal variance models finds that some significantly overpredict the temperature variance in the plume and most underpredict the thermal growth rate of the jet. The models yield very similar fluctuating temperature intensities in jets from straight pipes and smooth contraction nozzles, in contrast to data that indicate the latter should have noticeably higher values. For the particular low subsonic heated jet cases examined, changes in the turbulent Prandtl number had no effect on the centerline velocity decay.

computational fluid dynamics↗

NASA NTRS2019

On the Use of Optimization Techniques for Turbulence Model Calibration

This article discusses the use of numerical optimization procedures to aid in the calibration of turbulence model coefficients. Such methods would increase the rigor and repeatability of the calibration procedure by requiring clearly defined and objective optimization metrics, and could be used to identify unique combinations of coefficient values for specific flow problems. The approach is applied to the re-calibration of an explicit algebraic Reynolds stress model for the incompressible planar mixing layer using the Nelder-Mead simplex algorithm and a micro-genetic algorithm with minimally imposed constraints. Three composite fitness functions, each based upon the error in the mixing layer growth rate and the normal and shear components of the Reynolds stresses, are investigated. The results demonstrate a significant improvement in the target objectives through the adjustment of three pressure-strain coefficients. Adjustments of additional coefficients provide little further benefit. Issues regarding the effectiveness of the fitness functions and the efficiency of the optimization algorithms are also discussed.

Yoder, Dennis A.↗

DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗

NASA NTRS1973

A new approach to the design of time-varying control systems with application to the Space Shuttle Boost

A general approach to the analysis and design of closed-loop control for slowly time-varying linear systems is presented. The approach uses the method of generalized multiple scales in which slow and fast dynamics are systematically separated by employing different 'clocks' which measure time at varying rates. The clocks, which are necessarily nonlinear functions of time, are chosen in such a manner that the system dynamics are asymptotically invariant with respect to the new time scale. A transfer function relating the output to the input for general linear slowly time-varying systems is developed and represents the actual system under certain conditions. The clock function is shown to satisfy an algebraic characteristic equation and can be determined in terms of the coefficients. The approach can be viewed as an extension of time-invariant linear feedback control theory.

Ramanth, R. V.↗

DOE OSTI2024

Electromagnetic and two-photon transition form factors of the pseudoscalar mesons: An algebraic model computation

We compute electromagnetic and two-photon transition form factors of ground-state pseudoscalar mesons: π , K , η c , η b . To this end, we employ an algebraic model based upon the coupled formalism of Schwinger-Dyson and Bethe-Salpeter equations. Within this approach, the dressed quark propagator and the relevant Bethe-Salpeter amplitude encode the internal structure of the corresponding meson. Electromagnetic properties of the meson are probed via the quark-photon interaction. The algebraic model employed by us unifies the treatment of all ground-state pseudoscalar mesons. Its parameters are carefully fitted performing a global analysis of existing experimental data including the knowledge of the charge radii of the mesons studied. We then compute and predict electromagnetic and two-photon transition form factors for a wide range of probing photon momentum-squared which is of direct relevance to the experimental observations carried out thus far or planned at different hadron physics facilities such as the Thomas Jefferson National Accelerator Facility (JLab) and the forthcoming Electron-Ion Collider. We also present comparisons with other theoretical models and approaches and lattice quantum chromodynamics. Published by the American Physical Society 2024

Higuera-Angulo, I. M. (ORCID:0000000256008875)↗

DOE OSTI2023

Tomographic optical emission spectroscopy of an atmospheric pressure plasma jet and surface ionization waves on planar and structured surfaces

Here in this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N$_2^+$, and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N$_2^+$ and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N$_2^+$ that spatially converged to form the 'plasma bullet' or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.

plasma interactions with complex surfaces↗

DOE OSTI2022

A fast and accurate modeling approach for water and steam thermodynamics with practical applications in district heating system simulation

In U.S. district heating (DH) systems, steam is the most common heat transport medium. Industry demand for new advanced modeling capabilities of complete steam DH systems is increasing; however, the existing models for water/steam thermodynamics are too slow for large system simulations because of computationally expensive algebraic loops that require the solution to nonlinear systems of equations. For practical applications, this work presents a novel split-medium approach that implements numerically efficient liquid water models alongside various water/steam models, breaking costly algebraic loops by decoupling mass and energy balance equations. New component models for steam DH systems are also presented. We implemented the models in the equation based Modelica language and evaluated accuracy and computing speed across multiple scales: from fundamental thermodynamic properties to complete districts featuring 10 to 200 buildings. Compared to district models with the IF97 water/steam model and equipment models from the Modelica Standard Library, the new implementation improves the scaling rate for large districts from cubic to quadratic with negligible compromise to accuracy. Additionally, for an annual simulation with 180 buildings, this translates to a computing time reduction from 33 to 1-1.5 h. These results are critically important for industry practitioners to simulate steam DH systems at large scales.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1996

Generalized Bezout's Theorem and its applications in coding theory

This paper presents a generalized Bezout theorem which can be used to determine a tighter lower bound of the number of distinct points of intersection of two or more curves for a large class of plane curves. A new approach to determine a lower bound on the minimum distance (and also the generalized Hamming weights) for algebraic-geometric codes defined from a class of plane curves is introduced, based on the generalized Bezout theorem. Examples of more efficient linear codes are constructed using the generalized Bezout theorem and the new approach. For d = 4, the linear codes constructed by the new construction are better than or equal to the known linear codes. For d greater than 5, these new codes are better than the known codes. The Klein code over GF(2(sup 3)) is also constructed.

Berg, Gene A.↗

NASA NTRS1992

An analytical approach to grid sensitivity analysis

Sensitivity analysis in Computational Fluid Dynamics with emphasis on grids and surface parameterization is described. An interactive algebraic grid-generation technique is employed to generate C-type grids around NACA four-digit wing sections. An analytical procedure is developed for calculating grid sensitivity with respect to design parameters of a wing section. A comparison of the sensitivity with that obtained using a finite-difference approach is made. Grid sensitivity with respect to grid parameters, such as grid-stretching coefficients, are also investigated. Using the resultant grid sensitivity, aerodynamic sensitivity is obtained using the compressible two-dimensional thin-layer Navier-Stokes equations.

Sadrehaghighi, Ideen↗

DOE OSTI2021

Physics of the Inverted Harmonic Oscillator: From the lowest Landau level to event horizons

In this work, we present the inverted harmonic oscillator (IHO) Hamiltonian as a paradigm to understand the quantum mechanics of scattering and time-decay in a diverse set of physical systems. As one of the generators of area preserving transformations, the IHO Hamiltonian can be studied as a dilatation generator, squeeze generator, a Lorentz boost generator, or a scattering potential. In establishing these different forms, we demonstrate the physics of the IHO that underlies phenomena as disparate as the Hawking–Unruh effect and scattering in the lowest Landau level (LLL) in quantum Hall systems. We derive the emergence of the IHO Hamiltonian in the LLL in a gauge invariant way and show its exact parallels with the Rindler Hamiltonian that describes quantum mechanics near event horizons. This approach of studying distinct physical systems with symmetries described by isomorphic Lie algebras through the emergent IHO Hamiltonian enables us to reinterpret geometric response in the lowest Landau level in terms of relativistic effects such as Wigner rotation. Further, the analytic scattering matrix of the IHO points to the existence of quasinormal modes (QNMs) in the spectrum, which have quantized time-decay rates. We present a way to access these QNMs through wave packet scattering, thus proposing a novel effect in quantum Hall point contact geometries that parallels those found in black holes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗