SEARCH · Engineering Papers
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.
We present initial results of the studies of localized line-like sources for a scalar field in anti-de Sitter space. Such configurations reveal some subtleties and require novel numerical analysis. Dual field theory configurations involve simple models of lattices and defects.
In the present paper, we endeavor to perform a numerical analysis in connection with the nonlinear radiative stagnation-point flow and heat transfer to Sisko fluid past a stretching cylinder in the presence of convective boundary conditions. The influence of thermal radiation using nonlinear Rosseland approximation is explored. The numerical solutions of transformed governing equations are calculated through forth order Runge-Kutta method using shooting technique. With the help of graphs and tables, the influence of non-dimensional parameters on velocity and temperature along with the local skin friction and Nusselt number is discussed. The results reveal that the temperature increases however, heat transfer from the surface of cylinder decreases with the increasing values of thermal radiation and temperature ratio parameters. Moreover, the authenticity of numerical solutions is validated by finding their good agreement with the HAM solutions.
The computer modeling methods of a compact isochronous cyclotron are described. The main stages of analysis of accelerator facilities systems are considered. The described methods are based on theoretical fundamentals of cyclotron physics and mention highlights of creation of the physical project of a compact cyclotron. The main attention is paid to the analysis of the beam dynamics, formation of a magnetic field, stability of the movement, and a realistic assessment of intensity of the generated bunch of particles. In the article, the stages of development of the accelerator computer model, analytical ways of assessment of the accelerator parameters, and the basic technique of the numerical analysis of dynamics of the particles are described.
A numerical analysis of fusion welding of steel- and aluminum-based foam materials is carried out. The schemes of the structured and stochastic pore distribution are considered. The research results were used to conduct the experiments which confirmed the reliability of the numerical calculations.
We consider a mathematical model for the determination of the temperature field in a layer that emits, absorbs, and scatters infrared radiation both in the volume and on the surface. Within the framework of the model, we formulate nonlinear direct and inverse problems of emission infrared tomography of the temperature field in the layer according to known space-and-angular distributions of the flows of infrared radiation emitted into the ambient medium. The iterative-variational methods are proposed for the solution of the indicated direct and inverse problems. As a specific example, we perform the numerical analysis of the developed algorithms.
Bulk acoustic wave devices are typically operated in a resonant state to achieve enhanced acoustic amplitudes and high acoustofluidic forces for the manipulation of microparticles. Among other loss mechanisms related to the structural parts of acoustofluidic devices, damping in the fluidic cavity is a crucial factor that limits the attainable acoustic amplitudes. In the analytical part of this study, we quantify all relevant loss mechanisms related to the fluid inside acoustofluidic micro-devices. Subsequently, a numerical analysis of the time-harmonic visco-acoustic and thermo-visco-acoustic equations is carried out to verify the analytical results for 2D and 3D examples. The damping results are fitted into the framework of classical linear acoustics to set up a numerically efficient device model. For this purpose, all damping effects are combined into an acoustofluidic loss factor. Since some components of the acoustofluidic loss factor depend on the acoustic mode shape in the fluid cavity, we propose a two-step simulation procedure. In the first step, the loss factors are deduced from the simulated mode shape. Subsequently, a second simulation is invoked, taking all losses into account. Owing to its computational efficiency, the presented numerical device model is of great relevance for the simulation of acoustofluidic particle manipulation by means of acoustic radiation forces or acoustic streaming. For the first time, accurate 3D simulations of realistic micro-devices for the quantitative prediction of pressure amplitudes and the related acoustofluidic forces become feasible.
In the present work, a numerical analysis is performed to investigate the comparative contribution of the mechanisms responsible for electron gain and losses in laser spark ignition and plasma formation of H{sub 2}. The analysis considered H{sub 2} over pressure range 150 -3000 torr irradiated by a Nd:YAG laser radiation at wavelengths 1064 and 532 nm with pulse length 5.5 ns. The study based on a modified electron cascade model by one of the authors which solves numerically the time dependent Boltzmann equation as well as a set of rate equations that describe the rate of change of the excited states population. The model includes most of the physical processes that might take place during the interaction. Computations of The threshold intensity are performed for the combined and separate contribution of each of the gain and loss processes. Reasonable agreement with the measured values over the tested pressure range is obtained only for the case of the combined contribution. Basing on the calculation of the electron energy distribution function, the determined relations of the time evolution of the electrons density for selected values of the tested gas pressure region revealed that photo-ionization of the excited states could determine the time of electron generation and hence spark ignition. Collisional ionization contributes to this phenomenon only at the high pressure regime. Loss processes due to electron diffusion, vibrational excitation are found to have significant effect over examined pressure values for the two applied laser wavelengths.
Various techniques for interpolation of data, moment matching in stochastic applications and various methods in numerical analysis can be described using Vandermonde matrices. For this reason the properties of the determinant of the Vandermonde matrix and related matrices are interesting. Here the extreme points of the Vandermonde determinant, and related determinants, on some simple surfaces such as the unit sphere are analyzed, both numerically and analytically. Some results are also visualized in various dimensions. The extreme points of the Vandermonde determinant are also related to the roots of certain orthogonal polynomials such as the Hermite polynomials.
The waveguide considered coincides with a strip having two narrows of width ε. An electron wave function satisfies the Dirichlet boundary value problem for the Helmholtz equation. The part of the waveguide between the narrows serves as a resonator, and conditions for the electron resonant tunneling may occur. In the paper, asymptotic formulas as ε → 0 for characteristics of the resonant tunneling are used. The asymptotic results are compared with the numerical ones obtained by approximate calculation of the scattering matrix for energies in the interval between the second and third thresholds. The comparison allows us to state an interval of ε, where the asymptotic and numerical approaches agree. The suggested methods can be applied to more complicated models than that considered in the paper. In particular, the same approach can be used for asymptotic and numerical analysis of the tunneling in three-dimensional quantum waveguides of variable cross-sections. Bibliography: 3 titles.
Numerical analysis of gas bearing stability based on bearing film step-jump responses, and use of Duhamel theorem
We discuss diagrammatic modifications to the coupled cluster doubles (CCD) equations, wherein different groups of terms out of rings, ladders, crossed-rings, and mosaics can be removed to form approximations to the coupled cluster method, of interest due to their similarity with various types of random phase approximations. The finite uniform electron gas (UEG) is benchmarked for 14- and 54-electron systems at the complete basis set limit over a wide density range and performance of different flavours of CCD is determined. These results confirm that rings generally overcorrelate and ladders generally undercorrelate; mosaics-only CCD yields a result surprisingly close to CCD. We use a recently developed numerical analysis [J. J. Shepherd and A. Grüneis, Phys. Rev. Lett. 110, 226401 (2013)] to study the behaviours of these methods in the thermodynamic limit. We determine that the mosaics, on forming the Brueckner one-body Hamiltonian, open a gap in the effective one-particle eigenvalues at the Fermi energy. Numerical evidence is presented which shows that methods based on this renormalisation have convergent energies in the thermodynamic limit including mosaic-only CCD, which is just a renormalised MP2. All other methods including only a single channel, namely, ladder-only CCD, ring-only CCD, and crossed-ring-only CCD, appear to yield divergent energies; incorporation of mosaic terms prevents this from happening.
The effect of the current rise time in a relativistic electron beam pulse on the tracking force exerted by the low-conductivity Ohmic plasma channel on the beam has been analyzed using the “hard” beam model. It is shown using numerical analysis that an increase in the beam current rise time substantially reduces the tracking force considered here.
We show that the Higgs portal interactions involving extra dark Higgs field can save generically the original Higgs inflation of the standard model (SM) from the problem of a deep non-SM vacuum in the SM Higgs potential. Specifically, we show that such interactions disconnect the top quark pole mass from inflationary observables and allow multi-dimensional parameter space to save the Higgs inflation, thanks to the additional parameters (the dark Higgs boson mass m {sub φ}, the mixing angle α between the SM Higgs H and dark Higgs Φ, and the mixed quartic coupling) affecting RG-running of the Higgs quartic coupling. The effect of Higgs portal interactions may lead to a larger tensor-to-scalar ratio, 0.08 ∼< r ∼< 0.1, by adjusting relevant parameters in wide ranges of α and m {sub φ}, some region of which can be probed at future colliders. Performing a numerical analysis we find an allowed region of parameters, matching the latest Planck data.
The interaction of magnetron-sputtered metal iron with titanium-oxide films upon isothermal vacuum annealing is studied by X-ray phase analysis, secondary-ion mass spectrometry, atomic-force microscopy, and mathematical simulation. A mechanism for the formation of complex oxides at grain boundaries is suggested. The mechanism is based on the reaction diffusion of metal iron into titanium oxide. A quantitative model of reaction interdiffusion in two-layer polycrystalline metal–oxide film systems with limited component solubility is developed. From numerical analysis of the experimental distributions of the metal concentrations in the Fe–TiO{sub 2} film system, the individual diffusion coefficients are determined. It is found that, under the conditions of vacuum annealing at 1073 K, the diffusion coefficients of iron and titanium are 8.0 × 10{sup –13} and 3.0 × 10{sup –15} cm{sup 2} s{sup –1}, respectively.
Highlights: • It is studied emission correlations between particle emission events in black holes. • Particles emitted have highly correlated moments before the BH stops radiating. • Rainbow Gravity allows maintaining the rest of it stored in its remnant. In this work, we study emission correlations in Rainbow Gravity (RG) black holes known to have black hole remnants. Non-thermal corrections are responsible for generating non-zero correlations between particle emission events during the black hole evaporation process. With this in mind, we calculate the temperatures considering back-reaction effects. Then, we obtain the correlations between emission events for a RG metric proposed by Magueijo and Smolin. We also discuss through a numerical analysis the emission correlations behavior in the last stages of evaporation of the black hole.
The stochastic protein kinetic equations can be stiff for certain parameters, which makes their numerical simulation rely on very small time step sizes, resulting in large computational cost and accumulated round-off errors. For such situation, we provide a method of reducing stiffness of the stochastic protein kinetic equation by means of a kind of variable transformation. Theoretical and numerical analysis show effectiveness of this method. Its generalization to a more general class of stochastic differential equation models is also discussed.
We present results of analytical and numerical investigation of explosive evaporation and burning scenarios of hydrogen droplets in hydrogen/oxygen aerosols. The following two scenarios have been elucidated. The first scenario, corresponding to sufficiently large droplets, is characterized by three stages: (i) an essentially homogeneous heating of a droplet to a near-critical temperature by IR radiation from the hot gas; (ii) explosive evaporation; and (iii) burning of hydrogen cloud formed by evaporation. The second scenario, corresponding to small droplets, differs in that a droplet is heated mainly by thermal conduction from the hot gas. The heating is accompanied by evaporation which can become explosive at the final stage of evaporation. The crossover droplet size separating the two scenarios is calculated. Conservative finite-difference numerical analysis is used to explore the predicted scenarios and verify analytical estimates.
The performance of commercial nuclear spent fuel dry storage casks are typically evaluated through detailed numerical analysis of the system's thermal performance. These modeling efforts are performed by the vendor to demonstrate the performance and regulatory compliance and are independently verified by the Nuclear Regulatory Commission (NRC). Carefully measured data sets generated from testing of full sized casks or smaller cask analogs are widely recognized as vital for validating these models. Numerous studies have been previously conducted. Recent advances in dry storage cask designs have moved the storage location from above ground to below ground and significantly increased the maximum thermal load allowed in a cask in part by increasing the canister helium pressure. Previous cask performance validation testing did not capture these parameters. The purpose of the investigation described in this report is to produce a data set that can be used to test the validity of the assumptions associated with the calculations presently used to determine steady-state cladding temperatures in modern dry casks. These modern cask designs utilize elevated helium pressure in the sealed canister or are intended for subsurface storage. The BWR cask simulator (BCS) has been designed in detail for both the above ground and below ground venting configurations. The pressure vessel representing the canister has been designed, fabricated, and pressure tested for a maximum allowable pressure (MAWP) rating of 24 bar at 400 C. An existing electrically heated but otherwise prototypic BWR Incoloy-clad test assembly is being deployed inside of a representative storage basket and cylindrical pressure vessel that represents the canister. The symmetric single assembly geometry with well-controlled boundary conditions simplifies interpretation of results. Various configurations of outer concentric ducting will be used to mimic conditions for above and below ground storage configurations of vertical, dry cask systems with canisters. Radial and axial temperature profiles will be measured for a wide range of decay power and helium cask pressures. Of particular interest is the evaluation of the effect of increased helium pressure on heat load and the effect of simulated wind on a simplified below ground vent configuration.