Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 595 records · Page 33

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Entanglement Renormalization for Quantum Field Theories with Discrete Wavelet Transforms

We propose an adaptation of Entanglement Renormalization for quantum field theories that, through the use of discrete wavelet transforms, strongly parallels the tensor network architecture of the Multiscale Entanglement Renormalization Ansatz (a.k.a. MERA). Our approach, called wMERA, has several advantages of over previous attempts to adapt MERA to continuum systems. In particular, (i) wMERA is formulated directly in position space, hence preserving the quasi-locality and sparsity of entanglers; and (ii) it enables a built-in RG flow in the implementation of real-time evolution and in computations of correlation functions, which is key for efficient numerical implementations. As examples, we describe in detail two concrete implementations of our wMERA algorithm for free scalar and fermionic theories in (1+1) spacetime dimensions. Possible avenues for constructing wMERAs for interacting field theories are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring regular and turbulent flow states in active nematic channel flow via Exact Coherent Structures and their invariant manifolds

This work is a unified study of stable and unstable steady states of 2D active nematic channel flow using the framework of Exact Coherent Structures (ECS). ECS are stationary, periodic, quasiperiodic, or traveling wave solutions of the governing equations that, together with their invariant manifolds, organize the dynamics of nonlinear continuum systems. We extend our earlier work on ECS in the preturbulent regime by performing a comprehensive study of stable and unstable ECS for a wide range of activity values spanning the preturbulent and turbulent regimes. In the weakly turbulent regime, we compute more than 200 unstable ECS that co-exist at a single set of parameters, and uncover the role of symmetries in organizing the phase space geometry. We provide conclusive numerical evidence that in the preturbulent regime, generic trajectories shadow a series of unstable ECS before settling onto an attractor. Lastly, our studies hint at shadowing of quasiperiodic type ECS in the turbulent regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A new rotational relaxation model for use in hypersonic computational fluid dynamics

The theoretical basis for the Landau-Teller equation commonly used to model rotational nonequilibrium is reviewed. Several assumptions underlying this model are indicated to be unrealistic for rarefied hypersonic flow. A new rotational nonequilibrium model based on recent measurements up to 2000K of quantum state transition rates is presented. The new model is applied to the continuum study of normal shock wave structure using both the Burnett equations and a simplified nonlinear constitutive relation. Comparisons are made between shock wave temperature profiles generated from the new model and those generated using the Landau-Teller model. Comparisons of shock reciprocal thicknesses between experimental data for nitrogen and continuum solutions using both rotational models are made. The new rotational model agrees well with experiment up to Mach 6, and under predicts shock thickness at higher Mach numbers. The Landau-Teller model agrees well with experimental shock thickness up to Mach 3, and over predicts the thickness at higher Mach numbers. A modification to the rotational collision number in the Landau-Teller model is found to give results which agree with experimental shock thicknesses at all Mach numbers up to 11.

Chapman, Dean E.↗

Computation of atmospheric cooling rates by exact and approximate methods

Infrared fluxes and cooling rates for several standard model atmospheres, with and without water vapor, carbon dioxide, and ozone, have been calculated using a line-by-line method at 0.01/cm resolution. The sensitivity of the results to the vertical integration scheme and to the model for water vapor continuum absorption is shown. Comparison with similar calculations performed at NOAA/GFDL shows agreement to within 0.5 W/sq m in fluxes at various levels and 0.05 K/d in cooling rates. Comparison with a fast, parameterized radiation code used in climate models reveals a worst case difference, when all gases are included, of 3.7 W/sq m in flux; cooling rate differences are 0.1 K/d or less when integrated over a substantial layer with point differences as large as 0.3 K/d.

Ridgway, William L.↗

Comprehensive Micromechanics-Analysis Code - Version 4.0

Version 4.0 of the Micromechanics Analysis Code With Generalized Method of Cells (MAC/GMC) has been developed as an improved means of computational simulation of advanced composite materials. The previous version of MAC/GMC was described in "Comprehensive Micromechanics-Analysis Code" (LEW-16870), NASA Tech Briefs, Vol. 24, No. 6 (June 2000), page 38. To recapitulate: MAC/GMC is a computer program that predicts the elastic and inelastic thermomechanical responses of continuous and discontinuous composite materials with arbitrary internal microstructures and reinforcement shapes. The predictive capability of MAC/GMC rests on a model known as the generalized method of cells (GMC) - a continuum-based model of micromechanics that provides closed-form expressions for the macroscopic response of a composite material in terms of the properties, sizes, shapes, and responses of the individual constituents or phases that make up the material. Enhancements in version 4.0 include a capability for modeling thermomechanically and electromagnetically coupled ("smart") materials; a more-accurate (high-fidelity) version of the GMC; a capability to simulate discontinuous plies within a laminate; additional constitutive models of materials; expanded yield-surface-analysis capabilities; and expanded failure-analysis and life-prediction capabilities on both the microscopic and macroscopic scales.

Arnold, S. M.↗

Implementation of Surface Tension on a Reacting Flow Solver, PeleLM: Preprint

In liquid rocket engines, the fuel is supplied to the combustion chamber in the liquid state though injectors. Such fuel undergoes atomization, vaporization, and combustion processes. To design reliable and efficient injectors, it is required to understand the full processes. This research is part of an effort to develop a full atomization-vaporization-combustion solver from first principles. As an initial step to tackle the atomization process, a multiphase flow solver is under development. For the development, a library of the volume of fluid scheme for multiphase, IRL is coupled with a reacting Navier-Stokes equation solver, PeleLM. Furthermore, as the surface tension has considerable effects on spray breakup. surface tension is implemented in the momentum equation using the continuum surface force model and the improved height function technique.

height function↗

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Co-design for Particle Applications at Exascale

Co-design across the Exascale Computing Project (ECP) has been critical for both enabling science applications and bringing disparate communities together. Developing and porting applications to the various high-performance computing (HPC) architectures on pre-exascale and exascale computers has been quite challenging due to the diversity of hardware features and software stacks. The Co-design Center for Particle Applications (CoPA) has developed and enhanced the Cabana and PROGRESS/BML libraries to facilitate the creation of new particle applications, make existing particle applications exascale capable, and allow teams to explore new capabilities. Particle methods from atomistic, mesoscale, continuum, through cosmological scales have been built with Cabana, along with new possibilities for application coupling. Similarly, the PROGRESS/BML library has enabled quantum particle applications with linear algebra solvers to use advanced hardware. Across these CoPA-developed libraries, the co-design abstraction layer combines performance portability with math library support to facilitate separation of concerns and directly support science runs.

97 MATHEMATICS AND COMPUTING↗

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗

Some experience with arc-heater simulation of outer planet entry radiation

An electric arc heater was operated at 800 amperes and 100,000 pa (1 atm) with hydrogen, helium, and two mixtures of hydrogen and helium. A VUV-scanning monochromator was used to record the spectra from an end view while a second spectrometer was used to determine the plasma temperature using hydrogen continuum radiation at 562 nm. Except for pure helium, the plasma temperature was found to be too low to produce significant helium radiation, and the measured spectra were primarily the hydrogen spectra with the highest intensity in the pure hydrogen case. A radiation computer code was used to compute the spectra for comparison to the measurements and to extend the study to simulation of outer planet entry radiation. Conductive cooling prevented ablation of phenolic carbon material samples mounted inside the arc heater during a cursory attempt to produce radiation absorption by ablation gases.

Wells, W. L.↗

Formation of the O I resonance triplet and intercombination doublet in the solar chromosphere

Spectrum synthesis calculations are presented for the O I resonance triplet at 1304 A and the intercombination doublet at 1358 A for the solar atmosphere and several variants, allowing for triplet fluorescence by the Ly-beta emission of H I. Profiles, synthesized from a seven-level plus continuum O I atom are compared to observations taken with the high-resolution spectrometer on OSO 8. It is found that the O I triplet emission is dominated by the Ly-beta fluorescence and that the agreement between observations and profiles computed with current chromospheric models is much improved over earlier studies.

Skelton, D. L.↗

Application of a global solar wind/planetary obstacle interaction computational model: Earth, Venus, Mars, Jupiter and Saturn studies

This work relates to a series of collaborative investigations involving the application of a computational model for the determination of the detailed plasma and magnetic field properties associated with the global interaction of the solar wind with various planetary obstacles throughout the solar system. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past planetary obstacles. The investigations undertaken relate to studies of various solar wind interaction phenomena with Venus, Earth, Jupiter, and Saturn.

Stahara, S. S.↗

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗

Turbomachinery

The discipline research in turbomachinery, which is directed toward building the tools needed to understand such a complex flow phenomenon, is based on the fact that flow in turbomachinery is fundamentally unsteady or time dependent. Success in building a reliable inventory of analytic and experimental tools will depend on how the time and time-averages are treated, as well as on who the space and space-averages are treated. The raw tools at disposal (both experimentally and computational) are truly powerful and their numbers are growing at a staggering pace. As a result of this power, a case can be made that a situation exists where information is outstripping understanding. The challenge is to develop a set of computational and experimental tools which genuinely increase understanding of the fluid flow and heat transfer in a turbomachine. Viewgraphs outline a philosophy based on working on a stairstep hierarchy of mathematical and experimental complexity to build a system of tools, which enable one to aggressively design the turbomachinery of the next century. Examples of the types of computational and experimental tools under current development at Lewis, with progress to date, are examined. The examples include work in both the time-resolved and time-averaged domains. Finally, an attempt is made to identify the proper place for Lewis in this continuum of research.

Simoneau, Robert J.↗

Accuracy of the Burnett equations for hypersonic real gas flows

Burnett solutions for the structure of normal shock waves for Mach numbers ranging from 1.1 to 50 have been obtained by using the time-dependent equation set in conjunction with modern computational fluid dynamic techniques. The current work seeks to provide corrected and improved high-resolution direct simulation Monte Carlo (DSMC) results for hypersonic shocks in monatomic gases. These DSMC results will serve as a reliable test base for assessing the accuracy of various continuum equations more advanced than Navier-Stokes, such as the Burnett equations. The main conclusion of Fisko and Chapman (1988) that the Burnett equations are more accurate than the Navier-Stokes equations for all gases and at all Mach numbers in matching DSMC shock wave results is reassesed by comparing the improved and corrected DSMC results with continuum Burnett results. Various Burnett results for shocks in a diatomic gas, nitrogen, which include improvements to the work of Lumpkin and Chapman (1989) are presented. It is concluded that the Burnett equations represent a significant improvement over Navier-Stokes equations for flows exhibiting translational nonequilibrium.

Lumpkin, Forrest E., III↗