Engineering Papers⌕ Search

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.

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

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

GenASiS Mathematics: Object-oriented manifolds, operations, and solvers for large-scale physics simulations (version 2)

We report GenASiS Mathematics provides modern Fortran classes furnishing extensible object-oriented functionality for the solution of fields governed by selected partial differential equations. The initial release included extensible object-oriented implementations of simple meshes and the evolution of generic conserved currents thereon. This revision - Version 2 of Mathematics - includes significant reorganization and streamlining of these classes, higher-order reconstruction by a different method, a Poisson solver, coarsening to avoid Courant time step limitations near coordinate singularities, and the offloading of computational kernels to GPUs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Role of pressure diffusion in non-homogeneous shear flows

A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.

Demuren, A. O.↗
NASA NTRS1974

Vehicle control for fuel optimization

The problem of predicting the minimum fuel trajectory for a six degree of freedom vehicle which has a motion characterized by the first order differential equations of translational and rotational dynamics is considered. The thrust direction and center of gravity of the vehicle are assumed to be fixed with respect to the vehicle. Thrust magnitude and the control moment are used as control variables and appear linearly in the equation of motion. Pontryagin's Maximum Principle is used to solve the variational problem. With this formulation, the extremal controls are bang-bang with the exception of the singular case.

Dannenberg, K. D.↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory by the method of potential gradient

A generalized solution of the hyperbolic wave equation has been derived. The method used has been further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals have been evaluated in closed form, while numerical integration methods are suggested for treating more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient because line integrations along the other two sides need only be considered. Furthermore, all the integrands vanish on the Mach cone and the need for determining the hyperbolic curves of intersection of the cone with the lifting surface is avoided. Generalized aerodynamic coefficients for three AGARD planforms have been calculated and compared with the available results.

Jones, W. P.↗
NASA NTRS1983

Solution of Poisson equations for 3-dimensional grid generations

A method for generating three dimensional, finite difference grids about complicated geometries by using Poisson equations is developed. The inhomogenous terms are automatically chosen such that orthogonality and spacing restrictions at the body surface are satisfied. Spherical variables are used to avoid the axis singularity, and an alternating-direction-implicit (ADI) solution scheme is used to accelerate the computations. Computed results are presented that show the capability of the method. Since most of the results presented have been used as grids for flow-field computations, this is indicative that the method is a useful tool for generating three-dimensional grids about complicated geometries.

Fujii, K.↗
NASA NTRS1992

Use of the maximum entropy method to retrieve the vertical atmospheric ozone profile and predict atmospheric ozone content

A method is developed for prediction of ozone levels in planetary atmospheres. This method is formulated in terms of error covariance matrices, and is associated with both direct measurements, a priori first guess profiles, and a weighting function matrix. This is described by the following linearized equation: y = A(matrix) x X + eta, where A is the weighting matrix and eta is noise. The problems to this approach are: (1) the A matrix is near singularity; (2) the number of unknowns in the profile exceeds the number of data points, therefore, the solution may not be unique; and (3) even if a unique solution exists, eta may cause the solution to be ill conditioned.

Turner, B. Curtis↗
NASA NTRS2001

Finite Element Methods for real-time Haptic Feedback of Soft-Tissue Models in Virtual Reality Simulators

We have applied the linear elastic finite element method to compute haptic force feedback and domain deformations of soft tissue models for use in virtual reality simulators. Our results show that, for virtual object models of high-resolution 3D data (>10,000 nodes), haptic real time computations (>500 Hz) are not currently possible using traditional methods. Current research efforts are focused in the following areas: 1) efficient implementation of fully adaptive multi-resolution methods and 2) multi-resolution methods with specialized basis functions to capture the singularity at the haptic interface (point loading). To achieve real time computations, we propose parallel processing of a Jacobi preconditioned conjugate gradient method applied to a reduced system of equations resulting from surface domain decomposition. This can effectively be achieved using reconfigurable computing systems such as field programmable gate arrays (FPGA), thereby providing a flexible solution that allows for new FPGA implementations as improved algorithms become available. The resulting soft tissue simulation system would meet NASA Virtual Glovebox requirements and, at the same time, provide a generalized simulation engine for any immersive environment application, such as biomedical/surgical procedures or interactive scientific applications.

Frank, Andreas O.↗
NASA NTRS2004

Hamilton's Equations with Euler Parameters for Rigid Body Dynamics Modeling

A combination of Euler parameter kinematics and Hamiltonian mechanics provides a rigid body dynamics model well suited for use in strongly nonlinear problems involving arbitrarily large rotations. The model is unconstrained, free of singularities, includes a general potential energy function and a minimum set of momentum variables, and takes an explicit state space form convenient for numerical implementation. The general formulation may be specialized to address particular applications, as illustrated in several three dimensional example problems.

Shivarama, Ravishankar↗
DOE OSTI2023

Bound-preserving finite element approximations of the Keller–Segel equations

We report this paper aims to develop numerical approximations of the Keller–Segel equations that mimic at the discrete level the lower bounds and the energy law of the continuous problem. We solve these equations for two unknowns: the organism (or cell) density, which is a positive variable, and the chemoattractant density, which is a non-negative variable. We propose two algorithms, which combine a stabilized finite element method and a semi-implicit time integration. The stabilization consists of a nonlinear artificial diffusion that employs a graph-Laplacian operator and a shock detector that localizes local extrema. As a result, both algorithms turn out to be nonlinear and can generate cell and chemoattractant numerical densities fulfilling lower bounds. However, the first algorithm requires a suitable constraint between the space and time discrete parameters, whereas the second one does not. We design the latter to attain a discrete energy law on acute meshes. We report some numerical experiments to validate the theoretical results on blowup and nonblowup phenomena. In the blowup setting, we identify a locking phenomenon that relates the L ∞ (Ω)-norm to the L 1 (Ω)-norm limiting the growth of the singularity when supported on a macroelement.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An ideal MHD δW stability analysis that bypasses the Newcomb equation

In the following work, we demonstrate the efficacy of a Riccati δW ideal MHD stability analysis that bypasses the numerically intractable integration of the Newcomb equation. By transforming the linear Newcomb equation into a quadratic Riccati equation, an accurate and equivalent δW analysis is performed that is shown to enjoy some numerical advantages. Additionally, we demonstrate that the Riccati approach is better conditioned than its Newcomb counterpart at the magnetic axis, and we apply dynamical systems insights to examine its behavior at singular points. We further discuss the constraints involved in solving for the extremal admissible perturbations of a δW analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Stresses in bonded materials with a crack perpendicular to the interface.

The problem of two elastic bonded half planes containing a crack perpendicular to the interface is considered. First the solution for the semi-infinite crack under concentrated wedge loading is given. Then the problem of a finite crack fully imbedded in one of the half planes of terminating at the interface is considered. The Mellin transform in conjunction with the dislocations is used to formulate the problem and to derive the integral equation. The integral equation is solved; the stress intensity factors, the crack surface displacement, and the stresses around the crack tip terminating at the interface are obtained. Noting that the power of the singularity for the interface tip of the crack is not -1/2, a tentative fracture criterion dealing with the initiation of fracture propagation is proposed.

Cook, T. S.↗
NASA NTRS1986

Theoretical three-and four-axis gimbal robot wrists

In high-performance flight simulations, a four-axis gimbal system allows all possible rotations with acceptable gimbal angle rates while it avoids the so-callled 'gimbal lock' that occurs when gimbal rotational axes are colinear. In this paper, pertinent equations (including quaternions) are assembled for a hypothetical robot wrist, functionally equivalent to this four-axis gimbal system, and also for a true three-axis gimbal robot wrist. These equations are used to simulate the rotation of a robot hand by the robot wrist in response to operator rotational velocity commands to the robot hand. Near gimbal lock (wrist singularity), excessive rotational rates occur. Scaling the rates, which is necessary for the three-gimbal robot wrist to prevent rate limiting, introduces an undesirable time delay in the robot hand rotation with respect to the commanded rotation. However, the merit of the four-gimbal robot wrist is that the fourth gimbal angle keeps the robot wrist away from the singularity so that the robot hand moves exactly as commanded. It appears that in a 'worst-type' maneuver of the robot hand, the fourth gimbal angle can be defined so that none of the gimbal angle rates exceed about twice the commanded rates.

Barker, L. K.↗
NASA NTRS1989

Blow-up of unsteady two-dimensional Euler and Navier-Stokes solutions having stagnation-point form

The time-dependent form of the classic, two-dimensional stagnation-point solution of the Navier-Stokes equations is considered. If the viscosity is zero, a class of solutions of the initial-value problem can be found in closed form using Lagrangian coordinates. These solutions exhibit singular behavior in finite time, because of the infinite domain and unbounded initial vorticity. Thus, the blow-up found by Stuart in three dimensions using the stagnation-point form, also occurs in two. The singularity vanishes under a discrete, finite-dimensional 'point vortex' approximation, but is recovered as the number of vortices tends to infinity. We find that a small positive viscosity does not arrest the breakdown, but does strongly alter its form. Similar results are summarized for certain Boussinesq stratified flows.

Childress, S.↗
NASA NTRS1998

Nonlinear Interaction of Detuned Instability Waves in Boundary-Layer Transition: Amplitude Equations

The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant-triads is presented. In this part of the analysis, the system of partial differential critical-layer equations derived in Part I is solved analytically to yield the amplitude equations which are analyzed using a combination of asymptotic and numerical methods. Numerical solutions of the inviscid non-equilibrium oblique-mode amplitude equations show that the frequency-detuned self-interaction enhances the growth of the lower-frequency oblique modes more than the higher-frequency ones. All amplitudes become singular at the same finite downstream position. The frequency detuning delays the occurrence of the singularity. The spanwise-periodic mean-flow distortion and low-frequency nonlinear modes are generated by the critical-layer interaction between frequency-detuned oblique modes. The nonlinear mean flow and higher harmonics as well as the primary instabilities become as large as the base mean flow in the inviscid wall layer in the downstream region where the distance from the singularity is of the order of the wavelength scale.

Lee, Sang Soo↗
DOE OSTI2024

Asymptotic behavior of the Fourier coefficients and the analytic structure of the QCD equation of state

In this paper we study the universal properties of the baryon chemical potential Fourier coefficients in quantum chromodynamics. We show that by following a well-defined strategy, the Fourier coefficients can be used to locate Yang-Lee edge singularities associated with chiral phase transition (and by extension with the Roberge-Weiss) in the complex chemical potential plane. We comment on the viability of performing this analysis using lattice QCD data. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Nonperturbative resolution of strong coupling singularities in 4D 𝒩 = 1 heterotic M-theory

We investigate the interior of the moduli space of four-dimensional 𝒩 = 1 theories of gravity arising from compactifications of the 𝐸 8 × 𝐸 8 heterotic string on Calabi-Yau threefolds. By studying the threshold corrections to the coupling of the heterotic gauge groups, we infer the existence of a strong coupling singularity for one of the perturbative heterotic gauge groups, which effectively yields an additional finite distance boundary of the classical scalar field space. In heterotic M-theory, this boundary maps to a domain wall solution for which the gauge coupling and the warp factor on one of the Horava-Witten 9-branes diverge, thus highlighting the gravitational origin of the classical strong coupling singularity. The divergence of the warp factor is, however, regulated once nonperturbative effects are taken into account, as we demonstrate by studying the instanton corrections to the 5D BPS domain wall equations. This regularization implies that the classical strong coupling boundary of the scalar field 4D 𝒩 = 1 heterotic M-theory is resolved, indicating that, at the quantum level, the field space can be extended beyond this classical boundary.

Cvetič, Mirjam↗
NASA NTRS1979

Stress-intensity factors for internal surface cracks in cylindrical pressure vessels

The stress intensity factors were calculated by a three dimensional finite element method. The finite element models employed singularity elements along the crack front and linear strain elements elsewhere. The models had about 6500 degrees of freedom. The stress intensity factors were evaluated from a nodal force method. An equation for the stress intensity factors was obtained form the results of the present analysis. The equation applies over a wide range of configuration parameters and was within about 5 percent of the present results. A comparison was made between the present results and other analyses of internal surface cracks in cylinders. The results from a boundary integral equation method were in agreement (+ or - 2 percent) and those from another finite element were in fair agreement (+ or - 8 percent) with the present results.

Newman, J. C., Jr.↗
NASA NTRS1982

Out-of-Core Solutions of Complex Sparse Linear Equations

ETCLIB is library of subroutines for obtaining out-of-core solutions of complex sparse linear equations. Routines apply to dense and sparse matrices too large to be stored in core. Useful for solving any set of linear equations, but particularly useful in cases where coefficient matrix has no special properties that guarantee convergence with any of interative processes. The only assumption made is that coefficient matrix is not singular.

Yip, E. L.↗