Search NASA⌕ Search

SEARCH · Search NASA

Results for “smoothing 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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Three-dimensional self-adaptive grid method for complex flows

A self-adaptive grid procedure for efficient computation of three-dimensional complex flow fields is described. The method is based on variational principles to minimize the energy of a spring system analogy which redistributes the grid points. Grid control parameters are determined by specifying maximum and minimum grid spacing. Multidirectional adaptation is achieved by splitting the procedure into a sequence of successive applications of a unidirectional adaptation. One-sided, two-directional constraints for orthogonality and smoothness are used to enhance the efficiency of the method. Feasibility of the scheme is demonstrated by application to a multinozzle, afterbody, plume flow field. Application of the algorithm for initial grid generation is illustrated by constructing a three-dimensional grid about a bump-like geometry.

Djomehri, M. Jahed↗
DOE OSTI2021

On the numerical accuracy in finite-volume methods to accurately capture turbulence in compressible flows

The goal of the present article is to understand the impact of numerical schemes for the reconstruction of data at cell faces in finite-volume methods, and to assess their interaction with the quadrature rule used to compute the average over the cell volume. Here, third-, fifth- and seventh-order WENO-Z schemes are investigated. On a problem with a smooth solution, the theoretical order of convergence rate for each method is retrieved, and changing the order of the reconstruction at cell faces does not impact the results, whereas for a shock-driven problem all the methods collapse to first-order. Here, study of the decay of compressible homogeneous isotropic turbulence reveals that using a high-order quadrature rule to compute the average over a finite-volume cell does not improve the spectral accuracy and that all methods present a second-order convergence rate. However the choice of the numerical method to reconstruct data at cell faces is found to be critical to correctly capture turbulent spectra. In the context of simulations with finite-volume methods of practical flows encountered in engineering applications, it becomes apparent that an efficient strategy is to perform the average integration with a low-order quadrature rule on a fine mesh resolution, whereas high-order schemes should be used to reconstruct data at cell faces.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1997

Nonlinear Rescaling and Proximal-Like Methods in Convex Optimization

The nonlinear rescaling principle (NRP) consists of transforming the objective function and/or the constraints of a given constrained optimization problem into another problem which is equivalent to the original one in the sense that their optimal set of solutions coincides. A nonlinear transformation parameterized by a positive scalar parameter and based on a smooth scaling function is used to transform the constraints. The methods based on NRP consist of sequential unconstrained minimization of the classical Lagrangian for the equivalent problem, followed by an explicit formula updating the Lagrange multipliers. We first show that the NRP leads naturally to proximal methods with an entropy-like kernel, which is defined by the conjugate of the scaling function, and establish that the two methods are dually equivalent for convex constrained minimization problems. We then study the convergence properties of the nonlinear rescaling algorithm and the corresponding entropy-like proximal methods for convex constrained optimization problems. Special cases of the nonlinear resealing algorithm are presented. In particular a new class of exponential penalty-modified barrier functions methods is introduced.

Polyak, Roman↗
DOE OSTI2019

Development of a P 2 element with optimal L 2 convergence for biharmonic equation

It is well known that convergence rate of finite element approximation is suboptimal in the L 2 norm for solving biharmonic equations when P 2 or Q 2 element is used. The goal of this paper is to derive a weak Galerkin (WG) P 2 element with the L 2 optimal convergence rate by assuming the exact solution sufficiently smooth. In addition, our new WG finite element method can be applied to general mesh such as hybrid mesh, polygonal mesh or mesh with hanging node. The numerical experiments have been conducted on different meshes including hybrid meshes with mixed of pentagon and rectangle and mixed of hexagon and triangle.

Mu, Lin↗
DOE OSTI2016

Discovery of a massive supercluster system at z ~ 0.47

Aims. Superclusters are the largest relatively isolated systems in the cosmic web. Using the SDSS BOSS survey, we search for the largest superclusters in the redshift range 0.43 < z < 0.71. Methods. We generate a luminosity-density field smoothed over 8 h -1 Mpc to detect the large-scale over-density regions. Each individual over-density region is defined as single supercluster in the survey. Here, we define the superclusters so that they are comparable to the superclusters found in the SDSS main survey. Results. We find a system that we call the BOSS Great Wall (BGW), which consists of two walls with diameters 186 and 173 h -1 Mpc and two other major superclusters with diameters of 64 and 91 h -1 Mpc. As a whole, this system consists of 830 galaxies with the mean redshift 0.47. We estimate the total mass to be approximately 2 × 10 17 h -1 M ⊙ .The morphology of the superclusters in the BGW system is similar to the morphology of the superclusters in the Sloan Great Wall region. Conclusions. The BGW is one of the most extended and massive systems of superclusters found so far in the Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Robustness of waves with a high phase velocity

Norman Rostoker pioneered research of (1) plasma-driven accelerators and (2) beam-driven fusion reactors. The collective acceleration, coined by Veksler, advocates to drive above-ionization plasma waves by an electron beam to accelerate ions. The research on this, among others, by the Rostoker group incubated the idea that eventually led to the birth of the laser wakefield acceleration (LWFA), by which a large and robust accelerating collective fields may be generated in plasma in which plasma remains robust and undisrupted. Besides the emergence of LWFA, the Rostoker research spawned our lessons learned on the importance of adiabatic acceleration of ions in collective accelerators, including the recent rebirth in laser-driven ion acceleration efforts in a smooth adiabatic fashion by a variety of ingenious methods. Following Rostoker’s research in (2), the beam-driven Field Reversed Configuration (FRC) has accomplished breakthroughs in recent years. The beam-driven kinetic plasma instabilities have been found to drive the reactivity of deuteron-deuteron fusion beyond the thermonuclear yield in C-2U plasma that Rostoker started. This remarkable result in FRCs as well as the above mentioned LWFA may be understood with the aid of the newly introduced idea of the “robustness hypothesis of waves with a high phase velocity”. It posits that when the wave driven by a particle beam (or laser pulse) has a high phase velocity, its amplitude is high without disrupting the supporting bulk plasma. This hypothesis may guide us into more robust and efficient fusion reactors and more compact accelerators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Breakup of a liquid rivulet falling over an inclined plate: Identification of a critical Weber number

Here, we have numerically investigated the breakup of a rivulet falling over a smooth inclined plate using the volume of fluid method. The breakup the rivulet is a complex phenomenon that is dictated by many factors, such as solvent properties, contact angle, inertia, plate inclination, etc. An extensive simulation campaign was conducted wherein these factors were systematically investigated. Regimes for a stable rivulet and an unstable rivulet that leads to the breakup and formation of a droplet are examined in terms of a critical value of the Weber number (We cr ) that delineates these regimes. The effect of plate inclination on the breakup of the rivulet shows that the critical Weber number decreases with increased inclination angle (θ) owing to higher liquid velocity. However, the effect is negligible beyond θ > 60°. The impact of solvent properties is characterized using the Kapitza number (Ka). Variation of We cr with Ka shows two trends depending on the Ka value of the solvent. Solvents with lower Ka values, corresponding to high viscosities and/or low surface tensions, show smaller values of the critical Weber number and the variation is linear. While solvents with higher Ka values exhibit higher values of the We cr and the variation in We cr is steep. This behavior is more pronounced with increasing contact angle. Higher contact angles promote rivulet breakup so that inertia must be higher to the breakup. As a result, a phenomenological scaling for a critical Weber number with the Kapitza number and contact angle is presented that can offer insight into rivulet breakup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Single shot high bandwidth laser plasma probe

Herein a novel scheme to probe a plasma's ion acoustic resonances in single shot high-bandwidth pump-probe experiments is proposed. Our design utilizes concepts similar to smoothing by spectral dispersion; however, we propose a method of applying the bandwidth to the probe that avoids the effects of angular dispersion. Theoretical modeling and fluid simulations of the interaction are presented and compared. The results show the possibility of recovering a plasma's full complex index of refraction in a single shot pump-probe experiment. This wavelength dependent refractive index is directly related to the plasma's particle distribution function and can be used to recover plasma parameters including electron and ion temperatures, electron density, and plasma flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analytic marginalization of N(z) uncertainties in tomographic galaxy surveys

In this paper, we present a new method to marginalize over uncertainties in redshift distributions, N(z), within tomographic cosmological analyses applicable to current and upcoming photometric galaxy surveys. We allow for arbitrary deviations from the best-guess N(z) governed by a general covariance matrix describing the uncertainty in our knowledge of redshift distributions. In principle, this is marginalization over hundreds or thousands of new parameters describing potential deviations as a function of redshift and tomographic bin. However, by linearly expanding the theory predictions around a fiducial model, this marginalization can be performed analytically, resulting in a modified data covariance matrix that effectively downweights the modes of the data vector that are more sensitive to redshift distribution variations. We showcase this method by applying it to the galaxy clustering measurements from the Hyper Suprime-Cam first data release. We illustrate how to marginalize over sample variance of the calibration sample and a large general systematic uncertainty in photometric estimation methods, and explore the impact of priors imposing smoothness in the redshift distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Doped NiO: The mottness of a charge transfer insulator

The evolution of the electronic structures of strongly correlated insulators with doping has long been a central fundamental question in condensed matter physics; it is also of great practical relevance for applications. We have studied the evolution of NiO under hole and electron doping at low doping levels such that the system remains insulating using high-quality thin film and a wide range of experimental and theoretical methods. The evolution is in both cases very smooth with dopant concentration. The band gap is asymmetric under electron and hole doping, consistent with a charge-transfer insulator picture, and is reduced faster under hole doping than under electron doping. For both electron and hole doping, occupied states are introduced at the top of the valence band. The formation of deep donor levels under electron doping and the inability to pin otherwise empty states near the conduction-band edge are indicative that local electron addition and removal energies are dominated by a Mott-like Hubbard U interaction even though the global band gap is predominantly a charge-transfer-type gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Deploy threading in Nalu solver stack

The goal of the ExaWind project is to enable predictive simulations of wind farms composed of many MW-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary code in the ExaWind project is Nalu, which is an unstructured-grid solver for the acousticallyincompressible Navier-Stokes equations, and mass continuity is maintained through pressure projection. The model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as re-initialization of matrices and re-computation of preconditioners is required at every time step In this Milestone, we examine the effect of threading on the solver stack performance against flat-MPI results obtained from previous milestones using Haswell performance data full-turbine simulations. Whereas the momentum equations are solved only with the Trilinos solvers, we investigate two algebraic-multigrid preconditioners for the continuity equations: Trilinos/Muelu and HYPRE/BoomerAMG. These two packages embody smoothed-aggregation and classical Ruge-Stiiben AMG methods, respectively. In our FY18 Q2 report, we described our efforts to improve setup and solve of the continuity equations under flat-MPI parallelism. While significant improvement was demonstrated in the solve phase, setup times remained larger than expected. Starting with the optimized settings described in the Q2 report, we explore here simulation performance where OpenMP threading is employed in the solver stack. For Trilinos, threading is acheived through the Kokkos abstraction where, whereas HYPRE/BoomerAMG employs straight OpenMP. We examined results for our mid-resolution baseline turbine simulation configuration (229M DOF). Simulations on 2048 Haswell cores explored the effect of decreasing the number of MPI ranks while increasing the number of threads. Both HYPRE and Trilinos exhibited similar overal solution times, and both showed dramatic increases in simulation time in the shift from MPI ranks to OpenMP threads. This increase is attributed to the large amount of work per MPI rank starting at the single-thread configuration. Decreasing MPI ranks, while increasing threads, may be increasing simulation time due to thread synchronization and start-up overhead contributing to the latency and serial time in the model. These result showed that an MPI+OpenMP parallel decomposition will be more effective as the amount per MPI rank computation per MPI rank decreases and the communication latency increases. This idea was demonstrated in a strong scaling study of our low-resolution baseline model (29M DOF) with the Trilinos-HYPRE configuration. While MPI-only results showed scaling improvement out to about 1536 cores, engaging threading carried scaling improvements out to 4128 cores — roughly 7000 DOF per core. This is an important result as improved strong scaling is needed for simulations to be executed over sufficiently long simulated durations (i.e., for many timesteps). In addition to threading work described above, the team examined solver-performance improvements by exploring communication-overhead in the HYPRE-GMRES implementation through a communicationoptimal- GMRE algorithm (CO-GMRES), and offloading compute-intensive solver actions to GPUs. To those ends, a HYPRE mini-app was allow us to easily test different solver approaches and HYPRE parameter settings without running the entire Nalu code. With GPU acceleration on the Summitdev supercomputer, a 20x speedup was achieved for the overall preconditioner and solver execution time for the mini-app. A study on Haswell processors showed that CO-GMRES provides benefits as one increases MPI ranks.

17 WIND ENERGY↗
NASA NTRS1973

An efficient technique for determining apparent temperature distributions from antenna temperature measurements

A method by which the apparent microwave temperature characteristic of a flat scene is estimated from dual polarized measurements is derived and interpreted. Approximate linear relationships between antenna and apparent temperatures are established by weighting emission components in spherical bands under the assumption that the surface is isotropic. The weighting factors are formed by integrating the antenna pattern functions over these bands. The vector aspect of the formulation is retained to account for the difference between the definition of the antenna polarizations and the polarizations of the emitted fields. The method is largely applicable to the measurement of smooth temperature distributions by an antenna having good spatial resolution of the distributions and is considered efficient for inverting large volumes of measurements. Sample cases are presented and the implications of these cases on remote radiometer observations are discussed. It is shown that cross-coupling occurs between the polarizations of the emitted fields and the polarizations of the antenna. For this reason and because practical antennas have cross-polarized patterns associated with them, it is necessary to conduct measurements at both horizontal and vertical polarizations to realize the inversion. It is also made evident that thorough inversions require that the apparent temperatures be sampled at a sufficient number of points between nadir and zenith.

Claassen, J. P.↗
NASA NTRS1976

Design of active controls for the NASA F-8 digital fly-by-wire airplane

The design of a set of control laws for the NASA F-8 digital fly by wire research airplane is described. These control laws implement several active controls functions: maneuver load control, ride smoothing and departure boundary limiting. The criteria and methods which were used in the design of the control laws are also included. Results of linear analyses and nonlinear simulation are summarized.

Gera, J.↗
NASA NTRS1979

Gravitational spectra from direct measurements

A simple rapid method is described for determining the spectrum of a surface field (in spherical harmonics) from harmonic analysis of direct (in situ) measurements along great circle arcs. The method is shown to give excellent overall trends (smoothed spectra) to very high degree from even a few short arcs of satellite data. Three examples are taken with perfect measurements of satellite tracking over a planet made up of hundreds of point masses using (1) altimetric heights from a low-orbiting spacecraft, (2) velocity (range rate) residuals between a low and a high satellite in circular orbits, and (3) range rate data between a station at infinity and a satellite in a highly eccentric orbit. In particular, the smoothed spectrum of the earth's gravitational field is determined to about degree 400(50-km half wavelength) from 1 x 1 deg gravimetry and the equivalent of 11 revolutions of GEOS 3 and Skylab altimetry. This measurement shows that there is about 46 cm of geoid height (rms worldwide) remaining in the field beyond degree 180.

Wagner, C. A.↗
NASA NTRS1983

More-Uniform Heat Curing for Structural Repairs

Copper lining helps to distribute heat under heating blanket. Bond during setup employs heating blanket over copper sheet within vacuum bag. Copper sheet smoothes out temperature distribution under hotspots in heater. Method applicable to uniformly heating such low-thermal-conductivity metals as titanium or stainless steel.

Bauer, P. E.↗
NASA NTRS1981

A function space approach to smoothing with applications to model error estimation for flexible spacecraft control

A function space approach to smoothing is used to obtain a set of model error estimates inherent in a reduced-order model. By establishing knowledge of inevitable deficiencies in the truncated model, the error estimates provide a foundation for updating the model and thereby improving system performance. The function space smoothing solution leads to a specification of a method for computation of the model error estimates and development of model error analysis techniques for comparison between actual and estimated errors. The paper summarizes the model error estimation approach as well as an application arising in the area of modeling for spacecraft attitude control.

Rodriguez, G.↗