Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Jacobi Method”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Application of convolutional neural networks for stellar spectral classification

ABSTRACT Due to the ever-expanding volume of observed spectroscopic data from surveys such as SDSS and LAMOST, it has become important to apply artificial intelligence (AI) techniques for analysing stellar spectra to solve spectral classification and regression problems like the determination of stellar atmospheric parameters Teff, $\rm {\log g}$, and [Fe/H]. We propose an automated approach for the classification of stellar spectra in the optical region using convolutional neural networks (CNNs). Traditional machine learning (ML) methods with ‘shallow’ architecture (usually up to two hidden layers) have been trained for these purposes in the past. However, deep learning methods with a larger number of hidden layers allow the use of finer details in the spectrum which results in improved accuracy and better generalization. Studying finer spectral signatures also enables us to determine accurate differential stellar parameters and find rare objects. We examine various machine and deep learning algorithms like artificial neural networks, Random Forest, and CNN to classify stellar spectra using the Jacoby Atlas, ELODIE, and MILES spectral libraries as training samples. We test the performance of the trained networks on the Indo-U.S. Library of Coudé Feed Stellar Spectra (CFLIB). We show that using CNNs, we are able to lower the error up to 1.23 spectral subclasses as compared to that of two subclasses achieved in the past studies with ML approach. We further apply the trained model to classify stellar spectra retrieved from the SDSS data base with SNR > 20.

Sharma, Kaushal↗
NASA NTRS1988

The Jacobi Matrix Technique in Computational Fluid Dynamics

The Jacobi matrix technique was applied to the direct calculation of inviscid supersonic flow about two dimensional airfoils of varying thickness, angle of attack and camber and axisymmetric bodies of varying thickness and taper; and the design (inverse) calculation of inviscid supersonic flow past airfoils described by a given family of pressure distributions and axisymmetric bodies described by a given family of pressure distributions. The method was also applied to subsonic potential flow about two dimensional airfoils by modifying Jameson's FLO36. Results are discussed.

Sharp, H. Thomas↗
NASA NTRS2004

Fracture Mechanics Analysis of an Annular Crack in a Three-concentric-cylinder Composite Model

A boundary-value problem governing a three-phase concentric-cylinder model was analytically modeled to analyze annular interfacial crack problems with Love s strain functions in order to find the stress intensity factors (SIFs) and strain energy release rates (SERRs) at the tips of an interface crack in a nonhomogeneous medium. The complex form of a singular integral equation (SIE) of the second kind was formulated using Bessel s functions in the Fourier domain, and the SIF and total SERR were calculated using Jacoby polynomials. For the validity of the SIF equations to be established, the SIE of the three-concentric-cylinder model was reduced to the SIE for a two-concentric-cylinder model, and the results were compared with the previous results of Erdogan. A preliminary set of parametric studies was carried out to show the effect of interphase properties on the SERR. The method presented here provides insight about the effect of interphase properties on the crack driving force.

Kuguoglu, Latife H.↗
DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Two-dimensional computer simulation of EMVJ and grating solar cells under AMO illumination

A computer program, SCAP2D (Solar Cell Analysis Program in 2-Dimensions), is used to evaluate the Etched Multiple Vertical Junction (EMVJ) and grating solar cells. The aim is to demonstrate how SCAP2D can be used to evaluate cell designs. The cell designs studied are by no means optimal designs. The SCAP2D program solves the three coupled, nonlinear partial differential equations, Poisson's Equation and the hole and electron continuity equations, simultaneously in two-dimensions using finite differences to discretize the equations and Newton's Method to linearize them. The variables solved for are the electrostatic potential and the hole and electron concentrations. Each linear system of equations is solved directly by Gaussian Elimination. Convergence of the Newton Iteration is assumed when the largest correction to the electrostatic potential or hole or electron quasi-potential is less than some predetermined error. A typical problem involves 2000 nodes with a Jacobi matrix of order 6000 and a bandwidth of 243.

Gray, J. L.↗
DOE OSTI2023

Enhanced internal condensation of R1233zd(E) on micro- and nanostructured copper and aluminum surfaces

In-tube condensation of refrigerants is an important process which affects thermal efficiency in many applications, ranging from refrigeration and air conditioning to electronics thermal management. In-tube heat transfer and pressure drop are important to heat exchanger sizing and design. Here, in this work, micro- and nanostructured surfaces are applied to the internal wetted areas of copper and aluminum mini-channels to enhance the condensation heat transfer coefficient of hydrofluorocarbon R1233zd(E) refrigerant. To achieve scalable nanomanufacturing, surfaces were uniformly structured by relying on hydrochloric acid etching of aluminum and chemical oxidation of copper. The etched aluminum surfaces exhibited a 150% increase in heat transfer coefficient compared to smooth aluminum channels at specific qualities, with a 66% heat transfer coefficient improvement for complete phase change from saturated vapor to saturated liquid. Copper oxide structures showed no discernable difference in thermal-hydraulic performance when compared to smooth copper channels. Critical dimensionless parameters governing the heat transfer enhancement were identified by varying the tube internal diameter (2.3 mm to 4.7 mm), refrigerant mass flux (50 to 300 kg/(m 2 ·s)), and refrigerant quality (0 to 1). The dimensionless parameters include the Bond number normalized to the condensate film thickness, and the Weber number modified by the vapor friction factor. The relatively small increase in pressure drop (< 10%) associated with these surface enhancements further supports the promise of this method. The scalable and cost-effective techniques used to create these aluminum microstructures may reduce manufacturing cost when compared with current enhancement approaches such as extrusion, drawing, and welding.

42 ENGINEERING↗
DOE OSTI2012

Linearization Techniques for Controlled Piecewise Deterministic Markov Processes; Application to Zubov's Method

We aim at characterizing domains of attraction for controlled piecewise deterministic processes using an occupational measure formulation and Zubov's approach. Firstly, we provide linear programming (primal and dual) formulations of discounted, infinite horizon control problems for PDMPs. These formulations involve an infinite-dimensional set of probability measures and are obtained using viscosity solutions theory. Secondly, these tools allow to construct stabilizing measures and to avoid the assumption of stability under concatenation for controls. The domain of controllability is then characterized as some level set of a convenient solution of the associated Hamilton-Jacobi integral-differential equation. The theoretical results are applied to PDMPs associated to stochastic gene networks. Explicit computations are given for Cook's model for gene expression.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Polarizabilities of two-electron positive ions and Rydberg levels of lithium

Second-order sums are computed for the two-electron positive ions Li(+), Be(2+), and B(3+) by representing the intermediate states as discrete pseudostates. The Hamiltonian for the system is given and translated to semi-Jacobi variables for its derivation in reduced Rydberg units. The interaction potential is developed in its multipole form and employed in second-order perturbation theory. The perturbation is always a multipole-potential term, and the initial and intermediate states are constructed from two-particle Hylleraas basis sets. The resulting pseudostates are shown to represent the continuum with good convergence, and the asymptotic-optical-potential method is applied to the two-body quantities to determine the energy levels of three-electron systems in high Rydberg states. The method is shown to give accurate values for the polarizabilities of the two-electron isoelectronic systems.

Bhatia, A. K.↗
DOE OSTI2022

A high accuracy/resolution spectral element/Fourier–Galerkin method for the simulation of shoaling non-linear internal waves and turbulence in long domains with variable bathymetry

A high-order hybrid continuous-Galerkin numerical method, designed for the simulation of non-linear, non -hydrostatic internal waves and turbulence in long computational domains with complex bathymetry, is presented. The spatial discretization in the non-periodic wave-propagating directions, utilizes the nodal spectral element method. Such a high-order element-based discretization allows the highly accurate representation of complex domain geometry along with the flexibility of concentrating resolution in areas of interest. Under the assumption of the normal-to-isobath propagation of non-linear internal waves, a third periodic direction is incorporated via a Fourier-Galerkin discretization. The distinct non-hydrostatic nature of non-linear internal waves and, any instabilities and turbulence therein, necessitates the numerically challenging solution of the pressure Poisson problem. A defining feature of this work is the application of a domain decomposition approach, combined with block-Jacobi/deflation-based preconditioning to the pressure Poisson problem. Such a combined approach is particularly suitable for the long high aspect-ratio complex domains of interest and enables the efficient high-accuracy reproduction of the non-hydrostatic dynamics of non-linear internal waves. Implementation details are also described in the context of the stability of the solver and its parallelization strategy. A series of benchmarks of increasing complexity demonstrate the robustness of the flow solver. The benchmarks culminate with the three-dimensional simulation of a convectively breaking mode-one non-linear internal wave over a realistic South-China-Sea bathymetric transect and background current/stratification profiles.

Deflation↗
DOE OSTI2015

A Semi-linear Backward Parabolic Cauchy Problem with Unbounded Coefficients of Hamilton–Jacobi–Bellman Type and Applications to Optimal Control

We obtain weighted uniform estimates for the gradient of the solutions to a class of linear parabolic Cauchy problems with unbounded coefficients. Such estimates are then used to prove existence and uniqueness of the mild solution to a semi-linear backward parabolic Cauchy problem, where the differential equation is the Hamilton–Jacobi–Bellman equation of a suitable optimal control problem. Via backward stochastic differential equations, we show that the mild solution is indeed the value function of the controlled equation and that the feedback law is verified.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Ultra-compact and efficient heat exchanger for solar desalination with unprecedented scaling resistance

The objective is to design, prototype, and test a heat exchanger technology with an efficiency improvement of more than 150% compared to state-of-the-art heat exchanger technology (UA/volume), with significantly improved fouling/scaling resistance at elevated temperatures (>120°C) when compared to existing technology where scaling is a major performance-limiting issue, as well as at lower operating-temperatures (~63-110°C). We now report our overall progress on the project. An ultra-compact, efficient heat exchanger geometry has been developed that can maximize surface-area-to-volume ratio so that proposed target (UA/volume). Advanced manufacturing techniques, such additive manufacturing, have made it possible to create complex structures cost effectively. In addition to traditional manufacturing methods, additive manufacturing has been successfully used for the development of high-performance geometries. These methods allow the ability to create custom surface textures, which are deployed to dramatically control wettability and create energy-gradient surfaces. These techniques allow fouling/scaling nucleation to be spatially controlled, such that as the fouling layer grows it will be broken and removed by flow forces.

42 ENGINEERING↗
DOE OSTI2023

Control of Fractional Diffusion Problems via Dynamic Programming Equations

In this study, we explore the approximation of feedback control of integro-differential equations containing a fractional Laplacian term. To obtain feedback control for the state variable of this nonlocal equation, we use the Hamilton–Jacobi–Bellman equation. It is well known that this approach suffers from the curse of dimensionality, and to mitigate this problem we couple semi-Lagrangian schemes for the discretization of the dynamic programming principle with the use of Shepard approximation. This coupling enables approximation of high-dimensional problems. Numerical convergence toward the solution of the continuous problem is provided together with linear and nonlinear examples. The robustness of the method with respect to disturbances of the system is illustrated by comparisons with an open-loop control approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

The theory of stochastic cosmological lensing

On the scale of the light beams subtended by small sources, e.g. supernovae, matter cannot be accurately described as a fluid, which questions the applicability of standard cosmic lensing to those cases. In this article, we propose a new formalism to deal with small-scale lensing as a diffusion process: the Sachs and Jacobi equations governing the propagation of narrow light beams are treated as Langevin equations. We derive the associated Fokker-Planck-Kolmogorov equations, and use them to deduce general analytical results on the mean and dispersion of the angular distance. This formalism is applied to random Einstein-Straus Swiss-cheese models, allowing us to: (1) show an explicit example of the involved calculations; (2) check the validity of the method against both ray-tracing simulations and direct numerical integration of the Langevin equation. As a byproduct, we obtain a post-Kantowski-Dyer-Roeder approximation, accounting for the effect of tidal distortions on the angular distance, in excellent agreement with numerical results. Besides, the dispersion of the angular distance is correctly reproduced in some regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Measuring the x-ray resolving power of bent potassium acid phthalate diffraction crystals

This report presents the results from measuring the X-ray resolving power of a curved potassium acid phthalate (KAP(001)) spectrometer crystal using two independent methods. It is part of a continuing effort to measure the fundamental diffraction properties of bent crystals that are used to study various characteristics of high temperature plasmas. Bent crystals like KAP(001) do not usually have the same diffraction properties as corresponding flat crystals. Models that do exist to calculate the effect of bending the crystal on the diffraction properties have simplifying assumptions and their accuracy limits have not been adequately determined. The type of crystals that we measured is being used in a spectrometer on the Z machine at Sandia National Laboratories in Albuquerque, New Mexico. The first technique for measuring the crystal resolving power measures the X-ray spectral line width of the characteristic lines from several metal anodes. The second method uses a diode X-ray source and a double crystal diffractometer arrangement to measure the reflectivity curve of the KAP(001) crystal. The width of that curve is inversely proportional to the crystal resolving power. The measurement results are analyzed and discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1984

Secondary pattern computation of an arbitrarily shaped main reflector

The secondary pattern of a perfectly conducting offset main reflector being illuminated by a point feed at an arbitrary location was studied. The method of analysis is based upon the application of the Fast Fourier Transform (FFT) to the aperture fields obtained using geometrical optics (GO) and geometrical theory of diffraction (GTD). Key features of the reflector surface is completely arbitrary, the incident field from the feed is most general with arbitrary polarization and location, and the edge diffraction is calculated by either UAT or by UTD. Comparison of this technique for an offset parabolic reflector with the Jacobi-Bessel and Fourier-Bessel techniques shows good agreement. Near field, far field, and scan data of a large reflector are presented.

Lam, P. T. C.↗
NASA NTRS2002

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Bryson, Steve↗
DOE OSTI2011

Extremal trajectories in a nilpotent sub-Riemannian problem on the Engel group

On the Engel group a nilpotent sub-Riemannian problem is considered, a 4-dimensional optimal control problem with a 2-dimensional linear control and an integral cost functional. It arises as a nilpotent approximation to nonholonomic systems with 2-dimensional control in a 4-dimensional space (for example, a system describing the navigation of a mobile robot with trailer). A parametrization of extremal trajectories by Jacobi functions is obtained. A discrete symmetry group and its fixed points, which are Maxwell points, are described. An estimate for the cut time (the time of the loss of optimality) on extremal trajectories is derived on this basis. Bibliography: 25 titles.

97 MATHEMATICS AND COMPUTING↗