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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2020

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗

NASA NTRS1990

Fuselage design for a specified Mach-sliced area distribution

A procedure for designing a fuselage having a prescribed effective area distribution computed from -90 deg Mach slices is described. This type of calculation is an essential tool in designing a complete configuration with an effective area distribution that corresponds to a desired sonic boom signature shape. Sample calculations are given for M=2 and M=3 designs. The examples include fuselages constrained to have circular cross sections and fuselages having cross sections of arbitrary shape. It is found that, for a prescribed effective area distribution having sharp variations, the iterative procedure converges to a smoothed approximation to the prescribed distribution. For a smooth prescribed area distribution, the solution is not unique.

Barger, Raymond L.↗

DOE OSTI2020

Stochastic DC optimal power flow with reserve saturation

We propose an optimization framework for stochastic optimal power flow with uncertain loads and renewable generator capacity. Our model follows previous work in assuming that generator outputs respond to load imbalances according to an affine control policy, but introduces a model of saturation of generator reserves by assuming that when a generator’s target level hits its limit, it abandons the affine policy and produces at that limit. This is a particularly interesting feature in models where wind power plants, which have uncertain upper generation limits, are scheduled to provide reserves to balance load fluctuations. Here, the resulting model is a nonsmooth nonconvex two-stage stochastic program, and we use a stochastic approximation method to find stationary solutions to a smooth approximation. Computational results on 6-bus and 118-bus test instances demonstrate that by considering the effects of saturation, our model can yield solutions with lower expected generation costs (at the same target line violation probability level) than those obtained from a model that enforces the affine policy to stay within generator limits with high probability.

17 WIND ENERGY↗

NASA NTRS1973

A method for digital image registration using a mathematical programming technique

A new algorithm based on a nonlinear programming technique to correct the geometrical distortions of one digital image with respect to another is discussed. This algorithm promises to be superior to existing ones in that it is capable of treating localized differential scaling, translational and rotational errors over the whole image plane. A series of piece-wise 'rubber-sheet' approximations are used, constrained in such a manner that a smooth approximation over the entire image can be obtained. The theoretical derivation is included. The result of using the algorithm to register four channel S065 Apollo IX digitized photography over Imperial Valley, California, is discussed in detail.

Yao, S. S.↗

DOE OSTI2020

Solving Eco-Driving Problems Using Indirect Collocation Method and Smooth Representation

Here, we discuss the eco-driving problem, considering both electric and conventional powertrains, and presents a pathway to solving it numerically using an indirect collocation method. Despite the low-order system dynamics, the piecewise fuel/efficiency map, gear shifting, and real-world traffic/road situations bring system discontinuities/switchings and pure state constraints into the problem formulation, which make the problem highly nonlinear and nontrivial to solve. This letter introduces smooth approximations to convert the original problem to an unconstrained (and penalized) smooth boundary-value problem. This approach eliminates the discussion of the switching structure and leads to a lightweight Newton-method-based solution procedure.

42 ENGINEERING↗

DOE OSTI2024

Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method

Abstract An interface-modified reproducing kernel particle method (IM-RKPM) is introduced in this work to allow for a direct model construction from image pixels of heterogeneous polycrystalline Li-ion battery microstructures. The interface-modified reproducing kernel (IM-RK) approximation is constructed through scaling of a kernel function by a regularized distance function in conjunction with strategic placement of interface node locations. This leads to RK shape functions with either weak or strong discontinuities across material interfaces, suitable for modeling various interface mechanics. With the placement of a triple junction node and distance-based scaling of kernel functions, the resulting IM-RK shape function also possesses proper discontinuities at the triple junctions. This IM-RK approximation effectively remedies the well-known Gibb’s oscillation in the smooth approximation of discontinuities. Different from the conventional meshfree approaches for interface discontinuities, this IM-RK approach is done without additional degrees of freedom associated with the enrichment functions, and it is formulated with the standard procedures in the RK shape function construction. This work focuses on identifying the accuracy and convergence properties of IM-RKPM for modeling the coupled electro-chemo-mechanical system. A linear patch test is formulated and numerically tested for the electro-chemo-mechanical coupled problem with a Butler–Volmer boundary condition representing the physical conditions in Li-ion battery microstructures. This is followed by verification of the optimal rates of convergence of IM-RKPM for solving the coupled problem with higher order solutions. The image-based modeling of Li-ion battery microstructures in the numerical examples demonstrates the applicability of the proposed method to realistic Li-ion battery materials modeling.

25 ENERGY STORAGE↗

NASA NTRS1951

Concerning the flow about ring-shaped cowlings Part II : annular bodies of infinite length with circulation for smooth entrance

The investigations carried out in a previous report (NACA TM 1325) concerning the flow about ring-shaped cowlings were extended by taking a circulation about the cowling into consideration. The present second report treats bodies of infinite length with approximately smooth entrance. The circulation was caused by distributing vortex rings of constant density over a stream surface extending to infinity. Furthermore, the influence of a hub body on such cowlings was dealt with. The examples treated are meant to give the designer a basis for his design.

FLOW, INCOMPRESSIBLE↗

NASA NTRS1985

Multigrid method for nearly singular and slightly indefinite problems

This paper deals with nearly singular, possibly indefinite problems for which the usual multigrid solvers converge very slowly or even diverge. The main difficulty is related to some badly approximated smooth functions which correspond to eigenfunctions with nearly zero eigenvalues. A correction to the usual coarse-grid equations is derived, both in the correction scheme and in the full approximation scheme. The performance of the new algorithm using this correction is essentially as that of usual multigrid for definite problems.

Brandt, A.↗

NASA NTRS1986

Multigrid method for nearly singular and slightly indefinite problems

Brandt, A.↗

NASA NTRS1995

Observations of the 3-D distribution of interplanetary electrons and ions from solar wind plasma to low energy cosmic rays

The 3-D Plasma and Energetic Particle instrument on the GGS Wind spacecraft (launched November 1, 1994) is designed to make measurements of the full three-dimensional distribution of suprathermal electrons and ions from solar wind plasma to low energy cosmic rays, with high sensitivity, wide dynamic range, good energy and angular resolution, and high time resolution. Three pairs of double-ended telescopes, each with two or three closely sandwiched passivated ion implanted silicon detectors measure electrons and ions from approximately 20 keV to greater than or equal to 300 keV. Four top-hat symmetrical spherical section electrostatic analyzers with microchannel plate detectors, a large and a small geometric factor analyzer for electrons and a similar pair for ions, cover from approximately 3 eV to 30 keV. We present preliminary observations of the electron and ion distributions in the absence of obvious solar impulsive events and upstream particles. The quiet time electron energy spectrum shows a smooth approximately power law fall-off extending from the halo population at a few hundred eV to well above approximately 100 keV The quiet time ion energy spectrum also shows significant fluxes over this energy range. Detailed 3-D distributions and their temporal variations will be presented.

Lin, R. P.↗

DOE OSTI2020

Adaptive Activation Functions Accelerate Convergence in Deep and Physics-informed Neural Networks

We employ adaptive activation functions for regression in deep and physics-informed neural networks (PINNs) to approximate smooth and discontinuous functions as well as solutions of linear and nonlinear partial differential equations. In particular, we solve the nonlinear Klein-Gordon equation, which has smooth solutions, the nonlinear Burgers equation, which can admit high gradient solutions, and the Helmholtz equation. We introduce a scalable hyper-parameter in the activation function, which can be optimized to achieve best performance of the network as it changes dynamically the topology of the loss function involved in the optimization process. The adaptive activation function has better learning capabilities than the traditional one (fixed activation) as it improves greatly the convergence rate, especially at early training, as well as the solution accuracy. To better understand the learning process, we plot the neural network solution in the frequency domain to examine how the network captures successively different frequency bands present in the solution. We consider both forward problems, where the approximate solutions are obtained, as well as inverse problems, where parameters involved in the governing equation are identified. Our simulation results show that the proposed method is a very simple and effective approach to increase the efficiency, robustness and accuracy of the neural network approximation of nonlinear functions as well as solutions of partial differential equations, especially for forward problems. We theoretically prove that in the proposed method, gradient descent algorithms are not attracted to suboptimal critical points or local minima.

machine leaning, Bad minima, Inverse problems, Phy↗

NASA NTRS1978

Wave propagation through a random medium - The random slab problem

The first-order smoothing approximation yields integral equations for the mean and the two-point correlation function of a wave in a random medium. A method is presented for the approximate solution of these equations that combines features of the eiconal approximation and of the Born expansion. This method is applied to the problem of reflection and transmission of a plane wave by a slab of a random medium. Both the mean wave and the covariance are calculated to determine the reflected and transmitted amplitudes and intensities.

Acquista, C.↗

NASA NTRS1987

Numerical studies of the thermal design sensitivity calculation for a reaction-diffusion system with discontinuous derivatives

The aim of this study is to find a reliable numerical algorithm to calculate thermal design sensitivities of a transient problem with discontinuous derivatives. The thermal system of interest is a transient heat conduction problem related to the curing process of a composite laminate. A logical function which can smoothly approximate the discontinuity is introduced to modify the system equation. Two commonly used methods, the adjoint variable method and the direct differentiation method, are then applied to find the design derivatives of the modified system. The comparisons of numerical results obtained by these two methods demonstrate that the direct differentiation method is a better choice to be used in calculating thermal design sensitivity.

Hou, Jean W.↗

NASA NTRS1991

Computer-Aided Design Of Sheet-Material Parts

Computer-aided-design system partly automates tedious process of designing and guiding assembly of small pieces of flat sheet material into large surfaces that approximate smoothly curved surfaces having complicated three-dimensional shapes. Capability provides for flexibility enabling designer to assess quickly and easily effects of changes in design in making engineering compromises among various sizes and shapes. Saves time and money in both design and fabrication. Used in rocket-engine application and other applications requiring design of sheet-material parts.

Gilbert, Jeffrey L.↗

NASA NTRS1987

Structure and dynamics of amorphous water ice

Further insight into the structure and dynamics of amorphous water ice, at low temperatures, was obtained by trapping in it Ar, Ne, H2, and D2. Ballistic water-vapor deposition results in the growth of smooth, approximately 1 x 0.2 micrometer2, ice needles. The amorphous ice seems to exist in at least two separate forms, at T < 85 K and at 85 < T < 136.8 K, and transform irreversibly from one form to the other through a series of temperature-dependent metastable states. The channels formed by the water hexagons in the ice are wide enough to allow the free penetration of H2 and D2 into the ice matrix even in the relatively compact cubic ice, resulting in H2-(D2-) to-ice ratios (by number) as high as 0.63. The larger Ar atoms can penetrate only into the wider channels of amorphous ice, and Ne is an intermediate case. Dynamic percolation behavior explains the emergence of Ar and Ne (but not H2 and D2) for the ice, upon warming, in small and big gas jets. The big jets, each containing approximately 5 x 10(10) atoms, break and propel the ice needles. Dynamic percolation also explains the collapse of the ice matrix under bombardment by Ar , at a pressure exceeding 2.6 dyn cm-2, and the burial of huge amounts of gas inside the collapsed matrix, up to an Ar-to-ice of 3.3 (by number). The experimental results could be relevant to comets, icy satellites, and icy grain mantles in dense interstellar clouds.

NASA Program Exobiology↗

NASA NTRS1993

Experiment-specific cosmic microwave background calculations made easier - Approximation formula for smoothed delta T/T windows

Simple and easy to implement elementary function approximations are introduced to the spectral window functions needed in calculations of model predictions of the cosmic microwave backgrond (CMB) anisotropy. These approximations allow the investigator to obtain model delta T/T predictions in terms of single integrals over the power spectrum of cosmological perturbations and to avoid the necessity of performing the additional integrations. The high accuracy of these approximations is demonstrated here for the CDM theory-based calculations of the expected delta T/T signal in several experiments searching for the CMB anisotropy.

Gorski, Krzysztof M.↗

NASA NTRS1993

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully non-linear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel-dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel-dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate than that from ZA or TZA, (b) the error in the phase angle of Fourier components is worse than that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗

NASA NTRS1994

A test of the adhesion approximation for gravitational clustering

We quantitatively compare a particle implementation of the adhesion approximation to fully nonlinear, numerical 'N-body' simulations. Our primary tool, cross-correlation of N-body simulations with the adhesion approximation, indicates good agreement, better than that found by the same test performed with the Zel'dovich approximation (hereafter ZA). However, the cross-correlation is not as good as that of the truncated Zel'dovich approximation (TZA), obtained by applying the Zel'dovich approximation after smoothing the initial density field with a Gaussian filter. We confirm that the adhesion approximation produces an excessively filamentary distribution. Relative to the N-body results, we also find that: (a) the power spectrum obtained from the adhesion approximation is more accurate that that from ZA to TZA, (b) the error in the phase angle of Fourier components is worse that that from TZA, and (c) the mass distribution function is more accurate than that from ZA or TZA. It appears that adhesion performs well statistically, but that TZA is more accurate dynamically, in the sense of moving mass to the right place.

Melott, Adrian L.↗