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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

QCD Predictions for Physical Multimeson Scattering Amplitudes

We use lattice QCD calculations of the finite-volume spectra of systems of two and three mesons to determine, for the first time, three-particle scattering amplitudes with physical quark masses. Our results are for combinations of 𝜋 + and 𝐾 + , at a lattice spacing 𝑎 = 0.063 fm, and in the isospin-symmetric limit. We also obtain accurate results for maximal-isospin two-meson amplitudes, with those for 𝜋 + ⁢𝐾 + and 2⁢𝐾 + being the first determinations at the physical point. Dense lattice spectra are obtained using the stochastic Laplacian-Heaviside method, and the analysis leading to scattering amplitudes is done using the relativistic finite-volume formalism. Results are compared to chiral perturbation theory and to phenomenological fits to experimental data, finding good agreement.

hadron-hadron interactions↗
NASA NTRS1984

On the instability of a 3-dimensional attachment line boundary layer: Weakly nonlinear theory and a numerical approach

The instability of a three dimensional attachment line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time dependent Navier-Stokes equations for the attachment line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment line boundary layer is also investigated.

Hall, P.↗
NASA NTRS1989

Exploiting symmetries in the modeling and analysis of tires

A computational procedure is presented for reducing the size of the analysis models of tires having unsymmetric material, geometry and/or loading. The two key elements of the procedure when applied to anisotropic tires are: (1) decomposition of the stiffness matrix into the sum of an orthotropic and nonorthotropic parts; and (2) successive application of the finite-element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The proposed technique has high potential for handling practical tire problems with anisotropic materials, unsymmetric imperfections and asymmetric loading. It is also particularly useful for use with three-dimensional finite-element models of tires.

Noor, Ahmed K.↗
DOE OSTI2021

Quantum Simulation of Open Quantum Systems Using a Unitary Decomposition of Operators

Electron transport in realistic physical and chemical systems often involves the nontrivial exchange of energy with a large environment, requiring the definition and treatment of open quantum systems. Because the time evolution of an open quantum system employs a nonunitary operator, the simulation of open quantum systems presents a challenge for universal quantum computers constructed from only unitary operators or gates. Here, we present a general algorithm for implementing the action of any nonunitary operator on an arbitrary state on a quantum device. We show that any quantum operator can be exactly decomposed as a linear combination of at most four unitary operators. We demonstrate this method on a two-level system in both zero and finite temperature amplitude damping channels. Furthermore, the results are in agreement with classical calculations, showing promise in simulating nonunitary operations on intermediate-term and future quantum devices.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1986

On the instability of a three-dimensional attachment-line boundary layer - Weakly nonlinear theory and a numerical approach

The instability of a three-dimensional attachment-line boundary layer is considered in the nonlinear regime. Using weakly nonlinear theory, it is found that, apart from a small interval near the (linear) critical Reynolds number, finite-amplitude solutions bifurcate subcritically from the upper branch of the neutral curve. The time-dependent Navier-Stokes equations for the attachment-line flow have been solved using a Fourier-Chebyshev spectral method and the subcritical instability is found at wavenumbers that correspond to the upper branch. Both the theory and the numerical calculations show the existence of supercritical finite-amplitude (equilibrium) states near the lower branch which explains why the observed flow exhibits a preference for the lower branch modes. The effect of blowing and suction on nonlinear stability of the attachment-line boundary layer is also investigated.

Hall, P.↗
NASA NTRS1999

Drop Ejection From an Oscillating Rod

The dynamics of a drop of a Newtonian liquid that is pendant from or sessile on a solid rod that is forced to undergo time-periodic oscillations along its axis is studied theoretically. The free boundary problem governing the time evolution of the shape of the drop and the flow field inside it is solved by a method of lines using a finite element algorithm incorporating an adaptive mesh. When the forcing amplitude is small, the drop approaches a limit cycle at large times and undergoes steady oscillations thereafter. However, drop breakup is the consequence if the forcing amplitude exceeds a critical value. Over a wide range of amplitudes above this critical value, drop ejection from the rod occurs during the second oscillation period from the commencement of rod motion. Remarkably, the shape of the interface at breakup and the volume of the primary drop formed are insensitive to changes in forcing amplitude. The interface shape at times close to and at breakup is a multi-valued function of distance measured along the rod axis and hence cannot be described by recently popularized one-dimensional approximations. The computations show that drop ejection occurs without the formation of a long neck. Therefore, this method of drop formation holds promise of preventing formation of undesirable satellite droplets.

Wilkes, E. D.↗
NASA NTRS1976

Numerical simulation of low gravity draining

A boundary value problem was solved numerically for a liquid that is assumed to be inviscid and incompressible, having a motion that is irrotational and axisymmetric, and having a constant (5 degrees) solid-liquid contact angle. The avoidance of excessive mesh distortion, encountered with strictly Lagrangian or Eulerian kinematics, was achieved by introducing an auxiliary kinematic velocity field along the free surface in order to vary the trajectories used in integrating the ordinary differential equations simulating the moving boundary. The computation of the velocity potential was based upon a nonuniform triangular mesh which was automatically revised to varying depths to accommodate the motion of the free surface. These methods permitted calculation of draining induced axisymmetric slosh through the many (or fractional) finite amplitude oscillations that can occur depending upon the balance of draining, gravitational, and surface tension forces. Velocity fields, evolution of the free surface with time, and liquid residual volumes were computed for three and one half decades of Weber number and for two Bond numbers, tank fill levels, and drain radii. Comparisons with experimental data are very satisfactory.

Bizzell, G. D.↗
NASA NTRS1989

Numerical Study of Three-dimensional Spatial Instability of a Supersonic Flat Plate Boundary Layer

The behavior of spatially growing three-dimensional waves in a supersonic boundary layer was studied numerically by solving the complete Navier-Stokes equations. Satisfactory comparison with linear parallel and non-parallel stability theories, and experiment are obtained when a small amplitude inflow disturbance is used. The three-dimensional unsteady Navier-Stokes equations are solved by a finite difference method which is fourth-order and second-order accurate in the convection and viscous terms respectively, and second-order accurate in time. Spanwise periodicity is assumed. The inflow disturbance is composed of eigenfunctions from linear stability theory. By increasing the amplitude of the inflow disturbance, nonlinear effects in the form of a relaxation type oscillation of the time signal of rho(u) are observed.

Maestrello, Lucio↗
NASA NTRS1997

Wavelet Methods Developed to Detect and Control Compressor Stall

A "wavelet" is, by definition, an amplitude-varying, short waveform with a finite bandwidth (e.g., that shown in the first two graphs). Naturally, wavelets are more effective than the sinusoids of Fourier analysis for matching and reconstructing signal features. In wavelet transformation and inversion, all transient or periodic data features (as in compressor-inlet pressures) can be detected and reconstructed by stretching or contracting a single wavelet to generate the matching building blocks. Consequently, wavelet analysis provides many flexible and effective ways to reduce noise and extract signals which surpass classical techniques - making it very attractive for data analysis, modeling, and active control of stall and surge in high-speed turbojet compressors. Therefore, fast and practical wavelet methods are being developed in-house at the NASA Lewis Research Center to assist in these tasks. This includes establishing user-friendly links between some fundamental wavelet analysis ideas and the classical theories (or practices) of system identification, data analysis, and processing.

Le, Dzu K.↗
NASA NTRS1971

Focusing of Finite-Amplitude Cylindrical and Spherical Sound Waves in a Viscous and Heat-Conducting Medium

The focusing of acoustic pulses is studied analytically by considering the region of study in three parts: the converging, interaction and diverging regions. First, the linear problem of a pulse of infinitesimal amplitude is studied. For the spherical case, the expected phase change as a result of focusing is verified. The nonlinear case of finite-amplitude pulses leads to the development of M-waves, as determined by applying the method of matched-asymptotic expansions to Burges equation.

Chu, T.↗
NASA NTRS1985

Model-size reduction technique for the analysis of symmetric anisotropic structures

A two-step computational procedure is presented for reducing the size of the analysis model for an anisotropic symmetric structure to that of the corresponding orthotropic structure. The key elements of the procedure are: (1) decomposition of the stiffness matrix into the sum of an orthotropic and nonorthotropic (anisotropic) parts; and (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The global approximation vectors are selected to be the solution corresponding to zero nonorthotropic matrix and its various-order derivatives with respect to an anisotropic tracing parameter (identifying the nonorthotropic material coefficients). The size of the analysis model used in generating the global approximation vectors is identical to that of the corresponding orthotropic structure. The effectiveness of the proposed technique is demonstrated by means of numerical examples and its potential for solving other quasi-symmetric problems is discussed.

Noor, A. K.↗
NASA NTRS1995

Numerical thermal analyses of heat exchangers for the stirling engine application

The Regenerator, Cooler and Heater for the NASA Space Power Research Engine (SPRE) have been analyzed in detail for laminar, incompressible and oscillatory flow conditions. Each component has been analyzed independently and in detail with the regenerator being modeled as two-parallel-plates channel with a solid wall. The ends of the channel are exposed to two reservoir maintained at different temperature thus providing an axial temperature gradient along the channel. The cooler and heater components have been modeled as circular pipes with isothermal walls. Two different types of thermal boundary conditions have been investigated for the cooler and heater, namely, symmetric and asymmetric temperature inflow. In symmetric temperature inflow the flow enters the channel with the same temperature in throughout the velocity cycle whereas, in asymmetric temperature inflow the flow enters with a different temperature in each half cycle. The study was conducted over a wide range of Maximum Reynolds number (RE(max) varying from 75 to 60000, Valensi number (Va) from 2.5 to 800, and relative amplitude of fluid displacement (A(sub r) from 0.357 to 1.34. A two dimensional Finite volume method based on the SIMPLE algorithm was used to solve the governing partial differential equations. Post processing programs were developed to effectively describe the heat transfer mechanism under oscillatory flows. The computer code was validated by comparing with existing analytical solutions for oscillating flows. The thermal field have been studied with the help of temperature contour and three dimensional plots. The instantaneous friction factor, wall heat flux and heat transfer coefficient have been examined. It has been concluded that in general, the frictional factor and heat transfer coefficient are higher under oscillatory flow conditions when the Valensi number is high. Also, the thermal efficiency decreases for lower A(r) values. Further, the usual steady state definition for the heat transfer coefficient does not seem to be valid.

Kannapareddy, Mohan Raj↗
NASA NTRS1977

Computation of unsteady transonic flows by the indicial method

The indicial method is investigated for the computation of unsteady transonic force and moment coefficients for use in flutter analyses. This approach has the advantage that solutions for all reduced frequencies for a given mode of motion can be obtained from a single finite-difference flowfield computation. Comparisons of indicial and time-integration computations for oscillating airfoil and flap motions help define limits on the motion amplitude for the applicability of the indicial method to transonic flows. Within these limits, solutions for various motion modes can be superposed to obtain solutions for multiple-degree-of-freedom aeroelastic systems. Also, a simple aeroelastic problem is solved by an alternative approach in which the structural motion and flowfield equations are integrated simultaneously using a time-integration finite-difference procedure.

Ballhaus, W. F.↗
NASA NTRS1987

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗
NASA NTRS1987

Multiple-mode nonlinear free and forced vibrations of beams using finite element method

Effects of large deflection geometric nonlinearity and multiple modes for free and forced vibrations of beams are investigated. Inplane displacement and inplane inertia are included in the formulation. The finite element method is employed. A harmonic force matrix is derived for forced vibration analysis. The relations of amplitude, frequency ratio and forcing intensity for beams of different boundary conditions and various load conditions are presented.

Mei, Chuh↗
NASA NTRS1989

A technique for representinbg three-dimensional vertical circulation in baroclinic disturbances

A method is presented for producing three-dimensional representations of vertical circulations in frontal zones in cases where the circulations are not confined to the cross-front vertical plane. Various approaches for representing vertical circulations are reviewed. The numerical aspects of the channel-model simulation are summarized. A vector streamfunction, the psi vector, that describes the vertical velocity and the horizontal irrotational flow, is defined. The technique for determining the psi vector for a channel geometry is outlined. Results are presented from applying the psi vector method to four frontal systems in an f-plane primitive equation channel model of a finite-amplitude baroclinic wave.

Keyser, Daniel↗
NASA NTRS1990

Satellite formation in capillary jet breakup

Computations of finite-amplitude, spatially periodic wave growth on a cylindrical jet have been carried out using a boundary integral method. The initial wave growth is in agreement with Rayleigh's linear theory. When followed to completion these waves pinch off large drops separated by smaller satellite drops (spherules) that decrease in size with decreasing wavelength. The computed sizes of both drops and satellites agree with experiment. It is found that satellites will form for all unstable wave numbers. The small satellites that are computed at wave numbers near the critical wave number were not predicted by near-linear analysis but are observed in experimental photographs of jet breakup. Computation of the collapse of elongated satellites shows short waves propagating on their surfaces.

Mansour, Nagi N.↗