Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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,711 records · Page 95

Analytical methods for the development of Reynolds-stress closures in turbulence

The derivation of Reynolds-stress models for viscous incompressible turbulent flow on the basis of the Navier-Stokes and continuity equations is explored in an analytical review. The formulation of the basic equations is outlined, and particular attention is given to zero-equation and one-equation models based on eddy viscosity, two-equation (k-l, k-epsilon, and k-omega) models, and second-order closure models. The superior performance of the second-order models is demonstrated by numerical results for (1) the return to isotropy from anisotropic homogeneous turbulence, (2) homogeneous turbulent shear flow in a rotating frame, (3) and fully developed inhomogeneous turbulent channel flow in a rotating frame. The inherent limitations of the Reynolds-stress approach and areas for further improvement are discussed.

Speziale, Charles G.↗

Theoretical studies of lasers and converters

Discussed here is Doppler broadening and its effects upon the stimulated emission cross-section connecting an upper level with a lower level in an iodine laser. The stimulated emission cross section is given. The level transitions for the iodine laser are illustrated as are the relative intensities of these transitions. Also discussed is the Voigt profile, which considers the effect of both Doppler broadening and collision broadening upon the absorption line shape. The equations describing the Voigt profile were added to the continuous flow laser model laser simulation program. The results were compared with the standard absorption profile reported in an earlier study. There seems to be no advantage to using the Voigt profile as the laser power output is relatively insensitive to changes in the absorption cross section at the pressures being considered for a space laser. One disadvantage of using the Voigt profile is the excessive numerical computations required by the additional equations.

Heinbockel, John H.↗

Evolving the Reuse Process at the Flight Dynamics Division (FDD) Goddard Space Flight Center

This paper presents the interim results from the Software Engineering Laboratory's (SEL) Reuse Study. The team conducting this study has, over the past few months, been studying the Generalized Support Software (GSS) domain asset library and architecture, and the various processes associated with it. In particular, we have characterized the process used to configure GSS-based attitude ground support systems (AGSS) to support satellite missions at NASA's Goddard Space Flight Center. To do this, we built detailed models of the tasks involved, the people who perform these tasks, and the interdependencies and information flows among these people. These models were based on information gleaned from numerous interviews with people involved in this process at various levels. We also analyzed effort data in order to determine the cost savings in moving from actual development of AGSSs to support each mission (which was necessary before GSS was available) to configuring AGSS software from the domain asset library. While characterizing the GSS process, we became aware of several interesting factors which affect the successful continued use of GSS. Many of these issues fall under the subject of evolving technologies, which were not available at the inception of GSS, but are now. Some of these technologies could be incorporated into the GSS process, thus making the whole asset library more usable. Other technologies are being considered as an alternative to the GSS process altogether. In this paper, we outline some of issues we will be considering in our continued study of GSS and the impact of evolving technologies.

Condon, S.↗

Simulation of Delamination-Migration and Core Crushing in a CFRP Sandwich Structure

Following the onset of damage caused by an impact load on a composite laminate structure, delaminations often form propagating outwards from the point of impact and in some cases can migrate via matrix cracks between plies as they grow. The goal of the present study is to develop an accurate finite element modeling technique for simulation of the delamination-migration phenomena in laminate impact damage processes. An experiment was devised where, under a quasi-static indentation load, an embedded delamination in the facesheet of a laminate sandwich specimen migrates via a transverse matrix crack and then continues to grow on a new ply interface. The quasistatic nature of the indentation results in structural behavior equivalent to that seen in low-velocity impact and also allows for highly detailed real time damage characterization. Several finite element damage simulation methods were investigated. Comparing the experimental results with those of the different models reveals certain modeling features that are important to include in a numerical simulation of delamination-migration and some that may be neglected.

McElroy, M.↗

Flow and Noise From Septa Nozzles

This investigation pertains to distributed propulsion system. Basically, the jet exhaust is divided and channeled through multiple compartments or septa. In one version of the concept each septum is to be driven by an electric fan, the power being supplied by engines suitably mounted on the airframe. In the latter hybrid case the need for heavy power trains is eliminated and system studies indicate a net benefit in efficiency. From aeroacoustics point of view, however, there are concerns. An immediate question is whether the noise from the multiple jets would be equal to or greater than that from an equivalent single jet. When multiple jets are placed side by side sometimes they resonate and produce more noise. In order to address these issues a fundamental model-scale experiment is being carried out at NASA Glenn Research Center (GRC). Preliminary results from this study have been presented recently. One of the key findings is that the noise from the septa nozzle can actually be less than that from the equivalent single nozzle. Continued studies make it clear that the flow from the septa nozzle may be quite sensitive to upstream passage geometry in each septum. This can have a profound impact on the evolution of the flow downstream as well as noise. These issues are currently being researched. The proposed paper is to summarize the results from the mainly experimental study that is to be aided by limited numerical simulation.

nozzle flow↗

Case study of a low-reflectivity pulsating microburst - Numerical simulation of the Denver, 8 July 1989, storm

An intense, low-reflectivity, and pulsating microburst event, detected on July 8, 1989 by the Low-Level Windshear Alert System north of Denver Stapleton Airport is described using a Terminal Area Simulation System (TASS) model (Proctor, 1987). This 3d nonhydrostatic cloud model includes parametrizations for both liquid- and ice-phase microphysics. Simulation results show that the dissipating cloud induces a type-1 pulsating microburst event, with at least three distinct microburst pulses over a period of 20 min. Virga from the dissipating cloud continued to maintain hazardous windshear for a period much longer than the typical lifetime of a single microburst.

Proctor, Fred H.↗

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

42 ENGINEERING↗

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

Here we present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A method for screening of plant species for space use

A cost-effective methodology which monitors numerous dynamic aspects of carbon assimilation and allocation kinetics in live, intact plants is discussed. Analogous methods can apply to nitrogen uptake and allocation. This methodology capitalizes on the special properties of the short-lived, positron-gamma emitting isotope C-11 especially when applied as CO2-11 in a special extended square wave (ESW) pattern. The 20.4 minute half-life allows for repeated or continuous experiments on the same plant over periods of minutes, hours, days, or weeks. The steady-state isotope equilibrium approached during the ESW experiments, and the parameters which can be analyzed by this technique are also direct results of that short half-life. Additionally, the paired .511 MeV gamma rays penetrate any amount of tissue and their 180 deg opposite orientation provides good collimation and allows coincidence counting which nearly eliminates background.

Goeschl, J. D.↗

Comparison results for ignition in conjugate systems

The exothermic reaction of an arbitrarily shaped solid in a well-mixed continuous flow of constant density and velocity is investigated analytically, with a focus on a zeroth-order reaction in a pellet. Comparison theorems for the basic problem and the eigenvalue problem of ignition are derived, and the domain dependence of the bifurcation point is explored using the method of Joseph (1976). Some typical analytical results are presented in tables and compared with published values obtained using numerical methods.

Pushpavanam, S.↗

Analytical model for the density distribution in the Io plasma torus

An analytical model is developed for the diffusive equilibrium plasma density distribution in the Io plasma torus. The model has been employed successfully to follow the ray path of plasma waves in the multi-ion Jovian magnetosphere; it would also be valuable for other studies of the Io torus that require a smooth and continuous description of the plasma density and its gradients. Validity of the analytical treatment requires that the temperature of thermal electrons be much lower than the ion temperature and that superthermal electrons be much less abundant than the thermal electrons; these two conditions are satisfied in the warm outer region of the Io torus from L = 6 to L = 10. The analytical solutions agree well with exact numerical calculations for the most dense portion of the Io torus within 30 deg of the equator.

Mei, YI↗

Large Volume Airborne Contamination Monitoring To Support Nuclear Processes' Deactivation and Decommissioning

Current D and D operations at Hanford have demonstrated a flaw in the current state of the art capability of defining airborne contamination boundaries - Airborne particulate emissions of Pu-239 from CM2H operations on the Hanford Pu Finishing Plant were detected well beyond areas controlled for airborne Pu, putting numerous workers at risk for radiological assimilations - Current air samplers and CAM systems were surveying insufficient volumes of air to accurately predict where respiratory protection was required - Hanford is located in a high radon environment, making air monitoring for alpha emitting actinides challenging in high alpha radon induced backgrounds. - Inexpensive HEPA based home and industrial air purifiers filter significantly higher volumes of air than commercially available continuous-air monitoring (CAM) systems. - Inexpensive models capable of filtered air volumes exceeding 1000 times that of a CAM. - Detector system can be built to detect x-rays generated from actinide decay. - Branching ratios of x-rays are 4 orders of magnitude more intense from Plutonium decay than its gamma emissions. - Detector system can be a low resolution systems, as x-ray region of interest is not terribly congested. - Detection systems evaluated will be designed based on slabs of NaI or pixelated NaI or CsI panels. - Simple graphical user interface software will be developed to operate the detection system. - Concept is to deploy some of these air-purifying systems in nonradiological areas around SRNL to confirm radon can be rejected with confidence. - Generate some filters contaminated with plutonium via electroplated Pu or lab generated simulated particles. - then deploy units in known airborne radiological environments to establish systems ability to accurately measure actinide based hot particles Lab analyses will follow up the system analyses to ensure hot particles were correctly identified. Analysis by Scintillation: - Application of scintillation media to the surface of the filters is being explored. - Evaluating slabs of ZnS(Ag), application of powdered ZnS(Cu) and spray on Perkin Elmer Enhance. - Scintillation events would then be digitized with a digital camera and quantified. - Comparing against sensitivity of a PMT or SiPD readout. - Currently evaluating a Thorlabs 8 Megapixel Monochrome Scientific CCD Camera, hermetically sealed cooled package with a wide angle lens. - Wide angle lens allows complete view of HEPA filter from 7 inches away. - Images taken of glow in the dark paint, as well as plutonium-induced fluorescence. - Plutonium was flamed mounted on a 1 inch diameter stainless steel planchet, covered with a layer of mylar and a section of scintillating ZnS obtained from Eljon. - Currently working on reducing signal-to-noise levels to boost sensitivity. - Dark box to hold contaminated or electroplated source covered filters under fabrication. Analysis by x-ray Spectroscopy: - Branching ratios for x-ray emissions from Plutonium isotopes are orders of magnitude more intense than gamma-ray emissions. - The nuclear databases are incomplete on the branching ratios of some of the isotopes. - Evaluating some custom built SrI2 x-ray spectrometers vs a conventional windowed NaI detector. - SrI2 detectors are carbon-composite-windowed 51 mm x 51 mm. - MCNP calculations to establish self absorbance of filter media on 17 keV x-ray. Radon Rejection: - Activated charcoal pre-filter acts as a Radon trap. - Adding a time lapse feature to camera to aid in radon rejection. SRS Plutonium Fuel Form D and D Operations: Currently deploying air sampler to D and D operation of Pu-238 facility at SRS to generate some field samples to analyze.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Molten Salt Reactors: What about Materials? An Overview and Future Needs

ORNL Molten Salt Reactor Experiment (MSRE) – 60’s- early 70’s Developed INOR-8 (aka Hastelloy N); MSRE operated over the period Dec 1965 to 1969. (not continuous); Operated for >13,000 hr at T from ~550 to ~650°C using U, Li-7, Be and Zr fluoride salt fuel; Post-Irradiation Material Evaluation (corrosion + neutron irradiation) of surveillance specimens; Noted ”beneficial effect of Nb” in alloy (on intergranular Te embrittlement) More recent (post-1980) Materials R&D for Gen IV …. Alloys 316L, 316H, 304, Alloy 800H, Alloy 617 (!!), Hastelloy N (mod), plus numerous experimental alloys developed

36 - MATERIALS SCIENCE↗

Numerical, micro-mechanical prediction of crack growth resistance in a fibre-reinforced/brittle matrix composite

Micromechanics fracture models are incorporated into three distinct fracture process zones which contribute to the crack growth resistance of fibrous composites. The frontal process zone includes microcracking, fiber debonding, and some fiber failure. The elastic process zone is related only to the linear elastic creation of new matrix and fiber fracture surfaces. The wake process zone includes fiber bridging, fiber pullout, and fiber breakage. The R-curve predictions of the model compare well with empirical results for a unidirectional, continuous fiber C/C composite. Separating the contributions of each process zone reveals the wake region to contain the dominant crack growth resistance mechanisms. Fractography showed the effects of the micromechanisms on the macroscopic fracture behavior.

Jenkins, Michael G.↗

Numerical Modeling of Physical Vapor Transport in Contactless Crystal Growth Geometry

Growth from the vapor under conditions of limited contact with the walls of the growth ampoule is beneficial for the quality of the growing crystal due to reduced stress and contamination which may be caused by interactions with the growth container. The technique may be of a particular interest for studies on crystal growth under microgravity conditions: elimination of some factors affecting the crystal quality may make interpretation of space-conducted processes more conclusive and meaningful. For that reason, and as a part of our continuing studies on 'contactless' growth technique, we have developed a computational model of crystal growth process in such system. The theoretical model was built, and simulations were performed using the commercial computational fluid dynamics code, (CFD) ACE. The code uses an implicit finite volume formulation with a gray discrete ordinate method radiation model which accounts for the diffuse absorption and reflection of radiation throughout the furnace. The three-dimensional model computes the heat transfer through the crystal, quartz, and gas both inside and outside the ampoule, and mass transport from the source to the crystal and the sink. The heat transport mechanisms by conduction, natural convection, and radiation, and mass transport by diffusion and convection are modeled simultaneously and include the heat of the phase transition at the solid-vapor interfaces. As the thermal boundary condition, temperature profile along the walls of the furnace is used. For different thermal profiles and furnace and ampoule dimensions, the crystal growth rate and development of the crystal-vapor and source-vapor interfaces (change of the interface shape and location with time) are obtained. Super/under-saturation in the ampoule is determined and critical factors determining the 'contactless' growth conditions are identified and discussed. The relative importance of the ampoule dimensions and geometry, the furnace dimensions and its temperature, and the properties of the grown material are analyzed. The results of the simulations are compared with related experimental results on growth of CdTe, CdZnTe, ZnTe, PbTe, and PbSnTe crystals by this technique.

Palosz, W.↗

Lignin-derived electrochemical energy materials and systems

Electrochemical energy storage systems such as supercapacitors, rechargeable batteries and fuel cells have been proven the most effective technologies for energy conversion, storage, and management at different scales. Although a large number of electrochemical energy technologies have been developed in the past and they will continue to be optimized in terms of cost, lifetime, and performance, there is a substantial growing demand for advanced electrochemical energy systems. To deploy these advanced systems, the electrode and electrolyte materials with higher performance, longer life, and lower cost, must be developed. Lignin is the second most abundant natural polymer after cellulose, a byproduct from emerging cellulosic biorefineries, and a waste product from pulp and paper industries. Numerous researches have successfully demonstrated that lignin from different sources can be used as precursors or feedstocks for preparing high-performance electrochemical energy materials and components such as electrodes, electrolytes, membrane separators, and additives. Moreover, techno-economic analyses indicate that it is possible to prepare cost-effective carbons from lignin at engineering scales, compared to current carbon products. These facts suggest that scalable conversion of lignin into high-value energy materials will offer a promising pathway to not only promote the utilization and valorization of lignin but also boost the development of the advanced electrochemical energy systems. This review presents state of the arts of renewable energy materials derived from various lignin and their applications in electrochemical energy systems with emphasis on supercapacitors, rechargeable batteries, and fuel cells. Meanwhile, this article also aims to carve out the critical barriers for lignin-derived high-performance materials for energy applications, intending to identify viable approaches for synthesis of sustainable new energy materials.

09 BIOMASS FUELS↗