SEARCH · Engineering Papers

Results for “Smooth approximation”

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 NTRS1990

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, D. J.↗

NASA NTRS1993

A finite-volume high-order ENO scheme for two-dimensional hyperbolic systems

The finite-volume approach is presently used to obtain a 2D, high-order accurate and basically nonoscillatory shock-capture method whose high-order spatial accuracy is obtained by means of a piecewise polynomial approximation of the solution from cell averages. Attention is given to a high-order spatial operator that is able to both retain high-order accuracy in smooth regions and avoid the oscillations that are associated with interpolations across steep gradients. The operator is extended to hyperbolic systems of equations and curvilinear meshes.

Casper, Jay↗

NASA NTRS2019

A Year in the Life of GW 170817: The Rise and Fall of a Structured Jet from a Binary Neutron Star Merger

We present the results of our year-long afterglow monitoring of GW 170817, the first binary neutron star merger detected by Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) and Advanced Virgo. New observations with the Australian Telescope Compact Array and the Chandra X-ray Telescope were used to constrain its late-time behaviour. The broad-band emission, from radio to X-rays, is well-described by a simple power-law spectrum with index β ∼ 0.585 at all epochs. After an initial shallow rise ∝t0.9, the afterglow displayed a smooth turnover, reaching a peak X-ray luminosity of LX ≈ 5 × 1039 erg s−1 at 160 d, and has now entered a phase of rapid decline, approximately ∝t−2. The latest temporal trend challenges most models of choked jet/cocoon systems, and is instead consistent with the emergence of a relativistic structured jet seen at an angle of ≈22◦ from its axis. Within such model, the properties of the explosion (such as its blast wave energy EK ≈ 2 × 1050 erg, jet width θc ≈ 4◦, and ambient density n ≈ 3 × 10−3 cm−3) fit well within the range of properties of cosmological short gamma-ray bursts.

Acceleration of particles↗

DOE OSTI2026

Measurements of excitation functions and photoneutron cross sections of 96,98,99,104 Ru and 100 Mo

Photoneutron cross sections were extracted for isotopes of ruthenium and molybdenum on natural abundance targets via a reduced chi-squared analysis on excitation functions, or activation yields, as function of electron beam energy. Bremsstrahlung photons were produced via an electron linear accelerator and a tungsten radiator at the Idaho Accelerator Center. A total of 16 irradiations were performed with electron beam energies ranging from approximately 8 MeV to 23 MeV. Induced radioactivity was measured using a high purity germanium detector. The cross sections were extracted assuming a three parameter Lorentzian fit with a smooth truncated rise at reaction threshold.

07 - ISOTOPES AND RADIATION SOURCES↗

NASA NTRS1979

Precession of the epsilon ring of Uranus

It is noted that the outer and inner boundaries of the epsilon ring of Uranus can be fitted by aligned Keplerian ellipses. Four possible mechanisms for maintaining uniform precession in the epsilon ring are considered: the ring's self-gravity, precession due to a satellite, smooth pressure gradients, and shocklike phenomena. It is proposed that apse alignment is maintained by the self-gravity of the ring. In this case, a ring mass of approximately 5 x 10 to the 18th g and a mean surface density at quadrature of about 25 g/sq cm are estimated.

Goldreich, P.↗

NASA NTRS1981

General relaxation schemes in multigrid algorithms for higher order singularity methods

Relaxation schemes based on approximate and incomplete factorization technique (AF) are described. The AF schemes allow construction of a fast multigrid method for solving integral equations of the second and first kind. The smoothing factors for integral equations of the first kind, and comparison with similar results from the second kind of equations are a novel item. Application of the MD algorithm shows convergence to the level of truncation error of a second order accurate panel method.

Oskam, B.↗

NASA NTRS1989

Multigrid solution of the Navier-Stokes equations on triangular meshes

A new Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite-element Galerkin approximation, which can be shown to be equivalent to a finite-volume approximation for regular equilateral triangular meshes. Integration to steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. The convergence efficiency of the present method is shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, D. J.↗

NASA NTRS1921

Pressure Drop in Radiator Air Tubes

This report describes a method for measuring the drop in static pressure of air flowing through a radiator and shows (1) a reason for the discrepancy noted by various observers between head resistance and drop in pressure; (2) a difference in degree of contraction of the jet in entering a circular cell and a square cell; (3) the ratio of internal frictional resistance to total head resistance for two representative types; (4) the effect of smoothness of surface on pressure gradient; and (5) the effects of supplying heat to the radiator on pressure gradient. The fact that the pressure gradients are found to be approximately proportional to the square of the rate of flow of air appears to indicate turbulent flow, even in the short tubes of the radiator. It was found that the drop in the static pressure in the air stream through a cellular radiator and the pressure gradient in the air tubes are practically proportional to the square of the air flow in a given air density; that the difference between the head resistance per unit area and the fall of static pressure through the air tubes in radiators is apparent rather than real; and that radiators of different types differ widely in the amount of contraction of the jet at entrance. The frictional resistance was found to vary considerably, and in one case to be two-thirds of the head resistance in the type using circular cells and one-half of the head resistance of the radiator type using square cells of approximately the same dimensions.

Parsons, S R↗

NASA NTRS1995

Line-of-sight magnetic flux imbalances caused by electric currents

Several physical and observational effects contribute to the significant imbalances of magnetic flux that are often observed in active regions. We consider an effect not previously treated: the influence of electric currents in the photosphere. Electric currents can cause a line-of-sight flux imbalance because of the directionality of the magnetic field they produce. Currents associated with magnetic flux tubes produce larger imbalances than do smoothly-varying distributions of flux and current. We estimate the magnitude of this effect for current densities, total currents, and magnetic geometry consistent with observations. The expected imbalances lie approximately in the range 0-15%, depending on the character of the current-carying fields and the angle from which they are viewed. Observationally, current-induced flux imbalances could be indicated by a statistical dependence of the imbalance on angular distance from disk center. A general study of magnetic flux balance in active regions is needed to determine the relative importance of other- probably larger- effects such as dilute flux (too weak to measure or rendered invisible by radiative transfer effects), merging with weak background fields, and long-range connections between active regions.

Gary, G. Allen↗

NASA NTRS2014

Kilometer-Scale Topographic Roughness of Mercury: Correlation with Geologic Features and Units

We present maps of the topographic roughness of the northern circumpolar area of Mercury at kilometer scales. The maps are derived from range profiles obtained by the Mercury Laser Altimeter (MLA) instrument onboard the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) mission. As measures of roughness, we used the interquartile range of profile curvature at three baselines: 0.7 kilometers, 2.8 kilometers, and 11 kilometers. The maps provide a synoptic overview of variations of typical topographic textures. They show a dichotomy between the smooth northern plains and rougher, more heavily cratered terrains. Analysis of the scale dependence of roughness indicates that the regolith on Mercury is thicker than on the Moon by approximately a factor of three. Roughness contrasts within northern volcanic plains of Mercury indicate a younger unit inside Goethe basin and inside another unnamed stealth basin. These new data permit interplanetary comparisons of topographic roughness.

Kilometer-Scale↗

NASA NTRS2018

Summary of Ice Shape Geometric Fidelity Studies on an Iced Swept Wing

Understanding the aerodynamic impact of swept-wing ice accretions is a crucial component of the design of modern aircraft. Computer-simulation tools are commonly used to approximate ice shapes, so the necessary level of detail or fidelity of those simulated ice shapes must be understood relative to high-fidelity representations of the ice. Previous tests were performed in the NASA Icing Research Tunnel to acquire high-fidelity ice shapes. Some of those ice shapes are based on aircraft certification requirements. From this database, full-span artificial ice shapes were designed and manufactured for both an 8.9%-scale and 13.3%-scale semispan wing model of the CRM65 which has been established as the full-scale baseline for this swept-wing project. These models were tested in the Walter H. Beech wind tunnel at Wichita State University and at the ONERA (Office national d'etudes et de recherches aerospatiales) F1 facility, respectively. The data collected in the Wichita State University wind tunnel provided a low-Reynolds number baseline study while the pressurized F1 facility produced data over a wide range of Reynolds and Mach numbers with the highest Reynolds number studied being approximately Re = 11.9 by 10 (sup 6). Three different fidelity representations were created based on three different icing conditions. Lower-fidelity ice shapes were created by lofting a smooth ice shape between cross-section cuts of the high-fidelity ice shape. Grit roughness was attached to this smooth ice shape as another fidelity variant. The data indicates that the geometric fidelity of the ice shapes resulted in significant differences in lift and drag. These results were similar at both facilities over the wide range of test conditions utilized.

Aerodynamics↗

DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗

NASA NTRS2017

Looking at A 0535+26 at Low Luminosities with NuSTAR

We report on two NuSTAR observations of the high-mass X-ray binary A 0535+26 taken toward the end of its normal 2015 outburst at very low 3-50 keV luminosities of approximately 1.4 times x 10 (sup 36) ergs per second and approximately 5 times x 10 (sup 35) ergs per second, which are complemented by nine Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth exponential rollover at high energies seen in the first observation evolves to a much more abrupt steepening of the spectrum at 20-30 keV. The continuum evolution can be nicely described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from 0.75 plus or minus 0.04 to 0.56 plus 0.01 (sup) minus 0.04 (sub), the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy-resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at approximately 45 keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with the rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is, however, in disagreement with expectations.

Ballhausen, Ralf↗

DOE OSTI2025

Slicing Solutions for Wire Arc Additive Manufacturing

Both commercial and research applications of wire arc additive manufacturing (WAAM) have seen considerable growth in the additive manufacturing of metallic components. However, there remains a clear lack of a unified paradigm for toolpath generation when slicing parts for WAAM deposition. Existing toolpath generation options typically lack the appropriate features to account for all complexities of the WAAM process. This manuscript explores the key slicing challenges specific to toolpaths for WAAM geometry and pairs each consideration with multiple solutions to mitigate most negative effects on completed components. These challenges must be addressed to minimize voids, prevent bead collapse, and ensure deposited components accurately approximate the desired geometry. Slicing considerations are grouped into four general categories: geometric, process, thermal, and productivity. Geometric considerations are addressed with overhang compensation, corner-sharpening, and toolpath-smoothing features. Process considerations are addressed with start point configuration and controls for the bead lengths and end points. Thermal and productivity considerations are addressed with island optimization, multi-material printing, and connected insets. Finally, tools for the post-processing of generated G-code are explored. Overall, these solutions represent a critical set of slicing features used to improve generated toolpaths and the quality of the components deposited with those toolpaths.

36 MATERIALS SCIENCE↗

NASA NTRS1986

Aeorodynamic characteristics of an air-exchanger system for the 40- by 80-foot wind tunnel at Ames Research Center

A 1/50-scale model of the 40- by 80-Foot Wind Tunnel at Ames Research Center was used to study various air-exchange configurations. System components were tested throughout a range of parameters, and approximate analytical relationships were derived to explain the observed characteristics. It is found that the efficiency of the air exchanger could be increased (1) by adding a shaped wall to smoothly turn the incoming air downstream, (2) by changing to a contoured door at the inlet to control the flow rate, and (3) by increasing the size of the exhaust opening. The static pressures inside the circuit then remain within the design limits at the higher tunnel speeds if the air-exchange rate is about 5% or more. Since the model is much smaller than the full-scale facility, it is not possible to completely duplicate the tunnel, and it will be necessary to measure such characteristics as flow rate and tunnel pressures during implementation of the remodeled facility. The aerodynamic loads estimated for the inlet door and for nearby walls are also presented.

Rossow, V. J.↗

NASA NTRS1993

Application of a Reynolds stress model to separating boundary layers

Separating turbulent boundary layers occur in many practical engineering applications. Nonetheless, the physics of separation/reattachment of flows is poorly understood. During the past decade, various turbulence models were proposed and their ability to successfully predict some types of flows was shown. However. prediction of separating/reattaching flows is still a formidable task for model developers. The present study is concerned with the process of separation from a smooth surface. Features of turbulent separating boundary layers that are relevant to modeling include the following: the occurrence of zero wall shear stress, which causes breakdown of the boundary layer approximation; the law of the wall not being satisfied in the mean back flow region; high turbulence levels in the separated region; a significant low-frequency motion in the separation bubble; and the turbulence structure of the separated shear layer being quite different from that of either the mixing layers or the boundary layers. These special characteristics of separating boundary layers make it difficult for simple turbulence models to correctly predict their behavior.

Ko, Sung HO↗

NASA NTRS1999

Effects of Internal and External Hydrogen on Inconel 718

Internal hydrogen embrittlement (IHE) and hydrogen environment embrittlement (HEE) tensile and bend crack growth tests were performed on Inconel 718. For the IHE tests, the specimens were precharged to approximately 90 ppm hydrogen by exposure to 34.5 MPa H2 at 650 C. The HEE tests were performed in 34.5 MPa H2. Parameters evaluated were test temperature, strain rate for smooth and notch specimen geometries. The strain rate effect was very significant at ambient temperature for both IHE and HEE and decreased with increasing temperatures. For IHE, the strain rate effect was neglible at 260'C, and for HEE the strain rate effect was neglible at 400 C. At low temperatures, IHE was more severe than HEE, and at high temperatures HEE was more severe than IHE with a cross over temperature about 350 C. At 350 C, the equilibrium hydrogen concentration in Inconel 718 is about 50% lower than the hydrogen content of the precharged IHE specimens. Dislocation hydrogen sweeping of surface absorbed hydrogen was the likely transport mechanism for increasing the hydrogen concentration in the HEE tests sufficiently to produce the same degree of embrittlement as that of the more highly hydrogen charged IHE specimens. The main IHE fracture characteristic was formation of large, brittle flat facets, which decreased with increasing test temperature. The IHE fracture matrix surrounding the large facets ranged between brittle fine faceted to microvoid ductility depending upon strain rate, specimen geometry as well as temperature. The HEE fractures were characteristically fine featured, transgranular and brittle with a significant portion forming a "saw tooth" crystallographic pattern. Both IHE and HEE fractures were predominantly along the {1 1 1) slip and twin boundaries. With respect to embrittlement mechanism, it was postulated that dislocation hydrogen sweeping and hydrogen enhanced localized plasticity were active in HEE and IHE for concentrating hydrogen along (1 1 1) slip and twin planes. Final brittle failure occurred by hydrogen induced planer decohesion.

Walter, R. J.↗