SEARCH · Engineering Papers

Results for “embedded boundary”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Analytic and approximate boundary-fitted coordinate systems for fluid flow simulation

A method for constructing a direct functional relation between the computational and physical domains such that boundaries in the computational domain map into boundaries in the physical domain is presented. This method which is called the 'two-boundary technique' is based on generating functions that represent two disconnected boundaries in the physical domain and an independent function which connects one boundary with another to generate grids. The suggested connecting functions are simple linear and cubic parametric polynomials with parameters based on position and derivatives at the disconnected boundaries. The boundaries can be defined by analytical, approximate, and tension spline functions; the grid density in specified regions is controlled by embedding contraction functions in either the boundary representation or the connecting function.

Smith, R. E.↗

NASA NTRS2013

Recent Development in the CESE Method for the Solution of the Navier-Stokes Equations Using Unstructured Triangular or Tetrahedral Meshes With High Aspect Ratio

In the multidimensional CESE development, triangles and tetrahedra turn out to be the most natural building blocks for 2D and 3D spatial meshes. As such the CESE method is compatible with the simplest unstructured meshes and thus can be easily applied to solve problems with complex geometries. However, because the method uses space-time staggered stencils, solution decoupling may become a real nuisance in applications involving unstructured meshes. In this paper we will describe a simple and general remedy which, according to numerical experiments, has removed any possibility of solution decoupling. Moreover, in a real-world viscous flow simulation near a solid wall, one often encounters a case where a boundary with high curvature or sharp corner is surrounded by triangular/tetrahedral meshes of extremely high aspect ratio (up to 106). For such an extreme case, the spatial projection of a space-time compounded conservation element constructed using the original CESE design may become highly concave and thus its centroid (referred to as a spatial solution point) may lie far outside of the spatial projection. It could even be embedded beyond a solid wall boundary and causes serious numerical difficulties. In this paper we will also present a new procedure for constructing conservation elements and solution elements which effectively overcomes the difficulties associated with the original design. Another difficulty issue which was addressed more recently is the wellknown fact that accuracy of gradient computations involving triangular/tetrahedral grids deteriorates rapidly as the aspect ratio of grid cells increases. The root cause of this difficulty was clearly identified and several remedies to overcome it were found through a rigorous mathematical analysis. However, because of the length of the current paper and the complexity of mathematics involved, this new work will be presented in another paper.

Chang, Sin-Chung↗

NASA NTRS1985

Control plate for shock-boundary layer interaction

Paper describes tests and computations for a relatively unique technique to greatly reduce/eliminate the separation region for shock-boundary layer interactions. A number of studies have shown that the usual effects of such interactions include increased local heating and wall pressures, thickening of the boundary layer and a decrease in the momentum of the flow and, for stronger waves, flow separation. This flow situation is particularly prevalent in supersonic and hypersonic inlets where severe performance degradation can occur due to flow separation. High performance engine design generally requires a uniform entering flow field with little stagnation pressure loss. Previous approaches to the problem involved primarily active devices (e.g., suction or blowing); the present paper considers a passive device. The boundary layer separation control technique considered herein involves the placement of an embedded plate in the outer portion of the boundary layer and parallel to the wall. This control plate is situated such that the incident shock impinges upon and reflects from its surface, thus greatly lessening the pressure gradient in the low momentum near wall region.

Goodman, W. L.↗

DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗

DOE OSTI2026

Substrate-Directed Underlayer Growth of Bilayer MoS 2 Revealed by Mo Isotope Labeling

Direct control over the vertical formation sequence and stacking registry in van der Waals (vdW) bilayers is essential for device performance and moiré engineering yet difficult to resolve unambiguously with conventional probes. Here, we use Mo isotope labeling in a two-step chemical vapor deposition process to synthesize bilayer MoS 2 and trace its vertical formation on common substrates. By combining site-selective laser thinning, Raman spectroscopy, time-of-flight secondary ion mass spectrometry, and atomic-resolution scanning transimission electron microscopy (STEM), we find a clear substrate dependence: on SiO 2 /Si, the second layer nucleates and grows beneath the first (underlayer), whereas on sapphire, it forms on top (overlayer). Density functional theory indicates that a larger equilibrium interfacial separation and weaker MoS 2 –substrate interactions on amorphous SiO 2 permit confined interfacial diffusion and underlayer nucleation, whereas stronger interactions and smaller separations on sapphire favor overlayer growth. On SiO 2 , confined epitaxy templates commensurate 2H, 3R, and mixed bilayers, as confirmed by second harmonic generation spectroscopy and STEM. During underlayer coalescence, embedded mirror-twin grain boundaries stitch atomically sharp 2H|3R junctions via alternating 4|8 ring motifs. Molecular-dynamics simulations reveal that these alternating 4|8 motifs accommodate interlayer vdW coupling and locally modulate the stacking registry. These results provide mechanistic insight into confined epitaxial growth and establish isotope labeling as a powerful probe of two-dimensional materials synthesis.

MoS2↗

DOE OSTI2020

Interface learning of multiphysics and multiscale systems

Complex natural or engineered systems comprise multiple characteristic scales, multiple spatiotemporal domains, and even multiple physical closure laws. To address such challenges, we introduce an interface learning paradigm and put forth a data-driven closure approach based on memory embedding to provide physically correct boundary conditions at the interface. To enable the interface learning for hyperbolic systems by considering the domain of influence and wave structures into account, we put forth the concept of upwind learning toward a physics-informed domain decomposition. The promise of the proposed approach is shown for a set of canonical illustrative problems. Here, we highlight that high-performance computing environments can benefit from this methodology to reduce communication costs among processing units in emerging machine-learning-ready heterogeneous platforms toward exascale era.

42 ENGINEERING↗

NASA NTRS1982

Numerical simulation of viscous-inviscid interactions on indented nose tips

An implicit numerical algorithm to solve the unsteady thin-layer Navier-Stokes equations in a strong conservative form has been used to compute the viscous flow over indented nose tips placed in a supersonic free stream. Numerical solutions are presented for axisymmetric and three-dimensional indented configurations for laminar flow conditions. Results demonstrate the capability of the present numerical procedure to predict flow fields that contain strong viscous-inviscid interactions, including boundary-layer separation, reattachment, and embedded discontinuities. Comparisons with available experimental data for the shock shape and surface pressure are also given.

Rizk, Y. M.↗

NASA NTRS1982

A transonic rectangular grid embedded panel method

A method is presented that has the potential for solving transonic flow problems about the same complex aircraft configurations currently being analyzed by subsonic panel methods. This method does not require the generation of surface fitted grids. Instead it uses rectangular grids and subgrids together with embedded surface panels on which boundary conditions are imposed. Both the Euler and full potential equations are considered. The method of least squares is used to reduce the solution of these equations to the solution of a sequence of Poisson problems. The Poisson problems are solved using fast Fourier transforms and panel influence coefficient techniques. The overall method is still in its infancy but some two dimensional results are shown illustrating various key features.

Johnson, F. T.↗

NASA NTRS1986

Application of a computational model for vortex generators in subsonic internal flows

A model for the analysis of vortex generators in a fully viscous subsonic internal flow is evaluated. A vorticity source term is used in a modified form of the Parabolized Navier-Stokes equations to model the shed vortex. Computed results are compared with idealized flow vortex paths, and with experimental data for vortex generators embedded in a thick turbulent boundary layer. The analysis is also compared with experimental data for a separated diffusing S-duct and for a diffusing S-duct with vortex generators. Quantitative comparisons are shown for the latter three cases. Emphasis is placed on verifying the ability of the model to predict global distortions in the flow field.

Kunik, W. G.↗

NASA NTRS1986

Application of a computational model for vortex generators in subsonic internal flows

Kunik, W. G.↗

NASA NTRS1990

Numerical simulation of the incompressible internal flow through a tilting disk valve

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady and unsteady flow calculations are performed by solving the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Chang, I-Dee↗

NASA NTRS1991

Simulation of blood flow through an artificial heart

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady state and unsteady flow calculations are done by solving the incompressible Navier-Stokes equations in 3-D generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Kiris, Cetin↗

NASA NTRS1996

Spanwise Spacing Effects on the Initial Structure and Decay of Axial Vortices

The initial structure and axial decay of an array of streamwise vortices embedded in a turbulent pipe boundary layer is experimentally investigated. The vortices are shed in counter-rotating fashion from an array of equally-spaced symmetric airfoil vortex generators. Vortex structure is quantified in terms of crossplane circulation and peak streamwise vorticity. Flow conditions are subsonic and incompressible. The focus of this study is on the effect of the initial spacing between the parent vortex generators. Arrays with vortex generators spaced at 15 and 30 degrees apart are considered. When the spacing between vortex generators is decreased the circulation and peak vorticity of the shed vortices increases. Analysis indicates this strengthening results from regions of fluid acceleration in the vicinity of the vortex generator array. Decreased spacing between the constituent vortices also produces increased rates of circulation and peak vorticity decay.

Wendt, B. J.↗

NASA NTRS2017

Development of a Flow Field for Testing a Boundary-Layer-Ingesting Propulsor

The test section of the 8- by 6-Foot Supersonic Wind Tunnel at NASA Glenn Research Center was modified to produce the test conditions for a boundary-layer-ingesting propulsor. A test was conducted to measure the flow properties in the modified test section before the propulsor was installed. Measured boundary layer and freestream conditions were compared to results from computational fluid dynamics simulations of the external surface for the reference vehicle. Testing showed that the desired freestream conditions and boundary layer thickness could be achieved; however, some non-uniformity of the freestream conditions, particularly the total temperature, were observed.

embedded propulsion↗

NASA NTRS1976

Three-dimensional boundary layers approaching separation

The theory of semi-similar solutions of the laminar boundary layer equations is applied to several flows in which the boundary layer approaches a three-dimensional separation line. The solutions obtained are used to deduce the nature of three-dimensional separation. It is shown that in these cases separation is of the "ordinary" type. A solution is also presented for a case in which a vortex is embedded within the three-dimensional boundary layer.

Williams, J. C., III↗

DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗

DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1979

Off-axis impact of unidirectional composites with cracks: Dynamic stress intensification

The dynamic response of unidirectional composites under off axis (angle loading) impact is analyzed by assuming that the composite contains an initial flaw in the matrix material. The analytical method utilizes Fourier transform for the space variable and Laplace transform for the time variable. The off axis impact is separated into two parts, one being symmetric and the other skew-symmetric with reference to the crack plane. Transient boundary conditions of normal and shear tractions are applied to a crack embedded in the matrix of the unidirectional composite. The two boundary conditions are solved independently and the results superimposed. Mathematically, these conditions reduce the problem to a system of dual integral equations which are solved in the Laplace transform plane for the transformation of the dynamic stress intensity factor. The time inversion is carried out numerically for various combinations of the material properties of the composite and the results are displayed graphically.

Sih, G. C.↗