Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SERIES EXPANSION”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Constraining f(R) theories with cosmography

A method to set constraints on the parameters of extended theories of gravitation is presented. It is based on the comparison of two series expansions of any observable that depends on H(z). The first expansion is of the cosmographical type, while the second uses the dependence of H with z furnished by a given type of extended theory. When applied to f(R) theories together with the redshift drift, the method yields limits on the parameters of two examples (the theory of Hu and Sawicki [1], and the exponential gravity introduced by Linder [2]) that are compatible with or more stringent than the existing ones, as well as a limit for a previously unconstrained parameter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Using Hessian update formulae to construct modified Shepard interpolated potential energy surfaces: Application to vibrating surface atoms

Modified Shepard interpolation based on second order Taylor series expansions has proven to be a flexible tool for constructing potential energy surfaces in a range of situations. Extending this to gas–surface dynamics where surface atoms are allowed to move represents a substantial increase in the dimensionality of the problem, reflected in a dramatic increase in the computational cost of the required Hessian (matrix of second derivatives) evaluations. This work demonstrates that using approximate Hessians derived from well known Hessian update formulae and a single accurate Hessian can provide an effective way to avoid this expensive accurate Hessian determination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Dependence of the Energy Resolution of a Hemispherical Semiconductor Detector on the Bias Voltage

It is shown that the series expansion of the amplitude and variance of the hemispherical semiconductor detector signal in inverse bias voltage allows finding the Fano factor, the product of electron lifetime and mobility, the degree of inhomogeneity of the trap density in the semiconductor material, and the relative variance of the electronic channel gain. An important advantage of the proposed method is that it is independent of the electronic channel gain and noise.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Spectral softening in the X-RAY afterglow of GRB 130925A as predicted by the dust scattering model

Gamma-ray bursts (GRBs) usually occur in a dense star-forming region with a massive circumburst medium. The small-angle scattering of intense prompt X-ray emission off the surrounding dust grains will have observable consequences and sometimes can dominate the X-ray afterglow. In most of the previous studies, only the Rayleigh-Gans (RG) approximation is employed for describing the scattering process, which works accurately for the typical size of grains (with radius of a ≤ 0.1 μm) in the diffuse interstellar medium. When the size of the grains may significantly increase, as in a more dense region where GRBs would occur, the RG approximation may not be valid enough for modeling detailed observational data. In order to study the temporal and spectral properties of the scattered X-ray emission more accurately with potentially larger dust grains, we provide a practical approach using the series expansions of anomalous diffraction (AD) approximation based on the complicated Mie theory. We apply our calculations to understand the puzzling X-ray afterglow of recently observed GRB 130925A that showed a significant spectral softening. We find that the X-ray scattering scenarios with either AD or RG approximation adopted could well reproduce both the temporal and spectral profile simultaneously. Given the plateau present in the early X-ray light curve, a typical distribution of smaller grains as in the interstellar medium would be suggested for GRB 130925A.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Dynamic Simulation with on the fly Doppler Broadening treatment in RMC

In conventional transient analysis, thermal hydraulics calculations are generally coupled to deterministic neutronics solvers. Due to the limited accuracy of deterministic methods and the rapid development of computer power, Monte Carlo neutronics calculations have received increasing attention nowadays. One of challenges for Monte Carlo methods coupled to thermal hydraulics calculations is the temperature dependent cross sections. Taking the detailed temperature distribution feedback into account, each nuclide may be in a large range of temperatures, pre-generated libraries at every temperature point are no longer applicable owing to unaffordable memory footprint. Different solutions have been proposed to deal with this problem. TMS, accounting for Doppler broadening effect by the explicit treatment of thermal motion, was used to model fast transients in Serpent 2. Stochastic mixing, using two sets of libraries for each nuclide and mixing them to get the correct temperature, was used in transient analysis by Sjenitzer and Hoogenboom. Besides, there are other on the fly Doppler broadening methods being developed, including series expansion of cross sections, multipole representation, different integration algorithms based on Doppler broadening equation. In this paper, on the fly Doppler broadening treatment based on integration algorithms was applied to transient calculation in Reactor Monte Carlo code RMC . In the following section, the integration algorithms implemented and dynamic simulation method in RMC is briefly introduced. Finally, a pure physical transient fuel assembly model was chosen for numerical validation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

From liquid crystal models to the guiding-center theory of magnetized plasmas

Upon combining Northrop’s picture of charged particle motion with modern liquid crystal theories, this paper provides a new description of guiding center dynamics (to lowest order). This new perspective is based on a rotation gauge field (gyrogauge) that encodes rotations around the magnetic field. In liquid crystal theory, an analogue rotation field is used to encode the rotational state of rod-like molecules. Instead of resorting to sophisticated tools (e.g. Hamiltonian perturbation theory and Lie series expansions) that still remain essential in higher-order gyrokinetics, the present approach combines the WKB method with a simple kinematical ansatz, which is then replaced into the charged particle Lagrangian. The latter is eventually averaged over the gyrophase to produce the guiding-center equations. A crucial role is played by the vector potential for the gyrogauge field. A similar vector potential is related to liquid crystal defects and is known as wryness tensor in Eringen’s micropolar theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

ANALYTICAL CALCULATION OF STOKES PROFILES OF ROTATING STELLAR MAGNETIC DIPOLE

The observation of the polarization emerging from a rotating star at different phases opens up the possibility to map the magnetic field in the stellar surface thanks to the well-known Zeeman-Doppler imaging. When the magnetic field is sufficiently weak, the circular and linear polarization profiles locally in each point of the star are proportional to the first and second derivatives of the unperturbed intensity profile, respectively. We show that the weak-field approximation (for weak lines in the case of linear polarization) can be generalized to the case of a rotating star including the Doppler effect and taking into account the integration on the stellar surface. The Stokes profiles are written as a linear combination of wavelength-dependent terms expressed as series expansions in terms of Hermite polynomials. These terms contain the surface-integrated magnetic field and velocity components. The direct numerical evaluation of these quantities is limited to rotation velocities not larger than eight times the Doppler width of the local absorption profiles. Additionally, we demonstrate that in a rotating star, the circular polarization flux depends on the derivative of the intensity flux with respect to the wavelength and also on the profile itself. Likewise, the linear polarization depends on the profile and on its first and second derivatives with respect to the wavelength. We particularize the general expressions to a rotating dipole.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Index theorem and Majorana zero modes along a non-Abelian vortex in a color superconductor

Color superconductivity in high-density QCD exhibits the color-flavor-locked phase. To explore zero modes in the color-flavor-locked phase in the presence of a non-Abelian vortex with an SU(2) symmetry in the vortex core, we apply the index theorem to the Bogoliubov-de Gennes (BdG) Hamiltonian. From the calculation of the topological index, we find that triplet, doublet and singlet sectors of SU(2) have certain number of chiral Majorana zero modes in the limit of vanishing chemical potential. We also solve the BdG equation by the use of the series expansion to show that the number of zero modes and their chirality match the result of the index theorem. From particle-hole symmetry of the BdG Hamiltonian, we conclude that if and only if the index of a given sector is odd, one zero mode survives generically for a finite chemical potential. We argue that this result should hold nonperturbatively even in the high-density limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Determination of plasma frequency, damping constant, and size distribution from the complex dielectric function of noble metal nanoparticles

This paper develops a novel method for simultaneously determining the plasma frequency ω{sub P} and the damping constant γ{sub free} in the bulk damped oscillator Drude model, based on experimentally measured real and imaginary parts of the metal refractive index in the IR wavelength range, lifting the usual approximation that restricts frequency values to the UV-deep UV region. Our method was applied to gold, silver, and copper, improving the relative uncertainties in the final values for ω{sub p} (0.5%–1.6%) and for γ{sub free} (3%–8%), which are smaller than those reported in the literature. These small uncertainties in ω{sub p} and γ{sub free} determination yield a much better fit of the experimental complex dielectric function. For the case of nanoparticles (Nps), a series expansion of the Drude expression (which includes ω{sub p} and γ{sub free} determined using our method) enables size-dependent dielectric function to be written as the sum of three terms: the experimental bulk dielectric function plus two size corrective terms, one for free electron, and the other for bound-electron contributions. Finally, size distribution of nanometric and subnanometric gold Nps in colloidal suspension was determined through fitting its experimental optical extinction spectrum using Mie theory based on the previously determined dielectric function. Results are compared with size histogram obtained from Transmission Electron Microscopy (TEM)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Spectral Formulation for the Solution of Full-Wave Scattering from a Conducting Wedge Tipped with a Corrugated Cylinder

A spectral mode-matching technique is formulated to solve for the full-wave scattering of a corrugated cylinder-tipped wedge in the presence of an impressed electric or magnetic line source. Asymptotic approximations of large-order Bessel or Henkel functions for a fixed argument were introduced in order to overcome numerical difficulties in their regular series expansions. The corrugations on the conducting cylinder have the shape of annular sectors. The primary objective of this work is to investigate the impact of corrugations on the scattered field in the shadow region of the structure. An optimally designed corrugated cylinder placed at the tip of a conducting wedge can effectively suppress electromagnetic scattering in the shadow region. Obtained numerical results using the proposed approach prove the above concept. These results were validated against numerical data obtained using a nodal finite element method. The aim of this research is to utilize these corrugated tips in horn antenna design for the reduction of side-lobe level and the shaping of the respective E-plane radiation pattern.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Multi-layer hydrostatic equilibrium of planets and synchronous moons: theory and application to Ceres and to solar system moons

The hydrostatic equilibrium of multi-layer bodies lacks a satisfactory theoretical treatment despite its wide range of applicability. Here we show that by using the exact analytical potential of homogeneous ellipsoids we can obtain recursive analytical solutions and an exact numerical method for the hydrostatic equilibrium shape problem of multi-layer planets and synchronous moons. The recursive solutions rely on the series expansion of the potential in terms of the polar and equatorial shape eccentricities, while the numerical method uses the exact potential expression. These solutions can be used to infer the interior structure of planets and synchronous moons from their observed shape, rotation, and gravity. When applied to the dwarf planet Ceres, we show that it is most likely a differentiated body with an icy crust of equatorial thickness 30-90 km and a rocky core of density 2.4-3.1 g cm{sup –3}. For synchronous moons, we show that the J {sub 2}/C {sub 22} ≅ 10/3 and the (b – c)/(a – c) ≅ 1/4 ratios have significant corrections of order Ω{sup 2}/(πGρ), with important implications for how their gravitational coefficients are determined from fly-by radio science data and for how we assess their hydrostatic equilibrium state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Microscopic description of anisotropic low-density dipolar Bose gases in two dimensions

A microscopic description of the zero-energy two-body ground state and many-body static properties of anisotropic homogeneous gases of bosonic dipoles in two dimensions at low densities is presented and discussed. By changing the polarization angle with respect to the plane, we study the impact of the anisotropy, present in the dipole-dipole interaction, on the energy per particle, comparing the results with mean-field predictions. We restrict the analysis to the regime where the interaction is always repulsive, although the strength of the repulsion depends on the orientation with respect to the polarization field. We present a series expansion of the solution of the zero-energy two-body problem, which allows us to find the scattering length of the interaction and to build a suitable Jastrow factor that we use as a trial wave function for both a variational and diffusion Monte Carlo simulation of the infinite system. We find that the anisotropy has an almost negligible impact on the ground-state properties of the many-body system in the universal regime where the scattering length governs the physics of the system. We also show that scaling in the gas parameter persists in the dipolar case up to values where other isotropic interactions with the same scattering length yield different predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Numerical methods for high-dimensional probability density function equations

In this paper we address the problem of computing the numerical solution to kinetic partial differential equations involving many phase variables. These types of equations arise naturally in many different areas of mathematical physics, e.g., in particle systems (Liouville and Boltzmann equations), stochastic dynamical systems (Fokker–Planck and Dostupov–Pugachev equations), random wave theory (Malakhov–Saichev equations) and coarse-grained stochastic systems (Mori–Zwanzig equations). We propose three different classes of new algorithms addressing high-dimensionality: The first one is based on separated series expansions resulting in a sequence of low-dimensional problems that can be solved recursively and in parallel by using alternating direction methods. The second class of algorithms relies on truncation of interaction in low-orders that resembles the Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) framework of kinetic gas theory and it yields a hierarchy of coupled probability density function equations. The third class of algorithms is based on high-dimensional model representations, e.g., the ANOVA method and probabilistic collocation methods. A common feature of all these approaches is that they are reducible to the problem of computing the solution to high-dimensional equations via a sequence of low-dimensional problems. The effectiveness of the new algorithms is demonstrated in numerical examples involving nonlinear stochastic dynamical systems and partial differential equations, with up to 120 variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Study of the grazing-incidence X-ray scattering of strongly disturbed fractal surfaces

The applicability of different approaches to the description of hard X-ray scattering from rough surfaces is generally limited by a maximum surface roughness height of no more than 1 nm. Meanwhile, this value is several times larger for the surfaces of different materials subjected to treatment, especially in the initial treatment stages. To control the roughness parameters in all stages of surface treatment, a new approach has been developed, which is based on a series expansion of wavefield over the plane eigenstate-function waves describing the small-angle scattering of incident X-rays in terms of plane q-waves propagating through the interface between two media with a random function of relief heights. To determine the amplitudes of reflected and transmitted plane q-waves, a system of two linked integral equations was derived. The solutions to these equations correspond (in zero order) to the well-known Fresnel expressions for a smooth plane interface. Based on these solutions, a statistical fractal model of an isotropic rough interface is built in terms of root-mean-square roughness σ, two-point correlation length l, and fractal surface index h. The model is used to interpret X-ray scattering data for polished surfaces of single-crystal cadmium telluride samples.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Resonance scattering and radiation force calculations for an elastic cylinder using the translational addition theorem for cylindrical wave functions

The standard Resonance Scattering Theory (RST) of plane waves is extended for the case of any two-dimensional (2D) arbitrarily-shaped monochromatic beam incident upon an elastic cylinder with arbitrary location using an exact methodology based on Graf’s translational addition theorem for the cylindrical wave functions. The analysis is exact as it does not require numerical integration procedures. The formulation is valid for any cylinder of finite size and material that is immersed in a nonviscous fluid. Partial-wave series expansions (PWSEs) for the incident, internal and scattered linear pressure fields are derived, and the analysis is further extended to obtain generalized expressions for the on-axis and off-axis acoustic radiation force components. The wave-fields are expressed using generalized PWSEs involving the beam-shape coefficients (BSCs) and the scattering coefficients of the cylinder. The off-axial BSCs are expressed analytically in terms of an infinite PWSE with emphasis on the translational offset distance d. Numerical computations are considered for a zeroth-order quasi-Gaussian beam chosen as an example to illustrate the analysis. Acoustic resonance scattering directivity diagrams are calculated by subtracting an appropriate background from the expression of the scattered pressure field. In addition, computations for the radiation force exerted on an elastic cylinder centered on the axis of wave propagation of the beam, and shifted off-axially are analyzed and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

LOGISTIC FUNCTION PROFILE FIT: A least-squares program for fitting interface profiles to an extended logistic function

The extended logistic function provides a physically reasonable description of interfaces such as depth profiles or line scans of surface topological or compositional features. It describes these interfaces with the minimum number of parameters, namely, position, width, and asymmetry. Logistic Function Profile Fit (LFPF) is a robust, least-squares fitting program in which the nonlinear extended logistic function is linearized by a Taylor series expansion (equivalent to a Newton-Raphson approach) with no apparent introduction of bias in the analysis. The program provides reliable confidence limits for the parameters when systematic errors are minimal and provides a display of the residuals from the fit for the detection of systematic errors. The program will aid researchers in applying ASTM E1636-10, 'Standard practice for analytically describing sputter-depth-profile and linescan-profile data by an extended logistic function,' and may also prove useful in applying ISO 18516: 2006, 'Surface chemical analysis-Auger electron spectroscopy and x-ray photoelectron spectroscopy-determination of lateral resolution.' Examples are given of LFPF fits to a secondary ion mass spectrometry depth profile, an Auger surface line scan, and synthetic data generated to exhibit known systematic errors for examining the significance of such errors to the extrapolation of partial profiles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Polariton spectra under the collective coupling regime. I. Efficient simulation of linear spectra and quantum dynamics

We outline two general theoretical techniques to simulate polariton quantum dynamics and optical spectra under the collective coupling regimes described by a Holstein–Tavis–Cummings (HTC) model Hamiltonian. The first one takes advantage of sparsity of the HTC Hamiltonian, which allows one to reduce the cost of acting polariton Hamiltonian onto a state vector to the linear order of the number of states, instead of the quadratic order. The second one is applying the well-known Chebyshev series expansion approach for quantum dynamics propagation and to simulate the polariton dynamics in the HTC system; this approach allows us to use a much larger time step for propagation and only requires a few recursive operations of the polariton Hamiltonian acting on state vectors. These two theoretical approaches are general and can be applied to any trajectory-based non-adiabatic quantum dynamics methods. We apply these two techniques with our previously developed Lindblad-partially linearized density matrix approach to simulate the linear absorption spectra of the HTC model system, with both inhomogeneous site energy disorders and dipolar orientational disorders. Our numerical results agree well with the previous analytic and numerical work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗