Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary Conditions”

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 1,099 records · Page 61

Development of Solar Wind Model Driven by Empirical Heat Flux and Pressure Terms

We are developing a time stationary self-consistent 2D MHD model of the solar corona and solar wind as suggested by Sittler et al. (2003). Sittler & Guhathakurta (1999) developed a semiempirical steady state model (SG model) of the solar wind in a multipole 3-streamer structure, with the model constrained by Skylab observations. Guhathakurta et al. (2006) presented a more recent version of their initial work. Sittler et al. (2003) modified the SG model by investigating time dependent MHD, ad hoc heating term with heat conduction and empirical heating solutions. Next step of development of 2D MHD models was performed by Sittler & Ofman (2006). They derived effective temperature and effective heat flux from the data-driven SG model and fit smooth analytical functions to be used in MHD calculations. Improvements of the Sittler & Ofman (2006) results now show a convergence of the 3-streamer topology into a single equatorial streamer at altitudes > 2 R(sub S). This is a new result and shows we are now able to reproduce observations of an equatorially confined streamer belt. In order to allow our solutions to be applied to more general applications, we extend that model by using magnetogram data and PFSS model as a boundary condition. Initial results were presented by Selwa et al. (2008). We choose solar minimum magnetogram data since during solar maximum the boundary conditions are more complex and the coronal magnetic field may not be described correctly by PFSS model. As the first step we studied the simplest 2D MHD case with variable heat conduction, and with empirical heat input combined with empirical momentum addition for the fast solar wind. We use realistic magnetic field data based on NSO/GONG data, and plan to extend the study to 3D. This study represents the first attempt of fully self-consistent realistic model based on real data and including semi-empirical heat flux and semi-empirical effective pressure terms.

Sittler, Edward C., Jr.↗

Flame propagation under partially-premixed conditions

This study concentrates on developing a better understanding of triple flames. We relax the assumption of zero heat release, address the issue of stabilization, and, in order to investigate the role that heat release plays in flame propagation in partially premixed combustion, we return to a simple flow field and investigate the behavior of flames in a laminar environment. We solve the compressible Navier-Stokes equations in a two-dimensional domain. At the boundaries, we use an inflow boundary condition on the left and nearly-perfect reflective boundary conditions, required to avoid pressure drift, at the outflow and sides. After the flow and flame are initialized, the mixture fraction is varied at the inlet from its uniform stoichiometric value to a tanh profile varying from zero to one. As the mixture fraction gradient reaches the flame surface only the centerline is exposed to the stoichionetric mixture fraction and locally maintains the planar flame speed and reaction rate. Above this point the mixture is fuel rich, and below fuel lean. As a result, these regions of non-unity equivalence ratio burn less, the reaction rate drops, and the local flame speed is reduced. The excess fuel and oxidizer then combine behind the premixed flame along the stoichiometric surface and burn in a trailing diffusion flame. Thus the 'triple' flame refers to the fuel-rich premixed flame, the fuel-lean premixed flame, and the trailing diffusion flame. Due to heat release, the normal velocity across the flame is increased, whereas the tangential component remains unchanged. Far-field flame speed, local flame speed, and their differences are shown as a function of the local mixing thickness. It was also determined that the lateral position of the flame affects stabilization, and the distribution of the reaction rate along the premixed wings of triple flames affects flame propagation.

Ruetsch, Gregory R.↗

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗

High-Fidelity CFD Modeling of Cryogenic Hydrogen Isotope Extrusion for Fusion Reactor Pellet Fueling

This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.

ANSYS-Polyflow↗

A finite-element analysis of fatigue crack closure

The paper is concerned with the application of a two-dimensional, nonlinear, finite-element analysis using an incremental theory of plasticity to predict crack-closure and crack-opening stresses during the crack-growth process under cyclic loading. A two-dimensional finite-element computer program is developed, which accounts for both elastic-plastic material behavior and changing boundary conditions associated with crack extension and intermittent contact of the crack surfaces under cyclic loading. An efficient technique to account for changing boundary conditions under cyclic loading is also incorporated into the nonlinear analysis program. This program is used to study crack extension and crack closure behavior in a center-cracked panel under constant-amplitude and two-level block loading. The calculated crack-opening stresses are found to be quantitatively consistent with experimental measurements.

Newman, J. C., Jr.↗

Nonlinear, dispersive, elliptically polarized Alfven wavaes

The derivative nonlinear Schroedinger (DNLS) equation is derived by an efficient means that employs Lagrangian variables. An expression for the stationary wave solutions of the DNLS that contains vanishing and nonvanishing and modulated and nonmodulated boundary conditions as subcases is then obtained. The solitary wave solutions for elliptically polarized quasiparallel Alfven waves in the magnetohydrodynamic limit (nonvanishing, unmodulated boundary conditions) are obtained. These converge to the Korteweg-de Vries and the modified Korteweg-de Vries solitons obtained previously for oblique propagation, but are more general. It is shown that there are no envelope solitary waves if the point at infinity is unstable to the modulational instability. The periodic solutions of the DNLS are characterized.

Kennel, C. F.↗

Interaction and merging of vortex filaments

The asymptotic solutions of Navier-Stokes equations for vortex filaments of finite strength with small effective vortical cores are summarized with special emphasis placed on the physical meaning and the practical limit to the applicability of the asymptotic solution. Finite-difference solutions of Navier-Stokes equations for the marging of the filament(s) are described with a focus on the development of the approximate boundary conditions for the computational domain. An efficiency study employing a model problem is used to assess the advantages of the present approximate boundary condition method over previously used techniques. Applications of the present method are presented for the motion and decay of a 3:1 elliptic vortex ring, and for the merging process of a pair of coaxial vortex rings. A numerical procedure for the problem of local merging of vortex filaments, which requires the asymptotic analysis as well as the numerical Navier-Stokes solver, is also presented.

Liu, C. H.↗

Evaluation of non-Fourier heat waves influenced by nonlinear/linear boundary effects employing an explicit architecture and controlled stabilization

The present paper is concerned with the problem of heat waves in solids, where, the heat transport due to conduction occurs as propagating thermal disturbances which are transmitted at finite but high speeds. Starting from the general heat flux model of the Jeffrey's-type, and subsequent formulations leading to the Cattaneo-type heat flux model, an evaluation of the heat transport behavior is described for models influenced by non-Fourier effects and subjected to general nonlinear/linear boundary conditions. An explicit time-integration architecture is employed which effectively provides not only accurate representations of the relaxation effects and general boundary conditions but also seeks to provide an understanding of the representative thermal behavior and heat transport mechanisms for a variety of physical situations.

Tamma, Kumar K.↗

Development of a CFD code for internal flows in liquid fueled engines

To support the design efforts of engines in which liquid propellants are used, one is often required to analyze incompressible, two-dimensional or axisymmetric flows within ducts and cavities with rotating walls and complicated geometries. The steady-state solution is of interest in most cases. This code is intended to provide a tool for efficient computational fluid dynamics (CFD) analysis of such flow problems, taking advantage of the artificial compressibility concept and the Beam-Warming numerical scheme modified for second order accurate, implicit boundary conditions. These concepts ensure a stable, robust, and accurate algorithm due to the reduced speed of sound and the accuracy of the boundary conditions. The code is dedicated only to two-dimensional or axisymmetric flows with or without swirl. Three-dimensional computation is excluded to increase efficiency and speed. This paper briefly presents the theory of the code, as well as several benchmark applications with comparison to well known analytical solutions. In all these test cases, the code produced accurate results.

Dakhoul, Youssef↗

Numerical simulations of small amplitude acoustic wave propagation in a converging-diverging nozzle

Numerical simulations of the small amplitude acoustic wave propagation in a converging-diverging nozzle proposed in Category 5 are performed and presented here. The quasi 1-D unsteady flow equations in conservative form are discretized by the DRP scheme with the artificial damping developed by Tam and Webb. Characteristic boundary conditions are used at both the inlet and outlet of the nozzle. The effect of different numerical implementations of the subsonic inflow boundary conditions on the convergence of the solution to the steady state is studied. In the case of a subsonic outflow in which a shock is formed in the nozzle, the interaction between acoustics and the shock is also investigated.

Dong, Thomas Z.↗

Numerical Examination of Silicon Avalanche Photodiodes Operated in Charge Storage Mode

The behavior of silicon-based avalanche photodiodes (APD's) operated in the charge storage mode is examined. In the charge storage mode, the diodes are periodically biased to a sub-breakdown voltage and then open-circuited. During this integration period, photo-excited and thermally generated carriers are accumulated within the structure. The dynamics of this accumulation and its effects upon the avalanching of the diode warrants a detailed, fully numerical analysis. The salient features of this investigation include device sensitivity to the input photo-current including the self-quenching effect of the diode and its limitations in sensing low light levels, the dependence of the response on the bulk lifetime and hence on the generation current within the device, the initial gain, transient response, dependence of the device uniformity upon performance, and the quantity of storable charge within the device. To achieve these tasks our device simulator, STEBS-2D, was utilized. A modified current-controlled boundary condition is employed which allows for the simulation of the isolated diode after the initial reset bias has been applied. With this boundary condition, it is possible to establish a steady-state voltage on the ohmic contact and then effectively remove the device from the external circuit while still including effects from surface recombination, trapped surface charge, and leakage current from the read-out electronics.

Parks, Joseph W., Jr.↗

Global aspects of the motion of the heliospheric termination shock: A gasdynamic model

The heliospheric termination shock is expected to move in response to variation in upstream solar wind conditions. Using numerical techniques, we extend an earlier strictly one-dimensional analytic gas dynamic model of shock motion to two dimensions, to investigate the qualitative features of global behavior of shock motion, and the consequences of latitudinal variation in dynamic pressure. The boundary conditions of the calculation are given by the solar wind parameters as a function of latitude and time on an inner spherical boundary, and a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. Density variations, specified at the inner boundary as a function of time, are convected into the termination shock. Immediately after the interaction, the shock moves with speeds given by the earlier analytic model. However, as the termination shock propagates outward (or inward), it begins to slow down. After about 2 to 10 years, depending on details of boundary conditions, the signal from the shock interaction has reached the outer boundary and propagates inward to the position of the termination shock, strongly affecting the behavior of the shock. Assuming no further disturbances in the solar wind, the termination shock will reach its new equilibrium after some tens of years. In reality, large-scale variations in solar wind dynamic pressure occur on time scales short in comparison with the eleven year solar cycle, so that one expects that the termination shock is never in an equilibrium position, but rather oscillates inward and outward; this oscillation will vary with heliographic latitude. The effects of a variety of types of solar wind disturbances are investigated and summarized.

Ratkiewicz, R.↗

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗

Initial Results from the GEM Challenge on the Spacecraft Surface Charging Environment

Spacecraft surface charging during geomagnetically disturbed times is one of the most important causes of satellite anomalies. Predicting the surface charging environment is one prevalent task of the geospace environment models. Therefore, the Geospace Environment Modeling (GEM) Focus Group "Inner Magnetosphere Crossenergy/Population Interactions" initiated a communitywide challenge study to assess the capability of several inner magnetosphere ring current models in determining surface charging environment for the Van Allen Probes orbits during the 17 March 2013 storm event. The integrated electron flux between 10 and 50 keV is used as the metrics. Various skill scores are applied to quantitatively measure the modeling performance against observations. Results indicate that no model consistently perform the best in all of the skill scores or for both satellites. We find that from these simulations the ring current model with observational flux boundary condition and Weimer electric potential driver generally reproduces the most realistic flux level around the spacecraft. A simple and weaker VollandStern electric field is not capable of effectively transporting the same plasma at the boundary toward the Earth. On the other hand, if the ring current model solves the electric field selfconsistently and obtains similar strength and pattern in the equatorial plane as the Weimer model, the boundary condition plays another crucial role in determining the electron flux level in the inner region. When the boundary flux spectra based on magnetohydrodynamics (MHD) model/empirical model deviate from the shape or magnitude of the observed distribution function, the simulation produces poor skill scores along Van Allen Probes orbits.

Yu, Yiqun↗

Antennas in plasmas - The problem of boundaries

This paper discusses some of the difficulties associated with definition of appropriate boundary conditions for use in analyses of antennas and other structures immersed in plasmas. Distinctions may be made according to whether cold plasma, microscopic, or macroscopic theory is to be employed; whether the plasma is to be considered uniform or inhomogeneous, and whether or not there is a static magnetic field. The properties and limitations of these models, and the commonly-used boundary conditions of specular reflection and velocity diffusion (for opaque boundaries), and transparency (for grids or wires) are discussed. As a specific illustrative example, a series of approximations to the impedance of a magnetoplasma capacitor are reviewed.

Crawford, F. W.↗

Geometries with mismatched branes

We study Randall-Sundrum two brane setups with mismatched brane tensions. For the vacuum solutions, boundary conditions demand that the induced metric on each of the branes is either de Sitter, Anti-de Sitter, or Minkowski. For incompatible boundary conditions, the bulk metric is necessarily time-dependent. This introduces a new class of time-dependent solutions with the potential to address cosmological issues and provide alternatives to conventional inflationary (or contracting) scenarios. We take a first step in this paper toward such solutions. One important finding is that the resulting solutions can be very succinctly described in terms of an effective action involving only the induced metric on either one of the branes and the radion field. But the full geometry cannot necessarily be simply described with a single coordinate patch. We concentrate here on the time- dependent solutions but argue that supplemented with a brane stabilization mechanism one can potentially construct interesting cosmological models this way. This is true both with and without a brane stabilization mechanism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. Here we discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states, and the multiscale entanglement renormalization ansatz are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled pair states are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

97 MATHEMATICS AND COMPUTING↗

Nonlinear finite element models of laminated plates and shells

In this expository paper a first-order shear deformation laminate theory that accounts for the von Karman strains is described and its finite element model is developed. The element also accounts for Hill's elastic-perfectly plastic material model. Sample applications of the present element to linear and nonlinear bending of laminated plates and shells are presented. Numerical results showing the effects of boundary conditions, lamination scheme, and loading on the deflections and stresses are also presented. It is found that the effect of shear deformation and boundary conditions on deflections are significant.

Redddy, J. N.↗