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 667 records · Page 37

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Advanced panel-type influence coefficient methods applied to subsonic flows

An advanced technique for solving the linear integral equations of three-dimensional subsonic potential flows (steady, inviscid, irrotational and incompressible) about arbitrary configurations is presented. It involves assembling select, logically consistent networks whose construction comprises four tasks, which are described in detail: surface geometry definition; singularity strength definition; control point and boundary condition specification; and calculation of induced potential or velocity. The technique is applied to seven wing examples approached by four network types: source/analysis, doublet/analysis, source/design, and doublet/design. The results demonstrate the forgiveness of the model to irregular paneling and the practicality of combined analysis/design boundary conditions. The appearance of doublet strength mismatch is a valuable indicator of locally inadequate paneling.

Johnson, F. T.↗

NASA NTRS1984

Methods for analysis of cracks in three-dimensional solids

Various analytical and numerical methods used to evaluate the stress intensity factors for cracks in three-dimensional (3-D) solids are reviewed. Classical exact solutions and many of the approximate methods used in 3-D analyses of cracks are reviewed. The exact solutions for embedded elliptic cracks in infinite solids are discussed. The approximate methods reviewed are the finite element methods, the boundary integral equation (BIE) method, the mixed methods (superposition of analytical and finite element method, stress difference method, discretization-error method, alternating method, finite element-alternating method), and the line-spring model. The finite element method with singularity elements is the most widely used method. The BIE method only needs modeling of the surfaces of the solid and so is gaining popularity. The line-spring model appears to be the quickest way to obtain good estimates of the stress intensity factors. The finite element-alternating method appears to yield the most accurate solution at the minimum cost.

Raju, I. S.↗

NASA NTRS1990

A fully Sinc-Galerkin method for Euler-Bernoulli beam models

A fully Sinc-Galerkin method in both space and time is presented for fourth-order time-dependent partial differential equations with fixed and cantilever boundary conditions. The Sinc discretizations for the second-order temporal problem and the fourth-order spatial problems are presented. Alternate formulations for variable parameter fourth-order problems are given which prove to be especially useful when applying the forward techniques to parameter recovery problems. The discrete system which corresponds to the time-dependent partial differential equations of interest are then formulated. Computational issues are discussed and a robust and efficient algorithm for solving the resulting matrix system is outlined. Numerical results which highlight the method are given for problems with both analytic and singular solutions as well as fixed and cantilever boundary conditions.

Smith, R. C.↗

NASA NTRS2002

High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations

We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].

Bryson, Steve↗

NASA NTRS1985

The origins of cosmic rays and quantum effects on gravity

The energy spectrum of primary cosmic rays is explained by particles emitted during a thermal expansion of explosive objects inside and near the galaxy, remnants of which may be supernova and/or active talaxies, or even stars or galaxies that disappeared from our sight after the explosion. A power law energy spectrum for cosmic rays, E to the (-alpha -1, is obtained from an expansion rate T is proportional to R to the alpha. Using the solution of the Einstein equation, we obtain a spectrum which agrees very well with experimental data. The implication of an inflationary early universe on the cosmic ray spectrum is also discussed. It is also suggested that the conflict between this model and the singularity theorem in classical general relativity may be eliminated by quantum effects.

Tomozawa, Y.↗

NASA NTRS1989

Integrated tools for control-system analysis

The basic functions embedded within a user friendly software package (MATRIXx) are used to provide a high level systems approach to the analysis of linear control systems. Various control system analysis configurations are assembled automatically to minimize the amount of work by the user. Interactive decision making is incorporated via menu options and at selected points, such as in the plotting section, by inputting data. There are five evaluations such as the singular value robustness test, singular value loop transfer frequency response, Bode frequency response, steady-state covariance analysis, and closed-loop eigenvalues. Another section describes time response simulations. A time response for random white noise disturbance is available. The configurations and key equations used for each type of analysis, the restrictions that apply, the type of data required, and an example problem are described. One approach for integrating the design and analysis tools is also presented.

Ostroff, Aaron J.↗

NASA NTRS1991

Forward and inverse kinematics of double universal joint robot wrists

A robot wrist consisting of two universal joints can eliminate the wrist singularity problem found on many individual robots. Forward and inverse position and velocity kinematics are presented for such a wrist having three degrees of freedom. Denavit-Hartenberg parameters are derived to find the transforms required for the kinematic equations. The Omni-Wrist, a commercial double universal joint robot wrist, is studied in detail. There are four levels of kinematic parameters identified for this wrist; three forward and three inverse maps are presented for both position and velocity. These equations relate the hand coordinate frame to the wrist base frame. They are sufficient for control of the wrist standing alone. When the wrist is attached to a manipulator arm; the offset between the two universal joints complicates the solution of the overall kinematics problem. All wrist coordinate frame origins are not coincident, which prevents decoupling of position and orientation for manipulator inverse kinematics.

Williams, Robert L., II↗

NASA NTRS1992

Kinetic treatment of nonlinear magnetized plasma motions - General geometry and parallel waves

The expansion of kinetic equations in the limit of a strong magnetic field is presented. This gives a natural description of the motions of magnetized plasmas, which are slow compared to the particle gyroperiods and gyroradii. Although the approach is 3D, this very general result is used only to focus on the parallel propagation of nonlinear Alfven waves. The derivative nonlinear Schroedinger-like equation is obtained. Two new terms occur compared to earlier treatments, a nonlinear term proportional to the heat flux along the magnetic field line and a higher-order dispersive term. It is shown that kinetic description avoids the singularities occurring in magnetohydrodynamic or multifluid approaches, which correspond to the degenerate case of sound speeds equal to the Alfven speed, and that parallel heat fluxes cannot be neglected, not even in the case of low parallel plasma beta. A truly stationary soliton solution is derived.

Khabibrakhmanov, I. KH.↗

NASA NTRS2009

Stability of Detached Solidification

Bridgman crystal growth can be conducted in the so-called "detached" solidification regime, where the growing crystal is detached from the crucible wall. A small gap between the growing crystal and the crucible wall, of the order of 100 micrometers or less, can be maintained during the process. A meniscus is formed at the bottom of the melt between the crystal and crucible wall. Under proper conditions, growth can proceed without collapsing the meniscus. The meniscus shape plays a key role in stabilizing the process. Thermal and other process parameters can also affect the geometrical steady-state stability conditions of solidification. The dynamic stability theory of the shaped crystal growth process has been developed by Tatarchenko. It consists of finding a simplified autonomous set of differential equations for the radius, height, and possibly other process parameters. The problem then reduces to analyzing a system of first order linear differential equations for stability. Here we apply a modified version of this theory for a particular case of detached solidification. Approximate analytical formulas as well as accurate numerical values for the capillary stability coefficients are presented. They display an unexpected singularity as a function of pressure differential. A novel approach to study the thermal field effects on the crystal shape stability has been proposed. In essence, it rectifies the unphysical assumption of the model that utilizes a perturbation of the crystal radius along the axis as being instantaneous. It consists of introducing time delay effects into the mathematical description and leads, in general, to stability over a broader parameter range. We believe that this novel treatment can be advantageously implemented in stability analyses of other crystal growth techniques such as Czochralski and float zone methods.

Mazuruk, K.↗

DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗

NASA NTRS1971

Nonplanar Method for Predicting Incompressible Aerodynamic Coefficients of Rectangular Wings with Circular-Arc Camber

The development of a nonplanar lifting surface method having a continuous distribution of singularities and satisfying the tangent flow boundary condition on the mean camber surface is given. The method predicts some incompressible longitudinal aerodynamic coefficients of rectangular wings which have circular-arc camber. The solution method is of the integral-equation type and the resulting surface integrals are evaluated by either using numerical or analytical techniques, as are appropriate. Applications are made and the results compared with those from an exact two-dimensional circular-arc camber solution, a three-dimensional flat-wing solution which represents the camber by a projected slope onto the flat surface, and a flat-wing experiment. From these comparisons, the present method is found to predict well the flat-wing experiment and limiting values, in addition to the center of pressure variation at an angle of attack of zero for any camber. For wings having camber ratios larger than about 1.25% and moderate to high aspect ratios, the results deterioriate due to the inadequacy of lifting pressure modes employed.

Lamar, J. E.↗

DOE OSTI2024

Effects of final state interactions on Landau singularities

In certain kinematic and particle mass configurations, triangle singularities may lead to line-shapes which mimic the effects of resonances. This well-known effect is scrutinized here in the presence of final-state rescattering. The goal is achieved first by utilizing general arguments provided by Landau equations, and second by applying a modern scattering formalism with explicit two- and three-body unitarity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Intrinsic Anomalous Hall Conductivity in a Nonuniform Electric Field

We study how the intrinsic anomalous Hall conductivity is modified in two-dimensional crystals with broken time-reversal symmetry due to weak inhomogeneity of the applied electric field. Focusing on a clean noninteracting two-band system without band crossings, we derive the general expression for the Hall conductivity at small finite wave vector q to order q 2 , which governs the Hall response to the second gradient of the electric field. Using the Kubo formula, we show that the answer can be expressed through the Berry curvature, Fubini-Study quantum metric, and the rank-3 symmetric tensor which is related to the quantum geometric connection and physically corresponds to the gauge-invariant part of the third cumulant of the position operator. We further compare our results with the predictions made within the semiclassical approach. By deriving the semiclassical equations of motion, we reproduce the result obtained from the Kubo formula in some limits. We also find, however, that the conventional semiclassical description in terms of the definite position and momentum of the electron is not fully consistent because of singular terms originating from the Heisenberg uncertainty principle. We thus present a clear example of a case when the semiclassical approach inherently suffers from the uncertainty principle, implying that it should be applied to systems in nonuniform fields with extra care.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1993

Finite-surface method for the Maxwell equations in generalized coordinates

A finite-surface technique for the Maxwell equations in generalized nonorthogonal coordinates is developed. It directly applies the integral Faraday's and Ampere's laws to faces of primary and secondary grid cells, respectively. The technique features an accurate treatment of matching conditions at a material interface, grid singularities, and radiation conditions at outer boundaries.

Vinokur, Marcel↗

NASA NTRS2012

Quaternion-Based Control Architecture for Determining Controllability/Maneuverability Limits

Dynamic inversion has often been used in the simulation environment to rapidly prototype controls for the full flight envelope, because of its capacity for assessing a vehicle s maneuver performance and proper sizing of control surfaces. Generally, the architectures involve either a direct inversion of the entire set of equations of motion or a sequential set of inversions exploiting time scale separation in the vehicle dynamics where faster parameters are considered as controls for slower varying parameters. The proposed architecture builds on the latter using a quaternion formulation that provides singularity free tracking. Of interest, the proposed architecture simplifies the sequential approach by exploiting a simpler kinematic inversion in place of a more difficult inversion typically used. This kinematic relationship accurately describes the angular rate required to drive some reference frame of interest to a desired attitude at some desired quaternion error rate. A simple PID control is used to define the desired quaternion error rate. The paper develops the theoretical framework for the approach, and shows results in tracking a desired trajectory.

Bacon, Barton J↗

NASA NTRS1971

Singular perturbation methods and the warm plasma model

The application of techniques of the singular perturbation theory to the analysis of warm plasma is discussed. Typically, the cold plasma model can be applied over wide ranges of parameters and only over narrow ranges forming so-called boundary layers is the warm plasma model used. Simplified equations can be used and the solutions matched on both sides of the layer's boundary. Simple examples to illustrate the solution are presented. The analysis confirms that some results are highly sensitive to the values of: (1) wire radius or gap size for an antenna, (2) temperature of the medium, and (3) incident angle of a plane wave.

Lee, S. W.↗

NASA NTRS1981

Solutions of the Helmholtz equation with boundary conditions for force-free magnetic fields

It is shown that the solution, with one ignorable coordinate, for the Taylor minimum energy state (resulting in a force-free magnetic field) in either a straight cylindrical or a toroidal geometry with arbitrary cross section can be reduced to the solution of either an inhomogeneous Helmholtz equation or a Grad-Shafranov equation with simple boundary conditions. Standard Green's function theory is, therefore, applicable. Detailed solutions are presented for the Taylor state in toroidal and cylindrical domains having a rectangular cross section. The focus is on solutions corresponding to the continuous eigenvalue spectra. Singular behavior at 90 deg corners is explored in detail.

Rasband, S. N.↗

NASA NTRS1989

Dynamic grid refinement for partial differential equations on parallel computers

The fast adaptive composite grid method (FAC) is an algorithm that uses various levels of uniform grids to provide adaptive resolution and fast solution of PDEs. An asynchronous version of FAC, called AFAC, that completely eliminates the bottleneck to parallelism is presented. This paper describes the advantage that this algorithm has in adaptive refinement for moving singularities on multiprocessor computers. This work is applicable to the parallel solution of two- and three-dimensional shock tracking problems.

Mccormick, S.↗