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 559 records · Page 31

Useful extremum principle for the variational calculation of matrix elements. II

Recent work (Gerjuoy et al., 1974) on variational principles for diagonal bound state matrix elements of arbitrary Hermitian operators is extended. In particular, it is shown that the previously derived minimum principle for the trial auxiliary function appearing in such variational principles can be constructed using a modified Hamiltonian possessing not heretofore recognized positive definite properties. Thus there is at least one alternative to the particular modified Hamiltonian on which the results of Gerjuoy et al. (1974) originally were based.

Gerjuoy, E.↗

Lossy radial diffusion of relativistic Jovian electrons

The radial diffusion equation with synchrotron losses was solved by the Laplace transform method for near-equatorially mirroring relativistic electrons. The evolution of a power law distribution function was found and the characteristics of synchrotron burn-off are stated in terms of explicit parameters for an arbitrary diffusion coefficient. Emissivity from the radiation belts of Jupiter was studied. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.

Barbosa, D. D.↗

The application of the principles of invariance to the radiative transfer equation in plant canopies

Solutions of the radiative transfer equation describing photon interactions with vegetation canopies are important in remote sensing since they provide the canopy reflectance distribution required in the interpretation of satellite acquired information. The general one-dimensional two-angle transport problem for a finite copy of arbitrary leaf angle distribution is considered. Analytical solutions are obtained in terms of generalized Chandrasekhar's X- and Y-functions by invoking the principles of invariance. A critical step in the formulation involves the decomposition of the integral of the scattering phase function into a product of known functions of the incident and scattered photon directions. Several simplified cases previously considered in the literature are derived from the generalized solution. Various symmetries obeyed by the scattering operator and reciprocity relations are formally proved.

Ganapol, B. D.↗

Optimal control problems with mixed control-phase variable equality and inequality constraints

In this paper, necessary conditions are obtained for optimal control problems containing equality constraints defined in terms of functions of the control and phase variables. The control system is assumed to be characterized by an ordinary differential equation, and more conventional constraints, including phase inequality constraints, are also assumed to be present. Because the first-mentioned equality constraint must be satisfied for all t (the independent variable of the differential equation) belonging to an arbitrary (prescribed) measurable set, this problem gives rise to infinite-dimensional equality constraints. To obtain the necessary conditions, which are in the form of a maximum principle, an implicit-function-type theorem in Banach spaces is derived.

Makowski, K.↗

A numerical method for determination of source time functions for general three-dimensional rupture propagation

A method to determine the displacement and the stress on the crack plane for a three-dimensional shear crack of arbitrary shape propagating in an infinite, homogeneous medium which is linearly elastic everywhere off the crack plane is presented. The main idea of the method is to use a representation theorem in which the displacement at any given point on the crack plane is written as an integral of the traction over the whole crack plane. As a test of the accuracy of the numerical technique, the results are compared with known solutions for two simple cases.

Das, S.↗

Extremely luminous far-infrared colliding galaxies - Computation of colliding-galaxy luminosity function

A computational model used to determine the luminosity function for gaseous galaxies in collision is described. It is based on the assumption of high-velocity encounters of two thin, gas-rich disks having identical radii and thickness, able to approach each other from arbitrary directions, with arbitrary orientations and with arbitrary degree of overlap. The distribution of expected luminosities generated in the dissipation of translational energy is then found to closely fit the observed distribution of luminosities of extremely luminous IRAS galaxies. Bursts of star formation resulting from such collisions would not have the correct luminosity function.

Harwit, Martin↗

Approximating accelerator impedances with resonator networks

It is common in the accelerator community to use the impedance of accelerator components to describe wake interactions in the frequency domain. However, it is often desirable to understand such wake interactions in the time domain in a general manner for excitations that are not necessarily Gaussian in nature. Here, the conventional method for doing this involves taking the inverse Fourier Transform of the component impedance, obtaining the Green's Function, and then convolving it with the desired excitation distribution. This method can prove numerically cumbersome, for a convolution integral must be evaluated for each individual point in time when the wake function is desired. An alternative to this method would be to compute the wake function analytically, which would sidestep the need for repetitive integration. Only a handful of cases, however, are simple enough for this method to be tenable. One of these cases is the case where the component in question is an RLC resonator, which has a closed-form analytical wake function solution. This means that a component which can be represented in terms of resonators can leverage this solution. As it happens, common network synthesis techniques may be used to map arbitrary impedance profiles to RLC resonator networks in a manner the accelerator community has yet to take advantage of. In this work, we will use Foster Canonical Resonator Networks and partial derivative descent optimization to develop a technique for synthesizing resonator networks that well approximate the impedances of real-world accelerator components. We will link this synthesis to the closed-form resonator wake function solution, giving rise to a powerful workflow that may be used to streamline beam dynamics simulations.

43 PARTICLE ACCELERATORS↗

Drell-Yan q T resummation of fiducial power corrections at N 3 LL

We consider Drell-Yan production pp → V*X → LX at small q T << Q, where q T and Q are the total transverse momentum and invariant mass of the leptonic final state L. Experimental measurements require fiducial cuts on L, which in general introduce enhanced, linear power corrections in q T /Q. We show that they can be unambiguously predicted from factorization, and resummed to the same order as the leading-power contribution. For the fiducial q T spectrum, they constitute the complete linear power corrections. We thus obtain predictions for the fiducial q T spectrum to N 3 LL and next-to-leading-power in q T /Q. Matching to full NNLO ($α$$^{2}_{s}$), we find that the linear power corrections are indeed the dominant ones, and once included by factorization, the remaining fixed-order corrections become almost negligible below q T ≲ 40 GeV. We also discuss the implications for more complicated observables, and provide predictions for the fiducial Φ* spectrum at N 3 LL+NNLO. We find excellent agreement with ATLAS and CMS measurements of q T and Φ*. We also consider the $p$$^{ℓ}_{T}$ spectrum. We show that it develops leptonic power corrections in q T /(Q – 2$p$$^{ℓ}_{T}$), which diverge near the Jacobian peak $p$$^{ℓ}_{T}$ ~ Q/2 and must be kept to all powers to obtain a meaningful result there. Doing so, we obtain for the first time an analytically resummed result for the $p$$^{ℓ}_{T}$ spectrum around the Jacobian peak at N 3 LL+NNLO. Our method is based on performing a complete tensor decomposition for hadronic and leptonic tensors. We show that in practice this is equivalent to often-used recoil prescriptions, for which our results now provide rigorous, formal justification. Our tensor decomposition yields nine Lorentz-scalar hadronic structure functions, which for Z/γ* → ℓℓ or W → ℓν directly map onto the commonly used angular coefficients, but also holds for arbitrary leptonic final states. In particular, for suitably defined Born-projected leptons it still yields a LO-like angular decomposition even when including QED final-state radiation. Finally, we also discuss the application to q T subtractions. Including the unambiguously predicted fiducial power corrections significantly improves their performance, and in particular makes them applicable near kinematic edges where they otherwise break down due to large leptonic power corrections.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Medial axis and local thickness computation using the Fast Sweeping Method

This report describes an efficient and robust voxel-based methodology for computing the medial axis, local thickness, and distance-to-skeleton of arbitrary three-dimensional geometries. It is assumed that the object can be represented by an exact or approximate signed distance function on a discrete grid. The gradient of such function is used to formulate a hyperbolic partial differential equation (PDE) that models the collapse of the position vector in space. By exploiting the causality property of the PDE, the Fast Sweeping Method is able to obtain the solution in a finite number of sweeps independent of the mesh resolution. The intersection of characteristic lines leads to the formation of shocks and a discrete bisector function is used to identify the medial axis. The same PDE approach is used to compute the local thickness inside the object and obtain the distance-to-skeleton field. Multiple examples are given in two and three dimensions along with a resolution study. The methodology has optimal complexity and yields subsecond computational times for geometries with over a million zones on a single core. The methodology is also capable of parallelization across shared and distributed memory architectures.

97 MATHEMATICS AND COMPUTING↗

Adiabatic Hamiltonian of charged particle motion in a dipole field

The Hamiltonian for a dipole field is developed, and the result is expressed by an analytic approximation accurate to within about 1%. This allows extension of results derived for equatorial particles to particles with arbitrary pitch angles; in particular, it makes available even in the presence of electric fields orthogonal to the magnetic field a function K that is preserved by the bounce-averaged motion. This function provides at once the equations of drift paths in (alpha, beta) or of their projections onto the equatorial plane; the derivation of a pacing function that times the progress of particles along such drift paths is also described.

Chen, A. J.↗

MMIC devices for active phased array antennas

The study of printed circuit discontinuities is necessary in order to design, for example, transitions between rectangular waveguides and printed circuits. New developments with respect to the analytical approaches to this problem are discussed. A summary of the progress in the experimental approach is presented. The accurate solution for the modes in various millimeter-wave waveguides is essential in the analysis of many integrated circuit components, such as filters and impedance transformers. Problems associated with the numerical computation of these modes in two frequently used waveguide forms, namely, the finline and microstrip, are presented. The spectral domain method of formulation, with a moment method solution, is considered. This approach can be readily extended to analyze an arbitrary configuration of dielectric and metallized regions in a shielded enclosure. Galerkin's method is used, where the testing and basic functions are the same. It is shown that the mode functions, or eigenfunctions, are more sensitive to errors than the phase constants, or eigenvalues. The approximate mode functions do not satisfy the orthogonality relationship well, resulting in difficulties when these modal solutions are used to form an approximate Green's function or are used in a mode matching analysis.

Mittra, R.↗

Whistler damping at oblique propagation - Laminar shock precursors

This paper addresses the collisionless damping of whistlers observed as precursors standing upstream of oblique, low-Mach number terrestrial bow shocks. The linear theory of electromagnetic waves in a homogeneous Vlasov plasma with Maxwellian distribution functions and a magnetic field is considered. Numerical solutions of the full dispersion equation are presented for whistlers propagating at an arbitrary angle with respect to the magnetic field. It is demonstrated that electron Landau damping attenuates oblique whistlers and that the parameter which determines this damping is beta-e. In a well-defined range of parameters, this theory provides damping lengths which are the same order of magnitude as those observed. Thus electron Landau damping is a plausible process in the dissipation of upstream whistlers. Nonlinear plasma processes which may contribute to precursor damping are also discussed, and criteria for distinguishing among these are described.

Gary, S. P.↗

Light scattering of large rough particles application to cometary grains

While the electromagnetic field scattered by a spherical particle is classically obtained by the Helmholtz equation, the general case of an arbitrary particle may be investigated in the general framework of the interaction of a wave with a scattering potential. A wave function then satisfies the Schroedinger equation. The general solution of the Schroedinger equation is given. The main disadvantage of this approach are its restriction to large particles and its scalar nature preventing the calculation of the polarization. However, Perrin and Lamy have shown how to avoid the second limitation and retrieve a vectorial description. They proved that in the case of large spheres when the ad hoc assumptions are satisfied, the expression of the scattering amplitude may be approximated by an expansion series in partial waves, i.e., on a discrete basis. The analogy may be generalized, and the ratio of the two components for a rough particle obtained by taking the ratio of the reflectivities for the two directions of polarization. These reflectivities involve the simple and double reflections calculated following the method developed by Wolff for rough surfaces. The theory is further detailed.

Lamy, Philippe L.↗

A 32-bit Ultrafast Parallel Correlator using Resonant Tunneling Devices

An ultrafast 32-bit pipeline correlator has been implemented using resonant tunneling diodes (RTD) and hetero-junction bipolar transistors (HBT). The negative differential resistance (NDR) characteristics of RTD's is the basis of logic gates with the self-latching property that eliminates pipeline area and delay overheads which limit throughput in conventional technologies. The circuit topology also allows threshold logic functions such as minority/majority to be implemented in a compact manner resulting in reduction of the overall complexity and delay of arbitrary logic circuits. The parallel correlator is an essential component in code division multi-access (CDMA) transceivers used for the continuous calculation of correlation between an incoming data stream and a PN sequence. Simulation results show that a nano-pipelined correlator can provide and effective throughput of one 32-bit correlation every 100 picoseconds, using minimal hardware, with a power dissipation of 1.5 watts. RTD plus HBT based logic gates have been fabricated and the RTD plus HBT based correlator is compared with state of the art complementary metal oxide semiconductor (CMOS) implementations.

Kulkarni, Shriram↗

A Communication Channel Density Estimating Generative Adversarial Network

Autoencoder-based communication systems use neural network channel models to backwardly propagate message reconstruction error gradients across an approximation of the physical communication channel. In this work, we develop and test a new generative adversarial network (GAN) architecture for the purpose of training a stochastic channel approximating neural network. In previous research, investigators have focused on additive white Gaussian noise (AWGN) channels and/or simplified Rayleigh fading channels, both of which are linear and have well defined analytic solutions. Given that training a neural network is computationally expensive, channel approximation networks— and more generally the autoencoder systems—should be evaluated in communication environments that are traditionally difficult. To that end, our investigation focuses on channels that contain a combination of non-linear amplifier distortion, pulse shape filtering, intersymbol interference, frequency-dependent group delay, multipath, and non-Gaussian statistics. Each of our models are trained without any prior knowledge of the channel. We show that the trained models have learned to generalize over an arbitrary amplifier drive level and constellation alphabet. We demonstrate the versatility of our GAN architecture by comparing the marginal probability density function of several channel simulations with that of their corresponding neural network approximations

Smith, Aaron↗

Equilibrium of magnetic fields with arbitrary interweaving of the lines of force. I - Discontinuities in the torsion

Consideration is given to the static force-free equilibrium of a magnetic field in which all of the lines of force connect without knotting between parallel planes. The field is formed by continuous deformation from an initial uniform field, and is conventiently described in terms of the scalar function psi, which is the stream function for the incompressible wrapping and interweaving of the lines of force. Local compression and expansion of the lines of force is described in terms of the scalar function Phi. Equilibrium in the field requires satisfaction of two independent equations which cannot be accomplished without the full freedom of both psi and Phi. It is shown that discontinuities in the torsional characteristics of the lines occur when psi is predetermined by an arbitrary pattern. Discontinuities in the winding pattern of the lines can lead to discontinuities in the associated current sheets.

Parker, E. N.↗

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗

Resolution of a Few Problems in the Application of Quasilinear Theory to Calculating Diffusion Coefficients in Heliophysics

A theory of quasilinear diffusion for obliquely propagating electromagnetic waves was developed in the 1960's and applied in the 1970's to model scattering of relativistic electrons by a prescribed distribution of waves. In the latter work, a transformation of variables, from wavevector space to the temporal frequency and tangent of the wave normal angle, was used so that simple Gaussian functions of frequency and tangent of the wave normal angle could be multiplied together to define the distribution of wave power, although arbitrary distributions of power in these two variables is also permitted. Finally, in 2005, previous work was consolidated and has been widely used in heliophysics studies that require computation of quasilinear diffusion coefficients. Here it is shown that this transformation is inppropriate when the precise wave vector distribution is known. The correct transformation is derived and used to produce diffusion coefficients that can differ by orders of magnitude from those computed using the inappropriate transformation. The differences are largest when the distribution of wave power extends to wave normal angles near the resonance cone. When the ratio of the plasma frequency to the gyrofrequency is large, only low energies (keV) are affected, but as the ratio decreases higher energies (MeV) also show differences. It is also shown that the derivation from the 1960's uses a notation that results in the diffusion coefficients depending on the distribution of wave power with respect to the wave azimuthal angle whereas there should be no such dependence.

79 ASTRONOMY AND ASTROPHYSICS↗