SEARCH · Engineering Papers

Results for “Normal form 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

DOE OSTI2021

Independent normalization for γ -ray strength functions: The shape method

Here, the shape method, a novel approach to obtain the functional form of the γ-ray strength function (γSF), is introduced. In connection with the Oslo method the slope of the nuclear level density (NLD) and γSF can be obtained simultaneously even in the absence of neutron resonance spacing data. The foundation of the shape method lies in the primary γ-ray transitions which preserve information on the functional form of the γSF. The shape method has been applied to 56 Fe, 92 Zr, and 164 Dy, which are representative cases for the variety of situations encountered in typical NLD and γSF studies. The comparisons of results from the shape method to those from the Oslo method demonstrate that the functional form of the γSF is retained regardless of nuclear structure details or J π values of the states fed by the primary transitions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1993

Structure and length scales of rotational discontinuities

Results are presented of a series of 1D hybrid (kinetic ions and fluid electrons) simulations, using a nondynamic method to form rotational discontinuities (RDs) at an angle of 60 deg between the normal direction and the upstream magnetic field. Ion kinetic effects are found to be important with regard to the stability of the current layer, its thickness and scaling, and its dependence on initial conditions. A detailed examination of the downstream wave train and gradient scale size showed that RDs with alpha1 = +/-180 deg are stable within the framework of the simulation. With proper initialization, electron-sense rotations of -270 deg were found that do not disintegrate immediately but rather appear to consist of a superposition of a solitary wave with -360-deg rotation and a 90-deg ion-sense rotation.

Krauss-Varban, D.↗

NASA NTRS1987

An L-U implicit multigrid algorithm for the three-dimensional Euler equations

An LU implicit multigrid scheme is developed for the calculation of three-dimensional transonic flow through rotating cascades. This numerical method solves the unsteady Euler equations of gas dynamics in a finite-volume form. The implicit scheme makes it possible to take a much larger time step than is normally permitted in most explicit schemes, while the multigrid method is incorporated to accelerate the convergence rate for steady state calculations. Using this method, computational storage requirements are comparable to those of explicit schemes, while operation counts are considerably less than those found in the more widely-used ADI schemes.

Yokota, Jeffrey W.↗

NASA NTRS1984

Modification to shock fitting program

A modified form of the Lepping - Argentiero single spacecraft, shock normal determination procedure is presented. The modified method incorporates a simple predictor-corrector algorithm which allows a faster convergence rate and the use of average values of the parameters for the starting vector.

Acuna, M. H.↗

NASA NTRS1976

A new method for solving problems of sound radiation from a duct for the case of low Mach number

A method for rapidly solving duct end-plane impedance problems is presented which is based on deriving a form of the Helmholtz integral formula expressing the normal derivative of the acoustic pressure at a boundary point. A set of integral expressions is obtained by satisfying the existing boundary conditions and then solved using a general collocation method. The results obtained for the Levine-Schwinger (1948) problem agree well with the exact values. In the case of a duct flow having a temperature mismatch with the surroundings, the results show an increase in the magnitude of the reflection coefficient for the case of a cold core, and a decrease in its value for the case of a hot core, as compared with the case of exhaust into a uniform medium. A phase change of between pi/2 and pi is nearly always achieved, which indicates the tendency towards maintaining a constant pressure at the exit plane. The method can easily be extended to arbitrary duct shapes as long as they are axisymmetric.

Wahbah, M. M.↗

DOE OSTI2024

Extracting Topological Orders of Generalized Pauli Stabilizer Codes in Two Dimensions

In this paper, we introduce an algorithm for extracting topological data from translation invariant generalized Pauli stabilizer codes in two-dimensional systems, focusing on the analysis of anyon excitations and string operators. The algorithm applies to Z d qudits, including instances where d is a nonprime number. This capability allows the identification of topological orders that differ from the Z d toric codes. It extends our understanding beyond the established theorem that Pauli stabilizer codes for Z p qudits (with p being a prime) are equivalent to finite copies of Z p toric codes and trivial stabilizers. The algorithm is designed to determine all anyons and their string operators, enabling the computation of their fusion rules, topological spins, and braiding statistics. The method converts the identification of topological orders into computational tasks, including Gaussian elimination, the Hermite normal form, and the Smith normal form of truncated Laurent polynomials. Furthermore, the algorithm provides a systematic approach for studying quantum error-correcting codes. We apply it to various codes, such as self-dual CSS quantum codes modified from the two-dimensional honeycomb color code and non-CSS quantum codes that contain the double semion topological order or the six-semion topological order. Published by the American Physical Society 2024

Physics↗

NASA NTRS1982

A study of the high frequency limitations of series resonant converters

A transformer induced oscillation in series resonant (SR) converters is studied. It may occur in the discontinuous current mode. The source of the oscillation is an unexpected resonant circuit formed by normal resonance components in series with the magnetizing inductance of the output transformers. The methods for achieving cyclic stability are: to use a half bridge SR converter where q0.5. Q should be as close to 1.0 as possible. If 0.5q1.0, the instability will be avoided if psi2/3q-1/3. The second objective was to investigate a power field effect transistor (FET) version of the SR converter capable of operating at frequencies above 100 KHz, to study component stress and losses at various frequencies.

Stuart, T. A.↗

NASA NTRS1984

Modification to shock fitting program

A modified form of the Lepping-Argentiero single spacecraft, shock normal determination procedure is presented. The modified method incorporates a simple predictor-corrector algorithm which allows a faster convergence rate and the use of average values of the parameters for the starting vector.

Acuna, M. H.↗

NASA NTRS1997

Improvements in Block-Krylov Ritz Vectors and the Boundary Flexibility Method of Component Synthesis

A method of dynamic substructuring is presented which utilizes a set of static Ritz vectors as a replacement for normal eigenvectors in component mode synthesis. This set of Ritz vectors is generated in a recurrence relationship, proposed by Wilson, which has the form of a block-Krylov subspace. The initial seed to the recurrence algorithm is based upon the boundary flexibility vectors of the component. Improvements have been made in the formulation of the initial seed to the Krylov sequence, through the use of block-filtering. A method to shift the Krylov sequence to create Ritz vectors that will represent the dynamic behavior of the component at target frequencies, the target frequency being determined by the applied forcing functions, has been developed. A method to terminate the Krylov sequence has also been developed. Various orthonormalization schemes have been developed and evaluated, including the Cholesky/QR method. Several auxiliary theorems and proofs which illustrate issues in component mode synthesis and loss of orthogonality in the Krylov sequence have also been presented. The resulting methodology is applicable to both fixed and free- interface boundary components, and results in a general component model appropriate for any type of dynamic analysis. The accuracy is found to be comparable to that of component synthesis based upon normal modes, using fewer generalized coordinates. In addition, the block-Krylov recurrence algorithm is a series of static solutions and so requires significantly less computation than solving the normal eigenspace problem. The requirement for less vectors to form the component, coupled with the lower computational expense of calculating these Ritz vectors, combine to create a method more efficient than traditional component mode synthesis.

Carney, Kelly Scott↗

NASA NTRS1981

Theory of strong turbulence by renormalization

The hydrodynamical equations of turbulent motions are inhomogeneous and nonlinear in their inertia and force terms and will generate a hierarchy. A kinetic method was developed to transform the hydrodynamic equations into a master equation governing the velocity distribution, as a function of the time, the position and the velocity as an independent variable. The master equation presents the advantage of being homogeneous and having fewer nonlinear terms and is therefore simpler for the investigation of closure. After the closure by means of a cascade scaling procedure, the kinetic equation is derived and possesses a memory which represents the nonMarkovian character of turbulence. The kinetic equation is transformed back to the hydrodynamical form to yield an energy balance in the cascade form. Normal and anomalous transports are analyzed. The theory is described for incompressible, compressible and plasma turbulence. Applications of the method to problems relating to sound generation and the propagation of light in a nonfrozen turbulence are considered.

Tchen, C. M.↗

NASA NTRS1983

The Riemann disks. II - Stability

The Riemann disks are the two-dimensional analogs of the Riemann S-type ellipsoids. They form a simple yet rich family of galactic bar models. A method is developed to find their normal modes to arbitrary order, and the stability analysis is carried out for all modes with order not greater than 5.

Weinberg, M. D.↗

DOE OSTI2020

Theoretical Development of Controller Transfer applied to Dynamical Systems

Data-driven control strategies suffer from the necessity of training and its associated computational time and cost. This motivates the continuous increase of popularity of transfer learning based methods. In this paper, we have investigated transfer learning in the context of controller for dynamical systems. First, we have derived conditions under which controller can be transferred between two linear dynamical systems. Furthermore, we have introduced a novel algorithm to design the transferable and non-transferable control components for linear systems. Secondly, we have derived conditions under which a stabilizing controller can be transferred between two different nonlinear dynamical systems. Furthermore, we have used the notion of zero dynamics to transform the nonlinear dynamics into global normal form, to derive the transferability conditions on the controller. Finally, we have numerically evaluated the performance of our transfer learning based methods for two sets of example linear and nonlinear dynamical systems.

Chakraborty, Indrasis↗

NASA NTRS1956

A reevaluaion of data on atmospheric turbulence and airplane gust loads for application in spectral calculations

The available information on the spectrum of atmospheric turbulence is first briefly reviewed. On the basis of these results, methods are developed for the conversion of available gust statistics normally given in terms of counts of gusts or acceleration peaks into a form appropriate for use in spectral calculations. The fundamental quantity for this purpose appears to be the probability distribution of the root-mean-square gust velocity. Estimates of this distribution are derived from data for a number of load histories of transport operations; also, estimates of the variation of this distribution with altitude and weather condition are derived from available data and the method of applying these results to the calculation of airplane gust-response histories in operations is also outlined. (author)

Press, Harry↗

NASA NTRS1991

An approximate viscous shock layer approach to calculating hypersonic flows about blunt-nosed bodies

An approximate axisymmetric method has been developed which can reliably calculate fully viscous hypersonic flows over blunt-nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. Since the method is fully viscous, the problems associated with coupling a boundary-layer solution with an inviscid-layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed. Surface heat transfer and pressure predictions are comparable to both VSL results and experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher-order procedures which require starting solutions.

Cheatwood, F. MCN.↗

NASA NTRS1976

Fully unsteady subsonic and supersonic potential aerodynamics of complex aircraft configurations for flutter applications

The paper presents a general theory of steady, oscillatory or fully unsteady, subsonic and supersonic aerodynamics for aircraft of arbitrary shape. The finite element method is used to obtain the time functional relationship between aerodynamic potential and its normal derivative (normal wash) in a form suitable for computational analysis. The matrix of the aerodynamic influence coefficients, as necessary for flutter calculations, is obtained as a result. The theory is implemented in a computer program, and some typical results obtained with this program are presented.

Tseng, K.↗

NASA NTRS1992

The recovery of asteroids after two observations

It is shown that a generalization of the use of 'Vaisala orbits', briefly mentioned at the Asteroids 2 Conference, can be very conveniently accomplished by means of an inversion of the 'GEM' form of the Gauss method. The procedure can also be applied to Apollo objects and to indeterminate cases of normal three-observation orbit computation, and there is also a simple extension to situations involving four or more observations.

Marsden, Brian G.↗

NASA NTRS1992

Kinetic structure of slow shocks - Effects of the electromagnetic ion/ion cyclotron instability

The structure of slow magnetosonic shocks in the low beta regime is analyzed with attention given to ion heating and the effects of waves upstream of the electromagnetic ion/ion cyclotron (EMIIC) instability. Shock formation is assessed by means of three methods - a relaxation method and two based on dynamic flow interactions - to determine the effects of initialization and boundary conditions on the formation. Good solutions are found with the piston method and the similar flow-flow method in which the plasma is injected from two boundaries to form two slow shocks. Plasma parameters and shock normal angle are found to be the key variables dictating the structure of the magnetosonic shocks. Four unique classes of resultant shock structures are described in which classical, steady, or nonsteady behavior is found. The analysis also yields insight into the relationship between EMIIC instability and ion dissipation.

Omidi, N.↗