Engineering Papers⌕ Search

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.

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

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

The creation of local clusters in arbitrarily given grids

A method is presented to smoothly insert pointwise clusters into any given grid regardless of its origin, its topology, or its dimensionality. The process amounts to a local movement of the given coordinate curves or surfaces to more highly resolve an object. The object about which clustering is created can be a point, a curve, a surface, or segments of a curve or surface. The basic clustering capability is established by forming a grid operator for a single cluster. With a view toward multiple clusters being created about various objects, the basic operator is seen as an elementary operator. An algorithm is presented to execute the general elementary operation in three dimensions. In FORTRAN, this assumes the form of a subroutine which is fully operational and is presented to serve as a basic model for any such elementary clustering operation.

Eiseman, Peter R.↗
NASA NTRS1988

Magnetic reversal spurts: Rain gauges for comet showers

Abrupt increases in the rate of magnetic reversals (magnetic reversal spurts) were first studied by many others. They hypothesized that spurts result from increased turbulence in the earth's core dynamo during episodes of intense bolide bombardment of the earth. Mechanisms for creating episodes of intense bombardment of the earth involve gravitational perturbation of the Oort cloud of comets, either by a hidden planet, a solar companion, or massive matter in the galactic plane. Herein, the time variation in reversal rate is analyzed using methods of statistical density estimation. A smooth, continuous estimate of reversal rate is obtained using an adaptive kernel method, in which the kernel width is adjusted as a function of reversal rate. The estimates near the ends of the data series (at 165 my ago and the present) are obtained by extending the data by reflection. The results show that the reversal spurts are not associated demonstrably with extinctions or well-dated impacts. If the spurts do record episodes of intense bombardment of the earth, then the mass extinctions do not, in general, occur at times of impacts. Furthermore, the large impact craters seen are not obviously related to the spurts, suggesting that the craters may have been caused by bolides of a different nature and with a different temporal pattern. However, the most simple explanation seems to be that the spurts do not record comet showers, either because the recording mechanism suggested by Muller and Morris is not effective or because comet showers are not triggered in the ways considered by Hut et al.

Lutz, T. M.↗
NASA NTRS1993

Near-equilibrium growth of thick, high quality beta-SiC by sublimation

A close spaced near-equilibrium growth technique was used to produce thick, high quality epitaxial layers of beta-silicon carbide. The process utilized a sublimation method to grow morphologically smooth layers. The beta silicon carbide growth layers varied from about 200 to 750 microns in thickness. Chemical vapor deposition grown, 2-10 microns, beta silicon carbide films were used as seeds at 1860 and 1910 C growth temperatures. The respective average growth rates were 20 and 30 microns per hour. The layers are p-type with a 3.1 x 10 exp 17/cu cm carrier concentration. Electrical measurements indicate considerable improvement in the breakdown voltage of Schottky barriers on growth samples. Breakdown values ranged from 25 to 60 V. These measurements represent the highest values reported for 3C-SiC.

Shields, Virgil B.↗
NASA NTRS1993

SPH modelling of energy partitioning during impacts on Venus

Impact cratering of the Venusian planetary surface by meteorites was investigated numerically using the Smoothed Particle Hydrodynamics (SPH) method. Venus presently has a dense atmosphere. Vigorous transfer of energy between impacting meteorites, the planetary surface, and the atmosphere is expected during impact events. The investigation concentrated on the effects of the atmosphere on energy partitioning and the flow of ejecta and gas. The SPH method is particularly suitable for studying complex motion, especially because of its ability to be extended to three dimensions. In our simulations, particles representing impactors and targets are initially set to a uniform density, and those of atmosphere are set to be in hydrostatic equilibrium. Target, impactor, and atmosphere are represented by 9800, 80, and 4200 particles, respectively. A Tillotson equation of state for granite is assumed for the target and impactor, and an ideal gas with constant specific heat ratio is used for the atmosphere. Two dimensional axisymmetric geometry was assumed and normal impacts of 10km diameter projectiles with velocities of 5, 10, 20, and 40 km/s, both with and without an atmosphere present were modeled.

Takata, T.↗
NASA NTRS1994

Improved solution accuracy for Landsat-4 (TDRSS-user) orbit determination

This paper presents the results of a study to compare the orbit determination accuracy for a Tracking and Data Relay Satellite System (TDRSS) user spacecraft, Landsat-4, obtained using a Prototype Filter Smoother (PFS), with the accuracy of an established batch-least-squares system, the Goddard Trajectory Determination System (GTDS). The results of Landsat-4 orbit determination will provide useful experience for the Earth Observing System (EOS) series of satellites. The Landsat-4 ephemerides were estimated for the January 17-23, 1991, timeframe, during which intensive TDRSS tracking data for Landsat-4 were available. Independent assessments were made of the consistencies (overlap comparisons for the batch case and convariances for the sequential case) of solutions produced by the batch and sequential methods. The filtered and smoothed PFS orbit solutions were compared with the definitive GTDS orbit solutions for Landsat-4; the solution differences were generally less than 15 meters.

Oza, D. H.↗
NASA NTRS1994

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗
NASA NTRS1994

Hydrodynamics of binary coalescence. 1: Polytropes with stiff equations of state

We performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations on the verge of dynamical instability. We consider synchronized equilbria only and concentrate on stiff equations of state, with adiabatic Gamma greater than 5/3. We assume that the polytropic constants (K identically equal to P/(rho(exp Gamma)) are the same for both components. These conditions apply well to models of neutron star binaries. Accordingly, we discuss our results in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO) project, and we calculate the emission of gravitational radiation in the quadruple approximation. The fully nonlinear development of the instability is followed using SPH until a new equilibrium configuration is reached by the system. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An axisymmetric merged configuration is always produced when the adiabatic exponent Gamma approximately less than 2.3. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, triaxial merged configurations are obtained when Gamma approximately greater than 2.3, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios q not equal to 1 typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. For a mass ratio q approximately less than 0.5, however, dynamical mass transfer can temporarily retard the coalescence by causing a rapid reexpansion of the binary into a new, slightly eccentric but dynamically stable orbit. The maximum amplitude h(sub max) and peak luminosity L(sub max) of the gravitational waves emitted during the final coalescence are nearly independent of Gamma, but depend sensitively on the mass ratio q. The approximate scalings we find are h(sub max) varies as q(exp 2) and L(sub max) varies as q(exp 6) for q close to unity.

Rasio, Frederic A.↗
NASA NTRS1994

Impact on comet Shoemaker-Levy 9 on Jupiter

Three-dimensional numerical simulations of the impact of Comet Shoemaker - Levy 9 on Jupiter and the resulting vapor plume expansion were conducted using the Smoothed Particle Hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 km to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disk. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km, but no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Ahrens, Thomas J.↗
NASA NTRS1994

Comet Shoemaker-Levy 9: Impact on Jupiter and plume evolution

The impact of fragments of Comet Shoemaker-Levy 9 on Jupiter and the resulting vapor plume expansion are investigated by conducting three-dimensional numerical simulations using the smoothed particle hydrodynamics (SPH) method. An icy body, representing the cometary fragments, with a velocity of 60 km/sec and a diameter of 2 km can penetrate to 350 km below the 1-bar pressure level in the atmosphere. Most of the initial kinetic energy of the fragment is transferred to the atmosphere between 50 km and 300 km below the 1-bar pressure level. The shock-heated atmospheric gas in the wake is totally dissociated and partially ionized. Scaling our SPH results to other sizes indicates that fragments larger than approximately 100 m in diameter can penetrate to below the visible cloud decks. The energy deposited in the atmosphere is explosively released in the upward expansion of the resulting plume. The plume preferentially expands upward rather than horizontally due to the density gradient of the ambient atmosphere. It rises greater than or equal to 10(exp 2) km in approximately 10(exp 2) sec. Eventually the total atmospheric mass ejected to above 1 bar is greater than or equal to 40 times the initial mass of the impactor. The plume temperature at a radius approximately 10(exp 3) km is greater than 10(exp 3) K for 10(exp 3) sec for a 2-km fragment. We predict that impact-induced plumes will be observable with the remote sensing instruments of the Galileo spacecraft. As the impact site rotates into the view of Earth some 20 min after the impact, the plume expansion will be observable using the Hubble Space Telescope (HST) and from visible and infrared instruments on groundbased telescopes. The rising plume reaches approximately 3000 km altitude in approximately 10 min and will be visible from Earth.

Takata, Toshiko↗
NASA NTRS1994

Impact of comet Shoemaker-Levy 9 on Jupiter

We have employed three-dimensional numerical simulations of the impact of Comet Shoemaker-Levy 9 (SL9) on Jupiter and the resulting vapor plume expansion using the smoothed particle hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the atmosphere between -100 to -250 km. This energy is converted to potential energy of the resulting gas plume. The unconfined plume expands vertically and has a peak radiative power approximately equal to the total radiation from Jupiter's disc. The plume rises a few tens of atmospheric scale heights in approximately 10(exp 2) seconds. The rising plume reaches the altitude of approximately 3000 km; however, no atmospheric gas is accelerated to the escape velocity (approximately 60 km/s).

Takata, Toshiko↗
NASA NTRS2014

Characterizing the Background Corona with SDO/AIA

Characterizing the nature of the solar coronal background would enable scientists to more accurately determine plasma parameters, and may lead to a better understanding of the coronal heating problem. Because scientists study the 3D structure of the Sun in 2D, any line-of-sight includes both foreground and background material, and thus, the issue of background subtraction arises. By investigating the intensity values in and around an active region, using multiple wavelengths collected from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) over an eight-hour period, this project aims to characterize the background as smooth or structured. Different methods were employed to measure the true coronal background and create minimum intensity images. These were then investigated for the presence of structure. The background images created were found to contain long-lived structures, including coronal loops, that were still present in all of the wavelengths, 131, 171, 193, 211, and 335 A. The intensity profiles across the active region indicate that the background is much more structured than previously thought.

corona↗
NASA NTRS2014

Characterizing the True Background Corona with SDO/AIA

Characterizing the nature of the solar coronal background would enable scientists to more accurately determine plasma parameters, and may lead to a better understanding of the coronal heating problem. Because scientists study the 3D structure of the Sun in 2D, any line of sight includes both foreground and background material, and thus, the issue of background subtraction arises. By investigating the intensity values in and around an active region, using multiple wavelengths collected from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) over an eight-hour period, this project aims to characterize the background as smooth or structured. Different methods were employed to measure the true coronal background and create minimum intensity images. These were then investigated for the presence of structure. The background images created were found to contain long-lived structures, including coronal loops, that were still present in all of the wavelengths, 193 Angstroms,171 Angstroms,131 Angstroms, and 211 Angstroms. The intensity profiles across the active region indicate that the background is much more structured than previously thought.

Napier, Kate↗
NASA NTRS2003

An Analysis of Measured Pressure Signatures From Two Theory-Validation Low-Boom Models

Two wing/fuselage/nacelle/fin concepts were designed to check the validity and the applicability of sonic-boom minimization theory, sonic-boom analysis methods, and low-boom design methodology in use at the end of the 1980is. Models of these concepts were built, and the pressure signatures they generated were measured in the wind-tunnel. The results of these measurements lead to three conclusions: (1) the existing methods could adequately predict sonic-boom characteristics of wing/fuselage/fin(s) configurations if the equivalent area distributions of each component were smooth and continuous; (2) these methods needed revision so the engine-nacelle volume and the nacelle-wing interference lift disturbances could be accurately predicted; and (3) current nacelle-configuration integration methods had to be updated. With these changes in place, the existing sonic-boom analysis and minimization methods could be effectively applied to supersonic-cruise concepts for acceptable/tolerable sonic-boom overpressures during cruise.

Mack, Robert J.↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
DOE OSTI2021

Meshfree simulation and experimental validation of extreme thermomechanical conditions in friction stir extrusion

Abstract Friction stir extrusion (FSE) is a novel solid-phase processing technique that consolidates and extrudes metal powders, flakes, chips, or billets into high-performance parts by plastic deformation, which has the potential to save substantial processing time and energy. Currently, most studies on FSE are experimental and only a few numerical models have been developed to explain and predict the complex physics of the process. In this work, a meshfree simulation framework based on smoothed particle hydrodynamics (SPH) was developed for FSE. Unlike traditional grid-based methods, SPH is a Lagrangian particle-based method that can handle severe material deformations, capture moving interfaces and surfaces, and monitor the field variable histories explicitly without complicated tracking schemes. These aspects of SPH make it attractive for the FSE process, where in situ evolution of field variables is difficult to observe experimentally. To this end, a 3-D, fully thermomechanically coupled SPH model was developed to simulate the FSE of aluminum wires. The developed model was thoroughly validated by comparing the numerically predicted material flow, strain, temperature history, and extrusion force with experimental results for a certain set of process parameters. The validated SPH model can serve as an effective tool to predict and better understand the extreme thermomechanical conditions during the FSE process.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1993

Conservative high-order-accurate finite-difference methods for curvilinear grids

Two fourth-order-accurate finite-difference methods for numerically solving hyperbolic systems of conservation equations on smooth curvilinear grids are presented. The first method uses the differential form of the conservation equations; the second method uses the integral form of the conservation equations. Modifications to these schemes, which are required near boundaries to maintain overall high-order accuracy, are discussed. An analysis that demonstrates the stability of the modified schemes is also provided. Modifications to one of the schemes to make it total variation diminishing (TVD) are also discussed. Results that demonstrate the high-order accuracy of both schemes are included in the paper. In particular, a Ringleb-flow computation demonstrates the high-order accuracy and the stability of the boundary and near-boundary procedures. A second computation of supersonic flow over a cylinder demonstrates the shock-capturing capability of the TVD methodology. An important contribution of this paper is the dear demonstration that higher order accuracy leads to increased computational efficiency.

Rai, Man M.↗
NASA NTRS1986

Partial spline models for the inclusion of tropopause and frontal boundary information in otherwise smooth two- and three-dimensional objective analysis

A new method, based on partial spline models, is developed for including specified discontinuities in otherwise smooth two- and three-dimensional objective analyses. The method is appropriate for including tropopause height information in two- and three-dimensinal temperature analyses, using the O'Sullivan-Wahba physical variational method for analysis of satellite radiance data, and may in principle be used in a combined variational analysis of observed, forecast, and climate information. A numerical method for its implementation is described and a prototype two-dimensional analysis based on simulated radiosonde and tropopause height data is shown. The method may also be appropriate for other geophysical problems, such as modeling the ocean thermocline, fronts, discontinuities, etc.

Shiau, Jyh-Jen↗