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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Equations of motion in the linear approximation.

Attempt to develop a gauge-invariant theory of the motion of singularities in an n-dimensional Riemannian space. A gauge-invariant theory is developed for a nongeodetic particle represented by a time-like world line in the linearized field theory. The solution of the basic algebraic and differential equations is to be constructed from the source line, the future null cones emanating from it, a scale factor, and nothing else, and the leading term of the solution is the Schwarzschild field. A class of possible solutions is examined, and it is shown that this class contains just one acceptable member - namely, a particle in constant acceleration.

Robinson, I.↗
NASA NTRS1973

Laminar film boiling on inclined isothermal flat plates.

Laminar film boiling from an inclined flat plate has been investigated analytically. Using the singular perturbation scheme, the complete set of Navier-Stokes equations is solved. The zeroth-order perturbation coinciding with the boundary-layer equations for vertical flat plates governs the problem. The higher-order perturbations become important near the leading edge and for large values of the inclination angle from the vertical. The assumption of zero interfacial velocity shows that, except for fluids having large (rho x mu) ratios, the results can be predicted using the vertical flat plate results by defining a modified Grashof parameter containing a cos phi term. When the interfacial shear is considered, the solutions indicate that for fluids having large (rho x mu) ratios, the heat transfer rates will be larger (approximately 15% maximum) than those predicted by the simplified model using zero interfacial velocity. In general, the inclination decreases the rate of heat transfer as well as the rate of evaporation.

Nagendra, H. R.↗
NASA NTRS1976

Solution of linear systems by a singular perturbation technique

An approximate solution is obtained for a singularly perturbed system of initial valued, time invariant, linear differential equations with multiple boundary layers. Conditions are stated under which the approximate solution converges uniformly to the exact solution as the perturbation parameter tends to zero. The solution is obtained by the method of matched asymptotic expansions. Use of the results for obtaining approximate solutions of general linear systems is discussed. An example is considered to illustrate the method and it is shown that the formulas derived give a readily computed uniform approximation.

Ardema, M. D.↗
NASA NTRS1977

Interacting turbulent boundary layer over a wavy wall

The two dimensional supersonic flow of a thick turbulent boundary layer over a train of relatively small wave-like protuberances is considered. The flow conditions and the geometry are such that there exists a strong interaction between the viscous and inviscid flow. The problem cannot be solved without inclusion of interaction effects due to the occurrence of the separation singularity in classical boundary layer methods. The interacting boundary layer equations are solved numerically using a time-like relaxation method with turbulence effects represented by the inclusion of the eddy viscosity model. Results are presented for flow over a train of up to six waves for Mach numbers of 10 and 32 million/meter, and wall temperature rations (T sub w/T sub 0) of 0.4 and 0.8. Limited comparisons with independent experimental and analytical results are also given. Detailed results on the influence of small protuberances on surface heating by boundary layers are presented.

Polak, A.↗
NASA NTRS1983

An implicit, transonic, full-potential code for cascade flow on H-grid topology

A transonic, full-potential code is developed for computing the flow through two-dimensional cascades using an H-type grid topology that employs an implicit approximate-factorization scheme. The body-conforming H-grid is generated numerically by solving Poisson's equation. The flow-solution algorithm at the coordinate mapping singularity associated with this grid is investigated using two different types of finite-difference schemes. The grid-geometry effect on these schemes is also studied by noting free-stream capturing properties. It is found that by implementing a consistent spatial differencing scheme, the mapping singularities can be resolved numerically, and the grid-geometry-induced error minimized. The code is verified by computing model cascade flow problems.

Kwak, D.↗
NASA NTRS1983

Measurements of the scattering of sound by a line vortex

This paper presents measurements of the phase and magnitude of the scattered field arising from the incidence of a monochromatic plane sound field as a steady vortex. The amplitude of the scattered field was found to vary linearly with the vortex strength, and with the incident wave amplitude and frequency as predicted by solutions based on the Born approximation. The scattered field was observed to be nonsingular in the incidence direction, and this was similar to predictions by the Parabolic Equation Method (PEM) rather than the Born approximation, which predicts singular behavior in the incidence direction.

Horne, W. C.↗
NASA NTRS1987

The quadrupole gravitational lens

An asymmetric lens made of a spherically symmetric galaxy and a quadrupole tidal part allows a very simple pictorial examination of the multiplicity and positions of the images, for any mass density profile. The lens equations can be reduced to one equation in one variable. In the special case of a singular isothermal sphere with a tide (SIST), the lens properties are described by five parameters. Two parameters determine the angular scale and orientation, while three determine the configuration of the images. A comprehensive investigation of the SIST lens is performed, including the variation of the configurations with the parameter values, the critical surfaces in the parameter space, and the extremal solutions. Analytic expressions are derived for the SIST parameters, and additional image positions, in the case when the intensity ratio and positions relative to the galaxy for two of the images are given as input. Quadrupole lenses with a point mass and a pseudoisothermal sphere as the main part are also described in terms of direction diagrams. A direct relation between deflection profile and rotational velocity is derived.

Kovner, Israel↗
NASA NTRS1998

Optimization of Supersonic Transport Trajectories

This paper develops a near-optimal guidance law for generating minimum fuel, time, or cost fixed-range trajectories for supersonic transport aircraft. The approach uses a choice of new state variables along with singular perturbation techniques to time-scale decouple the dynamic equations into multiple equations of single order (second order for the fast dynamics). Application of the maximum principle to each of the decoupled equations, as opposed to application to the original coupled equations, avoids the two point boundary value problem and transforms the problem from one of a functional optimization to one of multiple function optimizations. It is shown that such an approach produces well known aircraft performance results such as minimizing the Brequet factor for minimum fuel consumption and the energy climb path. Furthermore, the new state variables produce a consistent calculation of flight path angle along the trajectory, eliminating one of the deficiencies in the traditional energy state approximation. In addition, jumps in the energy climb path are smoothed out by integration of the original dynamic equations at constant load factor. Numerical results performed for a supersonic transport design show that a pushover dive followed by a pullout at nominal load factors are sufficient maneuvers to smooth the jump.

Ardema, Mark D.↗
NASA NTRS1998

Development of Advanced Methods of Structural and Trajectory Analysis for Transport Aircraft

This paper develops a near-optimal guidance law for generating minimum fuel, time, or cost fixed-range trajectories for supersonic transport aircraft. The approach uses a choice of new state variables along with singular perturbation techniques to time-scale decouple the dynamic equations into multiple equations of single order (second order for the fast dynamics). Application of the maximum principle to each of the decoupled equations, as opposed to application to the original coupled equations, avoids the two point boundary value problem and transforms the problem from one of a functional optimization to one of multiple function optimizations. It is shown that such an approach produces well known aircraft performance results such as minimizing the Brequet factor for minimum fuel consumption and the energy climb path. Furthermore, the new state variables produce a consistent calculation of flight path angle along the trajectory, eliminating one of the deficiencies in the traditional energy state approximation. In addition, jumps in the energy climb path are smoothed out by integration of the original dynamic equations at constant load factor. Numerical results performed for a supersonic transport design show that a pushover dive followed by a pullout at nominal load factors are sufficient maneuvers to smooth the jump.

Ardema, Mark D.↗
NASA NTRS1999

Weakly Nonlinear Description of Parametric Instabilities in Vibrating Flows

This project focuses on the effects of weak dissipation on vibrational flows in microgravity and in particular on (a) the generation of mean flows through viscous effects and their reaction on the flows themselves, and (b) the effects of finite group velocity and dispersion on the resulting dynamics in large domains. The basic mechanism responsible for the generation of such flows is nonlinear and was identified by Schlichting [21] and Longuet-Higgins. However, only recently has it become possible to describe such flows self-consistently in terms of amplitude equations for the parametrically excited waves coupled to a mean flow equation. The derivation of these equations is nontrivial because the limit of zero viscosity is singular. This project focuses on various aspects of this singular problem (i.e., the limit C equivalent to (nu)((g)(h(exp 3)))exp -1/2 << 1,where nu is the kinematic viscosity and h is the liquid depth) in the weakly nonlinear regime. A number of distinct cases is identified depending on the values of the Bond number, the size of the nonlinear terms, distance above threshold and the length scales of interest. The theory provides a quantitative explanation of a number of experiments on the vibration modes of liquid bridges and related experiments on parametric excitation of capillary waves in containers of both small and large aspect ratio. The following is a summary of results obtained thus far.

Knobloch, E.↗
NASA NTRS1983

Energy state revisited

Kaiser (1944) has introduced the concept of "resultant height' in connection with aircraft minimum-time climbs. Its use as a state variable in trajectory work is attractive because it is a "slower' variable than either altitude or velocity. Kelley (1972, 1973) has made an attempt to synthesize "slow' state variables in connection with singular-perturbation procedures. In the present investigation, attempts are made to synthesize both "fast' and "slow' variables for the minimum-time-climb problem along lines explored by Kelley. Attention is given to climb equations, energy-modeling simplifications, "slow'-variable choice, "fast'-variable-choice considerations, a singular-perturbation analysis, the choise of a "fast' variable, and the climb-dash problem.

Kelley, H. J.↗
NASA NTRS1983

Energy state revisited

Kaiser (1944) has introduced the concept of 'resultant height' in connection with aircraft minimum-time climbs. Its use as a state variable in trajectory work is attractive because it is a 'slower' variable than either altitude or velocity. Kelley (1972, 1973) has made an attempt to synthesize 'slow' state variables in connection with singular-perturbation procedures. In the present investigation, attempts are made to synthesize both 'fast' and 'slow' variables for the minimum-time-to-climb problem along lines explored by Kelley. Attention is given to climb equations, energy-modeling simplifications, 'slow'-variable choice, 'fast'-variable-choice considerations, a singular-perturbation analysis, the choice of a 'fast' variable, and the climb-dash problem.

Kelley, H. J.↗
NASA NTRS1995

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗
DOE OSTI2024

A Solution Method for the Filtered Lifting Line Theory

The filtered lifting line theory presents a continuous form of the inviscid momentum equations of flow over a lifting device, such as a wing or rotor blade, using body forces without mathematical singularities. This theory is also consistent with an actuator line representation of a lifting device. In this work, we present a reformulation of the equations in terms of the local flow angle along the line, which allows solving the stand-alone equations using multivariate root-finding algorithms. This approach can be used to obtain a fast, computationally inexpensive solution of the loading distribution along a wing without the need to perform computational fluid dynamic simulations. We study the requirements in terms of resolution in the spanwise direction and establish the criteria for spacing and minimum amount of points required along the blade to obtain converged solutions. The solutions are compared to results from large-eddy simulations, and we observed excellent agreement with less than a percent difference in quantities along the blade between the methods.

17 WIND ENERGY↗
NASA NTRS1975

Perturbation solution of the Navier-Stokes equations and its relation to the Lighthill-Curle solution of aerodynamic sound

The aerodynamic sound described by the Lighthill-Curle solution is reexamined using a method of matched asymptotic expansions. The governing Navier-Stokes equations written in nondimensional form are expanded for a small Mach number. First- and second-order solutions for the pressure field are obtained, and the singular nature of the expansion at large distances is indicated. The nearfield pressure is governed by the Poisson equation, whereas the farfield equations describe a linear wave system in a dissipative medium. The pseudosound is related to the incompressible Reynolds stresses associated with a solenoidal velocity field, the velocity, the pressure perturbation, and their derivatives on the boundaries. A uniformly valid first-order solution for the pressure is obtained. It is shown that viscosity, thermal conductivity, and entropy in the flow do not contribute to the first-order noise generation, while the viscous stress contributes to noise only from some boundaries. The application of the proposed perturbation method to a subsonically moving surface and a hot jet is discussed.

Pan, Y. S.↗
NASA NTRS1988

Three-dimensional singular points in aerodynamics

When three-dimensional separation occurs on a body immersed in a flow governed by the incompressible Navier-Stokes equations, the geometrical surfaces formed by the three vector fields (velocity, vorticity and the skin-friction) and a scalar field (pressure) become interrelated through topological maps containing their respective singular points and extremal points. A mathematically consistent description of these singular points becomes inevitable when we want to study the geometry of the separation. A separated stream surface requires, for example, the existence of a saddle-type singular point on the skin-friction surface. This singular point is actually, in the proper language of mathematics, a saddle of index two. The index is a measure of the dimension of the outset (set leaving the singular point). Hence, when a saddle of index two is specified, a two dimensional surface that becomes separated from the osculating plane of the saddle is implied. The three-dimensional singular point is interpreted mathematically and the most common aerodynamical singular points are discussed through this perspective.

Unal, Aynur↗
DOE OSTI2024

Jacobian-based Model Diagnostics and Application to Equation Oriented Modeling of a Carbon Capture System

Equation-oriented (EO) modeling has the potential to enable the effective design and optimization of the operation of advanced energy systems. However, advanced modeling of energy systems results in a large number of variables and non-linear equations, and it can be difficult to search through these to identify the culprit(s) responsible for convergence issues. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms so they can be rescaled. A further singular value decomposition can be per-formed to identify degenerate sets of equations and remaining scaling issues. This work presents an EO model of a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. The IDAES diagnostics tools were successfully applied to this flowsheet to identify problems to improve model robustness and enable the optimization of process design and operating conditions of a carbon capture system.

Allan, Douglas↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗