SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Long-Range Dispersion Governs the Behavior of Near-Critical Fluids: Universal Scaling and Implications for Accurate Molecular Simulation

Computer simulations of near-critical and supercritical fluids often deviate from experimental results, a discrepancy commonly attributed to force-field resolutions and inaccuracies. We demonstrate that these errors primarily stem from using finite cutoff schemes for dispersion interactions, which become effectively long-range as the correlation length grows near the critical point. By employing the smooth particle-mesh Ewald method to account for full long-range dispersion, we show that coarse-grained models can also achieve high quantitative accuracy for n-undecane and propylbenzene. This approach enables the precise determination of density fluctuations and correlation lengths, which exhibit universal critical scaling and depend on the molecular size and shape. Our findings resolve a longstanding debate in molecular simulations of near-critical fluids, establishing that capturing long-wavelength fluctuations is essential for bridging the gap between the microscopic behavior and macroscopic critical phenomena.

Chen, Guang↗

DOE OSTI2023

Surface Merging Technique to Design GA-Optimized Heat Sinks

High-power density inverters require efficient and small heat dissipation systems. For liquid cooling systems, an effective heat sink design is important to enable higher heat transfer and keep the pressure drop within reasonable limits. A heat sink geometry generated from an extruded, 1D fast Fourier transform can achieve more even temperature distribution and improved heat transfer than a conventional heat sink. Building on this concept, this paper proposes a new heat sink geometry creation algorithm. Instead of one plane profile being extruded, multiple plane profiles, each created from the fast Fourier transform method, are subsequently merged by smooth surfaces along the flow direction of the heat sink. Thus, a 3D heat sink geometry is created that enables 3D coolant flow and efficient heat transfer. The volume of this design was reduced by 50% compared with pin fin heat sinks. Compared with extruded 1D fast Fourier transform heat sinks, the proposed design showed an approximately 5% device temperature reduction and a 20% pressure drop reduction.

Barua, Himel↗

DOE OSTI2024

Machine-assisted discovery of integrable symplectic mappings

Integrable systems possess a hidden symmetry associated with the existence of conserved quantities known as integrals of motion. These systems play an important role in understanding general dynamics in accelerators and have potential for future designs. This work will cover two automated methods for finding integrable symplectic maps of the plane. The first algorithm is based on the observation that the evolution of an integrable system in phase space is confined to a lower-dimensional submanifold of a specific type. The second algorithm relies on an analysis of dynamical variables. Both methods rediscover some of the famous McMillan-Suris integrable mappings and ultra-discrete Painlev\'e equations. Over 100 new integrable families are presented and analyzed, some of which are isolated in the space of parameters, while others are families with one parameter (or the ratio of parameters) being either continuous or discrete. In addition, the newly discovered maps are related to a general 2D symplectic map through the use of discrete perturbation theory. A method is proposed for constructing smooth near-integrable dynamical systems based on mappings with polygon invariants.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Machine-Assisted Discovery of Integrable Symplectic Mappings

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1984

Calculating C-grids with fine and embedded mesh regions

A program for calculating a C-type mesh around airfoil like shapes is described. The Jameson/Caughey approach is used: a parabolic transformation coupled with a shearing transformation. The algebraic algorithm is capable of efficiently generating nearly orthogonal grids. A high degree of grid control is possible. The user may specify grid boundaries, number of grid lines, and location of (and spacing in) trailing edge and leading edge fine mesh areas. The capability of embedding fine mesh regions, for use with new adaptive grid techniques, is being developed. Grids generated by the program were used in Euler flow flow calculatons by W. Usab. Results superior to results calculated on previous O-type grids were obtained. Specifically, calculations converged faster using C-grids rather than 0-grids, total pressure loss spikes at the trailing edge of the airfoil were eliminated, and the Ni method converged with zero artificial smoothing for a subcritical case (resulting overall total pressure loss was then nearly zero). These improvements were attributed to higher grid orthogonality, especially at the trailing edge. The program itself is fairly straightforward. Roughly half of the 800 code lines are comment lines.

Loyd, B.↗

NASA NTRS1985

Microwave inversion of leaf area and inclination angle distributions from backscattered data

The backscattering coefficient from a slab of thin randomly oriented dielectric disks over a flat lossy ground is used to reconstruct the inclination angle and area distributions of the disks. The disks are employed to model a leafy agricultural crop, such as soybeans, in the L-band microwave region of the spectrum. The distorted Born approximation, along with a thin disk approximation, is used to obtain a relationship between the horizontal-like polarized backscattering coefficient and the joint probability density of disk inclination angle and disk radius. Assuming large skin depth reduces the relationship to a linear Fredholm integral equation of the first kind. Due to the ill-posed nature of this equation, a Phillips-Twomey regularization method with a second difference smoothing condition is used to find the inversion. Results are obtained in the presence of 1 and 10 percent noise for both leaf inclination angle and leaf radius densities.

Lang, R. H.↗

NASA NTRS1988

A Navier-Stokes solver for cascade flows

A computer code for solving the Reynolds averaged full Navier-Stokes equations has been developed and applied using sheared H-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady state calculations. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications in less than 100 multigrid cycles on the finest mesh.

Arnone, A.↗

NASA NTRS1991

Multigrid calculation of three-dimensional viscous cascade flows

A 3-D code for viscous cascade flow prediction was developed. The space discretization uses a cell-centered scheme with eigenvalue scaling to weigh the artificial dissipation terms. Computational efficiency of a four stage Runge-Kutta scheme is enhanced by using variable coefficients, implicit residual smoothing, and a full multigrid method. The Baldwin-Lomax eddy viscosity model is used for turbulence closure. A zonal, nonperiodic grid is used to minimize mesh distortion in and downstream of the throat region. Applications are presented for an annular vane with and without end wall contouring, and for a large scale linear cascade. The calculation is validated by comparing with experiments and by studying grid dependency.

Arnone, A.↗

NASA NTRS1993

A Navier-Stokes solver for turbomachinery applications

A computer code for solving the Reynolds-averaged full Navier-Stokes equations has been developed and applied using H- and C-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady-state calculations. A grid independence analysis is presented for a transonic rotor blade. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications.

Arnone, A.↗

NASA NTRS2003

Measurement and Correlation of Ice Accretion Roughness

Measurements were taken of the roughness characteristics of ice accreted on NACA 0012 airfoils in the NASA Glenn Icing Research Tunnel (IRT). Tests were conducted with size scaled, using models with chords of 26.7, 53.3, and 80.0 cm, and with liquid-water content scaled, both according to previously-tested scaling methods. The width of the smooth zone which forms on either side of the leading edge of the airfoil and the diameter of the roughness elements are presented in non-dimensional form as functions of the accumulation parameter. The smooth-zone width was found to decrease with increasing accumulation parameter. The roughness-element diameter increased with accumulation parameter until a plateau was reached. This maximum diameter was about 0.06 times twice the model leading-edge radius. Neither smooth-zone width nor element diameter were affected by a change in freezing fraction from 0.2 to 0.4. Both roughness characteristics appeared to scale with model size and with liquid-water content.

Anderson, David N.↗

NASA NTRS2004

Friction and Wear Characteristics of Cu-4Al Foil Bearing Coating at 25 and 650 degree C

The friction and wear performance of a Cu-4Al top foil coating has been investigated in Generation I foil air bearings. The copper alloy was applied by a novel deposition technique (ion diffusion) and the journal was coated with PS304, a plasma spray deposited high temperature composite solid lubricant coating. The ion diffusion coating process deposits a desirable smooth layer compared to other methods like cathodic arc deposition. The tribological performance of bearings with and without Cu-4Al foil coatings were evaluated through start-stop tests on an air bearing test rig at 25 and 650 C. The results indicate that the Cu-4Al assists during the initial break-in period, gives more stable friction performance with respect to temperature, and appears to prevent top foil wear at high temperature. The measured load capacity coefficient was 0.5, which was comparable to earlier testing of more advanced design Generation III bearings coated with standard cathodic arc deposited Cu-4Al. However, further studies are needed to determine if deeper penetration of the copper alloy into the foil would help make the transition in friction behavior from contact with the Cu-4Al coated foil to contact with the base foil material more gradual. Also, future work is recommended to assess the performance of ion diffusion coatings with different Cu-based alloy compositions and to investigate the effect the coating has on the elastic modulus of the foil material.

Stanford, Malcolm K.↗

NASA NTRS2012

Uniform Foam Crush Testing for Multi-Mission Earth Entry Vehicle Impact Attenuation

Multi-Mission Earth Entry Vehicles (MMEEVs) are blunt-body vehicles designed with the purpose of transporting payloads from outer space to the surface of the Earth. To achieve high-reliability and minimum weight, MMEEVs avoid use of limited-reliability systems, such as parachutes and retro-rockets, instead using built-in impact attenuators to absorb energy remaining at impact to meet landing loads requirements. The Multi-Mission Systems Analysis for Planetary Entry (M-SAPE) parametric design tool is used to facilitate the design of MMEEVs and develop the trade space. Testing was conducted to characterize the material properties of several candidate impact foam attenuators to enhance M-SAPE analysis. In the current effort, four different Rohacell foams are tested at three different, uniform, strain rates (approximately 0.17, approximately 100, approximately 13,600%/s). The primary data analysis method uses a global data smoothing technique in the frequency domain to remove noise and system natural frequencies. The results from the data indicate that the filter and smoothing technique are successful in identifying the foam crush event and removing aberrations. The effect of strain rate increases with increasing foam density. The 71-WF-HT foam may support Mars Sample Return requirements. Several recommendations to improve the drop tower test technique are identified.

Patterson, Byron W.↗

DOE OSTI2023

Control of nanostructure and activity by alloying and/or segregation

A method for synthesis of PtNi smooth surface core/shell particles or Nano cages and porous nanocages from segregated nanoparticles.

Wang, Rongyue↗

DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Universal Passive Synchronization Method for Grid-Forming Inverters Without Mode Transition

Power systems are transforming with increasing levels of inverter-based resources (IBRs). This transformation requires critical roles of grid-forming (GFM) inverters replacing synchronous generators for bulk power system stabilization and ancillary services, also allowing flexible power system operation, such as microgrid that is operated by multiple GFM IBRs to achieve system resilience against contingencies. To realize the resilient power systems allowing flexible in-and-out operation of GFM IBRs potentially programmed with different primary controls, a synchronization method universally applicable, i.e., independent of control types, would be beneficial to ease the integration process, but it has not been actively studied. To fill the gap, this paper proposes a universal synchronization method that achieves a passive synchronization to enable a smooth transition in a grid with off-nominal system parameters, i.e., voltage and frequency. The logic proposed requires no modification on the primary control, thus applicable to any type of GFMs with a voltage reference input. To validate the concept, a simulation of an IEEE 13-bus benchmark system modified with 3 GFM inverters is presented. It simulates an inverter-driven black start scenario in which GFM inverters autonomously turn on and connect to the grid under heavy loading, using the synchronization logic. The case study demonstrates that GFM inverters can tune their voltage reference to smoothly synchronize without severe transients, and contribute to a seamless black start of the grid under unbalanced load conditions. Two GFM methods - Droop and dispatchable virtual oscillator control - are used for the demo to validate feasibility and interoperability of the passive synchronization.

black start↗

NASA NTRS1973

Aerodynamic influence coefficient method using singularity splines.

A new numerical formulation with computed results, is presented. This formulation combines the adaptability to complex shapes offered by paneling schemes with the smoothness and accuracy of the loading function methods. The formulation employs a continuous distribution of singularity strength over a set of panels on a paneled wing. The basic distributions are independent, and each satisfies all of the continuity conditions required of the final solution. These distributions are overlapped both spanwise and chordwise (termed 'spline'). Boundary conditions are satisfied in a least square error sense over the surface using a finite summing technique to approximate the integral.

Mercer, J. E.↗

DOE OSTI2021

A primal–dual algorithm for risk minimization

In this paper, we develop an algorithm to efficiently solve risk-averse optimization problems posed in reflexive Banach space. Such problems often arise in many practical applications as, e.g., optimization problems constrained by partial differential equations with uncertain inputs. Unfortunately, for many popular risk models including the coherent risk measures, the resulting risk-averse objective function is nonsmooth. Here, this lack of differentiability complicates the numerical approximation of the objective function as well as the numerical solution of the optimization problem. To address these challenges, we propose a primal–dual algorithm for solving large-scale nonsmooth risk-averse optimization problems. This algorithm is motivated by the classical method of multipliers and by epigraphical regularization of risk measures. As a result, the algorithm solves a sequence of smooth optimization problems using derivative-based methods. We prove convergence of the algorithm even when the subproblems are solved inexactly and conclude with numerical examples demonstrating the efficiency of our method.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1987

Optimal application of Morrison's iterative noise removal for deconvolution. Appendices

Morrison's iterative method of noise removal, or Morrison's smoothing, is applied in a simulation to noise-added data sets of various noise levels to determine its optimum use. Morrison's smoothing is applied for noise removal alone, and for noise removal prior to deconvolution. For the latter, an accurate method is analyzed to provide confidence in the optimization. The method consists of convolving the data with an inverse filter calculated by taking the inverse discrete Fourier transform of the reciprocal of the transform of the response of the system. Various length filters are calculated for the narrow and wide Gaussian response functions used. Deconvolution of non-noisy data is performed, and the error in each deconvolution calculated. Plots are produced of error versus filter length; and from these plots the most accurate length filters determined. The statistical methodologies employed in the optimizations of Morrison's method are similar. A typical peak-type input is selected and convolved with the two response functions to produce the data sets to be analyzed. Both constant and ordinate-dependent Gaussian distributed noise is added to the data, where the noise levels of the data are characterized by their signal-to-noise ratios. The error measures employed in the optimizations are the L1 and L2 norms. Results of the optimizations for both Gaussians, both noise types, and both norms include figures of optimum iteration number and error improvement versus signal-to-noise ratio, and tables of results. The statistical variation of all quantities considered is also given.

Ioup, George E.↗