Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite amplitude method”

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 199 records · Page 11

Hamiltonian Truncation Effective Theory

Hamiltonian truncation is a non-perturbative numerical method for calculating observables of a quantum field theory. The starting point for this method is to truncate the interacting Hamiltonian to a finite-dimensional space of states spanned by the eigenvectors of the free Hamiltonian H_0 H 0 with eigenvalues below some energy cutoff E_\text{max} E max . In this work, we show how to treat Hamiltonian truncation systematically using effective field theory methodology. We define the finite-dimensional effective Hamiltonian by integrating out the states above E_\text{max} E max . The effective Hamiltonian can be computed by matching a transition amplitude to the full theory, and gives corrections order by order as an expansion in powers of 1/E_\text{max} 1 / E max . The effective Hamiltonian is non-local, with the non-locality controlled in an expansion in powers of H_0/E_\text{max} H 0 / E max . The effective Hamiltonian is also non-Hermitian, and we discuss whether this is a necessary feature or an artifact of our definition. We apply our formalism to 2D \lambda\phi^4 λ ϕ 4 theory, and compute the the leading 1/E_\text{max}^2 1 / E max 2 corrections to the effective Hamiltonian. We show that these corrections nontrivially satisfy the crucial property of separation of scales. Numerical diagonalization of the effective Hamiltonian gives residual errors of order 1/E_\text{max}^3 1 / E max 3 , as expected by our power counting. We also present the power counting for 3D \lambda \phi^4 λ ϕ 4 theory and perform calculations that demonstrate the separation of scales in this theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interaction between Tollmien-Schlichting waves and free-stream disturbances in boundary-layer flows

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small freestream disturbance and a small localized variation of the surface geometry has been carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the freestream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the Tollmien-Schlichting waves.

Bodonyi, R. J.↗

Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit

The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.

Balsa, Thomas F.↗

Surface deformation and geoid anomalies over single and double-layered convective systems

Using a primitive variable formulation of the finite-element method, the differences in the surface observables, such as topography and geoid, produced by single- and double-layered thermal convection, were compared. Both constant and depth-dependent viscosities have been considered. For the same Rayleigh number, larger surface perturbations are produced by single-cell convection. For the same Nusselt number, the magnitudes of the surface observables are greater for double-layered convection. For the same surface heat-flux, surface topographies have similar magnitudes, but the relative amplitudes of geoid anomalies depend greatly on the style of viscosity stratification. This difference in the geoid between the two systems increases with greater surface heat-flow, regardless of viscosity structure.

Koch, M.↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

Boundary-layer receptivity due to a wall suction and control of Tollmien-Schlichting waves

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized suction slot on an otherwise flat surface was carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of T-S waves generated by the interaction between the free-stream disturbance and the suction slot, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves and the demonstration of the possible active control of the growth of T-S waves.

Bodonyi, R. J.↗

Neural Network Prediction of Failure of Damaged Composite Pressure Vessels from Strain Field Data Acquired by a Computer Vision Method

This effort used a new and novel method of acquiring strains called Sub-pixel Digital Video Image Correlation (SDVIC) on impact damaged Kevlar/epoxy filament wound pressure vessels during a proof test. To predict the burst pressure, the hoop strain field distribution around the impact location from three vessels was used to train a neural network. The network was then tested on additional pressure vessels. Several variations on the network were tried. The best results were obtained using a single hidden layer. SDVIC is a fill-field non-contact computer vision technique which provides in-plane deformation and strain data over a load differential. This method was used to determine hoop and axial displacements, hoop and axial linear strains, the in-plane shear strains and rotations in the regions surrounding impact sites in filament wound pressure vessels (FWPV) during proof loading by internal pressurization. The relationship between these deformation measurement values and the remaining life of the pressure vessels, however, requires a complex theoretical model or numerical simulation. Both of these techniques are time consuming and complicated. Previous results using neural network methods had been successful in predicting the burst pressure for graphite/epoxy pressure vessels based upon acoustic emission (AE) measurements in similar tests. The neural network associates the character of the AE amplitude distribution, which depends upon the extent of impact damage, with the burst pressure. Similarly, higher amounts of impact damage are theorized to cause a higher amount of strain concentration in the damage effected zone at a given pressure and result in lower burst pressures. This relationship suggests that a neural network might be able to find an empirical relationship between the SDVIC strain field data and the burst pressure, analogous to the AE method, with greater speed and simplicity than theoretical or finite element modeling. The process of testing SDVIC neural network analysis and some encouraging preliminary results are presented in this paper. Details are given concerning the processing of SDVIC output data such that it may be used as back propagation neural network (BPNN) input data. The software written to perform this processing and the BPNN algorithm are also discussed. It will be shown that, with limited training, test results indicate an average error in burst pressure prediction of approximately six percent,

Russell, Samuel S.↗

Normalization of ZZ instanton amplitudes in type 0B minimal superstring theory

Abstract We study ZZ instanton corrections in the (2,4k)$$ \mathcal{N} $$ N = 1 minimal superstring theory with the type 0B GSO projection, which becomes the type 0B$$ \mathcal{N} $$ N = 1 super-JT gravity in thek→ ∞ limit. Each member of the (2,4k) family of theories has two phases distinguished by the sign of the Liouville bulk cosmological constant. The worldsheet method for computing the one-loop normalization constant multiplying the instanton corrections gives an ill-defined answer in both phases. We fix these divergences using insights from string field theory and find finite, unambiguous results. Each member of the (2,4k) family of theories is dual to a double-scaled one-matrix integral, where the double-scaling limit can be obtained starting either from a unitary matrix integral with a leading one-cut saddle point, or from a hermitian matrix integral with a leading two-cut saddle point. The matrix integral exhibits a gap-closing transition, which is the same as the double-scaled Gross-Witten-Wadia transition whenk= 1. We also compute instanton corrections in the double-scaled matrix integral for allkand in both phases, and find perfect agreement with the string theory results.

Physics↗

Experimental realization of an extended Fermi-Hubbard model using a 2D lattice of dopant-based quantum dots

The Hubbard model is an essential tool for understanding many-body physics in condensed matter systems. Artificial lattices of dopants in silicon are a promising method for the analog quantum simulation of extended Fermi-Hubbard Hamiltonians in the strong interaction regime. However, complex atom-based device fabrication requirements have meant emulating a tunable two-dimensional Fermi-Hubbard Hamiltonian in silicon has not been achieved. Here, we fabricate 3 × 3 arrays of single/few-dopant quantum dots with finite disorder and demonstrate tuning of the electron ensemble using gates and probe the many-body states using quantum transport measurements. By controlling the lattice constants, we tune the hopping amplitude and long-range interactions and observe the finite-size analogue of a transition from metallic to Mott insulating behavior. We simulate thermally activated hopping and Hubbard band formation using increased temperatures. As atomically precise fabrication continues to improve, these results enable a new class of engineered artificial lattices to simulate interactive fermionic models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A numerical study of the interaction between unsteady free-stream disturbances and localized variations in surface geometry

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite-difference and spectral methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves.

Bodonyi, R. J.↗

Three-dimensional elastic-plastic finite-element analysis of fatigue crack propagation

Fatigue cracks are a major problem in designing structures subjected to cyclic loading. Cracks frequently occur in structures such as aircraft and spacecraft. The inspection intervals of many aircraft structures are based on crack-propagation lives. Therefore, improved prediction of propagation lives under flight-load conditions (variable-amplitude loading) are needed to provide more realistic design criteria for these structures. The main thrust was to develop a three-dimensional, nonlinear, elastic-plastic, finite element program capable of extending a crack and changing boundary conditions for the model under consideration. The finite-element model is composed of 8-noded (linear-strain) isoparametric elements. In the analysis, the material is assumed to be elastic-perfectly plastic. The cycle stress-strain curve for the material is shown Zienkiewicz's initial-stress method, von Mises's yield criterion, and Drucker's normality condition under small-strain assumptions are used to account for plasticity. The three-dimensional analysis is capable of extending the crack and changing boundary conditions under cyclic loading.

Goglia, G. L.↗

Wave-interactions in three-dimensional boundary layers

A Fourier-Chebyshev spectral method is used to study the nonlinear stability and disturbance wave-interactions in a two-dimensional boundary layer on a curved wall and a three-dimensional boundary layer on a rotating disk. Goertler and crossflow instability, respectively, are the dominant mechanisms in these boundary layers. For the range of parameters considered, the present computational model does not support the Goertler/TS interaction of the type proposed by Nayfeh (1981). According to his model, oblique TS waves with spanwise wavelengths twice that of the finite-amplitude Goertler vortex are excited in the boundary layer. The present nonlinear calculations show that both the Goertler and crosssflow vortices develop strong instabilities in their harmonics with half the wavelength of the primary vortices.

Malik, M. R.↗

Numerical modeling of the vortex/airfoil interaction

A modeling of the vortex-airfoil interaction is presented in which the finite-area of the real vortices is taken into consideration. Two vortex models are used. In the first, a disturbed piece of vorticity layer is simulated by four rows of discrete vortices of small strength. In the second, a number of discrete vortices is arranged within a circle. The first model may simulate a shear layer or a wake, while the second, a well-formed vortex. The method was applied to the calculation of the pressure induced on the surface of the airfoil by the interacting vortex. Both models give similar results. It was found that for large distances of the vortex from the surface of the airfoil, the consideration or not of the finite-area of the vortex is not a significant factor in determining the induced pressure field. However, when the distance of the vortex from the surface is reduced, its shape is distorted and the induced pressure pulses have lower amplitude than the ones induced by an equivalent point vortex. In the limit, where the vortex impinges on the leading edge of the airfoil, it is split into two and the time dependent pressure coefficient takes even negative values at some time intervals.

Panaras, Argyris G.↗

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗

Comparison of Methods to Predict Lower Bound Buckling Loads of Cylinders Under Axial Compression

Results from a numerical study of the buckling response of two different orthogrid stiffened circular cylindrical shells with initial imperfections and subjected to axial compression are used to compare three different lower bound buckling load prediction techniques. These lower bound prediction techniques assume different imperfection types and include an imperfection based on a mode shape from an eigenvalue analysis, an imperfection caused by a lateral perturbation load, and an imperfection in the shape of a single stress-free dimple. The STAGS finite element code is used for the analyses. Responses of the cylinders for ranges of imperfection amplitudes are considered, and the effect of each imperfection is compared to the response of a geometrically perfect cylinder. Similar behavior was observed for shells that include a lateral perturbation load and a single dimple imperfection, and the results indicate that the predicted lower bounds are much less conservative than the corresponding results for the cylinders with the mode shape imperfection considered herein. In addition, the lateral perturbation technique and the single dimple imperfection produce response characteristics that are physically meaningful and can be validated via testing.

Haynie, Waddy T.↗

Linearized frequency domain Landau-Lifshitz-Gilbert equation formulation

We present a general finite element linearized Landau-Lifshitz-Gilbert equation (LLGE) solver for magnetic systems under weak time-harmonic excitation field. The linearized LLGE is obtained by assuming a small deviation around the equilibrium state of the magnetic system. Inserting such expansion into LLGE and keeping only first order terms gives the linearized LLGE, which gives a frequency domain solution for the complex magnetization amplitudes under an external time-harmonic applied field of a given frequency. We solve the linear system with an iterative solver using generalized minimal residual method. We construct a preconditioner matrix to effectively solve the linear system. The validity, effectiveness, speed, and scalability of the linear solver are demonstrated via numerical examples.

36 MATERIALS SCIENCE↗

Computational Study of Near-limit Propagation of Detonation in Hydrogen-air Mixtures

A computational investigation of the near-limit propagation of detonation in lean and rich hydrogen-air mixtures is presented. The calculations were carried out over an equivalence ratio range of 0.4 to 5.0, pressures ranging from 0.2 bar to 1.0 bar and ambient initial temperature. The computations involved solution of the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing (TVD) scheme, and a point implicit, first-order-accurate, time marching algorithm. The hydrogen-air combustion was modeled with a 9-species, 19-step reaction mechanism. A multi-level, dynamically adaptive grid was utilized in order to resolve the structure of the detonation. The results of the computations indicate that when hydrogen concentrations are reduced below certain levels, the detonation wave switches from a high-frequency, low amplitude oscillation mode to a low frequency mode exhibiting large fluctuations in the detonation wave speed; that is, a 'galloping' propagation mode is established.

Yungster, S.↗

Numerical Simulations of Buoyancy Effects in low Density Gas Jets

This paper deals with the computational analysis of buoyancy effects in the near field of an isothermal helium jet injected into quiescent ambient air environment. The transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum were solved using a staggered grid finite volume method. Laminar, axisymmetric, unsteady flow conditions were considered for the analysis. An orthogonal system with non-uniform grids was used to capture the instability phenomena. Computations were performed for Earth gravity and during transition from Earth to different gravitational levels. The flow physics was described by simultaneous visualizations of velocity and concentration fields at Earth and microgravity conditions. Computed results were validated by comparing with experimental data substantiating that buoyancy induced global flow oscillations present in Earth gravity are absent in microgravity. The dependence of oscillation frequency and amplitude on gravitational forcing was presented to further quantify the buoyancy effects.

Satti, R. P.↗