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 523 records · Page 29

Stochastic Modeling in a Multimaterial Continuum Mixture Shock Physics Code

Stochastic modelling approaches are presented to capture random effects at multiple time and length scales. Random processes that occur at the microscale produce nondeterministic effects at the macroscale. Here we present three stochastic modeling approaches that describe random processes at microscopic length scales and map these processes to the macroscopic length scale. The first stochastic modeling approach is based upon a particle based numerical technique to solve a Stochastic Differential Equation (SDE) using an arbitrary diffusion process to capture random processes at the microstructural level. The second approach prescribes a Probability Density Function (PDF) for the drift and diffusion of the random variable derived using the forward and backward Kolmogorov equations. This method requires mean and drift evolution PDF transport equations. The third approach is the coupling of multiple random variables which are dependent on each other. The relationship of the PDFs and a coupling function, known as a copula, produces a Joint Probability Density Function (JPDF). These stochastic modeling approaches are implemented into a Multiple Component (MC) shock physics computational code and used to model statistical fracture and reactive flow applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Lossy radial diffusion of relativistic Jovian electrons

The radial diffusion equation with synchrotron losses is solved by the Laplace-transform method for near equatorially mirroring relativistic electrons. The evolution of a power-law distribution function is found, and the characteristics of synchrotron burnoff are stated in terms of explicit parameters for an arbitrary diffusion coefficient of a specific form. The peaking of the 10.4-cm volume emissivity from Jupiter at an L shell of about 1.8 provides an estimate of the diffusion coefficient in the radiation belts; one value is suggested as the appropriate modification, for an equatorial field strength of 4.2 G, of the Birmingham et al. (1974) result. Nonsynchrotron losses are included phenomenologically; from the phase-space densities reported by McIlwain and Fillius (1975), the particle lifetime is estimated. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.

Barbosa, D. D.↗

A mathematical model of post-instability in fluid mechanics

Postinstability of fluids is eliminated in numerical models by introducing multivalued velocity fields after discarding the principle of impenetrability. Smooth functions are shown to be incapable of keeping the derivatives from going towards infinity when iterating solutions for the governing equations such as those defined by Navier-Stokes. Enlarging the class of functions is shown to be necessary to eliminate the appearance of imaginary characteristic roots in the systems of arbitrary partial differential equations, a condition which leads to physically impossible motions. The enlarging is demonstrated to be achievable by allowing several individual particles with different velocities to appear at the same point of space, and the subsequent multivaluedness of the solutions is purely a mathematical concern, rather than one of actual physical existence. Applications are provided for an inviscid fluid and for turbulence.

Zak, M. A.↗

Radiative transfer in a sphere illuminated by a parallel beam - An integral equation approach

The problem of multiple scattering of nonpolarized light in a planetary body of arbitrary shape illuminated by a parallel beam is formulated using the integral equation approach. There exists a simple functional whose stationarity condition is equivalent to solving the equation of radiative transfer and whose value at the stationary point is proportional to the differential cross section. The analysis reveals a direct relation between the microscopic symmetry of the phase function for each scattering event and the macroscopic symmetry of the differential cross section for the entire planetary body, and the interconnection of these symmetry relations and the variational principle. The case of a homogeneous sphere containing isotropic scatterers is investigated in detail. It is shown that the solution can be expanded in a multipole series such that the general spherical problem is reduced to solving a set of decoupled integral equations in one dimension. Computations have been performed for a range of parameters of interest, and illustrative examples of applications to planetary problems as provided.

Shia, R.-L.↗

Time domain unsteady incompressible cascade airfoil theory for helicopter rotors in hover

A detailed derivation of a finite-time arbitrary-motion incompressible cascade theory is presented for both Laplace and frequency domains. The generalized cascade lift deficiency function is shown to be consistent with the generalized Theodorsen's lift deficiency function when the wake spacing approaches infinity or when the reduced frequency tends to infinity; it also yields a correct value for the zero reduced frequency limit. Efficient numerical procedures are presented for the evaluation of the cascade lift deficiency function. Numerical examples comparing the cascade lift deficiency function with Loewy's (1957) rotary-wing lift deficiency function in frequency domain are presented. Pade approximants of the cascade lift deficiency function are constructed using a Bode plot approach that allows for complex poles. A general-purpose optimization program is used to determine the coefficients of the approximant.

Dinyavari, M. A. H.↗

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗

Some Experiences with High Order Discretization Schemes for Compressible Fluid Flow

It is widely recognized that the structure and complexity of compressible fluid flow places high demands on numerical discretization techniques for the fluid flow equations. Fluid flows arising in external aerodynamics often contain both flow field discontinuities and fluid boundary-layers. Both must be accurately resolved to provide useful information to aerodynamic design and analysis engineers. These accuracy requirements motivated the present author to examine a class of finite-volume techniques on arbitrary triangulated domains based on linear or quadratic reconstruction of integral-averaged data followed by upwind flux function evaluation and small time evolution. More recently, we have considered some new upwind techniques which yield compact discretizations while maintaining higher order accuracy. In the mini-symposium talk we will discuss both of these techniques as well as demonstrate the relative merits of each method by computing a number of aerodynamic flows containing shock waves and boundary-layers.

Barth, Timothy J.↗

A Systematic Methodology for Constructing High-Order Energy-Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter (AIAA 2008-2876, 2008) was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables \energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗

High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables 'energy stable' modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗

A Systematic Methodology for Constructing High-Order Energy Stable WENO Schemes

A third-order Energy Stable Weighted Essentially Non{Oscillatory (ESWENO) finite difference scheme developed by Yamaleev and Carpenter [1] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables "energy stable" modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes; ESWENO schemes up to eighth order are presented in the appendix. New weight functions are also developed that provide (1) formal consistency, (2) much faster convergence for smooth solutions with an arbitrary number of vanishing derivatives, and (3) improved resolution near strong discontinuities.

Yamaleev, Nail K.↗

Design of a neutral thermal scattering (NeTS) module for hydrogen in light water

The accurate representation of thermal scattering law (TSL) data is integral to the design and characterization of many modern nuclear systems, particularly those using light water as a moderator/coolant. As a material-dependent distribution over energy-momentum phase space, the TSL may exhibit a variety of unique and relevant conditional (e.g., temperature, pressure) and compositional (e.g., porosity, stoichiometry, radiation damage) dependencies. Currently, there are various approaches to incorporate temperature dependence, which is especially important in coupled neutronic-thermal hydraulic simulations. In each approach, there is an inherent tradeoff between memory consumption and accuracy. Some techniques require tens to hundreds of MBs or more, while others fail to reproduce the underlying data to within 10% error over the considered input domain, despite having a reduced storage burden. This work aims to address both sides of the tradeoff simultaneously by implementing a novel deep learning (DL) approach to TSL representation. The neural thermal scattering (NeTS) concept, which is amenable to an arbitrary number of dependencies, is demonstrated via the inclusion of temperature dependence into a highly compact, highly accurate functional form of the multi- variate TSL (i.e., S(α, β, T)) for hydrogen in light water. Resulting storage requirements are on the order of 100 kB, and median and maximum percent deviations are on the order of 0.1% and 1%, respectively. These measures represent a step improvement over previous techniques. Notably, the developed neural network and feature methodology build on those employed in prior work on beryllium oxide. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Pick-and-Place Transfer of Arbitrary-Metal Electrodes for van der Waals Device Fabrication

Van der Waals electrode integration is a promising strategy to create nearly perfect interfaces between metals and 2D materials, with advantages such as eliminating Fermi-level pinning and reducing contact resistance. However, the lack of a simple, generalizable pick-and-place transfer technology has greatly hampered the wide use of this technique. Here, we demonstrate the pick-and-place transfer of prefabricated electrodes from reusable polished hydrogenated diamond substrates without the use of any sacrificial layers due to the inherent low-energy and dangling-bond-free nature of the hydrogenated diamond surface. The technique enables transfer of arbitrary-metal electrodes and an electrode array, as demonstrated by successful transfer of eight different elemental metals with work functions ranging from 4.22 to 5.65 eV. We also demonstrate the electrode array transfer for large-scale device fabrication. The mechanical transfer of metal electrodes from diamond to van der Waals materials creates atomically smooth interfaces with no interstitial impurities or disorder, as observed with cross-section high-resolution transmission electron microscopy and energy-dispersive X-ray spectroscopy. As a demonstration of its device application, we use the diamond transfer technique to create metal contacts to monolayer transition metal dichalcogenide semiconductors with high-work-function Pd, low-work-function Ti, and semimetal Bi to create n- and p-type field-effect transistors with low Schottky barrier heights. We also extend this technology to air-sensitive materials (trilayer 1T’ WTe 2 ) and other applications such as ambipolar transistors, Schottky diodes, and optoelectronics. This highly reliable and reproducible technology paves the way for new device architectures and high-performance devices.

2D materials↗

Chordwise and compressibility corrections for arbitrary planform slender wings

The Lomax and Sluder method for adapting slender-wing theory to delta or rectangular wings by making chordwise and compressibility corrections is extended to cover wings of any arbitrary planform in subsonic and supersonic flows. The numerical accuracy of the present work is better than that of the Lomax-Sluder results. Comparison of the results of this work with those of the vortex-lattice method and Kernel function method for a family of Gothic and arrowhead wings shows good agreement. A universal curve is proposed for the evaluation of the lift coefficient of a low aspect ratio wing of an arbitrary planform in subsonic flow. The location of the center of pressure can also be estimated.

Levin, D.↗

Sum rules and static local-field corrections of electron liquids in two and three dimensions

The dielectric functions of electron liquids which take into account short-range electron-electron correlations via the static local-field corrections are examined in the light of the frequency-moment sum rules. The formation is given for degenerate as well as classical electron liquids in arbitrary (d) spatial dimensions, which is suitable for comparison between the two- and three-dimensional cases. By using the virial equations of state it is shown that such dielectric functions cannot satisfy the compressibility sum rule and the third-frequency-moment sum rule simultaneously. In the degenerate case, the plasmon, single-pair, and multipair contributions to the sum rules are analyzed, and the reason for this incompatibility is discussed.

Iwamoto, N.↗

On hierarchical solutions to the BBGKY hierarchy

It is thought that the gravitational clustering of galaxies in the universe may approach a scale-invariant, hierarchical form in the small separation, large-clustering regime. Past attempts to solve the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy in this regime have assumed a certain separable hierarchical form for the higher order correlation functions of galaxies in phase space. It is shown here that such separable solutions to the BBGKY equations must satisfy the condition that the clustered component of the solution has cluster-cluster correlations equal to galaxy-galaxy correlations to all orders. The solutions also admit the presence of an arbitrary unclustered component, which plays no dyamical role in the large-clustering regime. These results are a particular property of the specific separable model assumed for the correlation functions in phase space, not an intrinsic property of spatially hierarchical solutions to the BBGKY hierarchy. The observed distribution of galaxies does not satisfy the required conditions. The disagreement between theory and observation may be traced, at least in part, to initial conditions which, if Gaussian, already have cluster correlations greater than galaxy correlations.

Hamilton, A. J. S.↗

Dynamic Multipoles in the RCS

The Rapid Cycling Synchrotron (RCS) of the Electron--Ion Collider (EIC) will accelerate polarized electrons up to 18 GeV at a repetition rate of 1 Hz. Dynamic multipoles driven by eddy currents in the vacuum chamber can affect beam dynamics and polarization preservation during the energy ramp. Building on the analytical treatment of S.Y. Lee for a ramped normal-conducting dipole magnet, we use conformal transformations to derive compact expressions for dynamic multipoles in higher-order ramped magnets with circular beam pipes. We further derive analytical frequency-domain transfer functions for both self-response and cross-response multipoles, providing compact dynamical models for their evolution during arbitrary ramp waveforms. The resulting analytical predictions are benchmarked against finite-element simulations, including the effects of practical pole geometry and beam-pipe misalignment. Finally, the expected dynamic multipoles in the RCS magnets are estimated and their implications for beam dynamics and polarization preservation are discussed.

43 PARTICLE ACCELERATORS↗