Engineering Papers⌕ Search

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.

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

At least 181 records · Page 10

Seismic Performance of Deeply Embedded Conceptual Advanced Reactors: Three-Dimensional Nonlinear Soil-Structure Interaction Analyses

The U.S. Nuclear Regulatory Commission’s functional containment concept provides advanced nuclear power plant designers with more flexibility in terms of the civil/structural design if the appropriate set of barriers for prevention of radioactive material release exist. Some of the conceptual advanced reactor structures, without the traditional pressure boundaries of large containment structures, are proposed to be deeply embedded or buried into soil. This approach is expected to provide (1) lesser seismic demands on the structures and safety-critical structures, (2) eased regulatory efforts and overall design against other external hazards such as aircraft impact, and (3) overall cost savings. Further, one of the important aspects of assessing the technical and economic viability of deeply embedding advanced reactor buildings is to assess the seismic performance with the understanding of effects with material and geometric nonlinearities. This study investigates the seismic response of deeply embedded or buried advanced reactors by conducting three-dimensional nonlinear soil-structure interaction analyses. Although the results indicate that there is a general trend of decreased seismic response with increased embedment depths, the change in the dynamic environment with different embedment depths and the nonlinear environment under high-intensity seismic inputs may result in increased peak response at increased embedment depths.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. 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 extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

User's manual for the Langley boundary layer noise propagation program (MRS-BLP)

A computer program, McAninch-Rawls-Spence Boundary Layer Propagation (MRS-BLP), is described. This program models the refractive and scattering effects on acoustic pressure waves propagating through a boundary layer encompassing an aircraft's fuselage. The noise source is assumed known and generated by a propeller. The fuselage is represented by an infinitely long cylinder embedded in a longitudinal flow. By matching a numerical solution inside the boundary layer with an analytical solution outside the boundary layer, the program calculates the acoustic pressure at the surface of the cylinder given the incident field at the top of the boundary layer. The boundary layer flow velocity and sound speed profiles, as well as the boundary layer thickness may be specified by the user. A detailed description of the input parameters and how to execute the program is given. Example executions of MRS-BLP showing results are also included.

Spence, Peter L.↗

Further studies of unsteady boundary layers with flow reversal

One set of calculations was performed using the first order, time dependent turbulent boundary layer equations, and extended earlier work by Nash and Patel to a wider range of flows. Another set of calculations was performed for laminar flow using the time dependent Navier-Stokes equations. The results of the calculations confirm previous conclusions concerning the existence of a regime of unseparated flow, containing an embedded region of reversal, which is accessible to first order boundary layer theory. However, certain doubts are cast on the precise nature of the events which accompany the eventual breakdown of the theory due to singularity onset. The earlier view that the singularity appears as the final event in a sequence involving rapid thickening of the boundary layer and the formation of a localized region of steep gradients is called into question by the present results. It appears that singularity onset is not necessarily preceded by rapid boundary layer thickening, or even necessarily produces immediate thickening.

Nash, J. F.↗

Boundary-fitted coordinate systems for numerical solution of partial differential equations - A review

A comprehensive review of methods of numerically generating curvilinear coordinate systems with coordinate lines coincident with all boundary segments is given. Some general mathematical framework and error analysis common to such coordinate systems is also included. The general categories of generating systems are those based on conformal mapping, orthogonal systems, nearly orthogonal systems, systems produced as the solution of elliptic and hyperbolic partial differential equations, and systems generated algebraically by interpolation among the boundaries. Also covered are the control of coordinate line spacing by functions embedded in the partial differential operators of the generating system and by subsequent stretching transformation. Dynamically adaptive coordinate systems, coupled with the physical solution, and time-dependent systems that follow moving boundaries are treated. References reporting experience using such coordinate systems are reviewed as well as those covering the system development.

Thompson, J. F.↗

Suppression of Collisionless Magnetic Reconnection in the High Ion β, Strong Guide Field Limit

Abstract In magnetic reconnection, the ion bulk outflow speed and ion heating have been shown to be set by the available reconnecting magnetic energy, i.e., the energy stored in the reconnecting magnetic field ( B r ). However, recent simulations, observations, and theoretical works have shown that the released magnetic energy is inhibited by upstream ion plasma beta β i —the relative ion thermal pressure normalized to magnetic pressure based on the reconnecting field—for antiparallel magnetic field configurations. Using kinetic theory and hybrid particle-in-cell simulations, we investigate the effects of β i on guide field reconnection. While previous works have suggested that guide field reconnection is uninfluenced by β i , we demonstrate that the reconnection process is modified and the outflow is reduced for sufficiently large β i > ( B r 2 + B g 2 ) / B r 2 . We develop a theoretical framework that shows that this reduction is consistent with an enhanced exhaust pressure gradient, which reduces the outflow speed as v out ∝ 1 / β i . These results apply to systems in which guide field reconnection is embedded in hot plasmas, such as reconnection at the boundary of eddies in fully developed turbulence like the solar wind or the magnetosheath as well as downstream of shocks such as the heliosheath or the mergers of galaxy clusters.

Giai, Carlos A.↗

Radiation and scattering from cylindrically conformal printed antennas

Microstrip patch antennas offer considerable advantages in terms of weight, aerodynamic drag, cost, flexibility, and observables over more conventional protruding antennas. These flat patch antennas were first proposed over thirty years ago by Deschamps in the United States and Gutton and Baisinot in France. Such antennas have been analyzed and developed for planar as well as curved platforms. However, the methods used in these designs employ gross approximations, suffer from extreme computational burden, or require expensive physical experiments. The goal of this thesis is to develop accurate and efficient numerical modeling techniques which represent actual antenna structures mounted on curved surfaces with a high degree of fidelity. In this thesis, the finite element method is extended to cavity-backed conformal antenna arrays embedded in a circular, metallic, infinite cylinder. Both the boundary integral and absorbing boundary mesh closure conditions will be used for terminating the mesh. These two approaches will be contrasted and used to study the scattering and radiation behavior of several useful antenna configurations. An important feature of this study will be to examine the effect of curvature and cavity size on the scattering and radiation properties of wraparound conformal antenna arrays.

Kempel, Leo C.↗

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗

The nature of AFGL 2591 and its associated molecular outflow: Infrared and millimeter-wave observations

The results of infrared photometry from 2 to 160 microns of AFGL and CO(12) observations of its associated molecular cloud and high velocity molecular outflow are presented and discussed. The observed solar luminosity is 6.7 x 10(4) at a distance of 2 kpc. The spectrum of AFGL 2591 is interpreted in the context of a model in which a single embedded object is the dominant source of the infrared luminosity. This object is determined to be surrounded by a compact, optically thick dust shell with a temperature in excess of several hundred degrees kelvin. The extinction to this source is estimated to be between 26 and 50 visual magnitudes. The absolute position of the infrared sources at 10 microns was determined to an accuracy of + or in. This indicates for the first time that the IR source and H2O source are not coincident. The CO(12) observations show the high-velocity molecular flow near AFGL 2591 to be extended, bipolar and roughly centered on the infrared emission. The observations suggest that the red-shifted flow component extends beyond the boundary of the ambient cloud within which AFGL 2591 is embedded. The CO(12) observations also show that AFGL 2591 is embedded in a molecular cloud with an LSR velocity of -5 km/s.

Lada, C. J.↗

The nature of AFGL 2591 and its associated molecular outflow Infrared and millimeter-wave observations

The results of infrared photometry from 2 to 160 microns of AFGL and CO(12) observations of its associated molecular cloud and high velocity molecular outflow are presented and discussed. The observed solar luminosity is 6.7 x 10(4) at a distance of 2 kpc. The spectrum of AFGL 2591 is interpreted in the context of a model in which a single embedded object is the dominant source of the infrared luminosity. This object is determined to be surrounded by a compact, optically thick dust shell with a temperature in excess of several hundred degrees kelvin. The extinction to this source is estimated to be between 26 and 50 visual magnitudes. The absolute position of the infrared sources at 10 microns was determined to an accuracy of + or - 1 in. This indicates for the first time that the IR source and H2O source are not coincident. The CO(12) observations show the high-velocity molecular flow near AFGL 2591 to be extended, bipolar and roughly centered on the infrared emission. The observations suggest that the red-shifted flow component extends beyond the boundary of the ambient cloud within which AFGL 2591 is embedded. The CO(12) observations also show that AFGL 2591 is embedded in a molecular cloud with an LSR velocity of -5 km/s.

Lada, C. J.↗

Nonlinear convergence in contact mechanics: Immersed boundary finite volume

In this report we present an immersed boundary finite volume (IBM) method for simulating quasistatic contact mechanics of linearly elastic domains at small strains. In IBM, all external boundaries and internal contacts of an object are represented by embedded surfaces inside a Cartesian mesh, which need not conform to the grid lines. The contact constraints consist of the non-penetrability condition and Coulomb’s friction law, which are discretized using special interpolation stencils and enforced via penalty parameters. The resulting nonlinear system depends on displacement unknowns only. To solve it, we use the Newton method but find that it diverges frequently. To understand the divergence pattern, we analyze a simplified 2-cell problem and show that the global convergence of Newton cannot be ensured for any choice of penalty parameters. We thus propose a modified Newton solver, which guarantees convergence for the 2-cell problem and is numerically verified to converge for all the challenging simulations considered herein. While both 1 st - and 2 nd -order variants of IBM, in displacement unknowns, are proposed, the modified Newton solver applies only to the 1 st -order variant.

42 ENGINEERING↗

Generation After Next Propulsor Research: Robust Design for Embedded Engine Systems

The National Aeronautics and Space Administration, United Technologies Research Center and Virginia Polytechnic and State University have contracted to pursue multi-disciplinary research into boundary layer ingesting (BLI) propulsors for generation after next environmentally responsible subsonic fixed wing aircraft. This Robust Design for Embedded Engine Systems project first conducted a high-level vehicle system study based on a large commercial transport class hybrid wing body aircraft, which determined that a 3 to 5 percent reduction in fuel burn could be achieved over a 7,500 nanometer mission. Both pylon-mounted baseline and BLI propulsion systems were based on a low-pressure-ratio fan (1.35) in an ultra-high-bypass ratio engine (16), consistent with the next generation of advanced commercial turbofans. An optimized, coupled BLI inlet and fan system was subsequently designed to achieve performance targets identified in the system study. The resulting system possesses an inlet with total pressure losses less than 0.5%, and a fan stage with an efficiency debit of less than 1.5 percent relative to the pylon-mounted, clean-inflow baseline. The subject research project has identified tools and methodologies necessary for the design of next-generation, highly-airframe-integrated propulsion systems. These tools will be validated in future large-scale testing of the BLI inlet / fan system in NASA's 8 foot x 6 foot transonic wind tunnel. In addition, fan unsteady response to screen-generated total pressure distortion is being characterized experimentally in a JT15D engine test rig. These data will document engine sensitivities to distortion magnitude and spatial distribution, providing early insight into key physical processes that will control BLI propulsor design.

Arend, David J.↗

Analysis of the heliospheric current sheet at Earth's orbit and model comparisons

IMP 8 magnetic field data for the first half of the year 1994, i.e., for about 6 solar rotations, are analyzed around regions of sector boundary crossings with the purpose of obtaining both gross- and fine-scale characteristics of the related heliospheric current sheets separating the observed sectors. For purposes of estimating the attitudes of the normals to the sector boundaries. analysis intervals (sometimes 30 min or more in length) allowing the field to fully complete an excursion of about 180 deg were used in the study, which consisted of variance analyses of the field within those intervals. The resulting boundary normals were analyzed and compared to known (generic) models of projected heliospheric current sheets and to a coronal field model for the same time period. One of the most outstanding features of the resulting ensemble of estimated boundary normals for this period is that they strongly prefer low inclinations, indicating that the observations do not support a 1 AU model that predicts a current sheet whose surface is approximately parallel with the sun's equator, such as the 'sombrero' model. They instead support a model that predicts a relatively high inclination current sheet at 1 AU. Also the normals assume a surprisingly large range of longitudes, somewhat favoring those consistent with a Parker model (45 deg and 225 deg) and/or radial alignment (0 deg and 180 deg). These boundary structures, as defined, are shown typically to be as broad as several hundred proton gyroradii, but having embedded within them very thin structures associated with stronger currents. Such thin structures have normals usually differing markedly from the gross boundary. For some crossings there are indications of a wave-like structure in the current sheet as it passed the spacecraft.

Lepping, R. P.↗

Helical spin structure in iron chains with hybridized boundaries

In this paper, we have compared the magnetic properties of well-controlled ultra-short (≤50 nm) atomic iron (Fe) chains embedded in Fe-phthalocyanine films with those in Fe–hydrogen (H2) phthalocyanine superlattices. Surprisingly, we found that the coercivity of the atomic chains with free boundary conditions is independent of the chain length, whereas the one subject to hybridization of the chain ends exhibits an unexpected length dependence. These findings suggest that ferromagnetism in the free-boundary condition system is caused by an intrinsic indirect exchange. On the other hand, controlled boundary conditions produce a helical spin structure due to an extrinsic indirect exchange, which arises from the interaction between iron atoms at the ends of the chain and the hydrogen in the H 2 phthalocyanine spacer. As a consequence, during magnetic reversal, ultra-short iron chains subject to boundary clamping develop a helical spin structure, leading to increased coercivity. These findings suggest unique insights and ideas for the design of atomic-scale ultra-dense magnetic storage nanodevices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A deep learning upscaling framework: Reactive transport and mineral precipitation in fracture-matrix systems

Pore-scale modeling has limited applicability at large scales due to its high computational cost. One common approach to upscale pore-scale models is the use of effective medium theories, which homogenize small-scale features in a porous structure and characterize the medium by macroscale properties (e.g., permeability) and equations (e.g., Darcy’s law). However, there are classes of physical processes for which effective medium approximations may become inaccurate, e.g., mineral precipitation and clogging during reactive transport. We have developed a deep learning upscaling framework, in which pore-scale modeling is directly employed in macroscale systems, without relying on effective medium approximations. The upscaling framework is first developed for general multiscale systems and then applied to modeling reactive transport with mineral precipitation in the altered layer in fracture-matrix structures. Solute transport from the fractures to the matrix is modeled as a wall boundary condition for the fractures, which, in turn, is predicted by recurrent neural networks using the concentration histories at the fracture-matrix boundary. Specifically, we consider a meter-scale fracture network embedded in sandstones, where the smallest feature is at the micron scale. Here the proposed framework allows us to span five orders of magnitude in length scales by capturing mineral precipitation in the altered layer of the rock matrix at the pore scale across the entire meter-scale fracture network.

42 ENGINEERING↗

Supersonic flow over ablated nosetips using an unsteady, implicit numerical procedure

The axisymmetric supersonic flow over passive, that is, nonablating, indented nosetips of reentry vehicles is determined using an unsteady implicit numerical algorithm which solves either the inviscid Euler equations or the 'thin-layer' Navier-Stokes equations. A nonorthogonal independent variable transformation is used to map the distorted physical region, containing multiple zones of embedded subsonic flow into a rectangular computational domain at whose boundaries the required permeable or impermeable boundary conditions are simulated. Use of the implicit algorithm results in faster convergence to the steady state because of a larger allowable time step over conventional explicit schemes. The numerical results obtained compare favorably with existing experimental data for very mildly and severely indented blunt nosetips.

Kutler, P.↗

Evaluation of Variable-Depth Liners with Slotted Cores

An initial evaluation of variable-depth, slotted-core liners fabricated with conventional, flightworthy com-ponents is presented. Liners with two distinct chamber lengths provide significant broadband absorption and their impedances are, to first order, predictable with the NASA impedance prediction code. These liners employ film caps embedded in alternating chambers to provide a near-rigid boundary condition. This enables improved control of the liner design process. A number of additively-manufactured, bent-chamber liners are also investigated to explore the e˙ects of the bend on the liner impedance. The liner impedance is negligibly a˙ected by changes to the partition width, but is significantly a˙ected by changes to the slot height. This o˙ers an option for achieving low-frequency sound absorption with relatively thin liners.

acoustic↗

Evaluation of Variable-Depth Liners with Slotted Cores

An initial evaluation of variable-depth, slotted-core liners fabricated with conventional, flightworthy components is presented. Liners with two distinct chamber lengths provide significant broadband absorption and their impedances are, to first order, predictable with the NASA impedance prediction code. These liners employ film caps embedded in alternating chambers to provide a near-rigid boundary condition. This enables improved control of the liner design process. A number of additively-manufactured, bent-chamber liners are also investigated to explore the effects of the bend on the liner impedance. The liner impedance is negligibly affected by changes to the partition width, but is significantly affected by changes to the slot height. This offers an option for achieving low-frequency sound absorption with relatively thin liners.

acoustic↗