SEARCH · Search NASA

Results for “Variational discretization”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Development of upwind schemes for the Euler equations

Described are many algorithmic and computational aspects of upwind schemes and their second-order accurate formulations based on Total-Variation-Diminishing (TVD) approaches. An operational unification of the underlying first-order scheme is first presented encompassing Godunov's, Roe's, Osher's, and Split-Flux methods. For higher order versions, the preprocessing and postprocessing approaches to constructing TVD discretizations are considered. TVD formulations can be used to construct relaxation methods for unfactored implicit upwind schemes, which in turn can be exploited to construct space-marching procedures for even the unsteady Euler equations. A major part of the report describes time- and space-marching procedures for solving the Euler equations in 2-D, 3-D, Cartesian, and curvilinear coordinates. Along with many illustrative examples, several results of efficient computations on 3-D supersonic flows with subsonic pockets are presented.

Chakravarthy, Sukumar R.↗

NASA NTRS2000

Hypersonic Boundary-Layer Trip Development for Hyper-X

Boundary layer trip devices for the Hper-X forebody have been experimentally examined in several wind tunnels. Five different trip configurations were compared in three hypersonic facilities, the LaRC 20-Inch Mach 6 Air Tunnel, the LaRC 31 -Inch Mach 10 Air Tunnel, and in the HYPULSE Reflected Shock Tunnel at GASL. Heat transfer distributions, utilizing the phosphor thermography and thin-film techniques, shock system details, and surface streamline patterns were measured on a 0.333-scale model of the Hyper-X forebody. Parametric variations include angles-of-attack of 0-deg, 2-deg, and 4-deg; Reynolds numbers based on model length of 1.2 to 15.4 million: and inlet cowl door simulated in both open and closed positions. Comparisons of transition due to discrete roughness elements have led to the selection of a trip configuration for the Hyper-X Mach 7 flight vehicle.

Berry, Scott A.↗

NASA NTRS1971

ELF propagation in the plasmasphere based on satellite observations of discrete and continuous forms

The propagation of electromagnetic waves in a nonhomogeneous anisotropic medium is examined from the point of view of geometrical optics. In particular, the propagation of ELF waves in the magnetosphere is described in terms of the electron and ion densities and the intensity and inclination of the earth's magnetic field. The analysis of the variations of wave normal angle along the ray path is extended to include the effects of ions. A comparison of the relative importance of each of the above parameters in controlling the orientation of the wave normals is made in the region of the magnetosphere where most of the ion whistlers have been detected.

Muzzio, J. L. R.↗

NASA NTRS1987

Variational approach to probabilistic finite elements

Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

NASA NTRS1991

Variational approach to probabilistic finite elements

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗

DOE OSTI2025

High-throughput micro-scale bandgap mapping for perovskite-inspired materials with complex composition space

Abstract To realize the full promise of high-throughput experimental workflows, the rate of sample synthesis must be matched by that of characterization. Of growing interest are contactless optical techniques that can rapidly measure material homogeneity and properties. Here, we present a hyperspectral imaging method to measure local optical bandgap distributions within samples, utilizing spatially-resolved reflectance spectra coupled with automated data analysis. We collect approximately one million optical bandgap data across the compositional space of Cs 3 (Bi x Sb 1-x ) 2 (Br y I 1-y ) 9 perovskite-inspired materials. Our results show non-monotonic bandgap variations (i.e., bandgap bowing) along six composition gradient sequences, in addition to identifying samples with multiple bandgaps in statistics. High-throughput transient absorption spectroscopy reveals that within these compositions, the depletion of the ground state carriers to excited states occurred at discrete energy levels with independent carrier dynamics, consistent with the bandgap observation and indicative of phase separation. This work demonstrates the potential for rapid optical measurements to assess material quality and homogeneity in a high-throughput experimental setting, supporting screening and recipe optimization of optoelectronic material candidates with desired carrier dynamics and optical properties.

Science & Technology - Other Topics↗

DOE OSTI2022

Investigations in the tape-to-tape contact resistance and contact composition in superconducting CORC ® wires

Conductor on Round Core (CORC®) wires and cables, constructed from multiple layers of helically wound REBa 2 Cu 3 O 7–δ tapes, are a promising cable technology for high field magnet applications. An important feature of high-temperature superconductor cables is the ability to share current between conductors, allowing current to bypass drops in I c and minimizing the risk of hot spot formation, which could lead to potential burnout in the superconductor. In contrast to stacked-tape cables, which have continuous contact between tapes, in CORC® the transfer points occur at discrete tape crossovers. The tape-to-tape contact resistance, R c , plays a critical role in the current sharing capabilities and current distribution in cables. For the work reported here, special CORC® wires were manufactured using different winding parameters to investigate variations in R c . Variations comprised inclusion of a lubricant, different lubricant conductivity, inclusion of pre-tinning, and heating briefly to melt the solder. Cables were first tested as straight lengths, followed by bending to a 10 cm diameter. In straight cables R c values ranged from 1 to over 1000 μΩ cm 2 , depending on cabling parameters, with the highest values being found for cables made by the present 'standard' process. Bending the cables to a 10 cm diameter decreased R c by a factor 2–5. Tinning with PbSn decreased R c by three orders of magnitude compared to standard CORC® wires, and heat treating wires with tinned conductor resulted in only a small further decrease in R c . Based on the measured R c at an electric field of 1 μV cm –1 the resulting current transfer length between layers can range from a few millimeters to a tens of centimeters. Examination of contacts with a laser confocal microscope showed plastic deformation of the copper at the edges of the contact overlap area, apparently caused by thicker plating at tape edges digging into the copper of neighboring layers. These images reveal that only a fraction of the total contact surface may actually be touching when there is nothing to compensate for height differential. Images of the PbSn coated tapes indicated that application of solder produces a much more uniform contact surface and higher contact area. Furthermore, imaging of CORC® cross-sections confirmed that in the non-tinned cables there are many regions where tapes are not in contact, while in contrast the PbSn cable shows significantly more contact between the tapes. These different imaging techniques reveal that tape surface morphology is a significant parameter in determining R c .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

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 NTRS1986

Comet Halley near-nucleus phenomena in 1986

Computer enhanced high resolution groundbased images taken during Halley's most active phases in 1985-86 exhibit all of the coma structures identified in 1910 photographs. The structures result mostly from discrete dust jet activity producing spirals, arcs, and jagged-edged envelopes. Multiple nucleus and antisunward jet phenomena are also observed. The curvature of the jets is consistent with the location of the Earth with respect to the nucleus equatorial plane. The variation in jet activity during the Vega and Giotto flybys is easily seen on the groundbased images. Images of CN jets show that they usually do not correlate with dust jets. Wide field, medium resolution images of faint outer envelopes are required to augment the smaller field, high resolution images of near-nucleus jets in modeling the jets.

Larson, S.↗

NASA NTRS1987

Irregular and diurnal variability in asynoptic measurements of stratospheric trace species

The consequences of irregular and diurnal variability on the interpretation of asynoptic measurements of stratospheric trace species are explored. In particular, the fidelity with which the continuous behavior may be recovered from discrete asynoptic measurements is examined. Dynamically created tracer variability in the form of an advected space-time random process is discussed. Diagnostics such as power spectra and time-mean fields are derived asynoptically and compared with the true behavior. Diurnal variations of a photochemically active species are considered, describing variability in terms of a solar waveform which propagates through a latitudinal envelope. Several diagnostics of the behavior are derived asynoptically and compared with the true variability.

Salby, Murry L.↗

NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

NASA NTRSDate not supplied

Wall-Modeled Large Eddy Simulation Method for Unstructured-Grid Navier-Stokes Solvers

This paper reports on the implementation and assessment of a Wall-Modeled Large-Eddy Simulation (WMLES) methodology in an unstructured-grid, node-centered flow solver, FUN3D that is developed and supported at the NASA Langley Research Center. Finite-volume (FV) and finite-element (FE) discretization schemes considered in the study provide formal second-order spatial accuracy. Large-Eddy Simulations (LES) resolve large-scale turbulent-flow features and filter out small-scale effects using the Vreman subgrid-scale model. At solid-wall boundaries, a shear-stress model is employed to provide a proper boundary-flux closure. The nonlinear equations are integrated in time using either an optimized backward difference formula or an implicit multistage Runge-Kutta temporal scheme. The implicit equations at each time step are solved by strong nonlinear iteration schemes. WMLES demonstrations are shown for two high-lift configurations, namely, the McDonnell Douglas 30P30N multielement airfoil and a NASA High-Lift Common Research Model. Results show that the WMLES approaches implemented in the FV and FE discretization methods produce consistent solutions and are capable of capturing key aerodynamic characteristics and flow structures for high-lift configurations at a wide range of angles of attack including maximum-lift conditions. In the 30P30N example, correct trends in the variations of integrated aerodynamic forces and moments, surface pressure distributions, and boundary-layer profiles are captured as the Reynolds number is increased.

CFD; turbulence modeling; High-Lift flow simulatio↗

NASA NTRS1991

Space Shuttle Redesigned Solid Rocket Motor nozzle natural frequency variations with burn time

The effects of erosion and thermal degradation on the Space Shuttle Redesigned Solid Rocket Motor (RSRM) nozzle's structural dynamic characteristics were analytically evaluated. Also considered was stiffening of the structure due to internal pressurization. A detailed NASTRAN finite element model of the nozzle was developed and used to evaluate the influence of these effects at several discrete times during motor burn. Methods were developed for treating erosion and thermal degradation, and a procedure was developed to account for internal pressure stiffening using differential stiffness matrix techniques. Results were verified using static firing test accelerometer data. Fast Fourier Transform and Maximum Entropy Method techniques were applied to the data to generate waterfall plots which track modal frequencies with burn time. Results indicate that the lower frequency nozzle 'vectoring' modes are only slightly affected by erosion, thermal effects and internal pressurization. The higher frequency shell modes of the nozzle are, however, significantly reduced.

Lui, C. Y.↗

DOE OSTI2025

Conditioned quantum-assisted deep generative surrogate for particle-calorimeter interactions

Particle collisions at accelerators like the Large Hadron Collider (LHC), recorded by experiments such as ATLAS and CMS, enable precise standard model measurements and searches for new phenomena. Simulating these collisions significantly influences experiment design and analysis but incurs immense computational costs, projected at millions of CPU-years annually during the high luminosity LHC (HL-LHC) phase. Currently, simulating a single event with Geant4 consumes around 1000 CPU seconds, with calorimeter simulations especially demanding. To address this, we propose a conditioned quantum-assisted generative model, integrating a conditioned variational autoencoder (VAE) and a conditioned restricted Boltzmann machine (RBM). Our RBM architecture is tailored for D-Wave’s Pegasus-structured advantage quantum annealer for sampling, leveraging the flux bias for conditioning. This approach combines classical RBMs as universal approximators for discrete distributions with quantum annealing’s speed and scalability. We also introduce an adaptive method for efficiently estimating effective inverse temperature, and validate our framework on Dataset 2 of CaloChallenge.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1979

Albedo and color contrasts on asteroid surfaces

Asteroids in general display only small or negligible variations in spectrum or albedo during a rotational cycle. Color variations with rotation are described in the literature but are usually comparable to the noise in the measurements. Twenty-four asteroids have been systematically monitored for such color changes. Only 3 Juno, 4 Vesta, 6 Hebe, 71 Niobe, 349 Dembowska, and 944 Hidalgo display color variations larger than 0.03 mag. In each of these cases the asteroid appears redder near maximum brightness. Of seven asteroids monitored polarimetrically, only 4 Vesta shows a convincing variation, attributed to an albedo change with rotation. The lightcurve can be explained by albedo differences alone; Vesta apparently has a nearly spheroidal shape. Nothwithstanding the above results, the degree of uniformity of most asteroid surfaces is remarkable. If asteroids exist with large discrete domains of ferrosilicate, metallic, and/or carbonaceous material together on their surfaces, they have not yet been identified.

Degewij, J.↗

NASA NTRS1983

A sampling procedure to guide the collection of narrow-band, high-resolution spatially and spectrally representative reflectance data

A multistage sampling procedure using image processing, geographical information systems, and analytical photogrammetry is presented which can be used to guide the collection of representative, high-resolution spectra and discrete reflectance targets for future satellite sensors. The procedure is general and can be adapted to characterize areas as small as minor watersheds and as large as multistate regions. Beginning with a user-determined study area, successive reductions in size and spectral variation are performed using image analysis techniques on data from the Multispectral Scanner, orbital and simulated Thematic Mapper, low altitude photography synchronized with the simulator, and associated digital data. An integrated image-based geographical information system supports processing requirements.

Brand, R. R.↗

NASA NTRS1993

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗

NASA NTRS1996

Multigrid Acceleration of Time-Accurate Navier-Stokes Calculations

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is fully implicit in time and is unconditionally stable (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability. This is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. A multigrid-multiblock, steady-state, three-dimensional Navier-Stokes solver, TLNS3D, was modified to iteratively invert the equations at each physical time step. The implementation of this procedure in TLNS3D is discussed. The implications of applying several popular turbulence models to unsteady flow are also considered. Numerical results are presented to show the application of the scheme to various two-dimensional turbulent flows. The results of a three-dimensional laminar flow calculation are also given.

Melson, N. Duane↗