Search NASA⌕ Search

SEARCH · Search NASA

Results for “discretization methods”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A tensor train-based isogeometric solver for large-scale 3D poisson problems

We introduce a three-dimensional (3D), fully tensor train (TT) assembled isogeometric analysis (IGA) framework, TT-IGA, for solving partial differential equations (PDEs). Our method reformulates IGA discrete operators into TT format, enabling efficient compression and computation. Geometry evaluations use the original NURBS description at sampling points and TT approximation is applied to geometry-derived coefficient fields and discrete operators. We demonstrate the effectiveness of the proposed TT-IGA framework on the three-dimensional Poisson equation, achieving substantial reductions in memory and computational cost without compromising solution quality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Experimental measurements and mathematical modeling of cold plate for aviation thermal management

Herein, this study, which has been motivated by the recent applications of the cold plate device in aviation thermal management, reports on physics-based mathematical models derived from the conservation laws of mass, momentum, and energy, and empiricism-based models. One of the objectives of the present work is to report on an elaborate and successful experimental work carried out on an additively manufactured device for the purpose of rigorously validating the numerical predictions. The excellent agreement between the numerical predictions and measured performance provides much needed confidence in the implementation, in the software package, of the offset-strip fin passage correlations, as well as in the software implementation of user-defined wavy fin correlations for aerospace heat exchangers and cold plates operating with ram air at Reynolds numbers below 8000. The contributions of this work can also be found in the development of a new and accurate thermal-hydraulic analysis procedure, referred to in this paper as plate-fin analogy. Results from this procedure are compared with those from thermal resistance network. The comparative study in this paper of the bulk and discrete enthalpy flux method is also new, as is the relative assessment of four off-set strip fin thermal-hydraulic models.

42 ENGINEERING↗
DOE OSTI2019

A FFT-based numerical implementation of mesoscale field dislocation mechanics: Application to two-phase laminates

In this study, we present an enhanced crystal plasticity elasto-viscoplastic fast Fourier transform (EVPFFT) formulation coupled with a phenomenological Mesoscale Field Dislocation Mechanics (MFDM) theory here named MFDM-EVPFFT formulation. In contrast with classic CP-EVPFFT, the model is able to tackle plastic flow and hardening due to polar dislocation density distributions or geometrically necessary dislocations (GNDs) in addition to statistically stored dislocations (SSDs). The model also considers GND mobility through a GND density evolution law numerically solved with a recently developed filtered spectral approach, which is here coupled with stress equilibrium. The discrete Fourier transform method combined with finite differences is applied to solve both lattice incompatibility and Lippmann-Schwinger equations in an augmented Lagrangian numerical scheme. Numerical results are presented for two-phase laminate composites with plastic channels and elastic second phase. It is shown that both GND densities and slip constraint at phase boundaries influence the overall and local hardening behavior. In contrast with the CP-EVPFFT formulation, a channel size effect is predicted on the shear flow stress with the present MFDM-EVPFFT formulation. The size effect originates from the progressive formation of continuous screw GND pile-ups from phase boundaries to the channel center. The effect of GND mean free path on local and global responses is also examined for the two-phase composite.

36 MATERIALS SCIENCE↗
DOE OSTI2015

An algebraic method for constructing stable and consistent autoregressive filters

In this paper, we introduce an algebraic method to construct stable and consistent univariate autoregressive (AR) models of low order for filtering and predicting nonlinear turbulent signals with memory depth. By stable, we refer to the classical stability condition for the AR model. By consistent, we refer to the classical consistency constraints of Adams–Bashforth methods of order-two. One attractive feature of this algebraic method is that the model parameters can be obtained without directly knowing any training data set as opposed to many standard, regression-based parameterization methods. It takes only long-time average statistics as inputs. The proposed method provides a discretization time step interval which guarantees the existence of stable and consistent AR model and simultaneously produces the parameters for the AR models. In our numerical examples with two chaotic time series with different characteristics of decaying time scales, we find that the proposed AR models produce significantly more accurate short-term predictive skill and comparable filtering skill relative to the linear regression-based AR models. These encouraging results are robust across wide ranges of discretization times, observation times, and observation noise variances. Finally, we also find that the proposed model produces an improved short-time prediction relative to the linear regression-based AR-models in forecasting a data set that characterizes the variability of the Madden–Julian Oscillation, a dominant tropical atmospheric wave pattern.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

POD/DEIM reduced-order strategies for efficient four dimensional variational data assimilation

This work studies reduced order modeling (ROM) approaches to speed up the solution of variational data assimilation problems with large scale nonlinear dynamical models. It is shown that a key requirement for a successful reduced order solution is that reduced order Karush–Kuhn–Tucker conditions accurately represent their full order counterparts. In particular, accurate reduced order approximations are needed for the forward and adjoint dynamical models, as well as for the reduced gradient. New strategies to construct reduced order based are developed for proper orthogonal decomposition (POD) ROM data assimilation using both Galerkin and Petrov–Galerkin projections. For the first time POD, tensorial POD, and discrete empirical interpolation method (DEIM) are employed to develop reduced data assimilation systems for a geophysical flow model, namely, the two dimensional shallow water equations. Numerical experiments confirm the theoretical framework for Galerkin projection. In the case of Petrov–Galerkin projection, stabilization strategies must be considered for the reduced order models. The new reduced order shallow water data assimilation system provides analyses similar to those produced by the full resolution data assimilation system in one tenth of the computational time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

AP-Cloud: Adaptive Particle-in-Cloud method for optimal solutions to Vlasov–Poisson equation

We propose a new adaptive Particle-in-Cloud (AP-Cloud) method for obtaining optimal numerical solutions to the Vlasov–Poisson equation. Unlike the traditional particle-in-cell (PIC) method, which is commonly used for solving this problem, the AP-Cloud adaptively selects computational nodes or particles to deliver higher accuracy and efficiency when the particle distribution is highly non-uniform. Unlike other adaptive techniques for PIC, our method balances the errors in PDE discretization and Monte Carlo integration, and discretizes the differential operators using a generalized finite difference (GFD) method based on a weighted least square formulation. As a result, AP-Cloud is independent of the geometric shapes of computational domains and is free of artificial parameters. Efficient and robust implementation is achieved through an octree data structure with 2:1 balance. We analyze the accuracy and convergence order of AP-Cloud theoretically, and verify the method using an electrostatic problem of a particle beam with halo. Simulation results show that the AP-Cloud method is substantially more accurate and faster than the traditional PIC, and it is free of artificial forces that are typical for some adaptive PIC techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Treating network junctions in finite volume solution of transient gas flow models

A finite volume scheme for the numerical solution of a non-isothermal non-adiabatic compressible flow model for gas transportation networks on non-flat topography is introduced. Unlike standard Euler equations, the model takes into account wall friction, variable height and heat transfer between the pipe and the environment which are source terms. The case of one single pipe was considered in a previous reference by the authors, , where a finite volume method with upwind discretization of the flux and source terms has been proposed in order to get a well-balanced scheme. The main goal of the present paper is to go a step further by considering a network of pipes. The main issue is the treatment of junctions for which container-like 2D finite volumes are introduced. The couplings between pipes (1D) and containers (2D) are carefully described and the conservation properties are analyzed. Numerical tests including real gas networks are solved showing the performance of the proposed methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

4.0 MOOSE: Enabling massively parallel Multiphysics simulation

Approaching 18 years of existence, MOOSE—the Multiphysics Object-Oriented Simulation Environment—is being developed at a higher pace than ever before. With significant support from four research institutions across the globe, and dozens of new contributors, the capabilities of the framework are being expanded to meet modeling challenges in a wide variety of fields from nuclear system design, to geomechanics, to material science. This includes new development in equation discretization techniques, solver methods, meshing capabilities, application deployment, and user interface improvements. Applications built on MOOSE benefit from all these improvements.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Measurements of neutron capture cross section for {sup 207,208}Pb

The neutron capture cross sections for {sup 207,208}Pb have been measured in the neutron energy region from 10 to 110 keV. The γ-rays cascaded from a capture state to the ground state or low-lying states of {sup 208,209}Pb were observed for the first time, using an anti-Compton Nal(Tl) spectrometer and a TOF method. The observed discrete γ-ray energy spectra enabled us to determine neutron capture cross sections for {sup 207,208}Pb with small systematic errors, since we could distinguish γ-ray of {sup 207,208}Pb(n,γ) reactions from background γ-ray with use of the γ-ray spectra. The obtained cross sections include both contributions of resonance and direct capture components different from the previous TOF measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Theoretical analysis of selectivity mechanisms in molecular transport through channels and nanopores

Selectivity is one of the most fundamental concepts in natural sciences, and it is also critically important in various technological, industrial, and medical applications. Although there are many experimental methods that allow to separate molecules, frequently they are expensive and not efficient. Recently, a new method of separation of chemical mixtures based on utilization of channels and nanopores has been proposed and successfully tested in several systems. However, mechanisms of selectivity in the molecular transport during the translocation are still not well understood. Here, we develop a simple theoretical approach to explain the origin of selectivity in molecular fluxes through channels. Our method utilizes discrete-state stochastic models that take into account all relevant chemical transitions and can be solved analytically. More specifically, we analyze channels with one and two binding sites employed for separating mixtures of two types of molecules. The effects of the symmetry and the strength of the molecular-pore interactions are examined. It is found that for one-site binding channels, the differences in the strength of interactions for two species drive the separation. At the same time, in more realistic two-site systems, the symmetry of interaction potential becomes also important. The most efficient separation is predicted when the specific binding site is located near the entrance to the nanopore. In addition, the selectivity is higher for large entrance rates into the channel. It is also found that the molecular transport is more selective for repulsive interactions than for attractive interactions. The physical-chemical origin of the observed phenomena is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

On control of continuous dynamical polysystems in discrete times

This paper considers control systems with a finite number of control parameters, i.e. dynamical polysystems. It is assumed that control switchings can occur only at certain discrete times. Statistical methods in number theory are used to investigate the controllability of these systems. Existence of control switching times is established through analysis of Diophantine equations. The values of control switching times are found through suitable numerical methods. Certain assertions on the controllability of the polysystems of the class under consideration are proved. Some examples are considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Cavity topologies of mode-locked fibre lasers: possibilities and prospects

We report, for the first time, the results of a comparative analysis of cavity topologies in mode-locked fibre lasers. It is shown that the use of various topologies and advantages of fibre technologies allows one to realise the specified combinations of dispersion, nonlinear, polarisation and other properties of fibre laser cavities, difficult to implement or not implementable using the methods of traditional discrete laser optics. The results of the work make it possible to consider the topological engineering of mode-locked fibre lasers as a new, rapidly growing field of research and development. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Hints of correlation between broad-line and radio variations for 3C 120

In this paper, we investigate the correlation between broad-line and radio variations for the broad-line radio galaxy 3C 120. By the z-transformed discrete correlation function method and the model-independent flux randomization/random subset selection (FR/RSS) Monte Carlo method, we find that broad Hβ line variations lead the 15 GHz variations. The FR/RSS method shows that the Hβ line variations lead the radio variations by a factor of τ{sub ob} = 0.34 ± 0.01 yr. This time lag can be used to locate the position of the emitting region of radio outbursts in the jet, on the order of ∼5 lt-yr from the central engine. This distance is much larger than the size of the broad-line region. The large separation of the radio outburst emitting region from the broad-line region will observably influence the gamma-ray emission in 3C 120.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

LONG-TERM PERIODICITY VARIATIONS OF THE SOLAR RADIUS

In order to study the long-term periodicity variations of the solar radius, daily solar radius data from 1978 February to 2000 September at the Calern Observatory are used. Continuous observations of the solar radius are difficult due to the weather, seasonal effects, and instrument characteristics. Thus, to analyze these data, we first use the Dixon criterion to reject suspect values, then we measure the cross-correlation between the solar radius and sunspot numbers. The result indicates that the solar radius is in complete antiphase with the sunspot numbers and shows lead times of 74 months relative to the sunspot numbers. The Lomb-Scargle and date compensated discrete Fourier transform methods are also used to investigate the periodicity of the solar radius. Both methods yield similar significance periodicities around {approx}1 yr, {approx}2.6 yr, {approx}3.6 yr, and {approx}11 yr. Possible mechanisms for these periods are discussed. The possible physical cause of the {approx}11 yr period is the cyclic variation of the magnetic pressure of the concentrated flux tubes at the bottom of the solar convection zone.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Hybrid Entanglement Distribution between Remote Microwave Quantum Computers Empowered by Machine Learning

Superconducting microwave circuits with Josephson junctions, the major platform for quantum computing, can only reach the full capability when connected. This requires an efficient protocol to distribute microwave entanglement. While quantum computers typically use discrete-variable (DV) methods for information encoding, the entire continuous-variable (CV) degree of freedom in electromagnetic fields must be utilized to achieve the highest entanglement distribution rate. Here, we propose a hybrid protocol to resolve the incompatibility between DV microwave quantum computers and CV quantum communications. CV microwave entanglement is distributed using optical swapping of optical-microwave entanglement pairs. To interface with DV microwave quantum computers, we further design a hybrid circuit to simultaneously convert and distill high-quality DV entanglement from noisy CV entanglement. The hybrid circuit is trained with machine-learning algorithms, ensuring high entanglement fidelity and generation rate. In conclusion, our work not only provides a practical method to realize efficient quantum links for superconducting microwave quantum computers, but also opens avenues to bridge the gap between DV and CV quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A novel software and conceptual design of the hardware platform for intensity modulated radiation therapy

Purpose: To deliver high quality intensity modulated radiotherapy (IMRT) using a novel generalized sparse orthogonal collimators (SOCs), the authors introduce a novel direct aperture optimization (DAO) approach based on discrete rectangular representation. Methods: A total of seven patients—two glioblastoma multiforme, three head & neck (including one with three prescription doses), and two lung—were included. 20 noncoplanar beams were selected using a column generation and pricing optimization method. The SOC is a generalized conventional orthogonal collimators with N leaves in each collimator bank, where N = 1, 2, or 4. SOC degenerates to conventional jaws when N = 1. For SOC-based IMRT, rectangular aperture optimization (RAO) was performed to optimize the fluence maps using rectangular representation, producing fluence maps that can be directly converted into a set of deliverable rectangular apertures. In order to optimize the dose distribution and minimize the number of apertures used, the overall objective was formulated to incorporate an L2 penalty reflecting the difference between the prescription and the projected doses, and an L1 sparsity regularization term to encourage a low number of nonzero rectangular basis coefficients. The optimization problem was solved using the Chambolle–Pock algorithm, a first-order primal–dual algorithm. Performance of RAO was compared to conventional two-step IMRT optimization including fluence map optimization and direct stratification for multileaf collimator (MLC) segmentation (DMS) using the same number of segments. For the RAO plans, segment travel time for SOC delivery was evaluated for the N = 1, N = 2, and N = 4 SOC designs to characterize the improvement in delivery efficiency as a function of N. Results: Comparable PTV dose homogeneity and coverage were observed between the RAO and the DMS plans. The RAO plans were slightly superior to the DMS plans in sparing critical structures. On average, the maximum and mean critical organ doses were reduced by 1.94% and 1.44% of the prescription dose. The average number of delivery segments was 12.68 segments per beam for both the RAO and DMS plans. The N = 2 and N = 4 SOC designs were, on average, 1.56 and 1.80 times more efficient than the N = 1 SOC design to deliver. The mean aperture size produced by the RAO plans was 3.9 times larger than that of the DMS plans. Conclusions: The DAO and dose domain optimization approach enabled high quality IMRT plans using a low-complexity collimator setup. The dosimetric quality is comparable or slightly superior to conventional MLC-based IMRT plans using the same number of delivery segments. The SOC IMRT delivery efficiency can be significantly improved by increasing the leaf numbers, but the number is still significantly lower than the number of leaves in a typical MLC.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2018

Modeling underwater noise propagation from marine hydrokinetic power devices through a time-domain, velocity-pressure solution

Marine hydrokinetic (MHK) devices generate electricity from the motion of tidal and ocean currents, as well as ocean waves, to provide an additional source of renewable energy available to the United States. These devices are a source of anthropogenic noise in the marine ecosystem and must meet regulatory guidelines that mandate a maximum amount of noise that may be generated. In the absence of measured levels from in situ deployments, a model for predicting the propagation of sound from an array of MHK sources in a real environment is essential. A set of coupled, linearized velocity-pressure equations in the time-domain are derived and presented in this paper, which are an alternative solution to the Helmholtz and wave equation methods traditionally employed. Discretizing these equations on a three-dimensional (3D), finite-difference grid ultimately permits a finite number of complex sources and spatially varying sound speeds, bathymetry, and bed composition. The solution to this system of equations has been parallelized in an acoustic-wave propagation package developed at Sandia National Labs, called Paracousti. This work presents the broadband sound pressure levels from a single source in two-dimensional (2D) ideal and Pekeris wave-guides and in a 3D domain with a sloping boundary. Furthermore, the paper concludes with demonstration of Paracousti for an array of MHK sources in a simple wave-guide.

13 HYDRO ENERGY↗