Search NASA⌕ Search

SEARCH · Search NASA

Results for “internal and external domain boundaries”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Minkowski difference-based advancing front packing technique for generating convex noncircular particles in complex domains

In this work, a Minkowski difference-based advancing front approach is proposed to generate convex and non-circular particles in a predefined computational domain. Two specific algorithms are developed to handle the contact conformity of generated particles with the boundaries of the computational domain. The first, called the open form, is used to handle the smooth contact of generated particles with (external) boundaries, while the other, called the closed form, is proposed to handle the internal boundaries of a computational domain with a complex cavity. The Gilbert-Johnson-Keerthi (GJK) method is used to efficiently solve the contact detection between the newly generated particle at the front and existing particles. Furthermore, the problem of one-sided particle lifting, which can cause some defects in the packing structure in existing advancing front methods during packing generation, is highlighted and an effective solution is developed. Several examples of increasing complexity are used to demonstrate the efficiency and applicability of the proposed packing generation approach. The numerical results show that the generated packing is not only more uniform, but also achieves a higher packing density than existing advancing front methods.

42 ENGINEERING↗
NASA NTRS1997

Advances in Numerical Boundary Conditions for Computational Aeroacoustics

Advances in Computational Aeroacoustics (CAA) depend critically on the availability of accurate, nondispersive, least dissipative computation algorithm as well as high quality numerical boundary treatments. This paper focuses on the recent developments of numerical boundary conditions. In a typical CAA problem, one often encounters two types of boundaries. Because a finite computation domain is used, there are external boundaries. On the external boundaries, boundary conditions simulating the solution outside the computation domain are to be imposed. Inside the computation domain, there may be internal boundaries. On these internal boundaries, boundary conditions simulating the presence of an object or surface with specific acoustic characteristics are to be applied. Numerical boundary conditions, both external or internal, developed for simple model problems are reviewed and examined. Numerical boundary conditions for real aeroacoustic problems are also discussed through specific examples. The paper concludes with a description of some much needed research in numerical boundary conditions for CAA.

Tam, Christopher K. W.↗
NASA NTRS2024

Disentangling the Entangled Linkages of Magnetic Helicity

Magnetic helicity, H, measures magnetic linkages in a volume. The early theoretical development of helicity focused on magnetically closed systems in V bounded by ∂V. For magnetically closed systems, V ∈ R -3 = V + V∗, no magnetic flux threads the boundary, n̂ · B|∂V = 0. Berger & Field (1984) and Finn & Antonsen (1985) extended the definition of helicity to relative helicity, H, for magnetically open systems where magnetic flux may thread the boundary. Berger (1999, 2003) expressed this relative helicity as two gauge invariant terms that describe the self helicity of magnetic field that closes inside V and the mutual helicity between the magnetic field that threads the boundary ∂V and the magnetic field that closes inside V. The total magnetic field that permeates V entangles magnetic fields that are produced by current sources J in V with magnetic fields that are produced by current sources J∗ in V∗. Building on this fact, we extend Berger’s expressions for relative magnetic helicity to eight gauge invariant quantities that simultaneously characterize both of these self and mutual helicities and attribute their origins to currents J in V and/or J∗ in V∗, thereby disentangling the domain of origin for these entangled linkages. We arrange these eight terms into novel expressions for internal and external helicity (self) and internal-external helicity (mutual) based on their domain of origin. The implications of these linkages for interpreting magnetic energy is discussed and new boundary observables are proposed for tracking the evolution of the field that threads the boundary.

Analytical mathematics↗
NASA NTRS1992

Non-local sub-characteristic zones of influence in unsteady interactive boundary-layers

The properties of incompressible, unsteady, interactive, boundary layers are examined for a model hypersonic boundary layer and internal flow past humps or, equivalently, external flow past short-scaled humps. Using a linear high frequency analysis, it is shown that the domains of dependence within the viscous sublayer may be a strong function of position within the sublayer and may be strongly influenced by the pressure displacement interaction, or the prescribed displacement condition. Detailed calculations are presented for the hypersonic boundary layer. This effect is found to carry over directly to the fully viscous problem as well as the nonlinear problem. In the fully viscous problem, the non-local character of the domains of dependence manifests itself in the sub-characteristics. Potential implications of the domain of dependence structure on finite difference computations of unsteady boundary layers are briefly discussed.

Rothmayer, A. P.↗
NASA NTRS1974

Nonlinear panel response by a Monte Carlo approach

The vibration of clamped and simply supported elastic panels due to subsonic and supersonic turbulent boundary-layer flows is investigated by a Monte Carlo technique. The resulting generalized random forces are simulated numerically from boundary-layer turbulence spectra, and the response analysis is performed in the time domain. The mutual interaction between panel motion and external and/or internal airflow is included. Response studies are performed with respect to rms response, probability structure, peak distribution, threshold crossing and spectral density. The effect on the response statistics of in-plane loading, static pressure differential, and cavity pressure is investigated.

Vaicaitis, R.↗
NASA NTRS1988

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
NASA NTRS1992

Enhancements to the GRIDGEN structured grid generation system for internal and external flow applications

GRIDGEN is a government domain software package for interactive generation of multiple block grids around general configurations. Though it has been freely available since 1989, it has not been widely embraced by the internal flow community due to a misconception that it was designed for external flow use only. In reality GRIDGEN has always worked for internal flow applications, and GRIDGEN ongoing enhancements are increasing the quality of and efficiency with which grids for external and internal flow problems may be constructed. The software consists of four codes used to perform the four steps of the grid generation process. GRIDBLOCK is first used to decompose the flow domain into a collection of component blocks and then to establish interblock connections and flow solver boundary conditions. GRIDGEN2D is then used to generate surface grids on the outer shell of each component block. GRIDGEN3D generates grid points on the interior of each block, and finally GRIDVUE3D is used to inspect the resulting multiple block grid. Three of these codes (GRIDBLOCK, GRIDGEN2D, and GRIDVUE3D) are highly interactive and graphical in nature, and currently run on Silicon Graphics, Inc., and IBM RS/6000 workstations. The lone batch code (GRIDGEN3D) may be run on any of several Unix based platforms. Surface grid generation in GRIDGEN2D is being improved with the addition of higher order surface definitions (NURBS and parametric surfaces input in IGES format and bicubic surfaces input in PATRAN Neutral File format) and double precision mathematics. In addition, two types of automation have been added to GRIDGEN2D that reduce the learning curve slope for new users and eliminate work for experienced users. Volume grid generation using GRIDGEN3D has been improved via the addition of an advanced hybrid control function formulation that provides both orthogonality and clustering control at the block faces and clustering control on the block interior.

Steinbrenner, John P.↗
NASA NTRS2014

Simulations of Emerging Magnetic Flux. II. The Formation of Unstable Coronal Flux Ropes and the Initiation of Coronal Mass Ejections

We present results from three-dimensional magnetohydrodynamic simulations of the emergence of a twisted convection zone flux tube into a pre-existing coronal dipole field. As in previous simulations, following the partial emergence of the sub-surface flux into the corona, a combination of vortical motions and internal magnetic reconnection forms a coronal flux rope. Then, in the simulations presented here, external reconnection between the emerging field and the pre-existing dipole coronal field allows further expansion of the coronal flux rope into the corona. After sufficient expansion, internal reconnection occurs beneath the coronal flux rope axis, and the flux rope erupts up to the top boundary of the simulation domain (approximately 36 Mm above the surface).We find that the presence of a pre-existing field, orientated in a direction to facilitate reconnection with the emerging field, is vital to the fast rise of the coronal flux rope. The simulations shown in this paper are able to self-consistently create many of the surface and coronal signatures used by coronal mass ejection (CME) models. These signatures include surface shearing and rotational motions, quadrupolar geometry above the surface, central sheared arcades reconnecting with oppositely orientated overlying dipole fields, the formation of coronal flux ropes underlying potential coronal field, and internal reconnection which resembles the classical flare reconnection scenario. This suggests that proposed mechanisms for the initiation of a CME, such as "magnetic breakout," are operating during the emergence of new active regions.

corona↗
NASA NTRS1979

Boundary conditions for implicit solutions to the compressible Navier-Stokes equations in finite computational domains

The admissible boundary conditions for the unsteady Navier-Stokes equations are insufficient to determine completely the flow variables at grid points on the boundaries of a finite computational domain. A general method is presented in which the boundary conditions are supplemented by conservatively differenced forms of the flow equations. For implicit algorithms, the method provides a fully implicit coupling between boundary and interior points. Implicit boundary-point computation schemes are given for most types of flow boundary encountered in internal and external flow problems. Convergence to steady state is generally faster than with locally explicit boundary-point schemes.

Thomas, P. D.↗
NASA NTRS1990

Coupled BE/FE/BE approach for scattering from fluid-filled structures

NASHUA is a coupled finite element/boundary element capability built around NASTRAN for calculating the low frequency far-field acoustic pressure field radiated or scattered by an arbitrary, submerged, three-dimensional, elastic structure subjected to either internal time-harmonic mechanical loads or external time-harmonic incident loadings. Described here are the formulation and use of NASHUA for solving such structural acoustics problems when the structure is fluid-filled. NASTRAN is used to generate the structural finite element model and to perform most of the required matrix operations. Both fluid domains are modeled using the boundary element capability in NASHUA, whose matrix formulation (and the associated NASTRAN DMAP) for evacuated structures can be used with suitable interpretation of the matrix definitions. After computing surface pressures and normal velocities, far-field pressures are evaluated using an asymptotic form of the Helmholtz exterior integral equation. The proposed numerical approach is validated by comparing the acoustic field scattered from a submerged fluid-filled spherical thin shell to that obtained with a series solution, which is also derived here.

Everstine, Gordon C.↗
NASA NTRS1990

Elliptic grid generation with orthogonality and spacing control on an arbitrary number of boundaries

A procedure for the generation of two and quasi-three-dimensional grids with control of orthogonality and spacing with respect to any and/or all boundaries of the domain is described. The elliptic grid generation equations of Thompson are solved implicitly. Control of the grid behavior is achieved through the introduction of forcing functions terms in the manner of Steger and Sorenson or in a modification of the method of Hilgenstock. The forcing function terms are constructed on the boundaries and propagated into the domain using transfinite Lagrangian bivariate interpolation. An anisotropic transfinite stencil is introduced and is shown to produce excellent grid behavior particularly in the vicinity of corner singularities. Emphasis is placed on the generation of viscous grids and the method is shown to be suited for use in the generation of grids for internal as well as external flow geometries. A FORTRAN program named PISCES has been written to implement the algorithm. Examples of grids for internal and external flows are given that highlight the characteristics and behavior of the algorithm.

White, J. A.↗
NASA NTRS2012

Quasi One-Dimensional Unsteady Modeling of External Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of an axisymmetric external compression supersonic inlet is being developed. The model utilizes compressible flow computational fluid dynamics to model the internal inlet segment as well as the external inlet portion between the cowl lip and normal shock, and compressible flow relations with flow propagation delay to model the oblique shocks upstream of the normal shock. The external compression portion between the cowl-lip and the normal shock is also modeled with leaking fluxes crossing the sonic boundary, with a moving CFD domain at the normal shock boundary. This model has been verified in steady state against tunnel inlet test data and it s a first attempt towards developing a more comprehensive model for inlet dynamics.

Kopasakis, George↗
DOE OSTI2024

Visualizing the Internal Nanocrystallinity of Calcite Due to Nonclassical Crystallization by 3D Coherent X–Ray Diffraction Imaging

The internal crystallinity of calcite is investigated for samples synthesized using two approaches: precipitation from solution and the ammonium carbonate diffusion method. Scanning electron microscopy (SEM) analyses reveal that the calcite products precipitated using both approaches have a well–defined rhombohedron shape, consistent with the euhedral crystal habit of the mineral. The internal structure of these calcite crystals is characterized using Bragg coherent diffraction imaging (BCDI) to determine the 3D electron density and the atomic displacement field. BCDI reconstructions for crystals synthesized using the ammonium carbonate diffusion approach have the expected euhedral shape, with internal strain fields and few internal defects. In contrast, the crystals synthesized by precipitation from solution have very complex external shapes and defective internal structures, presenting null electron density regions and pronounced displacement field distributions. These heterogeneities are interpreted as multiple crystalline domains, created by a nonclassical crystallization mechanism, where smaller nanoparticles coalescence into the final euhedral particles. The combined use of SEM, X–ray diffraction (XRD), and BCDI allows for structurally differentiating calcite crystals grown with different approaches, opening new opportunities to understand how grain boundaries and internal defects alter calcite reactivity.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Tangle-Free Finite Element Mesh Motion for Ablation Problems

In numerical simulations involving boundaries that evolve in time, the primary challenge is updating the computational mesh to reflect the physical changes in the domain. In particular, the fundamental objective for any such \mesh motion" scheme is to maintain mesh quality and suppress unphysical geometric anamolies and artifacts. External to a physical process of interest, mesh motion is an added component that determines the specifics of how to move the mesh given certain limited information from the main system. This paper develops a set of boundary conditions designed to eliminate tangling and internal collision within the context of PDE-based mesh motion (linear elasticity). These boundary conditions are developed for two- and three-dimensional meshes. The paper presents detailed algorithms for commonly occuring topological scenarios and explains how to apply them appropriately. Notably, the techniques discussed herein make use of none of the specifics of any particular formulation of mesh motion and thus are more broadly applicable. The two-dimensional algorithms are validated by an extensive verification procedure. Finally, many examples of diverse geometries in both two- and three-dimensions are shown to showcase the capabilities of the tangle-free boundary conditions.

Droba, Justin↗
NASA NTRS1992

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

Fundamental equations of aerodynamic sensitivity analysis and approximate analysis for the two dimensional thin layer Navier-Stokes equations are reviewed, and special boundary condition considerations necessary to apply these equations to isolated lifting airfoils on 'C' and 'O' meshes are discussed in detail. An efficient strategy which is based on the finite element method and an elastic membrane representation of the computational domain is successfully tested, which circumvents the costly 'brute force' method of obtaining grid sensitivity derivatives, and is also useful in mesh regeneration. The issue of turbulence modeling is addressed in a preliminary study. Aerodynamic shape sensitivity derivatives are efficiently calculated, and their accuracy is validated on two viscous test problems, including: (1) internal flow through a double throat nozzle, and (2) external flow over a NACA 4-digit airfoil. An automated aerodynamic design optimization strategy is outlined which includes the use of a design optimization program, an aerodynamic flow analysis code, an aerodynamic sensitivity and approximate analysis code, and a mesh regeneration and grid sensitivity analysis code. Application of the optimization methodology to the two test problems in each case resulted in a new design having a significantly improved performance in the aerodynamic response of interest.

Taylor, Arthur C., III↗
NASA NTRS2003

Effects of Wing Leading Edge Penetration with Venting and Exhaust Flow from Wheel Well at Mach 24 in Flight

A baseline solution for CFD Point 1 (Mach 24) in the STS-107 accident investigation was modified to include effects of: (1) holes through the leading edge into a vented cavity; and (2) a scarfed, conical nozzle directed toward the centerline of the vehicle from the forward, inboard corner of the landing gear door. The simulations were generated relatively quickly and early in the investigation because simplifications were made to the leading edge cavity geometry and an existing utility to merge scarfed nozzle grid domains with structured baseline external domains was implemented. These simplifications in the breach simulations enabled: (1) a very quick grid generation procedure; and (2) high fidelity corroboration of jet physics with internal surface impingements ensuing from a breach through the leading edge, fully coupled to the external shock layer flow at flight conditions. These simulations provided early evidence that the flow through a two-inch diameter (or larger) breach enters the cavity with significant retention of external flow directionality. A normal jet directed into the cavity was not an appropriate model for these conditions at CFD Point 1 (Mach 24). The breach diameters were of the same order or larger than the local, external boundary-layer thickness. High impingement heating and pressures on the downstream lip of the breach were computed. It is likely that hole shape would evolve as a slot cut in the direction of the external streamlines. In the case of the six-inch diameter breach the boundary layer is fully ingested. The intent of externally directed jet simulations in the second scenario was to approximately model aerodynamic effects of a relatively large internal wing pressure, fueled by combusting aluminum, which deforms the corner of the landing gear door and directs a jet across the windside surface. These jet interactions, in and of themselves, were not sufficiently large to explain observed aerodynamic behavior.

Gnoffo, Peter A.↗
NASA NTRS1986

Static current-sheet models of quiescent prominences

A particular class of theoretical models idealize the prominence to be a discrete flat electric-current sheet suspended vertically in a potential magnetic field. The weight of the prominence is supported by the Lorentz force in the current sheet. These models can be extended to have curved electric-current sheets and to vary three-dimensionally. The equation for force balance is 1 over 4 pi (del times B) times Bdel p- p9 z=zero. Using Cartesian coordinates we take, for simplicity, a uniform gravity with constant acceleration g in the direction -z. If we are interested not in the detailed internal structure of the prominence, but in the global magnetic configuration around the prominence, we may take prominence plasma to be cold. Consideration is given to how such equilibrium states can be constructed. To simplify the mathematical problem, suppose there is no electric current in the atmosphere except for the discrete currents in the cold prominence sheet. Let us take the plane z =0 to be the base of the atmosphere and restrict our attention to the domain z greater than 0. The task we have is to solve for a magnetic field which is everywhere potential except on some free surface S, subject to suit able to boundary conditions. The surface S is determined by requiring that it possesses a discrete electric current density such that the Lorentz force on it is everywhere vertically upward to balance the weight of the material m(S). Since the magnetic field is potential in the external atmosphere, the latter is decoupled from the magnetic field and its plane parallel hydrostatic pressure and density can be prescribed.

Wu, F.↗