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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

On the 'delta-equations' for vortex sheet evolution

We use a set of equations, sometimes referred to as the 'delta-equations', to approximate the two-dimensional inviscid motion of an initially circular vortex sheet released from rest in a cross-flow. We present numerical solutions of these equations for the case with delta-square = 0 (for which the equations are exact) and for delta-square greater than 0. For small values of the smoothing parameter delta, a spectral filter must be used to eliminate spurious instabilities due to round-off error. Two singularities appear simultaneously in the vortex sheet when delta-square = 0 at a critical time t(c). After t(c), the solutions do not converge as the computational mesh is refined. With delta-square greater than 0, converged solutions were found for all values of delta-square when t is less than t(c), and for all but the two smallest values of delta-square used when t is greater than t(c). Our results show that, when delta-square is greater than 0, the vortex sheet deforms into two doubly branched spirals some time after t(c). The limiting solution as delta approaching 0 clearly exists and equals the delta = 0 solution when t is less than t(c).

Rottman, James W.↗
NASA NTRS1992

Confining potential in momentum space

A method is presented for the solution in momentum space of the bound state problem with a linear potential in r space. The potential is unbounded at large r leading to a singularity at small q. The singularity is integrable, when regulated by exponentially screening the r-space potential, and is removed by a subtraction technique. The limit of zero screening is taken analytically, and the numerical solution of the subtracted integral equation gives eigenvalues and wave functions in good agreement with position space calculations.

Norbury, John W.↗
NASA NTRS1994

A Galerkin method for linear PDE systems in circular geometries with structural acoustic applications

A Galerkin method for systems of PDE's in circular geometries is presented with motivating problems being drawn from structural, acoustic, and structural acoustic applications. Depending upon the application under consideration, piecewise splines or Legendre polynomials are used when approximating the system dynamics with modifications included to incorporate the analytic solution decay near the coordinate singularity. This provides an efficient method which retains its accuracy throughout the circular domain without degradation at singularity. Because the problems under consideration are linear or weakly nonlinear with constant or piecewise constant coefficients, transform methods for the problems are not investigated. While the specific method is developed for the two dimensional wave equations on a circular domain and the equation of transverse motion for a thin circular plate, examples demonstrating the extension of the techniques to a fully coupled structural acoustic system are used to illustrate the flexibility of the method when approximating the dynamics of more complex systems.

Smith, Ralph C.↗
NASA NTRS1994

Influence of cross-flow on nonlinear Tollmien-Schlichting/vortex interaction

The transition of an incompressible three-dimensional boundary layer with strong cross-flow is considered theoretically and computationally in the context of vortex/wave interactions. Specifically the work centers on two lower-branch Tollmien-Schlichting waves which mutually interact nonlinearly to induce a longitudinal vortex flow. The vortex motion in turn gives rise to significant wave modulation via wall-shear forcing. The characteristic Reynolds number is large and, as a consequence, the waves' and the vortex motion are governed primarily by triple deck theory. The nonlinear interaction is captured by a viscous partial-differential system for the vortex coupled with a pair of amplitude equations for each wave pressure. Following analysis and computation over a wide range of parameters, three distinct responses are found to emerge in the nonlinear behavior of the flow solution downstream: an algebraic finite-distance singularity, far-downstream saturation or far-downstream wave decay leaving pure vortex flow. These depend on the input conditions, the wave angles and the size of the cross flow.

Davis, D. A. R.↗
NASA NTRS1998

The Real-Time Wall Interference Correction System of the NASA Ames 12-Foot Pressure Wind Tunnel

An improved version of the Wall Signature Method was developed to compute wall interference effects in three-dimensional subsonic wind tunnel testing of aircraft models in real-time. The method may be applied to a full-span or a semispan model. A simplified singularity representation of the aircraft model is used. Fuselage, support system, propulsion simulator, and separation wake volume blockage effects are represented by point sources and sinks. Lifting effects are represented by semi-infinite line doublets. The singularity representation of the test article is combined with the measurement of wind tunnel test reference conditions, wall pressure, lift force, thrust force, pitching moment, rolling moment, and pre-computed solutions of the subsonic potential equation to determine first order wall interference corrections. Second order wall interference corrections for pitching and rolling moment coefficient are also determined. A new procedure is presented that estimates a rolling moment coefficient correction for wings with non-symmetric lift distribution. Experimental data obtained during the calibration of the Ames Bipod model support system and during tests of two semispan models mounted on an image plane in the NASA Ames 12 ft. Pressure Wind Tunnel are used to demonstrate the application of the wall interference correction method.

Ulbrich, Norbert↗
NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗
NASA NTRSDate not supplied

A New Method for the Accurate Calculation of Derivatives at Structured Multiblock Grid Singularities

The new Structured Mapping At Reduced dimension Topology elements (SMART) method for the accurate calculation of derivatives at multiblock structured grid block interfaces is presented. In the SMART method, a ’local’ curvilinear mapping is defined at each interface point, giving a consistent multidimensional definition of the derivatives. The resulting system of equations is overdetermined, and a least squares method is used to determine the value of the derivatives in the local coordinate system. In the special case of a nonsingular interface point, the SMART method automatically recovers the spatial differencing stencils used at the interior points in the neighboring grid blocks.

Computational Aeroacoustics↗
DOE OSTI2021

SU(n) and Quantum SU(n) Symmetries in Physical Systems [Slides]

Presence of SU(n) or other Lie group symmetry in a physical system is its powerful, usually underutilized property. In many cases it allows for finding analytical solutions to nonlinear differential equations describing this system. Power of the method is presented on diversified examples from mathematical physics: Lie-group symmetries in finding solutions of generalized, multidimensional theory of gravity; analytical Dirac–equation solutions for description of conducting polymers; stability of qubit states in quantum computers; spatial defects in condensed matter; reconstruction of 3D object from its 2D tomographic image; significant improvement of numerical solutions stability for Euler equations. The next question after obtaining such Lie group symmetric solution is: does a generalized solution with appropriate quantum group symmetry exists for the given physical system, and if yes what is the physical meaning of the deformation parameter q introduced by such solution. In many cases it can be identified. Any SU(n) solution is by its nature singular, assuming a perfect symmetry of the physical system discussed. Such solution gives a powerful insight to theoretical physics, yet the assumption may be too demanding for experimental applications. Deformation parameter q from a quantum group symmetry allows for a continuum of solutions, more applicable to experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1972

Lifting-surface theory for calculating the loading induced on a wing by a flap

A method is described for using lifting-surface theory to obtain the pressure distribution on a wing with a trailing-edge flap or control surface. The loading has a logarithmic singularity at the flap edges, which may be determined directly by the method of matched asymptotic expansions. Expressions are given for the singular flap loading for various flap hinge line and side edge geometries, both for steady and unsteady flap deflection. The regular part of the flap loading must be obtained by inverting the lifting-surface-theory integral equation relating the pressure and the downwash on the wing: procedures are described to accomplish this for a general wing and flap geometry. The method is applied to several example wings, and the results are compared with experimental data. Theory and test correlate well.

Johnson, W. A.↗
NASA NTRS1973

Optimal low thrust geocentric transfer

A computer code which will rapidly calculate time-optimal low thrust transfers is being developed as a mission analysis tool. The final program will apply to NEP or SEP missions and will include a variety of environmental effects. The current program assumes constant acceleration. The oblateness effect and shadowing may be included. Detailed state and costate equations are given for the thrust effect, oblateness effect, and shadowing. A simple but adequate model yields analytical formulas for power degradation due to the Van Allen radiation belts for SEP missions. The program avoids the classical singularities by the use of equinoctial orbital elements. Kryloff-Bogoliuboff averaging is used to facilitate rapid calculation. Results for selected cases using the current program are given.

Edelbaum, T. N.↗
NASA NTRS1990

Doppler radar imaging of spherical planetary surfaces

This paper describes a technique for using echo power-spectra for making radar images of spherical planetary targets. In developing the Doppler-radar imaging system, the target's reflectivity distribution is expanded in a truncated spherical harmonic series; the distribution of echo power in rotational phase and the Doppler frequency are obtained in form of a system of linear equations which express the echo spectra as functions of the series coefficients. To estimate the coefficients from an observed phase-Doppler distribution, the inversion is cast as a least-squares problem, and is solved using singular value decompositions, yielding a linear imaging system which can be fully characterized by its impulse response as a function of latitude. Simulations were designed and used to explore the capabilities of Doppler radar imaging, and the sensitivity of imaging to such factors as subradar altitude coverage and signal-to-noise ratio was investigated.

Hudson, R. Scott↗
NASA NTRS2004

Kalman Filter Estimation of Spinning Spacecraft Attitude using Markley Variables

There are several different ways to represent spacecraft attitude and its time rate of change. For spinning or momentum-biased spacecraft, one particular representation has been put forward as a superior parameterization for numerical integration. Markley has demonstrated that these new variables have fewer rapidly varying elements for spinning spacecraft than other commonly used representations and provide advantages when integrating the equations of motion. The current work demonstrates how a Kalman filter can be devised to estimate the attitude using these new variables. The seven Markley variables are subject to one constraint condition, making the error covariance matrix singular. The filter design presented here explicitly accounts for this constraint by using a six-component error state in the filter update step. The reduced dimension error state is unconstrained and its covariance matrix is nonsingular.

Sedlak, Joseph E.↗
NASA NTRS1989

A hybrid-perturbation-Galerkin technique which combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be sued to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗
NASA NTRS1990

A hybrid perturbation-Galerkin technique that combines multiple expansions

A two-step hybrid perturbation-Galerkin method for the solution of a variety of differential equations type problems is found to give better results when multiple perturbation expansions are employed. The method assumes that there is parameter in the problem formulation and that a perturbation method can be used to construct one or more expansions in this perturbation coefficient functions multiplied by computed amplitudes. In step one, regular and/or singular perturbation methods are used to determine the perturbation coefficient functions. The results of step one are in the form of one or more expansions each expressed as a sum of perturbation coefficient functions multiplied by a priori known gauge functions. In step two the classical Bubnov-Galerkin method uses the perturbation coefficient functions computed in step one to determine a set of amplitudes which replace and improve upon the gauge functions. The hybrid method has the potential of overcoming some of the drawbacks of the perturbation and Galerkin methods as applied separately, while combining some of their better features. The proposed method is applied, with two perturbation expansions in each case, to a variety of model ordinary differential equations problems including: a family of linear two-boundary-value problems, a nonlinear two-point boundary-value problem, a quantum mechanical eigenvalue problem and a nonlinear free oscillation problem. The results obtained from the hybrid methods are compared with approximate solutions obtained by other methods, and the applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗
NASA NTRS2016

New Approaches to Inferences for Steep-Sided Domes on Venus

New mathematical approaches for the relaxation and emplacement of viscous lava domes are presented and applied to steep-sided domes on Venus. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry for two distinct scenarios. In the first scenario, dome relaxation is explored assuming a constant volume of fluid (i.e. lava) has been rapidly emplaced onto the surface. Cooling of lava is represented by a time-variable viscosity and singularities inherent in previous models for dome relaxation have been eliminated. At the onset of relaxation, bulk dynamic viscosities lie in the range between 10(exp 10) - 10(exp 16) Pa s, consistent with basaltic-andesite to rhyolitic compositions. Plausible relaxation times range from 5 to 5000 years, depending on initial lava viscosity. The first scenario, however, is only valid during the final stages of dome relaxation and does not consider the time taken for lava to be extruded onto the surface. In the second scenario, emplacement and growth of a steep-sided dome is considered when the volume of lava on the surface increases over time (i.e. time-variable volume approach). The volumetric flow rate may depend on an arbitrary power of the dome thickness, thus embracing Newtonian as well as other rheologies for describing terrestrial and planetary mass flows. The approach can be used to distinguish between basic flow rate models for fluid emplacement. The formalism results in radial expansion of a dome proportional to t(sup 1/2), consistent with the diffusive nature of the governing equation. The flow at the front is shown to thicken as the front advances for a constant rate of lava supply. Emplacement times are intimately correlated with the bulk rheology. Comparison of the theoretical profiles with the shape of a typical dome on Venus indicates that a Newtonian bulk rheology is most appropriate, consistent with prior studies. However, results here suggest a bulk dynamic viscosity of 10(exp 12) - 10(exp 13) Pa s and emplacement times of approximately 2-16 years. Both scenarios investigated give emplacement times significantly less than prior estimates and compositions consistent with basaltic andesite.

steep-sided domes↗
NASA NTRS1978

A double shooting scheme for certain unstable and singular boundary value problems

A scheme is presented to obtain the unique bounded solution for an exponentially unstable linear system. The scheme consists of choosing random data at large initial values and integrating forwards and backwards until accurate regular boundary values are obtained. Proofs of convergence are given for the case that the homogeneous equation has an exponential dichotomy. Applications to other types of problems are discussed and numerical results are presented.

Bayliss, A.↗
DOE OSTI2022

Parametric model-order reduction for radiation transport using multi-resolution proper orthogonal decomposition

For parametric high-fidelity simulations, it is often desirable to utilize a reduced-order model (ROM) to emulate, at a reduced computational cost, parametric solutions of the governing partial differential equations (PDEs) for unseen parameter values. One commonly employed option is to utilize a data-driven, projection-based ROM supplemented with subspace identification via proper orthogonal decomposition (POD). POD discovers the ROM subspace by computing the singular value decomposition (SVD) of a set of training data from the full-order model (FOM). In streaming-dominated radiation transport simulations with localized sources, solutions often greatly vary over the spatial domain by many orders of magnitude. In such cases, machine-precision arithmetic can be insufficient to obtain an accurate SVD, resulting in a poorly performing ROM. We present a method called multiresolution POD (mrPOD) that mitigates these inaccuracies. The mrPOD method works by decomposing the spatial domain into regions and performing proper orthogonal decomposition on the training dataset separately in each region. In conclusion, mrPOD is tested on single energy group and multigroup atmospheric shielding transport problems and is shown to outperform classic POD.

42 ENGINEERING↗
DOE OSTI2021

Geometry and entanglement in the scattering matrix

A formulation of nucleon–nucleon scattering is developed in which the S-matrix, rather than an effective-field theory (EFT) action, is the fundamental object. Spacetime plays no role in this description: the S-matrix is a trajectory that moves between RG fixed points in a compact theory space defined by unitarity. This theory space has a natural operator definition, and a geometric embedding of the unitarity constraints in four-dimensional Euclidean space yields a flat torus, which serves as the stage on which the S-matrix propagates. Trajectories with vanishing entanglement are special geodesics between RG fixed points on the flat torus, while entanglement is driven by an external potential. The system of equations describing S-matrix trajectories is in general complicated, however the very-low-energy S-matrix –that appears at leading-order in the EFT description– possesses a UV/IR conformal invariance which renders the system of equations integrable, and completely determines the potential. In this geometric viewpoint, inelasticity is in correspondence with the radius of a three-dimensional hyperbolic space whose two-dimensional boundary is the flat torus. This space has a singularity at vanishing radius, corresponding to maximal violation of unitarity. The trajectory on the flat torus boundary can be explicitly constructed from a bulk trajectory with a quantifiable error, providing a simple example of a holographic quantum error correcting code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗