Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Elliptic problems”

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 451 records · Page 25

Bowles-Tatnall Wake Vortex Encounter Hazard Metric

Wake vortex spacing standards constrict the terminal area throughput and impose severe constraints on the overall capacity and efficiency of the National Airspace System. For more than two decades starting in the early 1990s, the National Aeronautics and Space Administration conducted extensive research on characterizing the formation and evolution of aircraft wakes. This multidisciplinary work included comprehensive field experiments (Pruis et al. 2016), flight tests (Vicroy et al. 1998), and wind tunnel tests (Rossow 1994; Chow et al. 1997). Parametric studies using large eddy simulations (Proctor 1998; Proctor et al. 2006) were conducted in order to develop fast-time models for the prediction of wake transport and decay (Ahmad et al. 2016). Substantial effort was spent on the formulation of acceptable vortex hazard metrics (Tatnall 1995; Hinton and Tatnall 1997). Several wake encounter severity metrics have been suggested in the past, which include the wake circulation strength, vortex-induced rolling moment coefficient (Clv), bank angle, and the roll control ratio (Tatnall 1995; Hinton and Tatnall 1997; Van der Geest 2012). The vortex-induced rolling moment coefficient introduced by Bowles and Tatnall (Tatnall 1995; Gloudemans et al. 2016) has been used extensively for risk and safety analysis of newly proposed air traffic management concepts and procedures. The original method of Bowles and Tatnall assumed a constant wing loading (the wing lift-curve slope, CL is constant), which resulted in an overestimation of the vortexinduced rolling moment coefficient. Bowles (2014) suggested a correction to the original method that provides more accurate values of Clv and which is also consistent with the underlying physics of the problem. The overestimation of Clv in the original method can be corrected by assuming an elliptical lift distribution. Figure 1.1 illustrates the correction in Clv achieved by the modified method.

Joel Malissa↗

Hypersonic bodies of maximum drag for a given lift-to-drag ratio.

The problem considered in this paper is concerned with the aerodynamic design of the forebody shape of reentry vehicles in the blunt, homothetic, elliptic transversal contour, power-law longitudinal contour, raked-off configurational set. In particular, the forebody shape which maximizes the ratio of the forebody pressure drag to the free-stream dynamic pressure for a given lift-to-drag ratio and given geometric properties is determined. This problem is considered because recent survey articles indicate that its solution will provide useful qualitative design information about manned vehicles entering the earth's atmosphere from any of the foreseeable planetary missions. Single-integral equations relating the lift and drag in Newtonian hypersonic flow to the forebody geometry are derived and used to formulate the optimization problem which is solved by a direct numerical method.

Mcmillan, W., III↗

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗

Upscaling of Mixed Finite Element Discretization Problems by the Spectral AMGe Method

Here, we propose two multilevel spectral techniques for constructing coarse discretization spaces for saddle-point problems corresponding to PDEs involving a divergence constraint, with a focus on mixed finite element discretizations of scalar self-adjoint second order elliptic equations on general unstructured grids. We use element agglomeration algebraic multigrid (AMGe), which employs coarse elements that can have nonstandard shape since they are agglomerates of fine-grid elements. The coarse basis associated with each agglomerated coarse element is constructed by solving local eigenvalue problems and local mixed finite element problems. This construction leads to stable upscaled coarse spaces and guarantees the inf-sup compatibility of the upscaled discretization. Also, the approximation properties of these upscaled spaces improve by adding more local eigenfunctions to the coarse spaces. The higher accuracy comes at the cost of additional computational effort, as the sparsity of the resulting upscaled coarse discretization (referred to as operator complexity) deteriorates when we introduce additional functions in the coarse space. We also provide an efficient solver for the coarse (upscaled) saddle-point system by employing hybridization, which leads to a symmetric positive definite (s.p.d.) reduced system for the Lagrange multipliers, and to solve the latter s.p.d. system, we use our previously developed spectral AMGe solver. Numerical experiments, in both two and three dimensions, are provided to illustrate the efficiency of the proposed upscaling technique.

97 MATHEMATICS AND COMPUTING↗

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

Discrete restricted four-body problem: Existence of proof of equilibria and reproducibility of periodic orbits

We propose the discrete-time restricted four-body problem (d-R4BP), which approximates the orbits of the restricted four-body problem (R4BP). The d-R4BP is given as a special case of the discrete-time chain regularization of the general N-body problem published in Minesaki. Moreover, we analytically prove that the d-R4BP yields the correct orbits corresponding to the elliptic relative equilibrium solutions of the R4BP when the three primaries form an equilateral triangle at any time. Such orbits include the orbit of a relative equilibrium solution already discovered by Baltagiannis and Papadakis. Until the proof in this work, there has been no discrete analog that preserves the orbits of elliptic relative equilibrium solutions in the R4BP. For a long time interval, the d-R4BP can precisely compute some stable periodic orbits in the Sun–Jupiter–Trojan asteroid–spacecraft system that cannot necessarily be reproduced by other generic integrators.

79 ASTRONOMY AND ASTROPHYSICS↗

Applying Physics-Informed Neural Networks to Solve Navier–Stokes Equations for Laminar Flow around a Particle

In recent years, Physics-Informed Neural Networks (PINNs) have drawn great interest among researchers as a tool to solve computational physics problems. Unlike conventional neural networks, which are black-box models that “blindly” establish a correlation between input and output variables using a large quantity of labeled data, PINNs directly embed physical laws (primarily partial differential equations) within the loss function of neural networks. By minimizing the loss function, this approach allows the output variables to automatically satisfy physical equations without the need for labeled data. The Navier–Stokes equation is one of the most classic governing equations in thermal fluid engineering. This study constructs a PINN to solve the Navier–Stokes equations for a 2D incompressible laminar flow problem. Flows passing around a 2D circular particle are chosen as the benchmark case, and an elliptical particle is also examined to enrich the research. The velocity and pressure fields are predicted by the PINNs, and the results are compared with those derived from Computational Fluid Dynamics (CFD). Additionally, the particle drag force coefficient is calculated to quantify the discrepancy in the results of the PINNs as compared to CFD outcomes. The drag coefficient maintained an error within 10% across all test scenarios.

Hu, Beichao (ORCID:0009000163215151)↗

On the cost of approximating and recognizing a noise perturbed straight line or a quadratic curve segment in the plane

Approximation of noisy data in the plane by straight lines or elliptic or single-branch hyperbolic curve segments arises in pattern recognition, data compaction, and other problems. The efficient search for and approximation of data by such curves were examined. Recursive least-squares linear curve-fitting was used, and ellipses and hyperbolas are parameterized as quadratic functions in x and y. The error minimized by the algorithm is interpreted, and central processing unit (CPU) times for estimating parameters for fitting straight lines and quadratic curves were determined and compared. CPU time for data search was also determined for the case of straight line fitting. Quadratic curve fitting is shown to require about six times as much CPU time as does straight line fitting, and curves relating CPU time and fitting error were determined for straight line fitting. Results are derived on early sequential determination of whether or not the underlying curve is a straight line.

Cooper, D. B.↗

Meteoroid capture into earth orbit by atmospheric drag

A probabilistic analysis, based on orbital mechanics, is performed to evaluate the chances of meteoroids entering into earth orbit and the potential population of such objects. The problem is addressed in terms of meteoroids on earth collision courses, slowed by atmosphere entry/exit, entering elliptical orbits with apogees above the atmosphere. The overall capture probability is derived by integrating the capture fraction dependence on velocity and size over the probability frequency distribution of the collision courses. Account is taken of eventual orbit decay, lunar infall and ejection after encounters with the moon. The results indicate that the probability that a natural 10-100 m diam object has achieved earth orbit is negligibly small.

Friedlander, A. L.↗

A method for determining the camber and twist of a surface to support a given distribution of lift

A graphical method is described for finding the shape (camber and twist) of an airfoil having an arbitrary distribution of lift. The method consists in replacing the lifting surface and its wake with an equivalent arrangement of vortices and in finding the associated vertical velocities. By a division of the vortex pattern into circular strips concentric about the downwash point instead of into the usual rectangular strips, the lifting surface is reduced for each downwash point to an equivalent loaded line for which the induced velocity is readily computed. The ratio of the vertical velocity to the stream velocity is the slope of the surface in the free-stream direction. As an illustration, the shape of the wing consistent with the pressure distribution derived from the two-dimensional theories is found for two wings: a straight elliptical wing and one with 30 degree sweepback. Application of the method to solve the reverse problem - finding the lift distribution over a given surface - is briefly discussed.

Cohen, Doris↗

The prediction of noise radiation from supersonic elliptic jets

This paper describes the prediction of noise radiation from supersonic elliptic jets. The noise is associated with the large scale structures in the jet mixing layer. These structures are described as instability waves. The local characteristics of the instability waves are determined from a compressible, linear, analysis. The jet mean velocity and density are described in elliptic cylindrical coordinates. The local eigensolution for the instability waves is determined from a finite difference solution of the non-separable boundary value problem. This inner solution which is formulated in terms of the method of multiple scales is matched with the radiated field using the method of matched asymptotic expansions. The form of the far-field directivity is derived. Predictions are presented for the noise radiation by the several modes of instability in the elliptic jet. The radiated field is not axisymmetric and certain modes radiate strongly in the directions of the major and minor axes of the jet. The extension of the present work to other geometries and flow fields is discussed.

Morris, Philip J.↗

LOOPUS Mob-D: System concept for a public mobile satellite system providing integrated digital services for the Northern Hemisphere from an elliptical orbit

A new concept for a satellite based public mobile communications system, LOOPUS Mob-D, is introduced, whereby most of the classical problems in mobile satellite systems are approached in a different way. The LOOPUS system will offer a total capacity of 6000 high rate channels in three service areas (Europe, Asia, and North America), covering the entire Northern Hemisphere with a set of group special mobile (GSM) compatible mobile services, eventually providing the 'office in the car'. Special characteristics of the LOOPUS orbit and the communications network architecture are highlighted.

Kuhlen, H.↗

A solution of the variational equations for elliptic orbits in rotating coordinates

For elliptic reference orbits, formulas are given for the perturbation state transition matrix of the two-body problem. The formulas relate perturbations expressed in a local vertical rotating coordinate system and are valid for motion in the linear neighborhood of reference orbits with e in the range of 0 to 1. The elements of the state transition matrix are expressed in terms of natural parameters (horizontal and radial velocity, radius, eccentricity, true anomaly, etc.) at the initial and final points. In addition to the general form, a simplified version, valid for small eccentricity orbits, is given.

Jones, J. B.↗

Small Gas Turbine Combustor Primary Zone Study

The combustion research program, small gas turbine combustor primary zone study is summarized. The basic elements of a design methodology program to obtain the maximum performance potential of small reverse-flow annular combustors is described. Three preferred combustion design approaches for internal flame stabilization patterns were selected. Design features are incorporated in the combustors to address the performance limiting problem areas associated with smaller annular combustors. Performance is predicted by using a 3-D aerodynamic/chemical kinetic elliptic flow analysis, initially developed by Garrett Corporation for the USARTL. It is shown that the analytical flow field predictive models provide a useful design tool for understanding the combustion performance of a small reverse flow annular combustor.

Sullivan, R. E.↗

Fluid flow and heat convection studies for actively cooled airframes

The work done during the progress report period from May-October 1992 is summarized. The effect of wall thermal boundary conditions on flows over a step or rib when repeated rib roughness is used for heating augmentation is examined. In numerical investigations of various such laminar and turbulent flows, the local heat transfer coefficients on a forward-facing step or on a rib were found to be very sensitive to the wall thermal boundary condition. For the computation of constant property laminar flow, the wall thermal boundary conditions were either a uniform heat flux or a uniform temperature. Results (Nusselt number and isotherms) of the studies are included. The second part of the work consisted of using PHOENICS to solve the conjugate heat transfer problem of flow over a rib in channel. Finally, the algebraic stress model in the TEAM (Turbulent Elliptic Algorithm-Manchester) code was tested for jet impingement flow, but there needs to be an addition of the energy equation to the code.

Mills, A. F.↗

Canonical forms of multidimensional steady inviscid flows

Canonical forms and canonical variables for inviscid flow problems are derived. In these forms the components of the system governed by different types of operators (elliptic and hyperbolic) are separated. Both the incompressible and compressible cases are analyzed, and their similarities and differences are discussed. The canonical forms obtained are block upper triangular operator form in which the elliptic and non-elliptic parts reside in different blocks. The full nonlinear equations are treated without using any linearization process. This form enables a better analysis of the equations as well as better numerical treatment. These forms are the analog of the decomposition of the one dimensional Euler equations into characteristic directions and Riemann invariants.

Taasan, Shlomo↗

Investigation of a Technique for Measuring Dynamic Ground Effect in a Subsonic Wind Tunnel

To better understand the ground effect encountered by slender wing supersonic transport aircraft, a test was conducted at NASA Langley Research Center's 14 x 22 foot Subsonic Wind Tunnel in October, 1997. Emphasis was placed on improving the accuracy of the ground effect data by using a "dynamic" technique in which the model's vertical motion was varied automatically during wind-on testing. This report describes and evaluates different aspects of the dynamic method utilized for obtaining ground effect data in this test. The method for acquiring and processing time data from a dynamic ground effect wind tunnel test is outlined with details of the overall data acquisition system and software used for the data analysis. The removal of inertial loads due to sting motion and the support dynamics in the balance force and moment data measurements of the aerodynamic forces on the model is described. An evaluation of the results identifies problem areas providing recommendations for future experiments. Test results are validated by comparing test data for an elliptical wing planform with an Elliptical wing planform section with a NACA 0012 airfoil to results found in current literature. Major aerodynamic forces acting on the model in terms of lift curves for determining ground effect are presented. Comparisons of flight and wind tunnel data for the TU-144 are presented.

Graves, Sharon S.↗

Numerical Study of Multigrid Methods with Various Smoothers for the Elliptical Grid Generation Equations

A robust solver for the elliptic grid generation equations is sought via a numerical study. The system of PDEs is discretized with finite differences, and multigrid methods are applied to the resulting nonlinear algebraic equations. Multigrid iterations are compared with respect to the robustness and efficiency. Different smoothers are tried to improve the convergence of iterations. The methods are applied to four 2D grid generation problems over a wide range of grid distortions. The results of the study help to select smoothing schemes and the overall multigrid procedures for elliptic grid generation.

Golik, W. L.↗