Moving and reactive boundary conditions in moving-mesh hydrodynamics
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
Abstract not provided.
numerical study of laminar-flow heat trransfer enhancements
Absorbing boundary conditions contain differential operators even for first-order systems. There is a fundamental difficulty with this, since the conditions are applied on the ingoing variables, and the approximations necessarily become weakly unstable. This difficulty is more pronounced for inhomogeneous boundary conditions which occur if there is a source outside the computational domain D, or if the initial data are nonzero outside D. In this paper this is further investigated and it is shown that reasonable estimates can still be obtained if the solution is smooth. However, it is demonstrated that the approximations are less robust. A cure for this is proposed and an implementation of high-order conditions for first-order systems is described.
DRACO (Diffusion ReACtiOn) is a diffusion and chemistry code designed to: 1) Operate on 3D with an unstructured grid defining an arbitrary geometry of interacting parts. 2) Generate its own meshes and use meshes created by other software. 3)Model the transport of any number of diffusing quantities: Concentrations, pressures, temperature, etc. 4) Allow diffusion coefficients to depend in an arbitrary way on concentration, temperature, position, time, etc. 5) Model general chemistry between concentrations with arbitrary reaction rates. 6) Allow arbitrary initial conditions, boundary conditions, and sources/sinks. 7) Allow all of the above to be specified by the user.
ESCHER is an interactive mesh generation and editing program designed to help the user create a finite-element mesh, create additional input for finite-element analysis, including initial conditions, boundary conditions, and slidelines, and generate a NEUTRAL FILE that can be postprocessed for input into several finite-element codes, including ADINA, ADINAT, DYNA, NIKE, TSAAS, and ABUQUS. Two important ESCHER capabilities, interactive geometry creation and mesh archival storge are described in detail. Also described is the interactive command language and the use of interactive graphics. The archival storage and restart file is a modular, entity-based mesh data file. Modules of this file correspond to separate editing modes in the mesh editor, with data definition syntax preserved between the interactive commands and the archival storage file. Because ESCHER was expected to be highly interactive, extensive user documentation was provided in the form of an interactive HELP package.
An approximate analysis for buckling of a rectangular specially-orthotropic plate with a central circular hole is applied to symmetrically-laminated angle-ply plates. Results obtained from finite element analyses and experiments indicate that the approximate analysis predicts accurately the buckling loads of (+/-theta sub m)s plates with integer values of m not below 6 and with hole diameters up to 50 percent of the plate width. Moreover, the results indicate that the approximate analysis can be used to predict the buckling trends of plates with hole diameters up to 70 percent of the plate width. Results of a parametric study indicate the influence of hole size, plate aspect ratio, loading conditions, boundary conditions, and orthotropy on the buckling load. Results are also presented that indicate the relationship of the bending stiffness and the prebuckling load distribution to the buckling load of a plate with a hole.
In engineering, simulation software is often used as a convenient means for carrying out experiments to evaluate physical systems. The benefit of using simulations as 'numerical' experiments is that the experimental conditions can be easily modified and repeated at much lower cost than the comparable physical experiment. The goal of these experiments is to 'improve' the process or result of the experiment. In most cases, the computational experiments employ the same trial and error approach as their physical counterparts. When using this approach for complex systems, the cause and effect relationship of the system may never be fully understood and efficient strategies for improvement never utilized. However, it is possible when running simulations to accurately and efficiently determine the sensitivity of the system results with respect to simulation to accurately and efficiently determine the sensitivity of the system results with respect to simulation parameters (e.g., initial conditions, boundary conditions, and material properties) by manipulating the underlying computations. This results in a better understanding of the system dynamics and gives us efficient means to improve processing conditions. We begin by discussing the steps involved in performing simulations. Then we consider how sensitivity information about simulation results can be obtained and ways this information may be used to improve the process or result of the experiment. Next, we discuss optimization and the efficient algorithms which use sensitivity information. We draw on all this information to propose a generalized approach for integrating simulation and optimization, with an emphasis on software programming issues. After discussing our approach to simulation and optimization we consider an application involving crystal growth. This application is interesting because it includes radiative heat transfer. We discuss the computation of radiative new factors and the impact this mode of heat transfer has on our approach. Finally, we will demonstrate the results of our optimization.
The development of nuclear fuel and materials requires a continuous effort to investigate the acute and prolonged effects of irradiation, thermal-material stress, chemical change, or other conceivable damage mechanics acquired during normal operation or accident scenarios throughout its lifetime. As in-core fuel property measurement techniques advance to support in the real-time, non-invasive, and enhanced accuracy realm, it is the future of fuel development to pursue to higher degrees of control, predictability of integral test behavior via separate effects test (SET), and shorter test time intervals. The fuel development life cycle from initial concept to commercial licensing is approximated to be 20 years and current literature suggests by optimizing fuel performance codes with SETs, the process could possibly be compressed to 5 to 10 years. Recently, a research group in the Idaho National Laboratory (INL) is testing reduced scale fuel rods and increased power density to accelerate evolution of fuel phenomena in metallic fuels. In support of nuclear fuel rod development, compressing fuel test process, and accelerating fuel phenomena, it was the purpose of this study to investigate nuclear fuel performance phenomena via literature review and effectively scale the initial conditions, boundary conditions, and geometric properties to describe the time-dependent response including to fuel burnup, thermal and mechanical stress, transmutations and inter-diffusion, and other relevant observed phenomena. The study was based on BISON simulations of historic EBR-II metallic fuel experiments and Dynamical System Scaling (DSS) method are utilized to assess effects of scaling activity including fuel phenomena acceleration and calculations of time-dependent distortions. The research successfully scaled metallic fuel phenomena, accelerated fuel testing, and assessed the distortions for each scaled case derived.
The proper specification of boundary conditions at artificial boundaries for the simulation of time-dependent fluid flows has long been a matter of controversy. A general theory of asymptotic boundary conditions for dissipative waves is applied to the design of simple, accurate conditions at downstream boundary for incompressible flows. For Reynolds numbers far enough below the critical value for linear stability, a scaling is introduced which greatly simplifies the construction of the asymptotic conditions. Numerical experiments with the nonlinear dynamics of vortical disturbances to plane Poiseuille flow are presented which illustrate the accuracy of our approach. The consequences of directly applying the scalings to the equations are also considered.
During hypersonic flight, high temperatures and high heat fluxes are generated on the surfaces of vehicles. The Flight Loads Laboratory (FLL) at Dryden Flight Research Center (DFRC) is equipped with a calibration furnace, capable of calibrating heat flux gages up to 1100kW per square meters, and temperature sensors up to 2600 C. One heating configuration of the calibration furnace is a cylindrical blackbody cavity. Throughout the blackbody there are temperature gradients due to various boundary conditions. These boundary conditions include resistance heating, radiant heat transfer, and conduction to water-cooled electrodes. Also, an inert gas is purged through the graphite blackbody to prevent it from oxidizing. Consequently, the various modes of heat transfer present during operation of the blackbody cavity must be well understood in order to produce accurate heat flux gage and temperature sensor calibrations for use in ground testing or flight testing of hypersonic vehicles. The first step towards understanding the heat transfer in the blackbody cavity was to perform experiments at 1100 C, with and without outer surface insulation, while taking detailed temperature measurements inside the blackbody cavity. Steady state thermal models of the blackbody cavity were then developed. These models included detailed thermal analysis using commercial thermal analysis software. Conduction, radiation, and convection were considered in the thermal models for two cases: one with the outside of the blackbody cavity insulated and the second without insulation. This paper describes the experimental and numerical efforts used to characterize the steady state operation of the blackbody cavity. It describes the analysis of the test measurements, the boundary conditions used in the numerical models, and how the models were calibrated to fit the experimental data. Effects of various uncertainties, such as material properties, and convection are discussed.Initial thermal models predicted temperatures in the deepest part of the blackbody cavity within 7 C of the measured value and produced trends comparable to the experimental data, throughout the models. Adjustment of the boundary conditions, which were included in the thermal models, produced good agreement with measured temperatures. Free and forced convection of the purge gas inside the blackbody was found to be insignificant.
This paper discusses numerical solutions of a hyperbolic initial boundary value problem that arises from acoustic wave propagation in the atmosphere. Field equations are derived from the atmospheric fluid flow governed by the Euler equations. The resulting original problem is nonlinear. A first order linearized version of the problem is used for computational purposes. The main difficulty in the problem as with any open boundary problem is in obtaining stable boundary conditions. Approximate boundary conditions are derived and shown to be stable. Numerical results are presented to verify the effectiveness of these boundary conditions.
This paper discusses numerical solutions of a hyperbolic initial boundary value problem that arises from acoustic wave propagation in the atmosphere. Field equations are derived from the atmospheric fluid flow governed by the Euler equations. The resulting original problem is nonlinear. A first order linearized version of the problem is used for computational purposes. The main difficulty in the problem as with any open boundary problem is in obtaining stable boundary conditions. Approximate boundary conditions are derived and shown to be stable. Numerical results are presented to verify the effectiveness of these boundary conditions.