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 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Efficient screening of rare large pit anomalies on polished surfaces using a minimalist sampling scheme

Lawrence Livermore National Laboratory (LLNL) has made significant strides in generating clean energy through its inertial confinement fusion (ICF) experiments. These experiments rely on high-density carbon (HDC) coated shells to encapsulate the fusion fuel. The success of these experiments is heavily dependent on the surface quality of these shells, as even minor imperfections, such as deep pits, can negatively impact fusion yield. Ensuring the required smoothness involves an extensive surface-finishing process that spans approximately 20 stages, making it both time-intensive and resource-demanding. A critical challenge in this process is the need for high-resolution scans to detect rare deep pits, which can be costly and impractical if performed on every shell. This highlights the necessity of developing more efficient scanning methods to optimize time and cost without compromising accuracy. To address these challenges, we introduce a novel approach that employs the multivariate Dvoretzky–Kiefer–Wolfowitz (DKW) inequality to provide a probabilistic upper bound on the error in estimating pit distribution characteristics via a Kernel Density Estimator (KDE). This error bound enables efficient and reliable estimation of pit distribution characteristics at a specified statistical confidence level using a minimal number of surface scans. The integrated DKW-KDE approach was validated through surface-finishing experiments across two batches of HDC-coated shells, demonstrating consistent and robust performance across multiple stages of the surface-finishing experiments. The validation studies suggest that the integrated DKW-KDE approach achieves comparable accuracy in estimating the risk of deleterious large pits with six scans, thus conserving time and resources. Further evaluations show that performance remains consistent across batches and over multiple polishing stages. In conclusion, based on these findings, one can leverage the minimal-scan insights to strategically improve the bottleneck inspection process, thus enhancing the productivity and quality of shell polishing and similar challenging manufacturing processes.

Inertial confinement fusion↗

NASA NTRS1984

Flux-vector splitting and Runge-Kutta methods for the Euler equations

Runge-Kutta schemes have been used as a method of solving the Euler equations exterior to an airfoil. In the past this has been coupled with central differences and an artificial vesocity in space. In this study the Runge-Kutta time-stepping scheme is coupled with an upwinded space approximation based on flux-vector splitting. Several acceleration techniques are also considered including a local time step, residual smoothing and multigrid.

Turkel, E.↗

NASA NTRS1988

Interpolation methods for shaped reflector analysis

The diffraction analysis of reflector surfaces which are described only at a discrete set of locations usually leads to the requirement of an interpolation to determine the surface characteristics over a continuum of locations. Two methods of interpolation, the global and the local methods, are presented. The global interpolation representation is a closed-form or series expression valid over the entire surface. The coefficients of a series expression are found by an integration of all of the raw data. Since the number of coefficients used to describe the surface is much smaller than the number of raw data points, the integration effectively provides a smoothing of the raw data. The local interpolation provides a closed-form expression for only a small area of the reflector surface. The subreflector is divided into sectors each of which has constant discretized data. Each area segment is then locally described by a two-dimensional quadratic surface. The second derivative data give the desired smoothed values.

Galindo-Israel, Victor↗

NASA NTRS1996

Enhancing Adhesion: Relative Merits of Different Approaches

Adhesive performance is improved mainly by manipulation of the bimaterials interface zone, which is only a few molecules thick. There are three approaches to enhancement of interfacial adhesion at the molecular level. They are 1) changing the nonchemically bonded interactions across the interface from weak ones to strong ones, 2) making the true interfacial area much larger than the simple geometric area, and 3) inducing chemical bonding between the two materials forming the interface. Our goal this summer was to question some of the built-in assumptions contained within these approaches and to determine the most promising approach, both theoretically and practically, for enhancing adhesion in NASA structures. Our computations revealed that all three of these approaches have, in theory, the potential to enhance molecular adhesion approximately ten-fold. Experiments, however, revealed that this excellent level of enhancement is not likely to be reached in practice. Each approach was found to be severely limited by practical problems. In addition, some of the built-in assumptions associated with these approaches were found to be insufficient or inadequate. The first approach, changing the nonchemically bonded interactions from weak to strong, Is an example of one containing inadequate assumptions. The extensive literature on intermolecular interactions, based on solution studies, shows that certain functional group pairs interact much more strongly than others. It has always been assumed that these data can be reliably extended to systems where only one member of the pair is in solution and the other Is contained in a solid surface. Our experiments this summer demonstrated that solution data do not adequately predict the strength of functional group interaction at the solid-liquid interface. Furthermore, the strong solvents needed to dissolve the monomers or polymers to which the functional groups of interest are attached compete successfully with the solid surface for the functional group. As a result, functional groups in solution cannot pair with the complementary groups in the solid surface, and the expected enhancement of nonchemically bonded interactions is not realized. The second approach, increasing the true interfacial area, is an example of one containing inadequate assumptions and suffering from numerous practical problems. First, practitioners have assumed that material removal, such as bead blasting or etching, increases true surface area (and therefore interfacial area) in a meaningful way. Our geometric analysis demonstrated that removal methods increase area by a factor of two at most. To increase interfacial area by an order of magnitude or more. a thin layer of high porosity must be added to the substrate surface prior to application of the adhesive phase. Consistent with this finding, we attempted to create a thin layer of rigid, highly porous glass on the surface of our smooth glass substrate by means of sol-gel technology. We were unable to surmount a wide variety of practical problems and obtained only collapsed, nonporous layers. Thus this approach, appealing in principle, would require long term development and is not promising in the near term. The third approach, inducing chemical bonding at the interface, is an example of one having neither inadequate assumptions nor insurmountable practical problems. When silicate glass is the substrate, there are only a few chemical reactions that can be successfully conducted to create these chemical bonds, and these reactions usually involve silicon-containing reagents. We compared the silazane reagents to the silane reagents and found through experiment that the silazanes react with the glass surface much more readily, and under milder conditions, than the silanes. The functional groups attached to the glass surface by silazane reactions were not able to be removed by solvent extraction, elevated temperature exposure, or mechanical action. This clearly indicates that the formation of chemical bonds at the interface is the most effective approach for enhancing molecular adhesion.

Penn, L. S.↗

NASA NTRS2014

A Computational and Experimental Study of Coflow Laminar Methane/Air Diffusion Flames: Effects of Fuel Dilution, Inlet Velocity, and Gravity

The influences of fuel dilution, inlet velocity, and gravity on the shape and structure of laminar coflow CH4-air diffusion flames were investigated computationally and experimentally. A series of nitrogen-diluted flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) on board the International Space Station was assessed numerically under microgravity (mu g) and normal gravity (1g) conditions with CH4 mole fraction ranging from 0.4 to 1.0 and average inlet velocity ranging from 23 to 90 cm/s. Computationally, the MC-Smooth vorticity-velocity formulation was employed to describe the reactive gaseous mixture, and soot evolution was modeled by sectional aerosol equations. The governing equations and boundary conditions were discretized on a two-dimensional computational domain by finite differences, and the resulting set of fully coupled, strongly nonlinear equations was solved simultaneously at all points using a damped, modified Newton's method. Experimentally, flame shape and soot temperature were determined by flame emission images recorded by a digital color camera. Very good agreement between computation and measurement was obtained, and the conclusions were as follows. (1) Buoyant and nonbuoyant luminous flame lengths are proportional to the mass flow rate of the fuel mixture; computed and measured nonbuoyant flames are noticeably longer than their 1g counterparts; the effect of fuel dilution on flame shape (i.e., flame length and flame radius) is negligible when the flame shape is normalized by the methane flow rate. (2) Buoyancy-induced reduction of the flame radius through radially inward convection near the flame front is demonstrated. (3) Buoyant and nonbuoyant flame structure is mainly controlled by the fuel mass flow rate, and the effects from fuel dilution and inlet velocity are secondary.

flames↗

DOE OSTI2022

Deploying Traffic Smoothing Cruise Controllers Learned from Trajectory Data

Autonomous vehicle-based traffic smoothing con- trollers are often not transferred to real-world use due to challenges in calibrating many-agent traffic simulators. We show a pipeline to sidestep such calibration issues by collecting trajectory data and learning controllers directly from trajectory data that are then deployed zero-shot onto the highway. We construct a dataset of 772.3 kilometers of recorded drives on the I-24. We then construct a simple simulator using the recorded drives as the lead vehicle in front of a simulated platoon consisting of one autonomous vehicle and five human followers. Using policy-gradient methods with an asymmetric critic to learn the controller, we show that we are able to improve average MPG by 11% in simulation on congested trajectories. We deploy this controller to a mixed platoon of 4 autonomous Toyota RAV-4’s and 7 human drivers in a validation experiment and demonstrate that the expected time-gap of the controller is maintained in the real world test. Finally, we release the driving dataset [1], the simulator, and the trained controller at https://github.com/nathanlct/trajectory-training-icra.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2025

The Cosmic Evolution of C IV Absorbers at 1.4 < z < 4.5: Insights from 100,000 Systems in DESI Quasars

We present the largest catalog to date of triply ionized carbon (C IV ) absorbers detected in quasar spectra from the Dark Energy Spectroscopic Instrument. Using an automated matched-kernel convolution method with adaptive signal-to-noise thresholds, we identify 101,487 C IV systems in the redshift range 1.4 < z < 4.5 from 300,637 quasar spectra. Completeness is estimated via Monte Carlo simulations, and the catalog is 50% complete at EW C IV ≥ 0.4 Å. The differential equivalent width frequency distribution declines exponentially and shows weak redshift evolution. The absorber incidence per unit comoving path increases by a factor of 2–5 from z ≈ 4.5 to z ≈ 1.4, with stronger redshift evolution for strong systems. Using column densities derived from the apparent optical depth method, we constrain the cosmic mass density of C IV , Ω C IV , which increases by a factor of ∼3.8 from (0.82 ± 0.05) × 10 −8 at z ≈ 4.5 to (3.16 ± 0.2) × 10 −8 at z ≈ 1.4. From Ω C IV , we estimate a lower limit on intergalactic medium metallicity ${\mathrm{log}}({Z}_{{\rm{IGM}}}/{Z}_{\odot })\gtrsim -3.25$ at z ∼ 2.3, with a smooth decline at higher redshifts. These trends trace the cosmic star formation history and He II photoheating rate, suggesting a link between C IV enrichment, star formation, and UV background over ∼3 Gyr. The catalog also provides a critical resource for future studies connecting circumgalactic metals to galaxy evolution, especially near cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1981

A multigrid method for the transonic full potential equation discretized with finite elements on an arbitrary body fitted mesh

A multigrid method for the acceleration of transonic potential flow calculations based on a Galerkin finite element approach is described. In order to allow the use of arbitrary body fitted meshes it is necessary to introduce nonuniform interpolation and residual weighting. Emphasis is put on the construction of these operators consistent with the finite element approximation, while standard successive line overrelaxation is used as a smoothing step. Substantial convergence acceleration is obtained and results are presented for different transonic flow configurations including shocks.

Deconinck, H.↗

NASA NTRS1996

Convergence Acceleration of the Navier-Stokes Equations Through Time-Derivative Preconditioning

Chorin's method of artificial compressibility is extended to both compressible and incompressible fluids by using physical arguments to define artificial fluid properties that make up a local preconditioning matrix. In particular, perturbation expansions are used to provide appropriate temporal derivatives for the equations of motion at both low speeds and low Reynolds numbers. These limiting forms are then combined into a single function that smoothly merges into the physical time derivatives at high speeds so that the equations are left unchanged at transonic, high Reynolds number conditions. The effectiveness of the resulting preconditioning procedures for the Navier-Stokes equations is demonstrated for a wide speed and Reynolds number ranges by means of stability results and computational solutions. Nevertheless, the preconditioned equations sometimes fail to provide a solution for applications for which the non-preconditioned equations converge. Often this is because the reduced dissipation in the preconditioned equations results in an unsteady solution while the more dissipative non-preconditioned equations result in a steady state. Problems of this type represent a computational challenge; it is important to distinguish between non-convergence of algorithms, and the non-existence of steady state solutions.

Merkle, Charles L.↗

NASA NTRS2016

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Tanks

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involve the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations and knowledge were mainly carried out by extensive experimental studies. A Volume-Of-Fluid (VOF) based CFD program developed at NASA MSFC was applied to extract slosh damping in a baffled tank from the first principle. First, experimental data using water with subscale smooth wall tank were used as the baseline validation. CFD simulation was demonstrated to be capable of accurately predicting natural frequency and very low damping value from the smooth wall tank at different fill levels. The damping due to a ring baffle at different liquid fill levels from barrel section and into the upper dome was then investigated to understand the slosh damping physics due to the presence of a ring baffle. Based on this study, the Root-Mean-Square error of our CFD simulation in estimating slosh damping was less than 4.8%, and the maximum error was less than 8.5%. Scalability of subscale baffled tank test using water was investigated using the validated CFD tool, and it was found that unlike the smooth wall case, slosh damping with baffle is almost independent of the working fluid and it is reasonable to apply water test data to the full scale LOX tank when the damping from baffle is dominant. On the other hand, for the smooth wall, the damping value must be scaled according to the Reynolds number. Comparison of experimental data, CFD, with the classical and modified Miles equations for upper dome was made, and the limitations of these semi-empirical equations were identified.

Yang, H. Q.↗

DOE OSTI2026

Unveiling horizons in quantum critical collapse

Critical gravitational collapse offers a unique window into regimes of arbitrarily high curvature, culminating in a naked singularity arising from smooth initial data — thus providing a dynamical counterexample to weak cosmic censorship. Near the critical regime, quantum effects from the collapsing matter are expected to intervene before full quantum gravity resolves the singularity. Despite its fundamental significance, a self-consistent treatment has so far remained elusive. In this work, we perform a one-loop semiclassical analysis using the robust anomaly-based method in the canonical setup of Einstein gravity minimally coupled to a free, massless scalar field. Focusing on explicitly solvable near-critical solutions in both 2 + 1 and 3 + 1 dimensions, we analytically solve the semiclassical Einstein equations and obtain controlled, quantitative results for several long-standing questions within the dominant s-wave sector. We find that regularity uniquely selects a Boulware-like quantum state, encoding genuine vacuum polarization effects from the collapsing matter. Remarkably, the resulting quantum corrections manifest as a growing mode. Horizon-tracing analyses, incorporating both classical and quantum modes, reveal the emergence of a finite mass gap, signaling a phase transition from classical Type II to quantum-modified Type I behavior, thereby providing a quantum enforcement of the weak cosmic censorship. The most non-trivial aspect of our analysis involves dealing with non-conformal matter fields in explicitly time-dependent critical spacetimes. Along the way, we uncover intriguing and previously underexplored features of quantum field theory in curved spacetime.

2D Gravity↗

DOE OSTI2026

Radiative strength functions from the energy-localized Brink-Axel hypothesis

Radiative strength functions (RSFs) model the bulk electromagnetic response of highly excited nuclei and are critical inputs for statistical reaction codes. In this paper, we present a definition of the RSF that is consistent with Hauser-Feshbach reaction codes and that can be efficiently computed with the shell model using the Lanczos strength-function (LSF) method. Here, we introduce a variant of the shell-model LSF method that exploits the energy-localized Brink-Axel hypothesis, which makes it possible to compute both electric and magnetic RSFs across all energies relevant to capture reactions. We verify agreement with the conventional definition of RSFs with benchmark calculations of 24 Mg and then present novel results for 56 Fe. For 56 Fe we find that (i) the 𝑀⁢1 RSF shape evolves smoothly with excitation energy, consistent with the energy-localized Brink-Axel hypothesis; (ii) both 𝑀⁢1 and 𝐸⁢1 transitions contribute significantly to the radiative strength below the photoabsorption threshold; and (iii) within the sdpf model space, the strength below 3 MeV observed in Oslo-type experiments cannot be fully reproduced. These results pave the way for a coherent microscopic description of RSFs and further motivate the use of energy-dependent RSFs in modern reaction codes.

Physics - Nuclear physics and radiation physics↗

NASA NTRS1981

Precision heat forming of tetrafluoroethylene tubing

An invention that provides a method of altering the size of tetrafluoroethylene tubing which is only available in limited combination of wall thicknesses and diameter are discussed. The method includes the steps of sliding the tetrafluoroethylene tubing onto an aluminum mandrel and clamping the ends of the tubing to the mandrel by means of clamps. The tetrafluorethylene tubing and mandrel are then placed in a supporting coil which with the mandrel and tetrafluorethylene tubing are then positioned in a insulated steel pipe which is normally covered with a fiber glass insulator to smooth out temperature distribution therein. The entire structure is then placed in an event which heats the tetrafluorethylene tubing which is then shrunk by the heat to the outer dimension of the aluminum mandrel. After cooling the aluminum mandrel is removed from the newly sized tetrafluorethylene tubing by a conventional chemical milling process.

Ruiz, W. V.↗

NASA NTRS1987

Method for machining holes in composite materials

A method for boring well defined holes in a composite material such as graphite/epoxy is discussed. A slurry of silicon carbide powder and water is projected onto a work area of the composite material in which a hole is to be bored with a conventional drill bit. The silicon carbide powder and water slurry allow the drill bit, while experiencing only normal wear, to bore smooth, cylindrical holes in the composite material.

Daniels, Julia G.↗

NASA NTRS1924

The Small Angular Oscillations of Airplanes in Steady Flight

This investigation was carried out by the National Advisory Committee for Aeronautics at the request of the Army Air Service to provide data concerning the small angular oscillations of several types of airplanes in steady flight under various atmospheric conditions. The data are of use in the design of bomb sights and other aircraft instruments. The method used consisted in flying the airplane steadily in one direction for at least one minute, while recording the angle of the airplane with the sun by means of a kymograph. The results show that the oscillations differ but little for airplanes of various types, but that the condition of the atmosphere is an important factor. The average angular excursion from the mean in smooth air is 0.8 degrees in pitch, 1.4 degrees in roll, and 0.9 degrees in yaw, without special instruments to aid the pilot in holding steady conditions. In bumpy air the values given above are increased about 50 per cent. (author)

Norton, F H↗

NASA NTRS1994

Numerical study of the effects of icing on viscous flow over wings

An improved hybrid method for computing unsteady compressible viscous flows is presented. This method divides the computational domain into two zones. In the outer zone, the unsteady full-potential equation (FPE) is solved. In the inner zone, the Navier-Stokes equations are solved using a diagonal form of an alternating-direction implicit (ADI) approximate factorization procedure. The two zones are tightly coupled so that steady and unsteady flows may be efficiently solved. Characteristic-based viscous/inviscid interface boundary conditions are employed to avoid spurious reflections at that interface. The resulting CPU times are less than 60 percent of that required for a full-blown Navier-Stokes analysis for steady flow applications and about 60 percent of the Navier-Stokes CPU times for unsteady flows in non-vector processing machines. Applications of the method are presented for a rectangular NACA 0012 wing in low subsonic steady flow at moderate and high angles of attack, and for an F-5 wing in steady and unsteady subsonic and transonic flows. Steady surface pressures are in very good agreement with experimental data and are essentially identical to Navier-Stokes predictions. Density contours show that shocks cross the viscous/inviscid interface smoothly, so that the accuracy of full Navier-Stokes equations can be retained with a significant savings in computational time.

Sankar, L. N.↗

NASA NTRS1995

A Method for Exploiting Redundancy to Accommodate Actuator Limits in Multivariable Systems

This paper introduces a new method for accommodating actuator saturation in a multivariable system with actuator redundancy. Actuator saturation can cause significant deterioration in control system performance because unmet demand may result in sluggish transients and oscillations in response to setpoint changes. To help compensate for this problem, a technique has been developed which takes advantage of redundancy in multivariable systems to redistribute the unmet control demand over the remaining useful effectors. This method is not a redesign procedure, rather it modifies commands to the unlimited effectors to compensate for those which are limited, thereby exploiting the built-in redundancy. The original commands are modified by the increments due to unmet demand, but when a saturated effector comes off its limit, the incremental commands disappear and the original unmodified controller remains intact. This scheme provides a smooth transition between saturated and unsaturated modes as it divides up the unmet requirement over any available actuators. This way, if there is sufficiently redundant control authority, performance can be maintained.

Litt, Jonathan↗

DOE OSTI2020

Illuminating dark matter halo density profiles without subhaloes

Cold dark matter haloes consist of a relatively smooth dark matter component as well as a system of bound subhaloes. It is the prevailing practice to include all mass, including mass in subhaloes, in studies of halo density profiles in simulations. However, often in observational studies satellites are treated as having their own distinct dark matter density profiles in addition to the profile of the host. This difference can make comparisons between theoretical and observed results difficult. In this work, we investigate density profiles of the smooth components of host haloes by excluding mass contained within subhaloes. We find that the density profiles of the smooth halo component (without subhaloes) differ substantially from the conventional halo density profile, declining more rapidly at large radii. We also find that concentrations derived from smooth density profiles exhibit less scatter at fixed mass and a weaker mass dependence than standard concentrations. Both smooth and standard halo profiles can be described by a generalized Einasto profile, an Einasto profile with a modified central slope, with smaller residuals than either a Navarro–Frenk–White or Einasto profile. Furthermore, these results hold for both Milky Way-mass and cluster-mass haloes. This new characterization of smooth halo profiles can be useful for many analyses, such as lensing and dark matter annihilation, in which the smooth and clumpy components of a halo should be accounted for separately.

79 ASTRONOMY AND ASTROPHYSICS↗