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 91 records · Page 5

A simulation of the collapse and fragmentation of cooling molecular clouds

The application of the Smoothed Particle Hydrodynamics method to the fragmentation of rotating cloud and disk systems is described, allowing for molecular cooling due to H2 and CO. A novel approach to solving Poisson's equation for disklike structures which exploits the multigrid algorithm is also described. Numerical studies are presented which investigate the evolution of both rotating clouds and Maclaurin disks, in each case with both an isothermal equation of state and with molecular cooling. The results establish the influence of molecular cooling on the fragmentation of molecular clouds. The isothermal sequences, if they fragment at all, do so into far fewer lumps than the cooling sequences. This is not due to a cooling instability as such, but rather to the reduced thermal support. One of the sequences shows a remarkable similarity to the W49A star-forming region.

Monaghan, Joe J.↗

The origin of the moon and the single impact hypothesis. IV

The smoothed-particle hydrodynamics method developed to simulate the conditions under which a major planetary collision may have led to the formation of the moon is presently used in 41 simulations exploring the relevant parameter space. Two possible scenarios emerge: one in which the moon is formed directly in the collision or only after dissipation and evolution of the disk. Attention is given to the Giant Blowoff hypothesis, in which not only the primitive atmosphere but large amounts of mantle material flow over the surface of the orbiting disk and blow off in the forward direction. Large amounts of angular momentum can also be lost.

Cameron, A. G. W.↗

Investigations on entropy layer along hypersonic hyperboloids using a defect boundary layer

A defect approach coupled with matched asymptotic expansions is used to derive a new set of boundary layer equations. This method ensures a smooth matching of the boundary layer with the inviscid solution. These equations are solved to calculate boundary layers over hypersonic blunt bodies involving the entropy gradient effect. Systematic comparisons are made for both axisymmetric and plane flows in several cases with different Mach and Reynolds numbers. After a brief survey of the entropy layer characteristics, the defect boundary layer results are compared with standard boundary layer and full Navier-Stokes solutions. The entropy gradient effects are found to be more important in the axisymmetric case than in the plane one. The wall temperature has a great influence on the results through the displacement effect. Good predictions can be obtained with the defect approach over a cold wall in the nose region, with a first order solution. However, the defect approach gives less accurate results far from the nose on axisymmetric bodies because of the thinning of the entropy layer.

Brazier, J. P.↗

An adaptive grid/Navier-Stokes methodology for the calculation of nozzle afterbody base flows with a supersonic freestream

This paper describes an adaptive grid method for base flows in a supersonic freestream. The method is based on the direct finite-difference statement of the equidistribution principle. The weighting factor is a combination of the Mach number, density, and velocity first-derivative gradients in the radial direction. Two key ideas of the method are to smooth the weighting factor by using a type of implicit smoothing and to allow boundary points to move in the grid adaptation process. An AGARD nozzle afterbody base flow configuration is used to demonstrate the performance of the adaptive grid methodology. Computed base pressures are compared to experimental data. The adapted grid solutions offer a dramatic improvement in base pressure prediction compared to solutions computed on a nonadapted grid. A total-variation-diminishing (TVD) Navier-Stokes scheme is used to solve the governing flow equations.

Williams, Morgan↗

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

A Survey and Synthesis of Solar Cycle Prediction Techniques

A number of techniques for predicting solar activity on a solar cycle time scale are identified, described, and tested with historical data. Some techniques, e.g. regression and curve-fitting, work well as solar activity approaches maximum and provide a complete description of future activity, while others, e.g. geomagnetic precursors, work well near solar minimum but only provide an estimate of the amplitude of the cycle. A synthesis of different techniques is shown to provide a more accurate and useful forecast of solar cycle activity levels. A combination of two uncorrelated geomagnetic precursor techniques provides the most accurate prediction for the amplitude of a solar activity cycle at a time well before activity minimum. This Combined Precursor Method gives a smoothed sunspot number maximum of 154+/-21 at the 95% level of confidence for the next cycle maximum. A mathematical function dependent upon the time of cycle initiation and the cycle amplitude then describes the level of solar activity for the next complete cycle. As the time of cycle maximum approaches a better estimate of the cycle activity is obtained by including the fit between previous activity levels and this function. This Combined Solar Cycle Activity Forecast gives, as of January 1999, a smoothed sunspot maximum of 146+/-20 at the 95% level of confidence for the next cycle maximum. The success of the geomagnetic precursors in predicting future solar activity suggests that solar magnetic phenomena at latitudes above the sunspot activity belts are linked to solar activity which occurs many years later in the lower latitudes.

Hathaway, David H.↗

A Synthesis of Solar Cycle Prediction Techniques

A number of techniques currently in use for predicting solar activity on a solar cycle timescale are tested with historical data. Some techniques, e.g., regression and curve fitting, work well as solar activity approaches maximum and provide a month-by-month description of future activity, while others, e.g., geomagnetic precursors, work well near solar minimum but only provide an estimate of the amplitude of the cycle. A synthesis of different techniques is shown to provide a more accurate and useful forecast of solar cycle activity levels. A combination of two uncorrelated geomagnetic precursor techniques provides a more accurate prediction for the amplitude of a solar activity cycle at a time well before activity minimum. This combined precursor method gives a smoothed sunspot number maximum of 154 plus or minus 21 at the 95% level of confidence for the next cycle maximum. A mathematical function dependent on the time of cycle initiation and the cycle amplitude is used to describe the level of solar activity month by month for the next cycle. As the time of cycle maximum approaches a better estimate of the cycle activity is obtained by including the fit between previous activity levels and this function. This Combined Solar Cycle Activity Forecast gives, as of January 1999, a smoothed sunspot maximum of 146 plus or minus 20 at the 95% level of confidence for the next cycle maximum.

Hathaway, David H.↗

Status of Cycle 23 Forecasts

A number of techniques for predicting solar activity on a solar cycle time scale are identified, described, and tested with historical data. Some techniques, e.g,, regression and curve-fitting, work well as solar activity approaches maximum and provide a month- by-month description of future activity, while others, e.g., geomagnetic precursors, work well near solar minimum but provide an estimate only of the amplitude of the cycle. A synthesis of different techniques is shown to provide a more accurate and useful forecast of solar cycle activity levels. A combination of two uncorrelated geomagnetic precursor techniques provides the most accurate prediction for the amplitude of a solar activity cycle at a time well before activity minimum. This precursor method gave a smoothed sunspot number maximum of 154+21 for cycle 23. A mathematical function dependent upon the time of cycle initiation and the cycle amplitude then describes the level of solar activity for the complete cycle. As the time of cycle maximum approaches a better estimate of the cycle activity is obtained by including the fit between recent activity levels and this function. This Combined Solar Cycle Activity Forecast now gives a smoothed sunspot maximum of 140+20 for cycle 23. The success of the geomagnetic precursors in predicting future solar activity suggests that solar magnetic phenomena at latitudes above the sunspot activity belts are linked to solar activity, which occurs many years later in the lower latitudes.

Hathaway, D. H.↗

Modeling the Compression of Merged Compact Toroids by Multiple Plasma Jets

A fusion propulsion scheme has been proposed that makes use of the merging of a spherical distribution of plasma jets to dynamically form a gaseous liner. The gaseous liner is used to implode a magnetized target to produce the fusion reaction in a standoff manner. In this paper, the merging of the plasma jets to form the gaseous liner is investigated numerically. The Los Alamos SPHINX code, based on the smoothed particle hydrodynamics method is used to model the interaction of the jets. 2-D and 3-D simulations have been performed to study the characteristics of the resulting flow when these jets collide. The results show that the jets merge to form a plasma liner that converge radially which may be used to compress the central plasma to fusion conditions. Details of the computational model and the SPH numerical methods will be presented together with the numerical results.

Thio, Y. C. Francis↗

A Spatially Coherent Global Soil Moisture Product with Improved Temporal Resolution

Global soil moisture products that are completely independent of any type of ancillary data and solely rely on satellite observations are presented. Additionally, we further develop an existing downscaling technique that enhances the spatial resolution of such products to approximately 11 km. These products are based on internal modules of the Land Parameter Retrieval Model (LPRM), an algorithm that uses the radiative transfer equation to link soil moisture, vegetation optical depth and land surface temperature to observed brightness temperatures. The soil moisture product that is independent of any type of ancillary data uses the internally calculated dielectric constant as a soil moisture proxy. This data product is not influenced by errors associated with coarse-scale global soil property maps or by any other type of forcing (e.g. re-analysis) data and is therefore solely based on satellite microwave observations. The second step builds upon recent developments to increase the spatial resolution of the LPRM retrievals using a smoothing filter downscaling method. With this method we can attain a spatial resolution that can be more useful at the scale of local and regional hydrological studies as well. The steps presented in this paper were applied to observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). The newly derived data sets were validated using ground-based observations from the International Soil Moisture Network (ISMN). The internally calculated dielectric constant product results in significantly more days with valid retrievals than the original soil moisture data products, in particular over arid regions. The dielectric constant product resulted in similar correlations with in situ data as the original soil moisture data product. Together, these findings demonstrate the usefulness of this new dielectric constant product for the hydrological modeling community and climate studies. A case study on the Australian Fitzroy catchment demonstrated that the downscaled data product has a more detailed spatial description of soil moisture, especially during wet and dry conditions with more pronounced dry and wet regions within the catchment. The increased resolution data products could therefore improve runoff predictions and this study demonstrated the potential added value of a transitioning from a spatial resolution of 56 km toward a higher resolution of 11 km. The hydrological implications of these newly developed data records are not only linked to AMSR-E satellite data, but also to the next generation Soil Moisture Active and Passive (SMAP) mission where a 9 km spatial resolution is the target resolution for satellite soil moisture products. The new data products will not replace the current LPRM products, but will be added to the existing array of data products and will become publicly available through our data portals.

spatial↗

Airfoil self-noise and prediction

A prediction method is developed for the self-generated noise of an airfoil blade encountering smooth flow. The prediction methods for the individual self-noise mechanisms are semiempirical and are based on previous theoretical studies and data obtained from tests of two- and three-dimensional airfoil blade sections. The self-noise mechanisms are due to specific boundary-layer phenomena, that is, the boundary-layer turbulence passing the trailing edge, separated-boundary-layer and stalled flow over an airfoil, vortex shedding due to laminar boundary layer instabilities, vortex shedding from blunt trailing edges, and the turbulent vortex flow existing near the tip of lifting blades. The predictions are compared successfully with published data from three self-noise studies of different airfoil shapes. An application of the prediction method is reported for a large scale-model helicopter rotor, and the predictions compared well with experimental broadband noise measurements. A computer code of the method is given.

Brooks, Thomas F.↗

Smooth indentation of a composite laminate

Finite difference method is used to analyze the problem of smooth contact between a rigid indenter and a laminated circular plate clamped at the edges. The plate consists of transversely isotropic layers. The interlaminar shear stresses in the contact region are found to be much different from the laminate plate theory solutions. It has been found that low-modulus interlayers cause significant reduction in the maximum interlaminar shear stresses.

Sankar, B. V.↗

Generating Smooth Motions For Robotic Manipulators

In improved method for generating trajectory of robotic manipulator, each straight-line segment of trajectory composed of constant-velocity main portion sandwiched between smooth acceleration at start and smooth deceleration at finish. Algorithm implementing method computes velocity in each accelerating portion as sinusoidal function of position along line. This motion chosen for two reasons: closely approximates motion of human hand along straight-line trajectory, and provides very smooth transitions between constant-velocity portion and accelerated and decelerational end portions.

Bejczy, Antal K.↗

Effect of nose shape on three-dimensional stagnation region streamlines and heating rates

A new method for calculating the three-dimensional inviscid surface streamlines and streamline metrics using Cartesian coordinates and time as the independent variable of integration has been developed. The technique calculates the streamline from a specified point on the body to a point near the stagnation point by using a prescribed pressure distribution in the Euler equations. The differential equations, which are singular at the stagnation point, are of the two point boundary value problem type. Laminar heating rates are calculated using the axisymmetric analog concept for three-dimensional boundary layers and approximate solutions to the axisymmetric boundary layer equations. Results for elliptic conic forebody geometries show that location of the point of maximum heating depends on the type of conic in the plane of symmetry and the angle of attack, and that this location is in general different from the stagnation point. The new method was found to give smooth predictions of heat transfer in the nose region where previous methods gave oscillatory results.

Hassan, Basil↗

Formation of on-site normal points

We describe our methods of fitting a smoothing function to observational range differences from a predicted orbit, by deducing corrections to the orbit in the radial and along-track directions. The method has been used on observations of a variety of satellites, and using predicted orbits computed both by numerical integration using IRV's as starting values and analytically from orbital elements. The along-track corrections to the predicted orbit have been successfully used in the form of time biases to improve subsequent predictions, and a statistical test has been devised to ensure that the range residuals may be used to form unbiased quick look normal points.

Appleby, G. M.↗

Estimation of divergence and vorticity using multidimensional smoothing splines

Laplacian smoothing splines, smoothing splines on the sphere, and smoothing pseudo splines on the sphere are presented. The method of generalized cross validation to choose the smoothing parameter is described. These methods are applied to estimate divergence and vorticity of the atmosphere from wind speed and wind direction.

Wendelberger, J. G.↗

Multilevel elliptic smoothing of large three-dimensional grids

Elliptic grid generation methods have been used for many years to smooth and improve grids generated by algebraic interpolation schemes. However, the elliptic system that must be solved is nonlinear and convergence is generally very slow for large grids. In an attempt to make elliptic methods practical for large three-dimensional grids, a two-stage implementation is developed where the overall grid point locations are set using a coarse grid generated by the elliptic system. The coarse grid is then interpolated to generate a finer grid which is smoothed using only a few iterations of the elliptic system.

Mastin, C. Wayne↗

GEOS-2 C-band radar system project. Spectral analysis as related to C-band radar data analysis

Work performed on spectral analysis of data from the C-band radars tracking GEOS-2 and on the development of a data compaction method for the GEOS-2 C-band radar data is described. The purposes of the spectral analysis study were to determine the optimum data recording and sampling rates for C-band radar data and to determine the optimum method of filtering and smoothing the data. The optimum data recording and sampling rate is defined as the rate which includes an optimum compromise between serial correlation and the effects of frequency folding. The goal in development of a data compaction method was to reduce to a minimum the amount of data stored, while maintaining all of the statistical information content of the non-compacted data. A digital computer program for computing estimates of the power spectral density function of sampled data was used to perform the spectral analysis study.

Source record↗