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 289 records · Page 16

On the prediction of multigrid efficiency through local mode analysis

A single grid local mode analysis is used to predict the smoothing properties of numerical schemes for solving the Navier-Stokes equations with factorization based on Stone's Strongly Implicit Method. Four difference approximations for the convection terms are considered, namely, hybrid, central, second-order upwind, and third-order upwind. Smoothing factors from the analysis are compared with practical convergence factors in a multigrid method for flow over a backward facing step and it is found that the local mode analysis correctly predicts the effects of Reynolds number and higher-order schemes.

Wilson, R. V.↗

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

High-Order Spectral Volume Method for 2D Euler Equations

The Spectral Volume (SV) method is extended to the 2D Euler equations. The focus of this paper is to study the performance of the SV method on multidimensional non-linear systems. Implementation details including total variation diminishing (TVD) and total variation bounded (TVB) limiters are presented. Solutions with both smooth features and discontinuities are utilized to demonstrate the overall capability of the SV method.

Wang, Z. J.↗

Optimization of bonded joints

A procedure for minimizing the elastic shear stress concentration in adhesive lap joints is presented. The proposed method is based upon tapering the adherends to achieve smooth stiffness transitions and uniform shear stresses. Both single and double lap splices are considered, but numerical examples are restricted to the case of double lap joints. Nonisotropic materials and nonoptimum design limitations, such as minimum and maximum thickness adherends, load-line eccentricity, and peel stresses are treated, and typical results are presented.

Ojalvo, I. U.↗

A dynamo theory prediction for solar cycle 22 - Sunspot number, radio flux, exospheric temperature, and total density at 400 km

Using the 'dynamo theory' method to predict solar activity, a value for the smoothed sunspot number of 109 + or - 20 is obtained for solar cycle 22. The predicted cycle is expected to peak near December, 1990 + or - 1 year. Concommitantly, F(10.7) radio flux is expected to reach a smoothed value of 158 + or - 18 flux units. Global mean exospheric temperature is expected to reach 1060 + or - 50 K and global total average total thermospheric density at 400 km is expected to reach 4.3 x 10 to the -15th gm/cu cm + or - 25 percent.

Schatten, K. H.↗

Semianalytic modeling of aerodynamic shapes

Equations for the semianalytic representation of a class of surfaces that vary smoothly in cross-sectional shape are presented. Some methods of fitting together and superimposing such surfaces are described. A brief discussion is also included of the application of the theory in various contexts such as computerized lofting of aerodynamic surfaces and grid generation.

Barger, R. L.↗

A dynamo theory prediction for solar cycle 22: Sunspot number, radio flux, exospheric temperature, and total density at 400 km

Using the dynamo theory method to predict solar activity, a value for the smoothed sunspot number of 109 + or - 20 is obtained for solar cycle 22. The predicted cycle is expected to peak near December, 1990 + or - 1 year. Concommitantly, F(10.7) radio flux is expected to reach a smoothed value of 158 + or - 18 flux units. Global mean exospheric temperature is expected to reach 1060 + or - 50 K and global total average total thermospheric density at 400 km is expected to reach 4.3 x 10 to the -15th gm/cu cm + or - 25 percent.

Schatten, K. H.↗

An explicit Runge-Kutta method for 3D turbulent incompressible flows

A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.

Sung, Chao-Ho↗

High frequency scattering from corrugated stratified cylinders

Interest in applying radar remote sensing for the study of forested areas led to the development of a model for scattering from corrugated stratified dielectric cylinders. The model is used to investigate the effect of bark and its roughness on scattering from tree trunks and branches. The outer layer of the cylinder (bark) is assumed to be a low-loss dielectric material and to have a regular (periodic) corrugation pattern. The inner layers are treated as lossy dielectrics with smooth boundaries. A hybrid solution based on the moment method and the physical optics approximation is obtained. In the solution, the corrugations are replaced with polarization currents that are identical to those of the local tangential periodic corrugated surface, and the stratified cylinder is replaced with equivalent surface currents. New expressions for the equivalent physical-optics currents are used which are more convenient than the standard ones. It is shown that the bark layer and its roughness both reduce the radar cross-section. It is also demonstrated that the corrugations can be replaced by an equivalent anisotropic layer.

Sarabandi, Kamal↗

Computation of near-minimum-time maneuvers of flexible structures by parameter optimization

Near-minimum-time attitude maneuvers of space structures as well as ground based test articles are considered. Switching nature of the controls for rigid body maneuvers are illustrated using a control-cube and a critical control axis of rotation. The presence of torque smoothing and where appropriate, gravitational effects and connections to other bodies are explicitly included in the mathematical models of the systems to be optimized. A maximum fuel consumption constraint is included besides the required terminal conditions on attitude and angular velocities. The switch times, maximum thrust magnitudes, and smoothing parameters are determined using the Sequential Quadratic Programming method for parameter optimization. Results indicating attitude and angular velocity histories, thruster forces, and structural vibrations are presented for three, four, and five switch maneuvers, as well as maneuvers that involve large coasting arcs.

Vadali, S. R.↗

Regional Deformation Studies with GRACE and GPS

GRACE data indicate large seasonal variations in gravity that have been shown to be to be related to climate-driven fluxes of surface water. Seasonal redistribution of surface mass deforms the Earth, and our previous study using GRACE data demonstrate that annual radial deformations of +/-13 mm in the region of Amazon River Basin were observed by both GRACE and ten GPS sites in the region. For the GRACE determinations, we estimate in a least-squares solution for each Stokes coefficient parameters that represent the amplitudes of the annual variation. We then filter these parameters based on a statistical test that uses the scatter of the postfit residuals. We demonstrate by comparison to the GPS amplitudes that this method is more accurate, for this region, than Gaussian smoothing. Our model for the temporal behavior of the gravity coefficients includes a rate term, and although the time series are noisy, the glacial isostatic adjustment signal over Hudson s Bay can be observed. .

Davis, J. L.↗

STXM-C-XANES, IR Microspectroscopy, and Raman Microspectroscopy Analysis of Organic Matter in the Chwichiya 002 C3.00-Ungrouped Chondrite

Chwichiya 002 is a C3.00-ungrouped carbonaceous chondrite discovered in Western Sahara in 2018. The presence of hydrated phases in Chwichiya 002 was not confirmed by XRD analysis, and only a faint water and Si-O band attributed to phyllosilicate were observed by infrared transmission spectroscopy, indicating a very low degree of aqueous alteration (Meteoritical Bulletin, 2020). The structural order of polycyclic aromatic hydrocarbons evaluated by Raman spectroscopy indicates that the effect of thermal metamorphism is less severe than in the Semarkona meteorite, an LL3.00 ordinary chondrite (Meteoritical Bulletin, 2020). Chiwichiya 002 is a very primitive carbonaceous chondrite and an important information source for the early stage of the evolution of organic matter in the solar system. In this study, we investigated the molecular structure and distribution of organic matter in Chiwichiya 002 by Carbon X-ray Absorption Near Edge Structure (C-XANES) analysis using a Scanning Transmission X-ray Microscope (STXM). The STXM-C-XANES analysis has a spatial resolution of 20~100 nm and can observe the spatial distribution of organic materials in microstructures. In addition, to investigate the degree of thermal metamorphism of the organic matter, micro-IR, and micro-Raman spectroscopic analyses were performed. In addition, the Murchison CM2 chondrite was also analyzed for comparison. To evaluate contamination during analytical procedures, antigorite, which was heated in air at 500°C for 5 hours to eliminate organic matter was also analyzed in the same manner. The samples were milled into 100 nm-thick, thin sections using a Focused Ion Beam (FIB), and C-XANES spectra were obtained using the STXM at the High Energy Accelerator Research Organization, Photon Factory (PF), BL-19B. The obtained C-XANES spectra were smoothed and then fitted with Gaussian functions using the method of Le Guillou et al. (2018). The C-XANES showed that the ratio of oxygenated functional group (ketone, phenolic) carbons in Chwichiya 002 was lower than in Murchison, whereas the ratio of aliphatic carbons was almost the same. This may reflect the low degree of aqueous alteration of Chwichiya 002. In addition, contamination peaks were detected at 287.3-287.7 eV and 288.3-288.7 eV in the C-XANES spectra of some of the antigorite samples. The source of these peaks is under investigation. Raman spectroscopic analysis showed that the intensity ratio of the D band to the G band of Chwichiya 002 was slightly higher than that of Murchison. IR absorption spectra showed aromatic absorption at 1630 cm -1 , but there was little absorption due to aliphatic C-H around 2900 cm -1 . These results suggest that the organic matter in Chwichiya 002 may be more heated than in Murchison.

Chwichiya 002↗

High-Order/Low-Dissipation Chain-Rule Flux Solution Reconstruction Schemes in FUN3D

In this paper, we report progress in the development of economically high-order flux-solution-reconstruction (FSR) schemes, which are second-order accurate on general unstructured grids but achieve high-order accuracy when a grid is regular, i.e., has the same stencil (with the same spacing) throughout the domain. Two variants of the FSR schemes are discussed: chain-rule-flux-solution reconstruction (CFSR) and quadratic-form-flux-solution reconstruction (QFSR), where the former is based on the chain rule and the latter on the flux reconstruction expressed as a function of solution variables. These schemes are tested for flows with shock waves with a limiter incorporated in the flux and solution reconstructions. Improved results, compared with second-order methods are demonstrated for inviscid and viscous flows with smooth grids.

high-order↗

Full configuration-interaction study of the ionic-neutral curve crossing in LiF

Full configuration-interaction (FCI) calculations are used to assess the relative ability of methods for truncating the n-particle expansion in describing the ionic-neutral curve crossing between the two lowest Sigma(+) states of LiF. While the FCI calculations yield a smooth dipole moment function, MRCI calculations based on CASSCF orbitals optimized for the lowest state at all r values yield a discontinuous dipole moment function. However, when the orbitals are optimized using a state-averaged CASSCF procedure, with equal weights for the ionic and neutral solutions, both the CASSCF and MRCI dipole moment functions are smooth and in reasonable agreement with the FCI. No single-reference-based method is found to work satisfactorily. Potential curves for the lowest two Sigma(+) states are determined in both the adiabatic and diabatic representations using a large atomic natural orbit Gaussian basis set and a state-averaged CASSCF/MRCI treatment of electron correlation.

Bauschlicher, Charles W., Jr.↗

Adjustment, interpolation, and smoothing of GEOS 3 altimeter data

A computational scheme that provides a rapid method for constructing a local or regional geoid from GEOS 3 altimeter data is presented. It is shown that for an oceanic area not larger than 30 deg longitude by 30 deg latitude, an adjustment of GEOS 3 track data for the erroneous bias and tilt can reduce the data inconsistencies at track intersections to smaller than 50 cm. This makes it possible to construct a map of the geoid with the contour interval of 1 m in any local or regional areas where the altimeter data exist with reasonable density and uniformity. In addition, a formula is developed for computing an expected value of geoid height by taking a weighted average of the adjusted altimeter data. Using the data's autocovariance as a weight, the formula can give a profile of geoid height along a given traverse. The geoid data can thus be compared with other marine geophysical data that are recorded continually along a ship's track.

Horai, K.-I.↗

A Smoothed Eclipse Model for Solar Electric Propulsion Trajectory Optimization

Solar electric propulsion (SEP) is the dominant design option for employing low-thrust propulsion on a space mission. Spacecraft solar arrays power the SEP system but are subject to blackout periods during solar eclipse conditions. Discontinuity in power available to the spacecraft must be accounted for in trajectory optimization, but gradient-based methods require a differentiable power model. This work presents a power model that smooths the eclipse transition from total eclipse to total sunlight with a logistic function. Example trajectories are computed with differential dynamic programming, a second-order gradient-based method.

trajectory design↗

A Smoothed Eclipse Model for Solar Electric Propulsion Trajectory Optimization

Solar electric propulsion (SEP) is the dominant design option for employing low-thrust propulsion on a space mission. Spacecraft solar arrays power the SEP system but are subject to blackout periods during solar eclipse conditions. Discontinuity in power available to the spacecraft must be accounted for in trajectory optimization, but gradient-based methods require a differentiable power model. This work presents a power model that smooths the eclipse transition from total eclipse to total sunlight with a logistic function. Example trajectories are computed with differential dynamic programming, a second-order gradient-based method.

optimization↗