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

Electromagnetic Scattering From a Polygonal Thin Metallic Plate Using Quadrilateral Meshing

The problem of electromagnetic (EM) scattering from irregularly shaped, thin, metallic flat plates in free space is solved using the electric field integral equation (EFIE) approach in conjunction with the method of moments (MoM) with quadrilateral meshing. An irregularly shaped thin plate is discretized into quadrilateral patches and the unknown electric surface current over the plate is expressed in terms of proper basis functions over these patches. The basis functions for the electric surface current density that satisfy the proper boundary conditions on these quadrilateral patches are derived. The unknown surface current density on these quadrilateral patches is determined by setting up and solving the electric field integral equation by the application of the MoM. From the knowledge of the surface current density, the EM scattering from various irregularly shaped plates is determined and compared with the earlier published results. The novelty in the present approach is the use of quadrilateral patches instead of well known and often used triangular patches. The numerical results obtained using the quadrilateral patches compare favorably with measured results.

Deshpande, Manohar D.↗

Updates to Implicit Edge-Based Gradient Methods

In this paper, we report updates to the implicit edge-based gradient methods originally introduced in [H. Nishikawa, AIAA Paper 2020-3048, 2020]. First, we clarify the relationship between gradient accuracy and truncation error and show that the quadratic method involves a free parameter. Then, we provide a complete description and a simplified matrix form of the implicit gradient systems including a consistent boundary treatment, and derive a set of parameters for achieving fourth-order gradient accuracy on regular tetrahedral grids. A stability analysis is performed for a relaxation scheme used to solve the implicit gradient systems, and the result serves as a guide for choosing parameters. Numerical results are shown for accuracy verification and also for realistic inviscid flow problems in three dimensions, including flows with shock waves.

Weighted Least-Squares↗

Goertler vortices in growing boundary layers: The leading edge receptivity problem, linear growth and the nonlinear breakdown stage

Goertler vortices are thought to be the cause of transition in many fluid flows of practical importance. A review of the different stages of vortex growth is given. In the linear regime, nonparallel effects completely govern this growth, and parallel flow theories do not capture the essential features of the development of the vortices. A detailed comparison between the parallel and nonparallel theories is given and it is shown that at small vortex wavelengths, the parallel flow theories have some validity; otherwise nonparallel effects are dominant. New results for the receptivity problem for Goertler vortices are given; in particular vortices induced by free stream perturbations impinging on the leading edge of the walls are considered. It is found that the most dangerous mode of this type can be isolated and it's neutral curve is determined. This curve agrees very closely with the available experimental data. A discussion of the different regimes of growth of nonlinear vortices is also given. Again it is shown that, unless the vortex wavelength is small, nonparallel effects are dominant. Some new results for nonlinear vortices of 0(1) wavelengths are given and compared to experimental observations.

Hall, Philip↗

Boundary-layer effects in composite laminates. I - Free-edge stress singularities. II - Free-edge stress solutions and basic characteristics

The fundamental nature of the boundary-layer effect in fiber-reinforced composite laminates is formulated in terms of the theory of anisotropic elasticity. The basic structure of the boundary-layer field solution is obtained by using Lekhnitskii's stress potentials (1963). The boundary-layer stress field is found to be singular at composite laminate edges, and the exact order or strength of the boundary layer stress singularity is determined using an eigenfunction expansion method. A complete solution to the boundary-layer problem is then derived, and the convergence and accuracy of the solution are analyzed, comparing results with existing approximate numerical solutions. The solution method is demonstrated for a symmetric graphite-epoxy composite.

Wang, S. S.↗

A shear deformable theory of laminated composite shallow shell-type panels and their response analysis. I - Free vibration and buckling

This paper deals with the substantiation of a shear deformable theory of cross-ply laminated composite shallow shells. While the developed theory preserves all the advantages of the first order transverse shear deformation theory it succeeds in eliminating some of its basic shortcomings. The theory is further employed in the analysis of the eigenvibration and static buckling problems of doubly curved shallow panels. In this context, the state space concept is used in conjunction with the Levy method, allowing one to analyze these problems in a unified manner, for a variety of boundary conditions. Numerical results are presented and some pertinent conclusions are formulated.

Librescu, L.↗

Predicting aerodynamic characteristics of vortical flows on three-dimensional configurations using a surface-singularity panel method

A general low-order surface-singularity panel method is used to predict the aerodynamic characteristics of a problem where a wing-tip vortex from one wing closely interacts with an aft mounted wing in a low Reynolds Number flow; i.e., 125,000. Nonlinear effects due to wake roll-up and the influence of the wings on the vortex path are included in the calculation by using a coupled iterative wake relaxation scheme. The interaction also affects the wing pressures and boundary layer characteristics: these effects are also considered using coupled integral boundary layer codes and preliminary calculations using free vortex sheet separation modelling are included. Calculated results are compared with water tunnel experimental data with generally remarkably good agreement.

Maskew, B.↗

Optimal discrete-time LQR problems for parabolic systems with unbounded input: Approximation and convergence

An abstract approximation and convergence theory for the closed-loop solution of discrete-time linear-quadratic regulator problems for parabolic systems with unbounded input is developed. Under relatively mild stabilizability and detectability assumptions, functional analytic, operator techniques are used to demonstrate the norm convergence of Galerkin-based approximations to the optimal feedback control gains. The application of the general theory to a class of abstract boundary control systems is considered. Two examples, one involving the Neumann boundary control of a one-dimensional heat equation, and the other, the vibration control of a cantilevered viscoelastic beam via shear input at the free end, are discussed.

Rosen, I. G.↗

The NASA Low-Pressure Turbine Flow Physics Program

An overview of the NASA Lewis Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The program was established in response to the aero-engine industry's need for improved LPT efficiency and designs. Modern jet engines have four to seven LPT stages, significantly contributing to engine weight. In addition, there is a significant efficiency degradation between takeoff and cruise conditions, of up to 2 points. Reducing the weight and part count of the LPT and minimizing the efficiency degradation will translate into fuel savings. Accurate prediction methods of LPT flows and losses are needed to accomplish those improvements. The flow in LPT passages is at low Reynolds number, and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The affecting parameters traditionally considered are Reynolds number, freestream turbulence intensity, wake frequency parameter, and the pressure distribution (loading). Three-dimensional effects and additional parameters, particularly turbulence characteristics like length scales, spectra and other statistics, as well as wake turbulence intensity and properties also play a role. The flow of most interest is on the suction surface, where large losses are generated as the flow tends to separate at the low Reynolds numbers. Ignoring wakes, a common flow scenario, there is laminar separation, followed by transition on the separation bubble and turbulent reattachment. If transition starts earlier the separation will be eliminated and the boundary layer will be attached leading to the well known bypass transition issues. In contrast, transition over a separation bubble is closer to free shear layer transition and was not investigated as well, particularly in the turbine environment. Unsteadiness created by wakes complicates the picture. Wakes induce earlier transition, and the calmed regions trailing the induced turbulent spots can delay or eliminate separation via shear stress modification. Three-dimensional flow physics and geometry will have strong effects. Altogether a very complex and challenging problem emerges. The objective of the program is to provide improved models and physical understanding of the complex flow, which are essential for accurate prediction of flow and losses in the LPT. Experimental, computational and analytical work as complementing and augmenting approaches are used. The program involves industry, universities and research institutes, and other government laboratories. It is characterized by strong interaction among participants, quick dissemination of results, and responsiveness to industry's needs. The presentation will describe the work elements. Highlighting some activities in progress are experiments on simulated blade suction surface in low-speed wind tunnels, on curved wall, and on a flat-plate, both with pressure gradient. In the area of computation, assessment of existing models is performed using RANS (Reynolds Averaged Navier Stokes) simulations. Laminar flow DNS was completed. Analytical studies of instability and receptivity in attached and separated flows were started. In the near future the program is moving to include wake effects and development of improved modeling. Experimental work in preparation stages are: (1) Addition of wakes to the curved tunnel experiment; (2) Low-speed rotating rig experiment on GE90 engine LPT; and (3) Transonic cascade. In the area of computation, it is expected to move from model assessment towards development of improved models. In addition, a new project of Large Eddy Simulation (LES) of LPT is to begin and will provide numerical data bases. It is planned to implement the emerging improved models in a multistage turbomachinery code and to validate against the GE90 engine LPT.

Ashpis, David E.↗

Investigation of a Smooth Local Correlation-based Transition Model in a Discrete-Adjoint Aerodynamic Shape Optimization Algorithm

A smooth local correlation-based transition model is fully coupled to a RANS-based Newton-Krylov flow solver and discrete-adjoint gradient-based optimization algorithm. The free-transition optimization framework is evaluated using lift-constrained drag minimizations of airfoils at design conditions ranging from light to single-aisle aircraft and an infinite swept wing at design conditions representative of a transonic strut-braced wing aircraft. The impact of the streamwise grid resolution on the ability of the optimization algorithm to delay boundary-layer transition is investigated, with the results demonstrating that streamwise grid resolution requirements increase as the transition length decreases with increasing Reynolds number. The optimization problem at the light aircraft design conditions is demonstrated to be multi-modal, with the optimization algorithm producing two distinct designs: one with a thin, reflexed trailing edge and steep pressure recovery regions, the other with increased aft loading, with the latter design outperforming the former. A drag minimization of an airfoil at transonic design conditions demonstrates that the optimization algorithm successfully trades a decrease in viscous drag by delaying boundary-layer transition with an increase in wave drag, while the drag minimization of an infinite swept wing demonstrates the capability of the optimizational gorithm to delay both Tollmien-Schlichting and stationary crossflow instabilities.

AATT↗

A New Higher-Order Composite Theory for Analysis and Design of High Speed Tilt-Rotor Blades

A higher-order theory is developed to model composite box beams with arbitrary wall thicknesses. The theory, based on a refined displacement field, represents a three-dimensional model which approximates the elasticity solution. Therefore, the cross-sectional properties are not reduced to one-dimensional beam parameters. Both inplane and out-of-plane warping are automatically included in the formulation. The model accurately captures the transverse shear stresses through the thickness of each wall while satisfying all stress-free boundary conditions. Several numerical results are presented to validate the present theory. The developed theory is then used to model the load carrying member of a tilt-rotor blade which has thick-walled sections. The composite structural analysis is coupled with an aerodynamic analysis to compute the aeroelastic stability of the blade. Finally, a multidisciplinary optimization procedure is developed to improve the aerodynamic, structural and aeroelastic performance of the tilt-rotor aircraft. The Kreisselmeier-Steinhauser function is used to formulate the multiobjective function problem and a hybrid approximate analysis is used to reduce the computational effort. The optimum results are compared with the baseline values and show significant improvements in the overall performance of the tilt-rotor blade.

McCarthy, Thomas Robert↗

A space-time tensor formulation for continuum mechanics in general curvilinear, moving, and deforming coordinate systems

Tensor methods are used to express the continuum equations of motion in general curvilinear, moving, and deforming coordinate systems. The space-time tensor formulation is applicable to situations in which, for example, the boundaries move and deform. Placing a coordinate surface on such a boundary simplifies the boundary condition treatment. The space-time tensor formulation is also applicable to coordinate systems with coordinate surfaces defined as surfaces of constant pressure, density, temperature, or any other scalar continuum field function. The vanishing of the function gradient components along the coordinate surfaces may simplify the set of governing equations. In numerical integration of the equations of motion, the freedom of motion of the coordinate surfaces provides a potential for enhanced resolution of the continuum field function. An example problem of an incompressible, inviscid fluid with a top free surface is considered, where the surfaces of constant pressure (including the top free surface) are coordinate surfaces.

Avis, L. M.↗

Two-dimensional hybrid continuum/particle approach for rarefied flows

A hybrid numerical technique previously developed for one-dimensional rarefied gas flows is generalized to two dimensions. The method is based on the fact that the flowfield that develops near a body in a rarefied gas typically contains local regions of continuum, transitional and free-molecular flow. By utilizing the solution technique most appropriate for each region and coupling the techniques in an interface region where both are applicable, more computationally efficient solutions can be obtained, or equivalently, more complex flowfields can be analyzed. The present method combines finite difference solution of the Navier-Stokes equations in the continuum regions, with Direct Simulation Monte Carlo in the more rarefied regions. The two schemes are coupled interactively via a general conservative flux boundary condition. The method is tested by application to the model problem of pressure-driven rarefied flow through a slit. Results show the hybrid scheme offers a speedup of a factor of nearly two for the nominal conditions considered, due to the decrease in the size of the Direct Simulation domain.

Wadsworth, Dean C.↗

Vertical Profiling of Air Pollution at RAPCD

The interaction between local and regional pollution levels occurs at the interface of the Planetary Boundary Layer and the Free Troposphere. Measuring the vertical distribution of ozone, aerosols, and winds with high temporal and vertical resolution is essential to diagnose the nature of this interchange and ultimately for accurately forecasting ozone and aerosol pollution levels. The Regional Atmospheric Profiling Center for Discovery, RAPCD, was built and instrumented to address this critical issue. The ozone W DIAL lidar, Nd:YAG aerosol lidar, and 2.1 micron Doppler wind lidar, along with balloon- borne ECC ozonesondes form the core of the W C D instrumentation for addressing this problem. Instrumentation in the associated Mobile Integrated Profiling (MIPS) laboratory includes 91 5Mhz profiler, sodar, and ceilometer. The collocated Applied particle Optics and Radiometry (ApOR) laboratory hosts an FTIR along with MOUDI and optical particle counters. With MODELS-3 analysis by colleagues in the National Space Science and Technology Center on the UAH campus and the co- located National Weather Service Forecasting Office in Huntsville, AL we are developing a unique facility for advancing the state of the science of pollution forecasting.

Newchurch, Michael J.↗

Vibration isolation technology: Sensitivity of selected classes of experiments to residual accelerations

The solution was sought of a 2-D axisymmetric moving boundary problem for the sensitivity of isothermal and nonisothermal liquid columns and the sensitivity of thermo-capillary flows to buoyancy driven convection caused by residual accelerations. The sensitivity of a variety of space experiments to residual accelerations are examined. In all the cases discussed, the sensitivity is related to the dynamic response of a fluid. In some cases the sensitivity can be defined by the magnitude of the response of the velocity field. This response may involve motion of the fluid associated with internal density gradients, or the motion of a free liquid surface. For fluids with internal density gradients, the type of acceleration to which the experiment is sensitive will depend on whether buoyancy driven convection must be small in comparison to other types of fluid motion (such as thermocapillary flow), or fluid motion must be suppressed or eliminated (such as in diffusion studies, or directional solidification experiments). The effect of the velocity on the composition and temperature field must be considered, particularly in the vicinity of the melt crystal interface. As far as the response to transient disturbances is concerned the sensitivity is determined by both the magnitude and frequency the acceleration and the characteristic momentum and solute diffusion times.

Alexander, J. Iwan D.↗

A numerical solution method for acoustic radiation from axisymmetric bodies

A new and very efficient numerical method for solving equations of the Helmholtz type is specialized for problems having axisymmetric geometry. It is then demonstrated by application to the classical problem of acoustic radiation from a vibrating piston set in a stationary infinite plane. The method utilizes 'Green's Function Discretization', to obtain an accurate resolution of the waves using only 2-3 points per wave. Locally valid free space Green's functions, used in the discretization step, are obtained by quadrature. Results are computed for a range of grid spacing/piston radius ratios at a frequency parameter, omega R/c(sub 0), of 2 pi. In this case, the minimum required grid resolution appears to be fixed by the need to resolve a step boundary condition at the piston edge rather than by the length scale imposed by the wave length of the acoustic radiation. It is also demonstrated that a local near-field radiation boundary procedure allows the domain to be truncated very near the radiating source with little effect on the solution.

Caruthers, John E.↗

Supercritical wing sections II: A handbook

The numerical aspect of theoretical work on transonics and supercritical wing sections are compiled. A model of the trailing edge is introduced which eliminates the loss of 15 to 20 percent experienced with heavily aft-loaded models, and it is indicated how drag creep can be reduced at off-design conditions. A rotated finite difference scheme is presented which can handle supersonic as well as subsonic free stream Mach numbers and leads to an effective three-dimensional program for the computation of transonic flow past an oblique wing. In the case of two-dimensional flow, the method is extended to take into account the displacement thickness computed by a semiempirical turbulent boundary layer correction. A series of supercritical wing sections is discussed together with comparisons between experimental and theoretical data. Computer programs and a brief manual for their operation are listed. It is shown that the programs furnish a physically adequate computer simulation of the compressible flows that arise in problems of transonic aerodynamics.

Bauer, F.↗

Entropy Splitting for High Order Numerical Simulation of Compressible Turbulence

A stable high-order numerical scheme for direct numerical simulation (DNS) of shock-free compressible turbulence is presented. The method is applicable to general geometries. It contains no upwinding, artificial dissipation, or filtering. Instead the method relies on the stabilizing mechanisms of an appropriate conditioning of the governing equations and the use of compatible spatial difference operators for the interior points (interior scheme) as well as the boundary points (boundary scheme). An entropy-splitting approach splits the inviscid flux derivatives into conservative and nonconservative portions. The spatial difference operators satisfy a summation-by-parts condition, leading to a stable scheme (combined interior and boundary schemes) for the initial boundary value problem using a generalized energy estimate. A Laplacian formulation of the viscous and heat conduction terms on the right hand side of the Navier–Stokes equations is used to ensure that any tendency to odd–even decoupling associated with central schemes can be countered by the fluid viscosity. The resulting methods are able to minimize the spurious high-frequency oscillations associated with pure central schemes, especially for long time integration applications such as DNS. For validation purposes, the methods are tested in a DNS of compressible turbulent plane channel flow at low values of friction Mach number, where reference turbulence data bases exist. It is demonstrated that the methods are robust in terms of grid resolution, and in good agreement with published channel data. Accurate turbulence statistics can be obtained with moderate grid sizes. Stability limits on the range of the splitting parameter are determined from numerical tests.

N.D. Sandham↗

Initialization of Formation Flying Using Primer Vector Theory

In this paper, we extend primer vector analysis to formation flying. Optimization of the classical rendezvous or free-time transfer problem between two orbits using primer vector theory has been extensively studied for one spacecraft. However, an increasing number of missions are now considering flying a set of spacecraft in close formation. Missions such as the Magnetospheric MultiScale (MMS) and Leonardo-BRDF (Bidirectional Reflectance Distribution Function) need to determine strategies to transfer each spacecraft from the common launch orbit to their respective operational orbit. In addition, all the spacecraft must synchronize their states so that they achieve the same desired formation geometry over each orbit. This periodicity requirement imposes constraints on the boundary conditions that can be used for the primer vector algorithm. In this work we explore the impact of the periodicity requirement in optimizing each spacecraft transfer trajectory using primer vector theory. We first present our adaptation of primer vector theory to formation flying. Using this method, we then compute the AV budget for each spacecraft subject to different formation endpoint constraints.

Mailhe, Laurie↗