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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Application of the Discontinuous Galerkin Method to Acoustic Scatter Problems

The application of the quadrature-free form of the discontinuous Galerkin method to two problems from Category 1 of the Second Computational Aeroacoustics Workshop on Benchmark problems is presented. The method and boundary conditions relevant to this work are described followed by two test problems, both of which involve the scattering of an acoustic wave off a cylinder. The numerical test performed to evaluate mesh-resolution requirements and boundary-condition effectiveness are also described.

Atkins, H. L.↗

NASA NTRS2003

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows: Minimization of Delta * B Numerical Error - II

The generalization of a class of low-dissipative high order filter finite difference schemes for long time wave propagation of shock/turbulence/combustion compressible viscous gas dynamic flows to compressible MHD equations for structured curvilinear grids has been developed. The new scheme consists of a divergence free preserving high order spatial base scheme with a filter approach which can be divergence-free preserving depending on the type of filter operator being used, the method of applying the filter step, and the type of flow problem to be considered. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (Delta * B) numerical error in the sense that no standard divergence cleaning is required. Performance evaluation of these variants, and the key role that the proper treatment of their corresponding numerical boundary conditions can play will be illustrated. Many levels of grid refinement and detailed comparison with several commonly used compressible MHD shock-capturing schemes will be sought. For certain MHD 2-D test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Sjoegreen, Bjoern↗

NASA NTRS1986

Boundary integral techniques for multi-connected domains

Attention is given to iterative methods applicable to such changing domain cases of elliptic problems in multiconnected domains as those of incompressible inviscid flow with free surfaces. One such formulation is developed and tested on circular and elliptic annuli, using interpolated quadrature points to maintain accuracy when annuli regions are thin. The technique is noted to be applicable to the study of accelerating thin fluid shells.

Baker, G. R.↗

NASA NTRS1959

Experimental and Theoretical Studies of Axisymmetric Free Jets

Some experimental and theoretical studies have been made of axisymmetric free jets exhausting from sonic and supersonic nozzles into still air and into supersonic streams with a view toward problems associated with propulsive jets and the investigation of these problems. For jets exhausting into still air, consideration is given to the effects of jet Mach number, nozzle divergence angle, and jet static pressure ratio upon jet structure, jet wavelength, and the shape and curvature of the jet boundary. Studies of the effects of the ratio of specific heats of the jets are included are observations pertaining to jet noise and jet simulation. For jets exhausting into supersonic streams, an attempt has been made to present primarily theoretical certain jet interference effects and in formulating experimental studies. The primary variables considered are jet Mach number, free stream Mach number, jet static pressure ratio, ratio of specific heats of the jet, nozzle exit angle, and boattail angle. The simulation problem and the case of a hypothetical hypersonic vehicle are examined, A few experimental observations are included.

Love, Eugene S.↗

NASA NTRS1991

A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading

An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.

Norwood, D. S.↗

NASA NTRS1994

On the transition towards slow manifold in shallow-water and 3D Euler equations in a rotating frame

The long-time, asymptotic state of rotating homogeneous shallow-water equations is investigated. Our analysis is based on long-time averaged rotating shallow-water equations describing interactions of large-scale, horizontal, two-dimensional motions with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. These equations are obtained in two steps: first by introducing a Poincare/Kelvin linear propagator directly into classical shallow-water equations, then by averaging. The averaged equations describe interaction of wave fields with large-scale motions on time scales long compared to the time scale 1/f(sub o) introduced by rotation (f(sub o)/2-angular velocity of background rotation). The present analysis is similar to the one presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, since three-wave interactions in rotating shallow-water equations are forbidden, the final equations describing the asymptotic state are simplified considerably. Special emphasis is given to a new conservation law found in the asymptotic state and decoupling of the dynamics of the divergence free part of the velocity field. The possible rising of a decoupled dynamics in the asymptotic state is also investigated for homogeneous turbulence subjected to a background rotation. In our analysis we use long-time expansion, where the velocity field is decomposed into the 'slow manifold' part (the manifold which is unaffected by the linear 'rapid' effects of rotation or the inertial waves) and a formal 3D disturbance. We derive the physical space version of the long-time averaged equations and consider an invariant, basis-free derivation. This formulation can be used to generalize Waleffe's (1991) helical decomposition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on the cylinder surface and homogeneous in the vertical direction).

Mahalov, A.↗

NASA NTRS1989

A free surface vortex method with weak viscous effects

The advantage of vortex methods for inviscid flows of several fluids separated by interfaces is that problems of this kind may be posed and solved numerically entirely on the interface surface. For free surface problems, where one of the fluids has negligible inertia, a small amount of viscosity produces a thin, weak rotational layer adjacent to the potential flow vortex sheet. Here, a method is described for including this weak viscous effect as a modification of the boundary conditions for the potential part of the flow, thus allowing this slightly viscous flow to be formulated as a vortex method on the free surface. This is important for unsupported free drops where the viscosity at the free surface is the only source of damping for oscillations, and for gas-filled bubbles, where it provides both drag and damping. Small applied shear stresses from aerodynamic, electrostatic, or other sources may be included provided they are of the same order as the potential flow shear stress or smaller.

Lundgren, T. S.↗

NASA NTRS1977

A numerical method for the exact calculation of airloads associated with impulsively started wings

A numerical method is developed to calculate three-dimensional potential flows due to the steady and impulsive motion of isolated wing and wing-wing interaction problems. The velocity potential is represented by a discrete set of constant-doublet quadrilaterals on wing and wake surfaces. The exact surface boundary condition is enforced, and the solution is obtained in a step-by-step fashion, configurations being impulsively started from rest. Free-wake geometries are generated for each time step with Rankine or Lamb viscous vortex segments used in wake-velocity calculations. Sample results include calculated performance to steady state for a thick wing and indicial lift of a wing-wing interaction problem.

Summa, J. M.↗

NASA NTRS1974

Free convection over a vertical porous plate with transpiration

The problem of free convection over an isothermal vertical porous plate with transpiration is studied both numerically and experimentally. Numerical solutions to the variable-property transpired free-convection boundary layer equations have been obtained using the finite difference procedure of Patankar and Spalding (1967). The effects of uniform transpiration on heat transfer and on temperature and velocity profiles are predicted. Interferometrically measured nondimensional temperature profiles for the uniform wall temperature and transpiration case agreed closely with these numerical predictions.

Parikh, P. G.↗

NASA NTRS1977

Three-dimensional low Reynolds number flows with a free surface

The two-dimensional leveling problem (Degani, Gutfinger, 1976) is extended to three dimensions in the case where the flow Re number is very low and attention is paid to the free surface boundary condition with surface tension effects included. The no-slip boundary condition on the wall is observed. The numerical solution falls back on the Marker and Cell (MAC) method (Harlow and Welch, 1965) with the computation region divided into a finite number of stationary rectangular cells (or boxes in the 3-D case) and fluid flow traverses the cells (or boxes).

Degani, D.↗

NASA NTRS1982

A variational approach to the study of capillary phenomena

The problem of determining the free surface of a liquid in a capillary tube, and of a liquid drop, sitting first on a horizontal plane and then on more general surfaces is considered. With some modifications, the method applies to the study of pendent drops and of rotating drops as well. The standard capillary problem, i.e. the determination of the free surface of a liquid in a thin tube of general cross section, which resuls from the simultaneous action of surface tension, boundary adhesion and gravity is discussed. It turns out that in this case the existence of the solution surface depends heavily on the validity of a simple geometric condition about the mean curvature of the boundary curve of the cross section of the capillary tube. Some particular examples of physical interest are also be discussed. Liquid drops sitting on or hanging from a fixed horizontal plane are discussed. The symmetry of the solutions (which can actually be proved, as consequence of a general symmetrization argument) now plays the chief role in deriving both the existence and the regularity of energy-minimizing configurations. When symmetry fails (this is the case, for example, when the contact angle between the drop and the plate is not constant, or when the supporting surface is not itself symmetric), then more sophisticated methods must be used. Extensions in this direction are outlined.

Emmer, M.↗

NASA NTRS1997

A Comparison of PETSC Library and HPF Implementations of an Archetypal PDE Computation

Two paradigms for distributed-memory parallel computation that free the application programmer from the details of message passing are compared for an archetypal structured scientific computation a nonlinear, structured-grid partial differential equation boundary value problem using the same algorithm on the same hardware. Both paradigms, parallel libraries represented by Argonne's PETSC, and parallel languages represented by the Portland Group's HPF, are found to be easy to use for this problem class, and both are reasonably effective in exploiting concurrency after a short learning curve. The level of involvement required by the application programmer under either paradigm includes specification of the data partitioning (corresponding to a geometrically simple decomposition of the domain of the PDE). Programming in SPAM style for the PETSC library requires writing the routines that discretize the PDE and its Jacobian, managing subdomain-to-processor mappings (affine global- to-local index mappings), and interfacing to library solver routines. Programming for HPF requires a complete sequential implementation of the same algorithm, introducing concurrency through subdomain blocking (an effort similar to the index mapping), and modest experimentation with rewriting loops to elucidate to the compiler the latent concurrency. Correctness and scalability are cross-validated on up to 32 nodes of an IBM SP2.

Hayder, M. Ehtesham↗

NASA NTRS1989

Element-specific modal formulations for large-displacement multibody dynamics

Large dispacement assumed-mode modeling techniques are examined in the context of multibody elastodynamics. The range of both general and element-specific approaches are studied with the aid of examples involving beams, plates, and shells. For systems undergoing primarily structural bending and twisting with little or no membrane distortion, it is found that fully-linear, element-specific, modal formulations provide the most accurate time history solutions at the least expense. When membrane effects become dominant in structural problems due to loading and boundary conditions, one must naturally resort to a formulation involving a nonlinear stress-strain relationship in addition to nonlinear terms associated with large overall system motion. Such nonlinear models were investigated using assumed modes and found to lead to modal convergence difficulties when standard free-free structural modes are employed. A constrained mode formulation aimed at addressing the convergence problem is proposed.

Ryan, R. R.↗

NASA NTRS1994

Physical description of boundary-layer transition: Experimental evidence

The problems of understanding the origins of turbulent flow and transition to turbulent flow are the most important unsolved problems of fluid mechanics and aerodynamics. It is well known that the stability, transition, and turbulent characteristics of bounded shear layers are fundamentally different from those of free shear layers. Likewise, the stability, transition, and turbulent characteristics of open systems are fundamentally different from those of closed systems. Because of the influence of indigenous disturbances, surface geometry and roughness, sound, heat transfer, and ablation, it is not possible to develop general prediction schemes for transition location and the nature of turbulent structures in boundary-layer flows. At the present time no mathematical model exists that can predict the transition Reynolds number on a flat plate. The recent progress in this area is encouraging, in that a number of distinct transition mechanisms have been found experimentally. The theoretical work finds them to be amplitude and Reynolds-number dependent. The theory remains rather incomplete with regard to predicting transition. Amplitude and spectral characteristics of the disturbances inside the laminar viscous layer strongly influence which type of transition occurs. The major need in this area is to understand how freestream disturbances are entrained into the boundary layer, i.e., to answer the question of receptivity. We refer receptivity to the mechanism(s) that cause freestream disturbances to enter the boundary layer and create the initial amplitudes for unstable waves.

Saric, William S.↗

DOE OSTI2022

Machine learning moment closure models for the radiative transfer equation I: Directly learning a gradient based closure

In this paper, we take a data-driven approach and apply machine learning to the moment closure problem for the radiative transfer equation in slab geometry. Instead of learning the unclosed high order moment, we propose to directly learn the gradient of the high order moment using neural networks. This new approach is consistent with the exact closure we derive for the free streaming limit and also provides a natural output normalization. Finally, a variety of benchmark tests, including the variable scattering problem, the Gaussian source problem with both periodic and reflecting boundaries, and the two-material problem, show both good accuracy and generalizability of our machine learning closure model.

Computer science↗

NASA NTRS1991

Observations on practical turbulence modeling for high-speed jet/plume flowfields

An overview is presented of the current status in terms of simulating basic high-speed/straightback free jet flowfields with two-equation turbulence models. A building-block approach is described which provides a systematic method of calibrating such turbulence models. Data to support this approach focuses on obtaining turbulence modeling capabilities to reproduce basic mean flow structure for simple flows before progressing to more complex problems. For transverse jets and wall jets, new building-blocks of data are necessary with turbulence modeling capabilities unifying the free jet building-blocks discussed, and fundamental boundary layer and shock/boundary layer interaction building-blocks resulting from other studies.

Dash, Sanford M.↗

NASA NTRS1985

A study of internal and distributed damping for vibrating turbomachiner blades

Internal and distributed damping as possible methods for reducing the vibration response of turbomachine blades and theoretical methods for analyzing damped vibration were studied. It is demonstrated how the Ritz-Galerkin methods may be used to straightforwardly to analyze forced vibrations with damping. This is done directly without requiring the free vibration eigenfunctions. The Galerkin method is an effective technique for these types of problems. The Ritz method has the further advantage of not needing to satisfy the force type boundary conditions, which is particularly important for plates and shells. But proper functionals representing the forcing and damping terms must be developed, and this is done. Two types of damping--viscous and material (hysteretic) are considered. Both distributed and concentrated exciting forces are treated. Numerical results are obtained for cantilevered beams and rectangular plates. Studies showing the rates of convergence of the solutions are made. In the case of the cantilever beam, approximate solutions from the present methods are compared with the exact solutions.

Leissa, A. W.↗

NASA NTRS2006

Airframe Noise Prediction by Acoustic Analogy: Revisited

The present work follows a recent survey of airframe noise prediction methodologies. In that survey, Lighthill s acoustic analogy was identified as the most prominent analytical basis for current approaches to airframe noise research. Within this approach, a problem is typically modeled with the Ffowcs Williams and Hawkings (FW-H) equation, for which a geometry-independent solution is obtained by means of the use of the free-space Green function (FSGF). Nonetheless, the aeroacoustic literature would suggest some interest in the use of tailored or exact Green s function (EGF) for aerodynamic noise problems involving solid boundaries, in particular, for trailing edge (TE) noise. A study of possible applications of EGF for prediction of broadband noise from turbulent flow over an airfoil surface and the TE is, therefore, the primary topic of the present work. Typically, the applications of EGF in the literature have been limited to TE noise prediction at low Mach numbers assuming that the normal derivative of the pressure vanishes on the airfoil surface. To extend the application of EGF to higher Mach numbers, the uniqueness of the solution of the wave equation when either the Dirichlet or the Neumann boundary condition (BC) is specified on a deformable surface in motion. The solution of Lighthill s equation with either the Dirichlet or the Neumann BC is given for such a surface using EGFs. These solutions involve both surface and volume integrals just like the solution of FW-H equation using FSGF. Insight drawn from this analysis is evoked to discuss the potential application of EGF to broadband noise prediction. It appears that the use of a EGF offers distinct advantages for predicting TE noise of an airfoil when the normal pressure gradient vanishes on the airfoil surface. It is argued that such an approach may also apply to an airfoil in motion. However, for the prediction of broadband noise not directly associated with a trailing edge, the use of EGF does not appear to offer any advantages over the use of FSGF at the present stage of development. It is suggested here that the applications of EGF for airframe noise analysis be continued. As an example pertinent to airframe noise prediction, the Fast Scattering Code of NASA Langley is utilized to obtain the EGF numerically on the surface of a three dimensional wing with a flap and leading edge slat in uniform rectilinear motion. The interpretation and use of these numerical Green functions are then discussed.

Farassat, F.↗