Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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 829 records · Page 46

Informed total-error-minimizing priors: Interpretable cosmological parameter constraints despite complex nuisance effects

While Bayesian inference techniques are standard in cosmological analyses, it is common to interpret resulting parameter constraints with a frequentist intuition. This intuition can fail, for example, when marginalizing high-dimensional parameter spaces onto subsets of parameters, because of what has come to be known as projection effects or prior volume effects. We present the method of informed total-error-minimizing (ITEM) priors to address this problem. An ITEM prior is a prior distribution on a set of nuisance parameters, such as those describing astrophysical or calibration systematics, intended to enforce the validity of a frequentist interpretation of the posterior constraints derived for a set of target parameters (e.g., cosmological parameters). Our method works as follows. For a set of plausible nuisance realizations, we generate target parameter posteriors using several different candidate priors for the nuisance parameters. We reject candidate priors that do not accomplish the minimum requirements of bias (of point estimates) and coverage (of confidence regions among a set of noisy realizations of the data) for the target parameters on one or more of the plausible nuisance realizations. Of the priors that survive this cut, we select the ITEM prior as the one that minimizes the total error of the marginalized posteriors of the target parameters. As a proof of concept, we applied our method to the density split statistics measured in Dark Energy Survey Year 1 data. We demonstrate that the ITEM priors substantially reduce prior volume effects that otherwise arise and that they allow for sharpened yet robust constraints on the parameters of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical analysis of three complex three-dimensional shock wave / turbulent boundary layer interaction flows

A set of Reynolds-averaged Navier-Stokes (RANS) computations with turbulence closure provided by the Spalart-Allmaras (SA) model have been carried out for prediction of shock-induced three-dimensional (3D) turbulent separated flows. The experimental data for different aerodynamic test configurations have been used for assessing credibility of the numerical method employed. Particularly, shock wave / turbulent boundary layer interactions (SWTBLI) in the vicinity of an asymmetric sharp double-fin (DF) with different (7° and 11°) deflection angles mounted on a flat plate and two conically sharp cylindrical bodies at varying interbody distances and nose cone angles 60° and 40° mounted over a flat plate at freestream Mach number 4, as well as a transonic fan stage operating in the near-stall regime are considered. The gas dynamic structure and topology of 3D separated flows, surface flow patterns, and pressure distributions as well as body aerodynamic force prediction are analyzed. A transonic fan stage operating in the near-stall regime and a possibility of applying flow control is investigated. High Performance Computing was employed to make high resolution computations of these flows possible, and advanced 3D visualization techniques were employed in order to improve understanding of the separating flow phenomena.

42 ENGINEERING↗

Oscillations of a sessile droplet in open air

The open system consisting of a sessile drop, a neutral gas, and a substrate is analyzed by numerical methods. The mode with constant contact angle is considered. The model takes into account evaporation from drop surface, free and forced convection in gas, buoyancy, and Marangoni effect in the liquid phase. It was established that every considered mechanical and thermodynamical disturbance of the system leads to the drop surface oscillations, and thus a drop in an open air oscillates almost inevitably. The displacement of the liquid-gas interface due to oscillations is analyzed in terms of its impact on the accuracy of measurement of the surface tension by sessile drop method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dielectric function of sol-gel prepared nano-granular zinc oxide by spectroscopic ellipsometry

ZnO thin films have been prepared by sol gel and deposited by spin coating. The dielectric function has been determined by spectroscopic ellipsometry. Ellipsometric spectra are inverted by a direct numerical method without using the standard fitting procedures. The obtained dielectric function presents a broad excitonic effect. The dielectric function is studied using Elliot excitonic theory including exciton plus band-to-band Coulomb interactions with standard Lorentzian broadening. A modification of this model dielectric function with independent bound and unbound exciton contributions is empirically proposed to improve modelling of the band gap excitonic peak.

36 MATERIALS SCIENCE↗

On the structural and impedance characteristics of Li- doped PEO, using n-butyl lithium in hexane as dopant

Nowadays polymer based solid state electrolytes for applications in rechargeable battery systems are highly sought after materials, pursued extensively by various research groups worldwide. Numerous methods are discussed in literature to improve the fundamental properties like electrical conductivity, mechanical stability and interfacial stability of polymer based electrolytes. The application of these electrolytes in Li-ion cells is still in the amateur state, due to low ionic conductivity, low lithium transport number and the processing difficulties. The present work is an attempt to study the effects of Li doping on the structural and transport properties of the polymer electrolyte, poly-ethelene oxide (PEO) (Molecular weight: 200,000). Li doped PEO was obtained by treating PEO with n-Butyllithium in hexane for different doping concentrations. Structural characterization of the samples was done by XRD and FTIR techniques. Impedance measurements were carried out to estimate the ionic conductivity of Li doped PEO samples. It is seen that, the crystallinity of the doped PEO decreases on increasing the doping concentration. XRD and FTIR studies support this observation. It is inferred that, ionic conductivity of the sample is increasing on increasing the doping concentration since less crystallinity permits more ionic transport. Impedance measurements confirm the results quantitatively.

36 MATERIALS SCIENCE↗

Unraveling multi-spin effects in rotational resonance nuclear magnetic resonance using effective reduced density matrix theory

A quantum-mechanical model integrating the concepts of reduced density matrix and effective Hamiltonians is proposed to explain the multi-spin effects observed in rotational resonance (R{sup 2}) nuclear magnetic resonance (NMR) experiments. Employing this approach, the spin system of interest is described in a reduced subspace inclusive of its coupling to the surroundings. Through suitable model systems, the utility of our theory is demonstrated and verified with simulations emerging from both analytic and numerical methods. The analytic results presented in this article provide an accurate description/interpretation of R{sup 2} experimental results and could serve as a test-bed for distinguishing coherent/incoherent effects in solid-state NMR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of a static electric field on two-color photoassociation between different atoms

We study non-perturbative effects of a static electric field on two-color photoassociation of different atoms. A static electric field induces anisotropy in scattering between two different atoms and hybridizes field-free rotational states of heteronuclear dimers or polar molecules. In a previous paper [D. Chakraborty et al., J. Phys. B 44, 095201 (2011)], the effects of a static electric field on one-color photoassociation between different atoms has been described through field-modified ground-state scattering states, neglecting electric field effects on heteronuclear diatomic bound states. To study the effects of a static electric field on heteronuclear bound states, and the resulting influence on Raman-type two-color photoassociation between different atoms in the presence of a static electric field, we develop a non-perturbative numerical method to calculate static electric field-dressed heteronuclear bound states. We show that the static electric field induced scattering anisotropy as well as hybridization of rotational states strongly influence two-color photoassociation spectra, leading to significant enhancement in PA rate and large shift. In particular, for static electric field strengths of a few hundred kV/cm, two-color PA rate involving high-lying bound states in electronic ground-state increases by several orders of magnitude even in the weak photoassociative coupling regime.

36 MATERIALS SCIENCE↗

Non-linear dynamic characteristics and optimal control of giant magnetostrictive film subjected to in-plane stochastic excitation

The non-linear dynamic characteristics and optimal control of a giant magnetostrictive film (GMF) subjected to in-plane stochastic excitation were studied. Non-linear differential items were introduced to interpret the hysteretic phenomena of the GMF, and the non-linear dynamic model of the GMF subjected to in-plane stochastic excitation was developed. The stochastic stability was analysed, and the probability density function was obtained. The condition of stochastic Hopf bifurcation and noise-induced chaotic response were determined, and the fractal boundary of the system's safe basin was provided. The reliability function was solved from the backward Kolmogorov equation, and an optimal control strategy was proposed in the stochastic dynamic programming method. Numerical simulation shows that the system stability varies with the parameters, and stochastic Hopf bifurcation and chaos appear in the process; the area of the safe basin decreases when the noise intensifies, and the boundary of the safe basin becomes fractal; the system reliability improved through stochastic optimal control. Finally, the theoretical and numerical results were proved by experiments. The results are helpful in the engineering applications of GMF.

36 MATERIALS SCIENCE↗

Stochastic bifurcation and fractal and chaos control of a giant magnetostrictive film-shape memory alloy composite cantilever plate subjected to in-plane harmonic and stochastic excitation

Stochastic bifurcation and fractal and chaos control of a giant magnetostrictive film–shape memory alloy (GMF–SMA) composite cantilever plate subjected to in-plane harmonic and stochastic excitation were studied. Van der Pol items were improved to interpret the hysteretic phenomena of both GMF and SMA, and the nonlinear dynamic model of a GMF–SMA composite cantilever plate subjected to in-plane harmonic and stochastic excitation was developed. The probability density function of the dynamic response of the system was obtained, and the conditions of stochastic Hopf bifurcation were analyzed. The conditions of noise-induced chaotic response were obtained in the stochastic Melnikov integral method, and the fractal boundary of the safe basin of the system was provided. Finally, the chaos control strategy was proposed in the stochastic dynamic programming method. Numerical simulation shows that stochastic Hopf bifurcation and chaos appear in the parameter variation process. The boundary of the safe basin of the system has fractal characteristics, and its area decreases when the noise intensifies. The system reliability was improved through stochastic optimal control, and the safe basin area of the system increased.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of positive direct current corona pulse in coaxial wire-cylinder gap

In this paper, a system is designed and developed to measure the positive corona current in coaxial wire-cylinder gaps. The characteristic parameters of corona current pulses, such as the amplitude, rise time, half-wave time, and repetition frequency, are statistically analyzed and a new set of empirical formulas are derived by numerical fitting. The influence of space charges on corona currents is tested by using three corona cages with different radii. A numerical method is used to solve a simplified ion-flow model to explain the influence of space charges. Based on the statistical results, a stochastic model is developed to simulate the corona pulse trains. And this model is verified by comparing the simulated frequency-domain responses with the measured ones.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasi-bound states and continuum absorption background of polar Al{sub 0.5}Ga{sub 0.5}N/GaN quantum dots

A theoretical interpretation of the photoluminescence excitation spectra of self-organized polar GaN/(Al,Ga)N quantum dots is presented. A numerical method assuming a realistic shape of the dots and including the built-in electric field effects is developed to calculate their energy structure and hence their optical absorption. The electron and hole spectra show the existence of a set of quasi-bound states that does not originate from the wetting layer and plays a crucial role in the observed absorption spectrum of the GaN/(Al,Ga)N dots. Transitions involving these quasi-bound states and wetting layer states give a sufficient explanation for the observed continuum absorption background. The properties of this absorption band, especially its extension, depend strongly on the dot's size. Our simulation provides a natural explanation of the experimental luminescence excitation spectra of ensembles of dots of different heights. Our theoretical model can be convenient for future optical studies including systems with more complicated potentials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Limitation of parallel flow in double diffusive convection: Two- and three-dimensional transitions in a horizontal porous domain

Two-dimensional (2D) and three-dimensional (3D) numerical simulations of double diffusion natural convection in an elongated enclosure filled with a binary fluid saturating a porous medium are carried out in the present work. The Boussinesq approximation is made in the formulation of the problem, and Neumann boundary conditions for temperature and concentration are adopted, respectively, on vertical and horizontal walls of the cavity. The used numerical method is based on the control volume approach, with the third order quadratic upstream interpolation scheme in approximating the advection terms. A semi implicit method algorithm is used to handle the velocity-pressure coupling. To avoid the excessively high computer time inherent to the solution of 3D natural convection problems, full approximation storage with full multigrid method is used to solve the problem. A wide range of the controlling parameters (Rayleigh-Darcy number Ra, lateral aspect ratio Ay, Lewis number Le, and the buoyancy ration N) is investigated. We clearly show that increasing the depth of the cavity (i.e., the lateral aspect ratio) has an important effect on the flow patterns. The 2D perfect parallel flows obtained for small lateral aspect ratio are drastically destabilized by increasing the cavity lateral dimension. This yields a 3D fluid motion with a much more complex flow pattern and the usually considered 2D parallel flow model cannot be applied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mode composition analysis on experimental results of a Gigawatt-class Ka-band overmoded Cerenkov oscillator

In this paper, a modified numerical method is used to investigate the mode composition of a Gigawatt-class Ka-band overmoded Cerenkov oscillator, which has been proposed and studied in our previous experiment. In the experiment, the measured angular distribution of radiation did not fit a single TM{sub 01} mode. So the particle in cell code calculations and the antenna radiation calculations are carried out, which show a consistent picture: the dominant modes are the TM{sub 01} mode and the TM{sub 03} mode, and their phase relationship is constant with time; therefore, a steady radiation pattern is produced, which matches the experimental data. As a conclusion, the comprehensive analysis shows that the existing modes of the output microwave in our experiment are the first five TM{sub 0n} modes (n = 1–5), with corresponding power ratios of 36.64%, 0.78%, 56.26%, 5.70%, and 0.52%, and relative phase differences of 0°, 146°, 54°, 169°, and 133°.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asymmetric modes decomposition in an overmoded relativistic backward wave oscillator

Most of the investigated overmoded relativistic backward wave oscillators (RBWOs) are azimuthally symmetric; thus, they are designed through two dimensional (2-D) particle-in-cell (PIC) simulations. However, 2-D PIC simulations cannot reveal the effect of asymmetric modes on beam-wave interaction. In order to investigate whether asymmetric mode competition needs to be considered in the design of overmoded RBWOs, a numerical method of determining the composition of both symmetric and asymmetric modes in three dimensional (3-D) PIC simulations is introduced in this paper. The 2-D and 3-D PIC simulation results of an X-band overmoded RBWO are analyzed. Our analysis indicates that the 2-D and 3-D PIC simulation results of our device are quite different due to asymmetric mode competition. In fact, asymmetric surface waves, especially EH{sub 11} mode, can lead to serious mode competition when electron beam propagates near the surface of slow wave structures (SWSs). Therefore, additional method of suppressing asymmetric mode competition, such as adjusting the reflections at both ends of SWSs to decrease the Q-factor of asymmetric modes, needs to be utilized in the design of overmoded RBWOs. Besides, 3-D PIC simulation and modes decomposition are essential for designing overmoded RBWOs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study on the thermal resistance in secondary particles chain of silica aerogel by molecular dynamics simulation

In this article, molecular dynamics simulation was performed to study the heat transport in secondary particles chain of silica aerogel. The two adjacent particles as the basic heat transport unit were modelled to characterize the heat transfer through the calculation of thermal resistance and vibrational density of states (VDOS). The total thermal resistance of two contact particles was predicted by non-equilibrium molecular dynamics simulations (NEMD). The defects were formed by deleting atoms in the system randomly first and performing heating and quenching process afterwards to achieve the DLCA (diffusive limited cluster-cluster aggregation) process. This kind of treatment showed a very reasonable prediction of thermal conductivity for the silica aerogels compared with the experimental values. The heat transport was great suppressed as the contact length increased or defect concentration increased. The constrain effect of heat transport was much significant when contact length fraction was in the small range (<0.5) or the defect concentration is in the high range (>0.5). Also, as the contact length increased, the role of joint thermal resistance played in the constraint of heat transport was increasing. However, the defect concentration did not affect the share of joint thermal resistance as the contact length did. VDOS of the system was calculated by numerical method to characterize the heat transport from atomic vibration view. The smaller contact length and greater defect concentration primarily affected the longitudinal acoustic modes, which ultimately influenced the heat transport between the adjacent particles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Communication: Newton homotopies for sampling stationary points of potential energy landscapes

One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) solutions but they all exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor ϕ{sup 4} model and atomic clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical matching of the sheath and presheath solutions for a spherical probe in radial-motion theory

The theory of positive-ion collection by a probe immersed in a low-pressure plasma was reviewed and extended by Allen, Boyd, and Reynolds [Proc. Phys. Soc. 70, 297 (1957)]. For a given value of the ion current, the boundary values of the matched “nonneutral” or “sheath” solution V{sup ~}{sub nn}{sup (m)}(r; r{sub m}) were obtained from the “quasineutral” or “presheath” solution V{sup ~}{sub qn}(r) by choosing the small potential and electric-field values corresponding to some large “matching radius” r{sub m}. Here, a straightforward but efficient numerical method is presented for systematically determining an optimal value of the matching radius at which the presheath and sheath solutions are joined to yield the “matched” potential profile. Some suitable initial matching radius r{sub m1} is chosen and the related potential and electric-field values of the quasineutral solution are calculated. Using these as boundary conditions, Poisson's equation is integrated to yield the matched nonneutral solution including the corresponding potential at the probe surface. This procedure is repeated for increasing values r{sub m2}, r{sub m3},…. until the resulting potential at the probe surface becomes practically constant. The corresponding value of r{sub m} is taken as the “optimal” matching radius r{sub mo} at which the presheath and sheath solutions are ultimately joined to yield the “optimal” matched potential profile in the entire plasma-probe transition region. It is also shown that the Bohm criterion is inapplicable in the present problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evolution of chirped laser pulses in a magnetized plasma channel

The propagation of intense, short, sinusoidal laser pulses in a magnetized plasma channel has been studied. The wave equation governing the evolution of the radiation field is set up and a variational technique is used to obtain the equations describing the evolution of the laser spot size, pulse length and chirp parameter. Numerical methods are used to analyze the simultaneous evolution of these parameters. The effect of the external magnetic field on initially chirped as well as unchirped laser pulses on the spot size, pulse length and chirping has been analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗