Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Inflection points of microcanonical entropy: Monte Carlo simulation of q state Potts model on a finite square lattice

Traditional definition of phase transition involves an infinitely large system in thermodynamic limit. Finite systems such as biological proteins exhibit cooperative behavior similar to phase transitions. We employ recently discovered analysis of inflection points of microcanonical entropy to estimate the transition temperature of the phase transition in q state Potts model on a finite two dimensional square lattice for q=3 (second order) and q=8 (first order). The difference of energy density of states (DOS) Δ ln g(E) = ln g(E+ ΔE) −ln g(E) exhibits a point of inflexion at a value corresponding to inverse transition temperature. This feature is common to systems exhibiting both first as well as second order transitions. While the difference of DOS registers a monotonic variation around the point of inflexion for systems exhibiting second order transition, it has an S-shape with a minimum and maximum around the point of inflexion for the case of first order transition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Average dynamics of a finite set of coupled phase oscillators

We study the solutions of a dynamical system describing the average activity of an infinitely large set of driven coupled excitable units. We compared their topological organization with that reconstructed from the numerical integration of finite sets. In this way, we present a strategy to establish the pertinence of approximating the dynamics of finite sets of coupled nonlinear units by the dynamics of its infinitely large surrogate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Electromagnetic gyrokinetic turbulence in finite-beta helical plasmas

A saturation mechanism for microturbulence in a regime of weak zonal flow generation is investigated by means of electromagnetic gyrokinetic simulations. The study identifies a new saturation process of the kinetic ballooning mode (KBM) turbulence originating from the spatial structure of the KBM instabilities in a finite-beta Large Helical Device (LHD) plasma. Specifically, the most unstable KBM in LHD has an inclined mode structure with respect to the mid-plane of a torus, i.e., it has a finite radial wave-number in flux tube coordinates, in contrast to KBMs in tokamaks as well as ion-temperature gradient modes in tokamaks and helical systems. The simulations reveal that the growth of KBMs in LHD is saturated by nonlinear interactions of oppositely inclined convection cells through mutual shearing as well as by the zonal flow. The saturation mechanism is quantitatively investigated by analysis of the nonlinear entropy transfer that shows not only the mutual shearing but also a self-interaction with an elongated mode structure along the magnetic field line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Finite-density effects in the Fredrickson-Andersen and Kob-Andersen kinetically-constrained models

We calculate the corrections to the thermodynamic limit of the critical density for jamming in the Kob-Andersen and Fredrickson-Andersen kinetically-constrained models, and find them to be finite-density corrections, and not finite-size corrections. We do this by introducing a new numerical algorithm, which requires negligible computer memory since contrary to alternative approaches, it generates at each point only the necessary data. The algorithm starts from a single unfrozen site and at each step randomly generates the neighbors of the unfrozen region and checks whether they are frozen or not. Our results correspond to systems of size greater than 10{sup 7} × 10{sup 7}, much larger than any simulated before, and are consistent with the rigorous bounds on the asymptotic corrections. We also find that the average number of sites that seed a critical droplet is greater than 1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Influence of finite geometrical asymmetry of the electrodes in capacitively coupled radio frequency plasma

Capacitively coupled radio frequency (CCRF) plasmas are widely studied in last decades due to the versatile applicability of energetic ions, chemically active species, radicals, and also energetic neutral species in many material processing fields including microelectronics, aerospace, and biology. A dc self-bias is known to generate naturally in geometrically asymmetric CCRF plasma because of the difference in electrode sizes known as geometrical asymmetry of the electrodes in order to compensate electron and ion flux to each electrode within one rf period. The plasma series resonance effect is also come into play due to the geometrical asymmetry and excited several harmonics of the fundamental in low pressure CCRF plasma. In this work, a 13.56 MHz CCRF plasma is studied on the based on the nonlinear global model of asymmetric CCRF discharge to understand the influences of finite geometrical asymmetry of the electrodes in terms of generation of dc self-bias and plasma heating. The nonlinear global model on asymmetric discharge has been modified by considering the sheath at the grounded electrode to taking account the finite geometrical asymmetry of the electrodes. The ion density inside both the sheaths has been taken into account by incorporating the steady-state fluid equations for ions considering that the applied rf frequency is higher than the typical ion plasma frequency. Details results on the influences of geometrical asymmetry on the generation of dc self-bias and plasma heating are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Energy spectra of finite temperature superfluid helium-4 turbulence

A mesoscopic model of finite temperature superfluid helium-4 based on coupled Langevin-Navier-Stokes dynamics is proposed. Drawing upon scaling arguments and available numerical results, a numerical method for designing well resolved, mesoscopic calculations of finite temperature superfluid turbulence is developed. The application of model and numerical method to the problem of fully developed turbulence decay in helium II, indicates that the spectral structure of normal-fluid and superfluid turbulence is significantly more complex than that of turbulence in simple-fluids. Analysis based on a forced flow of helium-4 at 1.3 K, where viscous dissipation in the normal-fluid is compensated by the Lundgren force, indicate three scaling regimes in the normal-fluid, that include the inertial, low wavenumber, Kolmogorov k{sup −5/3} regime, a sub-turbulence, low Reynolds number, fluctuating k{sup −2.2} regime, and an intermediate, viscous k{sup −6} range that connects the two. The k{sup −2.2} regime is due to normal-fluid forcing by superfluid vortices at high wavenumbers. There are also three scaling regimes in the superfluid, that include a k{sup −3} range that corresponds to the growth of superfluid vortex instabilities due to mutual-friction action, and an adjacent, low wavenumber, k{sup −5/3} regime that emerges during the termination of this growth, as superfluid vortices agglomerate between intense normal-fluid vorticity regions, and weakly polarized bundles are formed. There is also evidence of a high wavenumber k{sup −1} range that corresponds to the probing of individual-vortex velocity fields. The Kelvin waves cascade (the main dynamical effect in zero temperature superfluids) appears to be damped at the intervortex space scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Propulsive performance of a finite-temperature plasma flow in a magnetic nozzle with applied azimuthal current

The plasma flow in a finite-electron-temperature magnetic nozzle, under the influence of an applied azimuthal current at the throat, is modeled analytically to assess its propulsive performance. A correction to the nozzle throat boundary conditions is derived by modifying the radial equilibrium of a magnetized infinite two-population cylindrical plasma column with the insertion of an external azimuthal body force for the electrons. Inclusion of finite-temperature effects, which leads to a modification of the radial density profile, is necessary for calculating the propulsive performance, which is represented by nozzle divergence efficiency and thrust coefficient. The solutions show that the application of the azimuthal current enhances all the calculated performance parameters through the narrowing of the radial density profile at the throat, and that investing power in this beam focusing effect is more effective than using the same power to pre-heat the electrons. The results open the possibility for the design of a focusing stage between the plasma source and the nozzle that can significantly enhance the propulsive performance of electron-driven magnetic nozzles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Migration of finite sized particles in a laminar square channel flow from low to high Reynolds numbers

The migration of neutrally buoyant finite sized particles in a Newtonian square channel flow is investigated in the limit of very low solid volumetric concentration, within a wide range of channel Reynolds numbers Re = [0.07-120]. In situ microscope measurements of particle distributions, taken far from the channel inlet (at a distance several thousand times the channel height), revealed that particles are preferentially located near the channel walls at Re > 10 and near the channel center at Re < 1. Whereas the cross-streamline particle motion is governed by inertia-induced lift forces at high inertia, it seems to be controlled by shear-induced particle interactions at low (but finite) Reynolds numbers, despite the low solid volume fraction (<1%). The transition between both regimes is observed in the range Re = [1-10]. In order to exclude the effect of multi-body interactions, the trajectories of single freely moving particles are calculated thanks to numerical simulations based on the force coupling method. With the deployed numerical tool, the complete particle trajectories are accessible within a reasonable computational time only in the inertial regime (Re > 10). In this regime, we show that (i) the particle undergoes cross-streamline migration followed by a cross-lateral migration (parallel to the wall) in agreement with previous observations, and (ii) the stable equilibrium positions are located at the midline of the channel faces while the diagonal equilibrium positions are unstable. At low flow inertia, the first instants of the numerical simulations (carried at Re = O(1)) reveal that the cross-streamline migration of a single particle is oriented towards the channel wall, suggesting that the particle preferential positions around the channel center, observed in the experiments, are rather due to multi-body interactions.

42 ENGINEERING↗
DOE OSTI2014

Weakly nonlinear Rayleigh-Taylor instability of a finite-thickness fluid layer

A weakly nonlinear (WN) model has been developed for the Rayleigh-Taylor instability of a finite-thickness incompressible fluid layer (slab). We derive the coupling evolution equations for perturbations on the (upper) “linearly stable” and (lower) “linearly unstable” interfaces of the slab. Expressions of temporal evolutions of the amplitudes of the perturbation first three harmonics on the upper and lower interfaces are obtained. The classical feedthrough (interface coupling) solution obtained by Taylor [Proc. R. Soc. London A 201, 192 (1950)] is readily recovered by the first-order results. Our third-order model can depict the WN perturbation growth and the saturation of linear (exponential) growth of the perturbation fundamental mode on both interfaces. The dependence of the WN perturbation growth and the slab distortion on the normalized layer thickness (kd) is analytically investigated via the third-order solutions. Comparison is made with Jacobs-Catton's formula [J. W. Jacobs and I. Catton, J. Fluid Mech. 187, 329 (1988)] of the position of the “linearly unstable” interface. Using a reduced formula, the saturation amplitude of linear growth of the perturbation fundamental mode is studied. It is found that the finite-thickness effects play a dominant role in the WN evolution of the slab, especially when kd < 1. Thus, it should be included in applications where the interface coupling effects are important, such as inertial confinement fusion implosions and supernova explosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Formation of large-scale structures with sharp density gradient through Rayleigh-Taylor growth in a two-dimensional slab under the two-fluid and finite Larmor radius effects

Two-fluid and the finite Larmor effects on linear and nonlinear growth of the Rayleigh-Taylor instability in a two-dimensional slab are studied numerically with special attention to high-wave-number dynamics and nonlinear structure formation at a low β-value. The two effects stabilize the unstable high wave number modes for a certain range of the β-value. In nonlinear simulations, the absence of the high wave number modes in the linear stage leads to the formation of the density field structure much larger than that in the single-fluid magnetohydrodynamic simulation, together with a sharp density gradient as well as a large velocity difference. The formation of the sharp velocity difference leads to a subsequent Kelvin-Helmholtz-type instability only when both the two-fluid and finite Larmor radius terms are incorporated, whereas it is not observed otherwise. It is shown that the emergence of the secondary instability can modify the outline of the turbulent structures associated with the primary Rayleigh-Taylor instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Thermodynamics of the polaron master equation at finite bias

We study coherent transport through a double quantum dot. Its two electronic leads induce electronic matter and energy transport and a phonon reservoir contributes further energy exchanges. By treating the system-lead couplings perturbatively, whereas the coupling to vibrations is treated non-perturbatively in a polaron-transformed frame, we derive a thermodynamic consistent low-dimensional master equation. When the number of phonon modes is finite, a Markovian description is only possible when these couple symmetrically to both quantum dots. For a continuum of phonon modes however, also asymmetric couplings can be described with a Markovian master equation. We compute the electronic current and dephasing rate. The electronic current enables transport spectroscopy of the phonon frequency and displays signatures of Franck-Condon blockade. For infinite external bias but finite tunneling bandwidths, we find oscillations in the current as a function of the internal bias due to the electron-phonon coupling. Furthermore, we derive the full fluctuation theorem and show its identity to the entropy production in the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Linear vs. nonlinear acceleration in plasma turbulence. II. Hall–finite-Larmor-radius magnetohydrodynamics

The local k-space ratio of linear and nonlinear accelerations associated with a variety of initial conditions undergoing steady relaxation is investigated for the Hall–finite-Larmor-radius magnetohydrodynamics (MHD) system in the presence of a mean magnetic field. Building on a related study (Paper I) where it was shown that discrepancies exist between describing the global and local characterizations of the pure MHD system with mean magnetic field, we find regions of the Fourier space that are consistently dominated by linear acceleration and other regions that are consistently dominated by nonlinear acceleration, independent of the overall system's description as linear, weakly nonlinear, or turbulent. In general, dynamics within a certain angular range of the mean magnetic field direction are predominantly linear, while dynamics adjacent the Hall scales along the field-parallel direction and dynamics adjacent the finite Larmor radius scales in the field-perpendicular direction can become strongly nonlinear. The nonlinear influences are particularly significant as the plasma beta increases from unity to higher values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical analysis for finite-range multitype stochastic contact financial market dynamic systems

In an attempt to reproduce and study the dynamics of financial markets, a random agent-based financial price model is developed and investigated by the finite-range multitype contact dynamic system, in which the interaction and dispersal of different types of investment attitudes in a stock market are imitated by viruses spreading. With different parameters of birth rates and finite-range, the normalized return series are simulated by Monte Carlo simulation method and numerical studied by power-law distribution analysis and autocorrelation analysis. To better understand the nonlinear dynamics of the return series, a q-order autocorrelation function and a multi-autocorrelation function are also defined in this work. The comparisons of statistical behaviors of return series from the agent-based model and the daily historical market returns of Shanghai Composite Index and Shenzhen Component Index indicate that the proposed model is a reasonable qualitative explanation for the price formation process of stock market systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Switching times of nanoscale FePt: Finite size effects on the linear reversal mechanism

The linear reversal mechanism in FePt grains ranging from 2.316 nm to 5.404 nm has been simulated using atomistic spin dynamics, parametrized from ab-initio calculations. The Curie temperature and the critical temperature (T{sup *}), at which the linear reversal mechanism occurs, are observed to decrease with system size whilst the temperature window T{sup *}<T<T{sub C} increases. The reversal paths close to the Curie temperature have been calculated, showing that for decreasing system size the reversal path becomes more elliptic at lower temperatures, consistent with the decrease in the Curie temperature arising from finite size effects. Calculations of the minimum pulse duration show faster switching in small grains and are qualitatively described by the Landau-Lifshitz-Bloch equation with finite size atomistic parameterization, which suggests that multiscale modeling of FePt down to a grain size of ≈3.5 nm is possible.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Semiclassical modelling of finite-pulse effects on non-adiabatic photodynamics via initial condition filtering: The predissociation of NaI as a test case

Femtosecond-laser pulse driven non-adiabatic spectroscopy and dynamics in molecular and condensed phase systems continue to be a challenge for theoretical modelling. One of the main obstacles is the “curse of dimensionality” encountered in non-adiabatic, exact wavepacket propagation. A possible route towards treating complex molecular systems is via semiclassical surface-hopping schemes, in particular if they account not only for non-adiabatic post-excitation dynamics but also for the initial optical excitation. One such approach, based on initial condition filtering, will be put forward in what follows. As a simple test case which can be compared with exact wavepacket dynamics, we investigate the influence of the different parameters determining the shape of a laser pulse (e.g., its finite width and a possible chirp) on the predissociation dynamics of a NaI molecule, upon photoexcitation of the A(0{sup +}) state. The finite-pulse effects are mapped into the initial conditions for semiclassical surface-hopping simulations. The simulated surface-hopping diabatic populations are in qualitative agreement with the quantum mechanical results, especially concerning the subpicosend photoinduced dynamics, the main deviations being the relative delay of the non-adiabatic transitions in the semiclassical picture. Likewise, these differences in the time-dependent electronic populations calculated via the semiclassical and the quantum methods are found to have a mild influence on the overall probability density distribution. As a result, the branching ratios between the bound and the dissociative reaction channels and the time-evolution of the molecular wavepacket predicted by the semiclassical method agree with those computed using quantum wavepacket propagation. Implications for more challenging molecular systems are given.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Electronic chemical response indexes at finite temperature in the canonical ensemble

Assuming that the electronic energy is given by a smooth function of the number of electrons and within the extension of density functional theory to finite temperature, the first and second order chemical reactivity response functions of the Helmholtz free energy with respect to the temperature, the number of electrons, and the external potential are derived. It is found that in all cases related to the first or second derivatives with respect to the number of electrons or the external potential, there is a term given by the average of the corresponding derivative of the electronic energy of each state (ground and excited). For the second derivatives, including those related with the temperature, there is a thermal fluctuation contribution that is zero at zero temperature. Thus, all expressions reduce correctly to their corresponding chemical reactivity expressions at zero temperature and show that, at room temperature, the corrections are very small. When the assumption that the electronic energy is given by a smooth function of the number of electrons is replaced by the straight lines behavior connecting integer values, as required by the ensemble theorem, one needs to introduce directional derivatives in most cases, so that the temperature dependent expressions reduce correctly to their zero temperature counterparts. However, the main result holds, namely, at finite temperature the thermal corrections to the chemical reactivity response functions are very small. Consequently, the present work validates the usage of reactivity indexes calculated at zero temperature to infer chemical behavior at room and even higher temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Large-deviation joint statistics of the finite-time Lyapunov spectrum in isotropic turbulence

One of the hallmarks of turbulent flows is the chaotic behavior of fluid particle paths with exponentially growing separation among them while their distance does not exceed the viscous range. The maximal (positive) Lyapunov exponent represents the average strength of the exponential growth rate, while fluctuations in the rate of growth are characterized by the finite-time Lyapunov exponents (FTLEs). In the last decade or so, the notion of Lagrangian coherent structures (which are often computed using FTLEs) has gained attention as a tool for visualizing coherent trajectory patterns in a flow and distinguishing regions of the flow with different mixing properties. A quantitative statistical characterization of FTLEs can be accomplished using the statistical theory of large deviations, based on the so-called Cramér function. To obtain the Cramér function from data, we use both the method based on measuring moments and measuring histograms and introduce a finite-size correction to the histogram-based method. We generalize the existing univariate formalism to the joint distributions of the two FTLEs needed to fully specify the Lyapunov spectrum in 3D flows. The joint Cramér function of turbulence is measured from two direct numerical simulation datasets of isotropic turbulence. Results are compared with joint statistics of FTLEs computed using only the symmetric part of the velocity gradient tensor, as well as with joint statistics of instantaneous strain-rate eigenvalues. When using only the strain contribution of the velocity gradient, the maximal FTLE nearly doubles in magnitude, highlighting the role of rotation in de-correlating the fluid deformations along particle paths. We also extend the large-deviation theory to study the statistics of the ratio of FTLEs. The most likely ratio of the FTLEs λ{sub 1} : λ{sub 2} : λ{sub 3} is shown to be about 4:1:−5, compared to about 8:3:−11 when using only the strain-rate tensor for calculating fluid volume deformations. The results serve to characterize the fundamental statistical and geometric structure of turbulence at small scales including cumulative, time integrated effects. These are important for deformable particles such as droplets and polymers advected by turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Finite-temperature coupled-cluster, many-body perturbation, and restricted and unrestricted Hartree–Fock study on one-dimensional solids: Luttinger liquids, Peierls transitions, and spin- and charge-density waves

One-dimensional (1D) solids exhibit a number of striking electronic structures including charge-density wave (CDW) and spin-density wave (SDW). Also, the Peierls theorem states that at zero temperature, a 1D system predicted by simple band theory to be a metal will spontaneously dimerize and open a finite fundamental bandgap, while at higher temperatures, it will assume the equidistant geometry with zero bandgap (a Peierls transition). We computationally study these unique electronic structures and transition in polyyne and all-trans polyacetylene using finite-temperature generalizations of ab initio spin-unrestricted Hartree–Fock (UHF) and spin-restricted coupled-cluster doubles (CCD) theories, extending upon previous work [He et al., J. Chem. Phys. 140, 024702 (2014)] that is based on spin-restricted Hartree–Fock (RHF) and second-order many-body perturbation (MP2) theories. Unlike RHF, UHF can predict SDW as well as CDW and metallic states, and unlike MP2, CCD does not diverge even if the underlying RHF reference wave function is metallic. UHF predicts a gapped SDW state with no dimerization at low temperatures, which gradually becomes metallic as the temperature is raised. CCD, meanwhile, confirms that electron correlation lowers the Peierls transition temperature. Furthermore, we show that the results from all theories for both polymers are subject to a unified interpretation in terms of the UHF solutions to the Hubbard–Peierls model using different values of the electron-electron interaction strength, U/t, in its Hamiltonian. The CCD wave function is shown to encompass the form of the exact solution of the Tomonaga–Luttinger model and is thus expected to describe accurately the electronic structure of Luttinger liquids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗