Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 541 records · Page 30

A half plane and a strip with an arbitrarily located crack

A technique is presented for dealing with the problem of an elastic domain containing an arbitrarily oriented internal crack. The problem is formulated as a system of integral equations for a fictitious layer of body forces imbedded in the plane along a closed smooth curve encircling the original domain. The problems of a half plane with a crack in the neighborhood of its free boundary and of an infinite strip containing a symmetrically located internal crack with an arbitrary orientation are considered as examples. In each case the stress intensity factors are computed and are given as functions of the crack angle.

Erdogan, F.↗

Proton dose approximation in arbitrary convex geometry

An expansion is derived for the solution to the transport equation in two dimensions subject to boundary conditions given for an arbitrary convex region. Questions of high-energy transport are considered along with the properties of the dose response function. The expansion of the solution of the transport equation is presented in terms of a parameter which measures the lateral dispersion of an unidirectional beam. This parameter is usually small and the expansion is expected to converge rapidly. The dominant term in the expansion is related to fluence-to-dose conversion factors in a semiinfinite slab for normal incidence. A convenient parameterization of the conversion factors is provided along with numerical examples.

Wilson, J. W.↗

A half plane and a strip with an arbitrarily located crack

The paper introduces a technique to deal with the problem of an elastic domain containing an arbitrarily oriented internal crack. The problem is formulated as a system of integral equations for a fictitious layer of body forces imbedded in the plane along a closed smooth curve encircling the original domain. The problems of a half plane with a crack in the neighborhood of its free boundary and of an infinite strip containing a symmetrically located internal crack with an arbitrary orientation are considered as examples. In each case the stress intensity factors are computed and are given as functions of the crack angle.

Erdogan, F.↗

A new determination of the halo luminosity density of the galaxy

The distribution of tangential velocities for arbitrary halo velocity dispersions and rotation rates are explicitly computed, assuming a particular form of the halo velocity distribution function near the sun. It is found that for a specific model of the halo velocity dispersion tensor that Schmidt's (1975) estimate of the mass carried by halo stars of 0.25-0.70 solar masses should be revised by a factor of between 1.0 and 2.0, depending on the rotation rate of the halo. Although obtaining a mass density in all stars from the stellar density in that mass range is extremely uncertain, it is found, in agreement with earlier work, that even under rather extreme assumptions about the remnants the observed stars and expected remnants fail to provide the theoretically expected mass. Comparison of the luminosity to the theoretical mass yields a mass to light ratio consistent with that of rich clusters of galaxies.

Richstone, D. O.↗

An accurate method for evaluating the kernel of the integral equation relating lift to downwash in unsteady potential flow

The method is capable of generating approximations of arbitrary accuracy. It is based on approximating the algebraic part of the nonelementary integrals in the kernel by exponential functions and then integrating termwise. The exponent spacing in the approximation is a geometric sequence. The coefficients and exponent multiplier of the exponential approximation are computed by least squares so the method is completely automated. Exponential approximates generated in this manner are two orders of magnitude more accurate than the exponential approximation that is currently most often used for this purpose. The method can be used to generate approximations to attain any desired trade-off between accuracy and computing cost.

Desmarais, R. N.↗

NLTE model stellar atmospheres with line blanketing near the series limits

In this paper we study the influence of line-merging regions at the intermediate long-wavelength side of a continuum threshold on the computed model atmosphere structure and predicted spectrum. In order to model these regions sufficiently accurately, we have developed two concepts. First, we have extended the occupation probability formalism of Hummer and Mihalas to non-local thermodynamic equilibrium (non-LTE) plasmas. Second, in order to treat the very complicated opacity in the line merging region, we have generalized the concept of opacity distribution functions to treat non-LTE situations. All Rydberg states are consistently included within this framework, so that no arbitrary cutoff of high (LTE) levels is made. We have calculated several pure hydrogen models atmospheres for two effective temperatures, T(sub eff) = 20000 and 35000 K, and discussed the differences between models calculated with various treatments of the line merging. In particular, we have shown that the error in the predicted profiles of Balmer lines resulting from the neglect of line merging is typically of the order of 3-4%, while the errors in the far-UV portion of the Balmer continuum reaches 15-35%. The errors generally decrease with increasing effective temperature. At the same time, the internal accuracy of the models is shown to be about or below 0.5% for all predicted spectral features. We conclude that for interpreting current high-accuracy spectrophotometric observations models including the line merging are necessary, and that the formalism developed in this paper is capable of providing a sufficiently accurate and robust modeling technique.

Hubeny, I.↗

An Electrothermodynamic Model with Distributed Properties for Effective Permittivities of Sea Ice

This paper presents a model to calculate the temperature dependence of effective permittivities for sea ice, a heterogeneous medium containing multiphase scatterers. With the strong permittivity fluctuation approach the model accounts for the electrodynamic scattering effect together with the quasi-static characteristics of multiple species and subspecies of inhomogeneities with distributed orientations, sizes, and shapes. Because of a preferential direction in the orientation distribution, the medium is effectively anisotropic. The size distribution is described with a probability density function in terms of normalized volumetric sizes. Scatterer shapes are nonuniform and have a general ellipsoidal form characterized by arbitrary axial ratios of correlation lengths which are related to physical geometries of the scatterers. In this formulation, sea ice consisting of solid ice, liquid brine, and gaseous inclusions is modeled to derive effective permittivities with thermodynamic phase redistribution and structural metamorphism. Theoretical results are in good agreement with experimental data at the C band frequency of 4.8 GHz for saline ice undergoing warming and cooling cycles. A competitive effect between the increase of liquid brine and the shape rounding of ellipsoidal scatterers at increasing temperatures explains the trend observed in measured data. Sensitivities of effective permittivities to structural and physical parameters characterizing sea ice are also studied.

Nghiem, S. V.↗

The selection of approximating functions for tabulated numerical data

A computer program was developed that selects, from a list of candidate functions, the approximating functions and associated coefficients which result in the best curve fit of a given set of numerical data. The advantages of the approach used here are: (1) Multivariable approximations can be performed. (2) Flexibility with respect to the type of approximations used is available. (3) The program is designed to choose the best terms to be used in the approximation from an arbitrary list of possible terms so that little knowledge of the proper approximating form is required. (4) Recursion relations are used in determining the coefficients of the approximating functions, which reduces the computer execution time of the program.

Ingram, H. L.↗

Triangle Method for Dense ReLU Layers [SWR-25-72]

This software is an implementation of the methods for initializing and training neural networks to be more efficient per parameter, described more fully below and in the related publication: In theory, depth should make a ReLU network EXPONENTIALLY more efficient by enabling it to produce an exponential number of piecewise linear sections in its output. This reasoning is largely based on the work of mathematicians that have hand-constructed networks that make good use of depth. In practice however, even very deep ReLU networks that have been randomly initialized will behave identically to their shallow counterparts - missing an entire exponential dimension of efficiency. The triangle method is a first attempt at realizing the exponential potential of deep networks. Instead of randomly setting weights, we force pairs of neurons in each layer learn to build triangles (i.e. functions from [0,1] -> [0,1] that look like triangles). This is a very efficient pattern for generating lots of linear pieces because composing two triangular functions doubles the number of pieces with each composition. The triangle method is more than just a different initialization, it is a new paradigm of training. Instead of making direct updates to the matrix weights, we do an extra step of backpropagation to collect the derivatives of the loss function with respect to the shapes of the triangles, training them to tilt left or right. This process essentially holds the networks hand throughout the loss landscape and forces it to always use depth effectively by producing triangular shapes internally. This can produce several orders of magnitude of improvement on convex one-dimensional regression problems. Much more theoretical work is needed to realize its full potential beyond this context, but the implementation in this repository will still work in arbitrary numbers of dimensions. The file Triangle_Method.py is a generalized form of the method that will build each neuron its own custom 1-d convex activation function (with exponential efficiency). Example usage on one dimensional problems can be found in Example_Usage.ipynb and an example of using this in a real neural network can be found in Example_VGG16_CIFAR10.ipynb.

Milkert, Max [National Renewable Energy Laboratory↗

Rotational symmetry relation for efficient response function generation in the coarse mesh transport method COMET

The coarse mesh transport code COMET is a continuous energy hybrid stochastic-deterministic neutronics solver with high fidelity and formidable computation speed in solving reactor core problems. Its method is based on the incident flux expansion theory. In this work, we take advantage of the local geometric symmetry in many reactor cores lattices (e.g., fuel lattices and reflector blocks) to develop relations among the flux response expansion coefficients for symmetric surfaces to further improve the computational efficiency of the COMET response function generation tool (method). This is done by a rigorous derivation of the transformation matrices for the angular and spatial expansion moments resulting from a rotation of a coarse mesh by an arbitrary angle. The relations for the response coefficients for the symmetric surfaces can be then written as the Kronecker product of those transformation matrices. The method is implemented into COMET and tested on two advanced high temperature reactor (AHTR) full-length single assembly benchmark problems. The COMET results using the response function library based on the symmetry relations were compared to those using the library directly generated by continuous energy Monte Carlo for all surfaces. It was found that the eigenvalues and stripe-wise fission densities using the two libraries are in statistical agreement as expected. This indicates that the new method maintains the high fidelity of the original COMET method while improving the computational efficiency in the response function generation by 270% to 400%, depending on the local geometric symmetry. This method also reduces the size of the response function library by the same magnitude (270% to 400%). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Porous cage-derived nanomaterial inks for direct and internal three-dimensional printing

The convergence of 3D printing techniques and nanomaterials is generating a compelling opportunity space to create advanced materials with multiscale structural control and hierarchical functionalities. While most nanoparticles consist of a dense material, less attention has been payed to 3D printing of nanoparticles with intrinsic porosity. Here, we combine ultrasmall (about 10 nm) silica nanocages with digital light processing technique for the direct 3D printing of hierarchically porous parts with arbitrary shapes, as well as tunable internal structures and high surface area. Thanks to the versatile and orthogonal cage surface modifications, we show how this approach can be applied for the implementation and positioning of functionalities throughout 3D printed objects. Furthermore, taking advantage of the internal porosity of the printed parts, an internal printing approach is proposed for the localized deposition of a guest material within a host matrix, enabling complex 3D material designs.

36 MATERIALS SCIENCE↗

Three-dimensional analytical model for the spatial variation of the foreshock electron distribution function - Systematics and comparisons with ISEE observations

A model which is consistent with the solar wind and shock surface boundary conditions for the foreshock electron distribution in the absence of wave-particle effects is formulated for an arbitrary location behind the magnetic tangent to the earth's bow shock. Variations of the gyrophase-averaged velocity distribution are compared and contrasted with in situ ISEE observations. It is found that magnetic mirroring of solar wind electrons is the most important process by which nonmonotonic reduced electron distributions in the foreshock are produced. Leakage of particles from the magnetosheath is shown to be relatively unimportant in determining reduced distributions that are nonmonotonic. The two-dimensional distribution function off the magnetic field direction is the crucial contribution in producing reduced distributions which have beams. The time scale for modification of the electron velocity distribution in velocity space can be significantly influenced by steady state spatial gradients in the background imposed by the curved shock geometry.

Fitzenreiter, R. J.↗

A general kinetic framework for dislocation mobility derived from probabilistic cellular automaton simulations of the kink-pair mechanism

Dislocation mobility laws are essential components of dislocation-density-based crystal plasticity models. For dislocations governed by the kink-pair mechanism, however, existing formulations are typically limited to specific regimes due to the com plex interplay between stochastic kink-pair nucleation and lateral kink migration. In this work, we develop a general kinetic framework that expresses the average dislocation velocity as a function of mechanism-level variables: positive/negative kink pair nucleation rates, kink migration velocity, dislocation segment length, critical kink-pair width, and kink height. Probabilis tic cellular automaton simulations are used to capture the behavior of conceptual dislocation segments between the limiting conditions of migration outpacing nucleation on the one end and nucleation outpacing migration on the other. An elemen tary functional form that captures the system dynamics is suggested and fitted against the simulation results. This framework remains valid for arbitrary combinations of the six variables and is, therefore, compatible with any admissible constitutive re lations that describe their stress and temperature dependence. Comparisons with established approaches and experimental results confirm the robustness and physical consistency of the formulation, making it broadly applicable to material systems in which dislocation motion is governed by the kink-pair mechanism.

36 MATERIALS SCIENCE↗

Flow-driven spectral chaos (FSC) method for simulating long-time dynamics of arbitrary-order non-linear stochastic dynamical systems

Uncertainty quantification techniques such as the time-dependent generalized polynomial chaos (TD-gPC) use an adaptive orthogonal basis to better represent the stochastic part of the solution space (aka random function space) in time. However, because the random function space is constructed using tensor products, TD-gPC-based methods are known to suffer from the curse of dimensionality. Here, we introduce a new numerical method called the flow-driven spectral chaos (FSC) which overcomes this curse of dimensionality at the random-function-space level. The proposed method is not only computationally more efficient than existing TD-gPC-based methods but is also far more accurate. The FSC method uses the concept of enriched stochastic flow maps to track the evolution of a finite-dimensional random function space efficiently in time. To transfer the probability information from one random function space to another, two approaches are developed and studied herein. In the first approach, the probability information is transferred in the mean-square sense, whereas in the second approach the transfer is done exactly using a new theorem that was developed for this purpose. The FSC method can quantify uncertainties with high fidelity, especially for the long-time response of stochastic dynamical systems governed by ODEs of arbitrary order. Six representative numerical examples, including a nonlinear problem (the Van-der-Pol oscillator), are presented to demonstrate the performance of the FSC method and corroborate the claims of its superior numerical properties. Finally, a parametric, high-dimensional stochastic problem is used to demonstrate that when the FSC method is used in conjunction with Monte Carlo integration, the curse of dimensionality can be overcome altogether.

(nonlinear) stochastic dynamical systems↗

Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product Report

Objective analysis (OA) is a method of mapping unstructured data to a structured grid. In the context of scanning radar data, OA is used to interpolate data in antenna coordinates (range, azimuth, and elevation) onto a regularly spaced Cartesian grid (Trapp and Doswell 2000). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Mapped Moments to a Cartesian Grid (MMCG) Value-Added Product (VAP) uses the Python ARM Radar Toolkit (Py-ART), a data model-driven interactive architecture for working with weather radar data, to map the data to a Cartesian grid (Helmus and Collis 2016). MMCG, with Py-ART built in, has the ability to take radar data in antenna coordinates and map the gates to a Cartesian grid using inverse distance weight functions such as Cressman (square) and Barnes (exponential), but also can filter the data during the interpolation. MMCG also allows arbitrary formulations for the radius of influence, which are matched to particular radar scanning strategies. This creates a complex parameter space for optimizing the retention of storm structure detail while minimizing artifacts. MMCG takes data processed with ARM’s Corrected Precipitation Radar Moments in Antenna Coordinates (CMAC) VAP and maps it to a Cartesian grid as the output product. A variety of fields that have been mapped to the Cartesian grid are then saved to plots to complement each grid file.

54 ENVIRONMENTAL SCIENCES↗

The eXtended virtual element method for elliptic problems with weakly singular solutions

This paper introduces a novel eXtended virtual element method, an extension of the conforming virtual element method. The X-VEM is formulated by incorporating appropriate enrichment functions in the local spaces. The method is designed to handle highly generic enrichment functions, including singularities arising from fractured domains. By achieving consistency on the enrichment space, the method is proven to achieve arbitrary approximation orders even in the presence of singular solutions. The paper includes a complete convergence analysis under general assumptions on mesh regularity, and numerical experiments validating the method’s accuracy on various mesh families, demonstrating optimal convergence rates in the L 2 - and H 1 - norms on fractured or L-shaped domains.

97 MATHEMATICS AND COMPUTING↗

Propagation of waves of acoustic frequencies in curved ducts

The propagation of waves of acoustic frequencies in curved ducts is studied for the first four modes. The analysis makes use of Bessel functions to construct curves of wave number in the duct versus imposed wave number. The results apply to ducts of arbitrary width and arbitrary radii of curvature. The characteristics of motion in a bend are compared with propagation of waves in a straight duct, and important differences in the behavior of waves are noted.

Rostafinski, W.↗