SEARCH · Engineering Papers

Results for “HOMOGENIZATION METHODS”

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.

SourceResearch connectionFrom yearThrough year

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

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Analytic Development of a Reference Profile for the First Entry in a Skip Atmospheric Entry

This note shows that a feasible reference drag profile for the first entry portion of a skip entry can be generated as a polynomial expression of the velocity. The coefficients of that polynomial are found through the resolution of a system composed of m + 1 equations, where m is the degree of the drag polynomial. It has been shown that a minimum of five equations (m = 4) are required to establish the range and the initial and final conditions on velocity and flight path angle. It has been shown that at least one constraint on the trajectory can be imposed through the addition of one extra equation in the system, which must be accompanied by the increase in the degree of the drag polynomial. In order to simplify the resolution of the system of equations, the drag was considered as being a probability density function of the velocity, with the velocity as a distribution function of the drag. Combining this notion with the introduction of empirically derived constants, it has been shown that the system of equations required to generate the drag profile can be successfully reduced to a system of linear algebraic equations. For completeness, the resulting drag profiles have been flown using the feedback linearization method of differential geometric control as a guidance law with the error dynamics of a second order homogeneous equation in the form of a damped oscillator. Satisfactory results were achieved when the gains in the error dynamics were changed at a certain point along the trajectory that is dependent on the velocity and the curvature of the drag as a function of the velocity. Future work should study the capacity to update the drag profile in flight when dispersions are introduced. Also, future studies should attempt to link the first entry, as presented and controlled in this note, with a more standard control concept for the second entry, such as the Apollo entry guidance, to try to assess the overall skip entry performance. A guidance law that includes an integral feedback term, as is the case in the actual Space Shuttle entry guidance and as is proposed in Ref 29, could be tried in future studies to assess whether its use results in an improvement of the tracking performance, and to evaluate the design needs when determining the control gains.

Garcia-Llama, Eduardo↗

NASA NTRS2004

Structural Evaluation of Radially Expandable Cardiovascular Stents Encased in a Polyurethane Film

A method of encasing cardiovascular stents with an expandable polyurethane coating has been developed to provide a smooth homogeneous inner wall allowing for a confluent growth of endothelial cells. In this design, the metal wire stent structure is completely covered by the polyurethane film minimizing biocorrosion of the metal (stainless steel or nitinol), and providing a homogeneous surface for surface treatment and incorporation of various eluting drugs to prevent platelet aggregation while supporting endothelialization. The polyurethane surface was treated with a helium plasma for sterilization and promotes growth of cells. The paper details the performance of the coated film to expand with the metal stent up to 225 % during deployment. We present stress/strain behavior of polyurethane films, and subsequent plasma treatment of the surface and the adhesion of the coating to the stent structure upon expansion. A film of less than 25 tm was found to be sufficient for corrosion resistance and flexibility without producing any excess stress on the stent structure. Straining the film to 225 % and plasma modification did not affect the mechanical and surface properties while allowing for improved biocompatibility as determined by the critical surface tension, surface chemistry, and roughness.

Trigwell, Steve↗

DOE OSTI2020

Improvements to the Modeling of the TREAT Reactor and Experiments

This paper summarizes the latest improvements and lessons learned from the modeling and simulation of the transient test reactor at Idaho National Laboratory using the MAM- MOTH reactor physics application. MAMMOTH is a MOOSE-based, Finite Element Method application that specializes in the analysis of the spatial dynamics behavior of nuclear reactors. Since early 2018 several transient tests have been conducted at TREAT, thus providing the opportunity to apply and benchmark modern modeling and simulation tools. MAMMOTH was used to provide predictions of the power coupling factor between the core and the experiment for various experiments. Even though the power coupling factor predictions agree very well with the experimental data, within the bounds of the experimental uncertainty, one shortcoming was the underprediction of the total energy deposited in the core and in the experiment. Determination of the sources for this discrepancy are ongoing, but several key problems have been identified and resolved, thus providing valuable insights for future research. This paper discusses several of these lessons learned. First, the heat capacity data for the TREAT fuel has some significant problems due to limitations of the measurement techniques used circa 1960s. The sensitivity of the peak power and the total energy deposition to various representations of the heat capacity is approximately 5%. Second, the effects of the biological shield and thermal column on the modeling of the core are non-negligible, since they affect the mean generation time and the effective reflection of neutrons back into the core, which is suspected to be important during the core heat up. Matching the reactor period resolves the fact that the reduced spatial domain used in the MAMMOTH model underpredicts the mean generation time. The neutron reflection from these regions is marginally improved with the use of an albedo boundary condition. Third, modeling of the control rod movement with a multi-scheme method is introduced and its current limitations are exposed. Fourth, we explore the effects of using a homogenized model with Superhomogenization equivalence and how that differs from fully heterogeneous simulations. Finally, the energy condensation effects for this graphite core are significant. Solutions with 10 and 26 energy groups show the benefits of using a finer coarse group structure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1978

Computation of infrared cooling rates in the water vapor bands

A fast but accurate method for calculating the infrared radiative terms due to water vapor has been developed. It makes use of the far wing approximation to scale transmission along an inhomogeneous path to an equivalent homogeneous path. Rather than using standard conditions for scaling, the reference temperatures and pressures are chosen in this study to correspond to the regions where cooling is most significant. This greatly increased the accuracy of the new method. Compared to line by line calculations, the new method has errors up to 4% of the maximum cooling rate, while a commonly used method based upon the Goody band model (Rodgers and Walshaw, 1966) introduces errors up to 11%. The effect of temperature dependence of transmittance has also been evaluated; the cooling rate errors range up to 11% when the temperature dependence is ignored. In addition to being more accurate, the new method is much faster than those based upon the Goody band model.

Chou, M. D.↗

NASA NTRS1976

Condensation of water vapor in rarefaction waves. I - Homogeneous nucleation

A detailed theoretical investigation has been made of the condensation of water vapor/carrier gas mixtures in the nonstationary rarefaction wave generated in a shock tube. It is assumed that condensation takes place by homogeneous nucleation. The equations of motion together with the nucleation rate and the droplet growth equations were solved numerically by the method of characteristics and Lax's method of implicit artificial viscosity. It is found that, for the case considered, the condensation wave formed by the collapse of the metastable nonequilibrium state is followed by a shock wave generated by the intersection of characteristics of the same family. The expansion is practically isentropic up to the onset of condensation. The condensation front accelerates in the x,t plane. The results of the computations for a chosen case of water vapor/nitrogen mixture are presented by plotting variations of pressure, nucleation rate, number density of critical clusters, and condensate mass-fraction along three particle paths. Some consideration is given to homogeneous condensation experiments conducted in a shock tube. Although a direct comparison of the present theoretical work and these experiments is not possible, several worthwhile interpretative features have resulted nevertheless.

Sislian, J. P.↗

NASA NTRS1990

Development of Renormalization Group Analysis of Turbulence

The renormalization group (RG) procedure for nonlinear, dissipative systems is now quite standard, and its applications to the problem of hydrodynamic turbulence are becoming well known. In summary, the RG method isolates self similar behavior and provides a systematic procedure to describe scale invariant dynamics in terms of large scale variables only. The parameterization of the small scales in a self consistent manner has important implications for sub-grid modeling. This paper develops the homogeneous, isotropic turbulence and addresses the meaning and consequence of epsilon-expansion. The theory is then extended to include a weak mean flow and application of the RG method to a sequence of models is shown to converge to the Navier-Stokes equations.

L M Smith↗

DOE OSTI2026

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2022

Multiscale and multiphysics FEA simulation and materials optimization for laser ultrasound transducers

In this study, the relationship between the nanocomposite design and the laser ultrasound transducer (LUT) characteristics was investigated through simulations in multiple scale levels for material behavior, device response, and acoustic wave propagation in media. First, the effects of the nanoparticle size and concentration on the effective properties of composites were quantitatively investigated with the finite element analysis (FEA) method. Second, the effective properties of the nanocomposite were assigned to the layer, which is modeled as a homogeneous material, in the FEA for the LUT simulating the energy conversion from the incident laser to the acoustic wave. Finally, the ultrasound propagation in the water was calculated by a theoretical wave propagation model. The FEA-based prediction was compared with the experimental data in the literature and a theoretical analysis for LUT based on Thermal-Acoustic coupling. As a result, the ultrasound waves on the transducer surface and at a distance in the water could be predicted. Based on the hierarchically integrated prediction procedure, the optimal conditions of the photoacoustic nanocomposites were investigated through the parametric study with the particle size and concentration as variables. The results guide the material designs optimized for different device characteristics, such as high pressure and broad bandwidth.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Design and simulation of multilayer coatings for a multi-channel Wolter-like x-ray imager with large field of view and high resolution

X-ray diagnostics are key instruments for understanding the physics behind inertial confinement fusion experiments. We report on the multilayer design optimization for the Toroidal X-ray Imager (TXI), a hard x-rays microscope instrument designed by Commissariat à l’énergie atomique (CEA) and Laboratoire Charles Fabry (LCF) to be installed on the National Ignition Facility. TXI includes six channels designed for three different energy bands centered on 8.7, 13, and 17.5 keV. Each channel is made up of two toroidal mirrors arranged in a Wolter-like configuration. The required field of view is 800 × 400 µm2, and the resolution should be better than 5 µm. In addition, we seek to estimate the spatial distribution of the temperature, which requires no spectral overlap of the different energy bands and a good spectral homogeneity of the image produced. The development of the multilayer coatings was performed in a two-step method. First, the coatings were optimized to obtain proper energy bands. Then, an x-ray tracing code was used to calculate the integrated optical response of each channel and adjust the response of the mirror to fulfill the requirements. To fulfill all the specifications, we propose an original design using a combination of two aperiodic coatings, one with a narrow bandwidth and the other one with a larger bandwidth.

hard x-rays↗

DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Evaluating Effects of Shear Processing on 2D Crystalline Materials in 3D Metal Matrices: Atomistic Understanding of High Shear Deformation of Copper Graphene Composites

Composites made of copper and graphene demonstrate high strength, lubricity and enhanced electrical and thermal conductivity compared pure copper. However, manufacturing these composites at bulk volumes for industrial applications has been a big challenge. Shear deformation assisted processing is an effective method for manufacturing materials such as copper-graphene composites demonstrating ultra-fine grain structures and compositional homogeneity. Nevertheless, microstructural evolution of the composites and their property development under such conditions is not clearly understood currently. To rectify this gap in literature, high strain shear deformation of copper coated graphene foils was performed using a tribometer pin in this study. Changes in microstructure of the composite as well as the constituent components under shear deformation was correlated to process conditions. A sharp increase in the coefficient of friction attributed to rupture and smearing of graphene layer into copper substrate was observed during the shear processing. The coefficient of friction of the sheared copper/graphene composite was lower than that of pure copper, suggesting that partially worn graphene is effectively lubricious at the macroscale. A multimodal characterization of the processed region further revealed a shear deformation-induced ultrafine two-phase grain-structure consisting of copper and graphitic domains. Shear deformation reduced the copper grain size from around 50 – 100 µm to ~200 nm on an average and ~2 – 5 nm in some locations. The semicrystalline graphene films were observed to fracture into 10 – 50 µm size flakes. Oxygen enrichment was observed in the processed region. Finally, graphitic domain were identified for the first time in the copper matrix and not just at the grain boundaries providing evidence for a metastable composite microstructure as a result of solid phase processing at room temperature.

36 MATERIALS SCIENCE↗

NASA NTRS1980

A fast and accurate radiance algorithm for application to inversion of limb measurements

A method is devised to accurately calculate broad band limb radiance based on a precalculated data base consisting of homogeneous path emissivity as a function of mass path data for various temperatures and pressures. Accurately estimating limb radiance by splicing together interpolated information from the emissivity curves requires no temperature or pressure weighting, and only two interpolations per layer path. The emissivity growth approximation (EGA) technique is applied in a simulated inversion study to the inference of O3, H2O, NO2 and HNO3, and it is found to have comparable accuracies to more sophisticated many-degree band models. The EGA technique is also found to be two orders of magnitude faster, while requiring one second per inversion for one constituent at one-km resolution, allowing a completely nonlinear relaxation time inversion to be done at linear speeds. The technique is also expected to work well for nadir calculations, and prove adequate for estimating occultation radiances.

Gordley, L. L.↗

NASA NTRS1987

Exact image method for Gaussian beam problems involving a planar interface

Exact image method, recently introduced for the solution of electromagnetic field problems involving sources above a planar interface or two homogeneous media, is shown to be valid also for sources located in complex space, which makes its application possible for Gaussian beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular shift of a TE polarized beam are correctly given as asymptotic results by the exact reflection image theory. Also, the apparent image location giving the correct Gaussian beam transmitted through the interface is obtained as another asymptotic check. The present theory makes it possible to calculate the exact coupling from the Gaussian beam to the reflected and refracted beams, as well as to the surface wave.

Lindell, I. V.↗

NASA NTRS1987

Exact-image method for Gaussian-beam problems involving a planar interface

The exact-image method, recently introduced for the solution of electromagnetic field problems involving sources above a planar interface between two homogeneous media, is shown to be valid also for sources located in complex space, which makes its application possible for Gaussian-beam analysis. It is demonstrated that the Goos-Hanchen shift and the angular shift of a TE-polarized beam are correctly given as asymptotic results by the exact-reflection-image theory. Also, the apparent-image location giving the correct Gaussian beam transmitted through the interface is obtained as another asymptotic check. The theory described here makes it possible to calculate the exact coupling from the Gaussian beam to the reflected and refracted beams as well as to the surface wave.

Lindell, Ismo V.↗

NASA NTRS1992

Reactant conversion in homogeneous turbulence - Mathematical modeling, computational validations, and practical applications

The presently obtained closed-form analytical expressions, which predict the limiting rate of mean reactant conversion in homogeneous turbulent flows under the influence of a binary reaction, are derived via the single-point pdf method based on amplitude mapping closure. With this model, the maximum rate of the mean reactant's decay can be conveniently expressed in terms of definite integrals of the parabolic cylinder functions. The results obtained are shown to be in good agreement with data generated by direct numerical simulations.

Madnia, C. K.↗

NASA NTRS1993

Bifurcation analysis of a simple analytic model of self-propagating star formation

We investigate the structure and stability of rotationally symmetric nonhomogeneous time-independent solutions derived from a simple analytic model of self-propagating star formation. For this purpose we employ two methodologies: We use bifurcation theoretical methods to prove the existence of nonhomogeneous axisymmetric stationary solutions of an appropriate nonlinear evolution equation for the stellar density. We show that the nonhomogeneous solution branch bifurcates from the homogeneous one at a critical parameter value of the star formation rate. Further, the analytical theory allows us to show that the new solution set is stable in the weakly nonlinear regime near the bifurcation point. To follow the solution branch further, we use numerical methods. The numerical calculation shows the structure and stability of these solutions. We conclude that no periodic time-dependent solutions of this special model exist, and no further bifurcations can be found. The same results have been found in simulations of stochastic self-propagating star formation based on similar models. Therefore, our findings provide a natural explanation, why long-lived large-scale structure have not been found in those simulations.

Neukirch, Thomas↗

NASA NTRS1996

Transient Creep of a Composite Lower Crust: Constitutive Theory - 1

A composite model is proposed to describe the time-dependent response of the Earth's lower crust. The motivation for such it model is twofold: First, new observations of widespread postseismic deformation indicate that the deep continental crust responds viscoelastically, having both long-and short-term decay times. Second, by any number of observationally based rationales, the lower crust is compositionally and structurally heterogeneous over many length scales. For heterogeneities that have much smaller characteristic lengths than the minimum deformation wavelength of interest, the aggregate rheology can be described by composite media theory. For wavelengths of the order of the thickness of the lower crust (approx. = 25-40 km) and larger, composite theory may be applied to heterogeneities that are smaller than about several hundred meters, or equivalent to the vertical extent of a thick lower crustal mylonitic shear zone. The composite media theory developed here is constructed using both Eshelhy-Mori-Tanaka theory for aligned generalized spheroidal inclusions and a generalized self-consistent method. The inclusions and matrix are considered to be Maxwellian viscoelastic: a rheology that is consistent with past homogeneous models of postseismic stress relaxation. The composite theory presented here introduces a transient response to a suddenly imposed stress field which does not appear in homogeneous Maxwell models. Analytic expressions for the amplitude and duration of the transient and for the effective long-and short-term viscosities of the composite are given which describe the sensitivity to inclusion concentration (phi), to shape, and to ratio of inclusion-to-matrix viscosity (R).

Ivins, Erik R.↗

NASA NTRS2012

Adhesion of Mineral and Soot Aerosols can Strongly Affect their Scattering and Absorption Properties

We use the numerically exact superposition T-matrix method to compute the optical cross sections and the Stokes scattering matrix for polydisperse mineral aerosols (modeled as homogeneous spheres) covered with a large number of much smaller soot particles. These results are compared with the Lorenz-Mie results for a uniform external mixture of mineral and soot aerosols. We show that the effect of soot particles adhering to large mineral particles can be to change the extinction and scattering cross sections and the asymmetry parameter quite substantially. The effect on the phase function and degree of linear polarization can be equally significant.

Mishchenko, Michael I.↗