Search NASA⌕ Search

SEARCH · Search NASA

Results for “free boundary problem”

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 289 records · Page 16

Modeling and Rapid Prototyping of Integrated Transmission-Distribution OPF Formulations with PowerModelsITD.jl

Conventional electric power systems are composed of different unidirectional power flow stages of generation, transmission, and distribution, managed independently by transmission system and distribution system operators. However, as distribution systems increase in complexity due to the integration of distributed energy resources, coordination between transmission and distribution networks will be imperative for the optimal operation of the power grid. However, coupling models and formulations between transmission and distribution is non-trivial, in particular due to the common practice of modeling transmission systems as single-phase, and distribution systems as multi-conductor phase-unbalanced. To enable the rapid prototyping of power flow formulations, in particular in the modeling of the boundary conditions between these two seemingly incompatible data models, we introduce PowerModelsITD.jl, a free, open-source toolkit written in Julia for integrated transmission-distribution (ITD) optimization that leverages mature optimization libraries from the InfrastructureModels.jl-ecosystem. The primary objective of the proposed framework is to provide baseline implementations of steady-state ITD optimization problems, while providing a common platform for the evaluation of emerging formulations and optimization problems. In this work, we introduce the nonlinear formulations currently supported in PowerModelsITD.jl, which include AC-polar, AC-rectangular, current-voltage, and a linear network transportation model. Results are validated using combinations of IEEE transmission and distribution networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Reentry-Vehicle Shape Optimization Using a Cartesian Adjoint Method and CAD Geometry

A DJOINT solutions of the governing flow equations are becoming increasingly important for the development of efficient analysis and optimization algorithms. A well-known use of the adjoint method is gradient-based shape. Given an objective function that defines some measure of performance, such as the lift and drag functionals, its gradient is computed at a cost that is essentially independent of the number of design variables (e.g., geometric parameters that control the shape). Classic aerodynamic applications of gradient-based optimization include the design of cruise configurations for transonic and supersonic flow, as well as the design of high-lift systems. are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric computer-aided design (CAD). In previous work on Cartesian adjoint solvers, Melvin et al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the two-dimensional Euler equations using a ghost-cell method to enforce the wall boundary conditions. In Refs. 18 and 19, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm were the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The accuracy of the gradient computation was verified using several three-dimensional test cases, which included design variables such as the free stream parameters and the planform shape of an isolated wing. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Factors under consideration include the computation of mesh sensitivities that provide a reliable approximation of the objective function gradient, as well as the computation of surface shape sensitivities based on a direct-CAD interface. We present detailed gradient verification studies and then focus on a shape optimization problem for an Apollo-like reentry vehicle. The goal of the optimization is to enhance the lift-to-drag ratio of the capsule by modifying the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design.

Nemec, Marian↗

On the Solution of the Three-Dimensional Flowfield About a Flow-Through Nacelle

The solution of the three dimensional flow field for a flow through nacelle was studied. Both inviscid and viscous inviscid interacting solutions were examined. Inviscid solutions were obtained with two different computational procedures for solving the three dimensional Euler equations. The first procedure employs an alternating direction implicit numerical algorithm, and required the development of a complete computational model for the nacelle problem. The second computational technique employs a fourth order Runge-Kutta numerical algorithm which was modified to fit the nacelle problem. Viscous effects on the flow field were evaluated with a viscous inviscid interacting computational model. This model was constructed by coupling the explicit Euler solution procedure with a flag entrainment boundary layer solution procedure in a global iteration scheme. The computational techniques were used to compute the flow field for a long duct turbofan engine nacelle at free stream Mach numbers of 0.80 and 0.94 and angles of attack of 0 and 4 deg.

Compton, William Bernard↗

Steady and unsteady incompressible free-wake analysis

The flows around highly sweptback wings and bodies of revolution at high angle of attack are described, and inviscid model approximations and mathematical formulation of the problem are given to steady and unsteady incompressible flows. A general presentation of the methods of solution is given, with emphasis on current computational techniques. Detailed descriptions of the nonlinear vortex-lattice and vortex-panel techniques are presented to show how the boundary conditions are enforced using iteration. Typical numerical results are compared with the available experimental data.

Kandil, O. A.↗

Experimental Investigation of Separated and Transitional Boundary Layers Under Low-Pressure Turbine Airfoil Conditions

Modern low-pressure turbine airfoils are subject to increasingly stronger pressure gradients as designers impose higher loading in an effort to improve efficiency and to reduce part count. The adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation, particularly under cruise conditions. Separation bubbles, notably those which fail to reattach, can result in a significant degradation of engine efficiency. Accurate prediction of separation and reattachment is hence crucial to improved turbine design. This requires an improved understanding of the transition flow physics. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions, has a strong influence on subsequent reattachment, and may even eliminate separation. Further complicating the problem are the high free-stream turbulence levels in a real engine environment, the strong pressure gradients along the airfoils, the curvature of the airfoils, and the unsteadiness associated with wake passing from upstream stages. Because of the complicated flow situation, transition in these devices can take many paths that can coexist, vary in importance, and possibly also interact, at different locations and instances in time. The present work was carried out in an attempt to systematically sort out some of these issues. Detailed velocity measurements were made along a flat plate subject to the same nominal dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil ('Pak-B'). The Reynolds number based on wetted plate length and nominal exit velocity, Re, was varied from 50;000 to 300; 000, covering cruise to takeoff conditions. Low, 0.2%, and high, 7%, inlet free-stream turbulence intensities were set using passive grids. These turbulence levels correspond to about 0.2% and 2.5% turbulence intensity in the test section when normalized with the exit velocity. The Reynolds number and free-stream turbulence level do not have a significant effect on the location of boundary-layer separation unless they are high enough to induce transition upstream of separation. The location and extent of the transition zone, in contrast, depend strongly on Re and TI. The beginning of reattachment closely follows the onset of transition. Under low free-stream turbulence conditions the boundary layer is laminar at separation and then begins to exhibit fluctuations in a finite frequency band in the shear layer over the separation bubble. These fluctuations are due to instability waves. The fluctuations grow in magnitude, higher harmonics are generated, and finally lead to a breakdown to turbulence. Transition begins in the shear layer, but quickly spreads to the near wall region and causes the boundary layer to reattach. The transition is rapid and the resulting turbulence contains a full range of high and low frequencies. Under high free-stream turbulence conditions, slowly growing low-frequency fluctuations are induced in the pretransitional boundary layer by the free-stream. The separation bubbles are considerably thinner than in the low TI cases, resulting in thinner boundary layers at the end of the test wall. At Re=50,000 and 100,000, the pre-transitional boundary layer separates at about the same location as in the low TI cases. Transition occurs through a bypass mode, begins upstream of the corresponding low-TI location, and proceeds in a manner similar to that of an attached boundary layer. Under high TI at Re=200,000 and 300,000, transition begins before separation. The boundary layer may separate, but if it does the separation bubble is very short and does not significantly affect the downstream development of the boundary layer. A comparison is made to previous work in a simulated cascade.

Hultgren, Lennart S.↗

S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the 2 → 2 scattering matrix S 2→2 such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider 3 + 1 dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves k ℓ (s) that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves f ℓ (s), for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of modified aerodynamic strip theories on rotor blade aeroelastic stability

Various existing unsteady aerodynamic strip theories which have been developed in the past for both fixed and rotary wing aeroelastic analyses are modified in the paper so as to make them applicable to the coupled flap-lag-torsional aeroelastic problem of a rotor blade in hover. These corrections are primarily due to constant angle of attack, constant inflow and variable free stream velocity due to lead-lag motion. Next, the modified strip theories are incorporated in a coupled flap-lag-torsional aeroelastic analysis of the rotor blade in hover and the sensitivity of the aeroelastic stability boundaries to the aerodynamic assumptions is examined.

Friedmann, P.↗

Electrostatic disturbance forces on a 3-axis drag-free sensor

The electrostatic analysis of a multiple-capacitance 3-axis drag-free sensor is presented in this paper. The instrument consists of a proof-mass (a dense metallic ball) floating freely inside a spherical cavity enclosed by the sensor plates and the shield. Since the ball and the cavity are not necessarily concentric, the problem in three-dimensional potential theory for electrostatics is solved by the method of boundary perturbations and specifically in terms of spherical harmonics. The capacitance outputs of the instrument and the electrostatic forces acting on the system are derived as non-linear functions of the ball position, ball charge and the sensor plate potentials. The instrument sensitivity and cross-coupling effects are discussed. The analysis may also be useful for electrostatic gyros and suspensions.

Hagstrom, T.↗

A note on the distortion of a flat-plate boundary layer by free-stream vorticity normal to the plate

The purpose of this note is to construct a local solution that eliminates a residual velocity discontinuity in the inviscid portion of a solution obtained in a recent paper by Goldstein, Leib and Cowley (1992). This result is of importance because it shows that the solution obtained in that paper is entirely non-singular outside the viscous wall boundary layer and that any singularity in the problem will have to arise in the usual way through a breakdown in the viscous boundary layer.

Goldstein, M. E.↗

Two-dimensional radiative transfer. I - Planar geometry

Differential-equation methods for solving the transfer equation in two-dimensional planar geometries are developed. One method, which uses a Hermitian integration formula on ray segments through grid points, proves to be extremely well suited to velocity-dependent problems. An efficient elimination scheme is developed for which the computing time scales linearly with the number of angles and frequencies; problems with large velocity amplitudes can thus be treated accurately. A very accurate and efficient method for performing a formal solution is also presented. A discussion is given of several examples of periodic media and free-standing slabs, both in static cases and with velocity fields. For the free-standing slabs, two-dimensional transport effects are significant near boundaries, but no important effects were found in any of the periodic cases studied.

Mihalas, D.↗

On the role of constant-stress surfaces in the problem of minimizing elastic stress concentration

Cases involving antiplane shear deformation, axisymmetric torsion, and plane strain theory, with surfaces of constant stress magnitude optimal in terms of minimizing stress, are investigated. Results for the plane theory refer to exterior doubly connected domains. Stresses generated by torsion of an elastic solid lying within a radially convex region of revolution with plane ends, body force absent, and lateral surface traction-free, are examined. The unknown portion of the boundary of such domains may involve a hole, fillet, or notch.

Wheeler, L.↗

Smoothers for Matrix-Free Algebraic Multigrid Preconditioning of High-Order Finite Elements

We investigate smoothers for use in matrix-free algebraic multigrid (AMG) preconditioning of high-order finite element problems. These AMG preconditioners are matrix-free in the sense that they are built from a related low-order refined finite element problem whose system matrix can be much more rapidly assembled than the high-order problem. Our proposed smoother, which we call distributive relaxation, is more robust to the anisotropy present in many low-order refined meshes which feature a clustering of nodes near the boundaries between high-order finite elements. For solving the low-order refined problem, we show that this new distributive relaxation smoother possesses significantly improved performance compared to more traditional smoothers.

97 MATHEMATICS AND COMPUTING↗

An analytical study of reduced-gravity liquid reorientation using a simplified marker and cell technique

A computer program called HOPI was developed to predict reorientation flow dynamics, wherein liquids move from one end of a closed, partially filled, rigid container to the other end under the influence of container acceleration. The program uses the simplified marker and cell numerical technique and, using explicit finite-differencing, solves the Navier-Stokes equations for an incompressible viscous fluid. The effects of turbulence are also simulated in the program. HOPI can consider curved as well as straight walled boundaries. Both free-surface and confined flows can be calculated. The program was used to simulate five liquid reorientation cases. Three of these cases simulated actual NASA LeRC drop tower test conditions while two cases simulated full-scale Centaur tank conditions. It was concluded that while HOPI can be used to analytically determine the fluid motion in a typical settling problem, there is a current need to optimize HOPI. This includes both reducing the computer usage time and also reducing the core storage required for a given size problem.

Betts, W. S., Jr.↗

Time development of a perturbed-spherical nucleus in a pure supercooled liquid. I - Power-law growth of morphological instabilities. II - Nonlinear development

The linear growth stage of the morphological instabilities developing at the liquid-solid interface during crystal growth from an axisymmetric spherical nucleus is analyzed. The corresponding growth rate parameters are calculated numerically, and it is shown that morphological instabilities for free growth evolve according to a power law in agreement with WKB results and contrary to an exponential law found in a quasi-stationary approximation. The time evolution of an arbitrary perturbed liquid-solid interface from a linear into the nonlinear stage is studied. The initial perturbations include the eigenfunctions for the linear problem, localized perturbations of the grain-boundary type, and a stochastic noise. It is shown that the perturbations grow and spread in a wavelike manner. The formation of a predendritic growth stage is characterized by establishment of constant values of tip radius and velocity.

Pines, V.↗

Large-amplitude finite element flutter analysis of composite panels in hypersonic flow

Consideration is given to a finite-element approach for determining the nonlinear flutter characteristics of 3D thin laminated composite panels using the full third-order-piston transverse-loading aerodynamic theory. The unsteady hypersonic aerodynamic theory and the von Karman large-deflection-plate theory are used to formulate the aeroelasticity problem. Nonlinear flutter analyses are performed to assess the influence of the higher-order aerodynamic theory on the structure's limit-cycle amplitude and the dynamic pressure of the flow velocity. A solution procedure is presented to solve the nonlinear panel flutter and large-amplitude free-vibration finite-element equations. Nonlinear flutter analyses are performed for different boundary support-conditions and for various system parameters. Linear finite-element flutter for isotropic and composite panels and large-amplitude isotropic panel flutter results are compared with existing classical solutions. The large-amplitude panel flutter results using the full third-order piston aerodynamic theory are presented to assess the influence of the nonlinear aerodynamic theory.

Gray, Carl E., Jr.↗

Towards Perfectly Absorbing Boundary Conditions for Euler Equations

In this paper, we examine the effectiveness of absorbing layers as non-reflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordinate directions and introducing absorption coefficients in each split equation. This methodology is similar to that used by Berenger for the numerical solutions of Maxwell's equations. Specifically, we apply this methodology to three physical problems shock-vortex interactions, a plane free shear flow and an axisymmetric jet- with emphasis on acoustic wave propagation. Our numerical results indicate that the use of absorbing layers effectively minimizes numerical reflection in all three problems considered.

Hayder, M. Ehtesham↗

Interlaminar stresses in composite laminates: A perturbation analysis

A general method of solution for an elastic balanced symmetric composite laminate subject to a uniaxial extension was developed based upon a perturbation analysis of a limiting free body containing an interfacial plane. The solution satisfies more physical requirements and boundary conditions than previous investigations, and predicts smooth continuous interlaminar stresses with no instabilities. It determines the finite maximum intensity for the interlaminar normal stress in all laminates, provides mathematical evidences for the singular stresses in angle-ply laminates, suggests the need for the experimental determination of an important problem parameter, and introduces a viable means for solving related problems of practical interest.

Hsu, P. W.↗

Calculation of two-dimensional inlet flow fields by an implicit method including viscous effects: User's manual

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time marching finite difference procedure to solve the thin layer Navier-Stokes equations formulated in body fitted coordinates. The mathematical formulation of the problem and the solution algorithm are given. Numerical stability and accuracy as well as the initial and boundary conditions used are discussed. Instructions for program use and operation along with the overall program logic are also given.

Biringen, S.↗