SEARCH · Engineering Papers

Results for “singular equations”

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 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Numerical methods for nonlocal and fractional models

Partial differential equations (PDEs) are used, with huge success, to model phenomena arising across all scientific and engineering disciplines. However, across an equally wide swath, there exist situations in which PDE models fail to adequately model observed phenomena or are not the best available model for that purpose. On the other hand, in many situations, nonlocal models that account for interaction occurring at a distance have been shown to more faithfully and effectively model observed phenomena that involve possible singularities and other anomalies. In this article, we consider a generic nonlocal model, beginning with a short review of its definition, the properties of its solution, its mathematical analysis, and specific concrete examples. We then provide extensive discussions about numerical methods, including finite element, finite difference, and spectral methods, for determining approximate solutions of the nonlocal models considered. In that discussion, we pay particular attention to a special class of nonlocal models that are the most widely studied in the literature, namely those involving fractional derivatives. The article ends with brief considerations of several modeling and algorithmic extensions which serve to show the wide applicability of nonlocal modeling.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

New superstrata from three-dimensional supergravity

We find a two-parameter family of generalized superstrata that emerge as smooth, supersymmetric solutions in three-dimensional gauged supergravity coupled to additional scalar fields. This new family of generalized superstrata are smooth microstate geometries and may be thought of as supersymmetric Coulomb-branch extensions of the original superstrata in which the underlying supertube undergoes an elliptical deformation. These solutions had already been obtained numerically, and as series solutions, to the equations of motion, and some of them were conjectured to be supersymmetric. Here we prove the supersymmetry of an entire two-parameter family and we obtain a highly non-trivial analytic and smooth solution for a one-parameter limit in which the global symmetry of the metric is enhanced to SO(3). We also confirm that the other known families of microstrata are not supersymmetric. We conclude with a cursory analysis of some of the singular brane distributions that can be accessed from three-dimensional gauged supergravity while preserving the same supersymmetries as the superstratum, and therefore of the three-charge black hole.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1972

Mode coupling and wave particle interactions for unstable ion acoustic waves

A theory for the spatial development of linearly unstable, coupled waves is presented in which both quasi-linear and mode coupling effects are treated in a self-consistent manner. Steady state excitation of two waves is assumed at the boundary x = 0, the plasma being homogeneous in the y and z directions. Coupled equations are derived for the x dependence of the amplitudes of the primary waves and the secondary waves, correct through second order terms in the wave amplitude, but without usual approximation of small growth rates. This general formalism is then applied to the case of coupled ion acoustic waves driven unstable by an ion beam streaming in the direction of the x axis. If the modifications of the ion beam by the waves (quasi-linear effects) are ignored, explosive instabilities (singularities in all of the amplitudes at finite x) are found, even when all of the waves have positive energy. If these wave-particle interactions are included, the solutions are no longer singular, and all of the amplitudes have finite maxima.

Martin, P.↗

NASA NTRS1992

Formally biorthogonal polynomials and a look-ahead Levinson algorithm for general Toeplitz systems

Systems of linear equations with Toeplitz coefficient matrices arise in many important applications. The classical Levinson algorithm computes solutions of Toeplitz systems with only O(n(sub 2)) arithmetic operations, as compared to O(n(sub 3)) operations that are needed for solving general linear systems. However, the Levinson algorithm in its original form requires that all leading principal submatrices are nonsingular. An extension of the Levinson algorithm to general Toeplitz systems is presented. The algorithm uses look-ahead to skip over exactly singular, as well as ill-conditioned leading submatrices, and, at the same time, it still fully exploits the Toeplitz structure. In our derivation of this algorithm, we make use of the intimate connection of Toeplitz matrices with formally biorthogonal polynomials.

Freund, Roland W.↗

NASA NTRS1997

Aerodynamic Design Optimization on Unstructured Grids with a Continuous Adjoint Formulation

A continuous adjoint approach for obtaining sensitivity derivatives on unstructured grids is developed and analyzed. The derivation of the costate equations is presented, and a second-order accurate discretization method is described. The relationship between the continuous formulation and a discrete formulation is explored for inviscid, as well as for viscous flow. Several limitations in a strict adherence to the continuous approach are uncovered, and an approach that circumvents these difficulties is presented. The issue of grid sensitivities, which do not arise naturally in the continuous formulation, is investigated and is observed to be of importance when dealing with geometric singularities. A method is described for modifying inviscid and viscous meshes during the design cycle to accommodate changes in the surface shape. The accuracy of the sensitivity derivatives is established by comparing with finite-difference gradients and several design examples are presented.

Anderson, W. Kyle↗

NASA NTRS2004

Large-Eddy Simulation of Laminar-Turbulent Transition in a Swept-Wing Boundary Layer

Huai, X.↗

NASA NTRS1972

Mode-coupling and wave-particle interactions for unstable ion-acoustic waves.

A theory for the spatial development of linearly unstable, coupled waves is presented in which both quasilinear and mode-coupling effects are treated in a self-consistent manner. Steady-state excitation of two waves is assumed at the boundary x = 0, the plasma being homogeneous in the y and z directions. Coupled equations are derived for the x dependence of the amplitudes of the primary waves and the secondary waves, correct through terms of second order in the wave amplitude, but without the usual approximation of small growth rates. This general formalism is then applied to the case of coupled ion-acoustic waves driven unstable by an ion beam streaming in the direction of the x axis. If the modifications of the ion beam by the waves (quasilinear effects) are ignored, explosive instabilities (singularities in all of the amplitudes at finite x) are found even when all of the waves have positive energy. If these wave-particle interactions are included, the solutions are no longer singular, and all of the amplitudes have finite maxima.

Martin, P.↗

NASA NTRS1983

Rain estimation from satellites: An examination of the Griffith-Woodley technique

The Griffith-Woodley Technique (GWT) is an approach to estimating precipitation using infrared observations of clouds from geosynchronous satellites. It is examined in three ways: an analysis of the terms in the GWT equations; a case study of infrared imagery portraying convective development over Florida; and the comparison of a simplified equation set and resultant rain map to results using the GWT. The objective is to determine the dominant factors in the calculation of GWT rain estimates. Analysis of a single day's convection over Florida produced a number of significant insights into various terms in the GWT rainfall equations. Due to the definition of clouds by a threshold isotherm the majority of clouds on this day did not go through an idealized life cycle before losing their identity through merger, splitting, etc. As a result, 85% of the clouds had a defined life of 0.5 or 1 h. For these clouds the terms in the GWT which are dependent on cloud life history become essentially constant. The empirically derived ratio of radar echo area to cloud area is given a singular value (0.02) for 43% of the sample, while the rainrate term is 20.7 mmh-1 for 61% of the sample. For 55% of the sampled clouds the temperature weighting term is identically 1.0. Cloud area itself is highly correlated (r=0.88) with GWT computed rain volume. An important, discriminating parameter in the GWT is the temperature defining the coldest 10% cloud area. The analysis further shows that the two dominant parameters in rainfall estimation are the existence of cold cloud and the duration of cloud over a point.

Negri, A. J.↗

NASA NTRS2015

Grid Convergence for Turbulent Flows(Invited)

A detailed grid convergence study has been conducted to establish accurate reference solutions corresponding to the one-equation linear eddy-viscosity Spalart-Allmaras turbulence model for two dimensional turbulent flows around the NACA 0012 airfoil and a flat plate. The study involved three widely used codes, CFL3D (NASA), FUN3D (NASA), and TAU (DLR), and families of uniformly refined structured grids that differ in the grid density patterns. Solutions computed by different codes on different grid families appear to converge to the same continuous limit, but exhibit different convergence characteristics. The grid resolution in the vicinity of geometric singularities, such as a sharp trailing edge, is found to be the major factor affecting accuracy and convergence of discrete solutions, more prominent than differences in discretization schemes and/or grid elements. The results reported for these relatively simple turbulent flows demonstrate that CFL3D, FUN3D, and TAU solutions are very accurate on the finest grids used in the study, but even those grids are not sufficient to conclusively establish an asymptotic convergence order.

Diskin, Boris↗

DOE OSTI2021

Quantum algorithms for open lattice field theory

Certain aspects of some unitary quantum systems are well described by evolution via a non-Hermitian effective Hamiltonian, as in the Wigner-Weisskopf theory for spontaneous decay. Conversely, any non-Hermitian Hamiltonian evolution can be accommodated in a corresponding unitary system + environment model via a generalization of Wigner-Weisskopf theory. This demonstrates the physical relevance of novel features such as exceptional points in quantum dynamics, and opens up avenues for studying many-body systems in the complex plane of coupling constants. In the case of lattice field theory, sparsity lends these channels the promise of efficient simulation on standardized quantum hardware. We thus consider quantum operations that correspond to Suzuki-Lie-Trotter approximation of lattice field theories undergoing nonunitary time evolution, with potential applicability to studies of spin or gauge models at finite chemical potential, with topological terms, to quantum phase transitions—a range of models with sign problems. We develop non-Hermitian quantum circuits and explore their promise on a benchmark, the quantum one-dimensional Ising model with complex longitudinal magnetic field, showing that observables can probe the Lee-Yang edge singularity. The development of attractors past critical points in the space of complex couplings indicates a potential for study on near-term noisy hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1982

Topological analysis of computed three-dimensional viscous flow fields

Computed solutions of the time-dependent, Reynolds-averaged Navier-Stokes equations for three dimensional flows having thin shear layers are analyzed using topological concepts. Specific examples include the transonic flow over a body of revolution with conical afterbody at moderate angles of incidence to the free stream. Experimental flow-visualization techniques are simulated graphically to visualize the computed flow. Scalar and vector fluid dynamics properties such as pressure, shear stress, and vorticity on the body surface are presented as topological maps, and their relationship to one another in terms of orientation and singular points is discussed. The extrapolation from these surface topologies toward the understanding of external flow-field behavior is and demonstrated.

Deiwert, G. S.↗

NASA NTRS1984

Rain estimation from satellites - An examination of the Griffith-Woodley Technique

The Griffith-Woodley Technique (GWT) is an approach to estimating precipitation using infrared observations of clouds from geosynchronous satellites. It is examined in three ways: an analysis of the terms in the GWT equations; a case study of infrared imagery portraying convective development over Florida; and the comparison of a simplified equation set and resultant rain map to results using the GWT. The objective is to determine the dominant factors in the calculation of GWT rain estimates. Analysis of a single day's convection over Florida produced a number of significant insights into various terms in the GWT rainfall equations. Due to the definition of clouds by a threshold isotherm the majority of clouds on this day did not go through an idealized life cycle before losing their identity through merger, splitting, etc. As a result, 85 percent of the clouds had a defined life of 0.5 or 1 h. For these clouds, the terms in the GWT which are dependent on cloud life history become essentially constant. The empirically derived ratio of radar echo area to cloud area is given a singular value (0.02) for 43 percent of the sample, while the rainrate term is 20.7 mmh-1 for 61 percent of the sample. For 55 percent of the sampled clouds, the temperature weighting term is identically 1.0. Cloud area itself is highly correlated (r = 0.88) with GWT computed rain volume. An important, discriminating parameter in the GWT is the temperature defining the coldest 10 percent cloud area. The analysis further shows that the two dominant parameters in rainfall estimation are the existence of cold cloud and the duration of cloud over a point. Previously announced in STAR as N84-13735

Negri, A. J.↗

DOE OSTI2024

Smoothing and flattening the universe through slow contraction versus inflation

In a systematic study, we use an equivalent pair of improved numerical relativity codes based on a tetrad-formulation of the classical Einstein-scalar field equations to examine whether slow contraction or inflation (or both) can resolve the homogeneity, isotropy and flatness problems. Our finding, based on a set of gauge/frame invariant diagnostics and the models considered, is that slow contraction robustly and rapidly smooths and flattens spacetime beginning from initial conditions that are outside the perturbative regime of the flat Friedmann-Robertson-Walker metric, whereas inflation fails these tests. We present new numerical evidence supporting the conjecture that the combination of ultralocal evolution and an effective equation-of-state with pressure much greater than energy density is the key to having robust and rapid smoothing. The opposite of ultralocality occurs in expanding spacetimes, which is the leading obstruction to smoothing following a big bang.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1982

Topological analysis of computed three-dimensional viscous flow fields

Computed solutions of the time-dependent, Reynolds-averaged Navier-Stokes equations for three-dimensional flows having thin shear layers are analyzed, using topological concepts. Specific examples include the transonic flow over a body of revolution with conical afterbody at moderate angles of incidence to the free stream. Experimental flow-visualization techniques are simulated graphically to visualize the computed flow. Scalar and vector fluid dynamic properties, such as pressure, shear stress, and vorticity on the body surface, are presented as topological maps, and their relationship to one another in terms of orientation and singular points is discussed. The extrapolation from these surface topologies toward the understanding of external flow-field behavior is discussed and demonstrated.

Deiwert, G. S.↗

NASA NTRS1985

Taylor-Goertler instabilities of Tollmien-Schlichting waves and other flows governed by the interactive boundary layer equations

The Taylor-Gortler vortex instability equations are formulated for steady and unsteady interacting boundary layer flows of the type which arise in triple-deck theory. The effective Gortler number is shown to be a function of the all shape in the boundary layer and the possibility of both steady and unsteady Taylor-Gortler modes exists. As an example the steady flow in a symmetrically constricted channel is considered and it is shown that unstable Gortler vortices exist before the boundary layers at the wall develop the Goldstein singularity. As an example of an unsteady spatially varying basic state the instability of high frequency large amplitude Tollmien-Schlichting waves in a curved channel were considered. It is shown that they are unstable in the first Stokes layer stage of the hierarchy of nonlinear states. The Tollmien-Schlichting waves are shown to be unstable in the presence of both convex and concave curvature.

Hall, P.↗

DOE OSTI2020

A State-Space Model of an Inverter-Based Microgrid for Multivariable Feedback Control Analysis and Design

In this work, a synchronous model for grid-connected and islanded microgrids is presented. The grid-connected model is based on the premise that the reference frame is synchronized with the AC bus. The quadrature component of the AC bus voltage can be cancelled, which allows to express output power as a linear equation for nominal values in the AC bus amplitude voltage. The model for the islanded microgrid is developed by integrating all the inverter dynamics using a state-space model for the load currents. This model is presented in a comprehensive way such that it could be scalable to any number of inverter-based generators using inductor–capacitor–inductor (LCL) output filters. The use of these models allows designers to assess microgrid stability and robustness using modern control methods such as eigenvalue analysis and singular value diagrams. Both models were tested and validated in an experimental setup to demonstrate their accuracy in describing microgrid dynamics. In addition, three scenarios are presented: non-controlled model, Linear-Quadratic Integrator (LQI) power control, and Power-Voltage (PQ/Vdq) droop–boost controller. Experimental results demonstrate the effectiveness of the control strategies and the accuracy of the models to describe microgrid dynamics.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1984

The collapse of the cores of slowly rotating isothermal clouds

A generalized model which accounts for the effects of initially uniform and slow rotation is defined for the spherical collapse of a singular isothermal sphere such as protosolar and binary nebulae. An initial unstable equilibrium state is described for a sound speed of 0.35 km/sec and a rotation rate of 10 to the -14th/sec for the molecular cloud surrounding the accreting core. The total angular momentum and mass of the inner cloud is set equal to solar system values. The evolution of the collapse is traced by applying a perturbation analysis to the similarity solution for a nonrotating condition, and matched asymptotic expansions solve the hydrodynamic equations. The model is concluded a valid tool for studying star and nebular disk formation.

Terebey, S.↗

NASA NTRS2004

A Study of Supersonic Surface Sources: The Ffowcs Williams-Hawkings Equation and the Kirchhoff Formula

In this paper we address the mathematical problem of noise generation from high speed moving surfaces. The problem we are solving is the linear wave equation with sources on a moving surface. The Ffowcs Williams-Hawkings (FW-H) equation as well as the govern- ing equation for deriving the Kirchhoff formula for moving surfaces are both this type of partial differential equation. We give a new exact solution of this problem here in closed form which is valid for subsonic and supersonic motion of the surface but it is particularly suitable for supersonically moving surfaces. This new solution is the simplest of all high speed formulations of Langley and is denoted formulation 4 following the tradition of numbering of our major results for the prediction of the noise of rotating blades. We show that for a smooth surface moving at supersonic speed, our solution has only removable singularities. Thus it can be used for numerical work.

Farassat, F.↗