Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Structural Concepts Study of Non-circular Fuselage Configurations

A preliminary study of structural concepts for noncircular fuselage configurations is presented. For an unconventional flying-wing type aircraft, in which the fuselage is inside the wing, multiple fuselage bays with non-circular sections need to be considered. In a conventional circular fuselage section, internal pressure is carried efficiently by a thin skin via hoop tension. If the section is non-circular, internal pressure loads also induce large bending stresses. The structure must also withstand additional bending and compression loads from aerodynamic and gravitational forces. Flat and vaulted shell structural configurations for such an unconventional, non-circular pressurized fuselage of a large flying-wing were studied. A deep honeycomb sandwich-shell and a ribbed double-wall shell construction were considered. Combinations of these structural concepts were analyzed using both analytical and simple finite element models of isolated sections for a comparative conceptual study. Weight, stress, and deflection results were compared to identify a suitable configuration for detailed analyses. The flat sandwich-shell concept was found preferable to the vaulted shell concept due to its superior buckling stiffness. Vaulted double-skin ribbed shell configurations were found to be superior due to their weight savings, load diffusion, and fail-safe features. The vaulted double-skin ribbed shell structure concept was also analyzed for an integrated wing-fuselage finite element model. Additional problem areas such as wing-fuselage junction and pressure-bearing spar were identified.

Mukhopadhyay, Vivel↗
DOE OSTI2016

Position-Dependency of Fuel Assembly Homogenization in a Pressurized Water Reactor

2-step nuclear reactor core analysis has long been established as the conventional calculational procedure for the pressurized water reactor (PWR) due to its efficiency. This is because the procedure significantly avoids the potentially huge computational burden arising from an explicit geometrical modeling and consideration of over ∼ 500 energy groups. It instead approaches the problem with spatial homogenization and group condensation: a 2-D lattice transport calculation to homogenize the fuel assembly cross-sections and a 3-D whole core nodal diffusion. One notes that the few-group constant generation using deterministic method could be complicated due to its limited flexibilities in geometrical modeling and inherent drawbacks of the multi-group approximation. Furthermore, accuracy of the deterministic few-group constant generation can especially be compromised during the resonance treatment calculation. As such, an alternative lattice calculation using the Monte Carlo (MC) approach, in place of the deterministic fuel assembly homogenization (few-group constant generation), has previously been proposed. During the assembly homogenization, directions and magnitudes of the current at boundaries are unknowns. Therefore, zero net current boundary conditions (infinite lattice) are instead used as the assembly node. The assumption is quite reasonable considering heterogeneities inside the node are more dominant than heterogeneities due to nodal location in the core. In fact, this approximation has successfully been applied to many LWR core analyses to reasonable accuracy. It should nonetheless be noted that the net current zero boundary assumption effectively ignores position-dependency of the fuel assembly homogenization in the general 2-step method. As such, there are a number of well-known researches on the extended assembly homogenization techniques, which introduce correction factors to compensate effects of the assembly environment, such as Generalized Equivalence Theory (GET), Super-homogenization method (SPH), and color set calculation method. Even though the assembly environment effect has long been acknowledged, its effect has not properly been investigated due to the limited computing power in the past and due to the reasonable accuracy of the extended assembly homogenization techniques. Nonetheless, quantification of position-dependency of the assembly homogenization in the core could still potentially offer additional value in terms of calculation accuracy to the 2-step reactor analysis. It is upon this motivation that the paper is prepared. In this research, position-dependent few-group constants were calculated and compared against those of infinite lattices. All calculations were completed using the MC Serpent2 code with ENDF/B7.0 nuclear data library. In this study, position-dependency of the fuel assembly homogenized cross sections in a PWR core has been investigated through comparison of the conventional 2-step reactor analysis with the reference high fidelity Monte Carlo simulation. It was found that the homogenized two-group constants are clearly position-dependent as the discrepancies between infinite lattice and whole-core calculations are quite significant especially in a rodded assembly. As such, the corresponding correction required due to the assembly environment effect is not negligible compared to the correction required due to assembly heterogeneity itself. It is therefore quite clear that position-dependency of the group cross sections should be considered in the 2-step procedure for a more accurate PWR core analysis. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Hyperbolic Reformulation Approach to Enable Efficient Simulation of Groundwater Flow and Reactive Transport

We apply Cattaneo's relaxation approach to the one-dimensional coupled Boussinesq groundwater flow and advection-diffusion-reaction equations, commonly used in engineering applications to simulate contaminant transport in the subsurface. The diffusion-type governing equations are reformulated as a hyperbolic system, augmented by an equation that can be interpreted as a momentum balance. The hyperbolization enables an efficient unified computation of the primary variable and its gradients, for example piezometric head and unit discharge in the Boussinesq equation. An augmented Roe scheme is used to solve the hyperbolic system. The hyperbolized system of equations is studied in a set of steady state and transient test cases with idealized geometry. These test cases confirm the equivalence of the hyperbolic system to its original formulation. The larger time step size of the hyperbolic equation is verified theoretically by means of a stability analysis and numerically in the test cases. Finally, a reach-scale application of flow and transport across a river meander is considered. This application case shows that the performance of the hyperbolic relaxation approach holds for more realistic groundwater flow and transport problems, relevant to water resources management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

THE TELEGRAPH APPROXIMATION FOR FOCUSED COSMIC-RAY TRANSPORT IN THE PRESENCE OF BOUNDARIES

Diffusive cosmic-ray transport in nonuniform large-scale magnetic fields in the presence of boundaries is considered. Reflecting and absorbing boundary conditions are derived for a modified telegraph equation with a convective term. Analytical and numerical solutions of illustrative boundary problems are presented. The applicability and accuracy of the telegraph approximation for focused cosmic-ray transport in the presence of boundaries are discussed, and potential applications to modeling cosmic-ray transport are noted.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1972

Application of finite-element method in the computation of temperature with emphasis on radiative exchanges.

Analyses pertaining to the solution of heat transfer problems in combined modes based on the finite-element method are presented. Two elements - a triangular element employing two spatial variables and a multi-faceted bar element employing one spatial variable - with nonlinear radiation on the boundaries are detailed. The radiative effects considered on the diffuse-gray surface elements include: (1) directional radiant fluxes from distant sources, (2) radiative exchanges including surfaces of prescribed temperatures, and (3) radiative exchanges including elements whose temperatures are not known a priori. The nonlinear part of the boundary condition is treated in two different ways: (1) a consistent linearization method, and (2) a direct energy distribution method. Applications of these elements together with solution algorithms to three sample problems exhibit solution capability and obtainable accuracy.

Lee, H.-P.↗
DOE OSTI2016

A variational treatment of material configurations with application to interface motion and microstructural evolution

We present a unified variational treatment of evolving configurations in crystalline solids with microstructure. The crux of our treatment lies in the introduction of a vector configurational field. This field lies in the material, or configurational, manifold, in contrast with the traditional displacement field, which we regard as lying in the spatial manifold. We identify two distinct cases which describe (a) problems in which the configurational field's evolution is localized to a mathematically sharp interface, and (b) those in which the configurational field's evolution can extend throughout the volume. The first case is suitable for describing incoherent phase interfaces in polycrystalline solids, and the latter is useful for describing smooth changes in crystal structure and naturally incorporates coherent (diffuse) phase interfaces. These descriptions also lead to parameterizations of the free energies for the two cases, from which variational treatments can be developed and equilibrium conditions obtained. For sharp interfaces that are out-of-equilibrium, the second law of thermodynamics furnishes restrictions on the kinetic law for the interface velocity. The class of problems in which the material undergoes configurational changes between distinct, stable crystal structures are characterized by free energy density functions that are non-convex with respect to configurational strain. For physically meaningful solutions and mathematical well-posedness, it becomes necessary to incorporate interfacial energy. This we have done by introducing a configurational strain gradient dependence in the free energy density function following ideas laid out by Toupin (Arch. Rat. Mech. Anal., 11, 1962, 385-414). The variational treatment leads to a system of partial differential equations governing the configuration that is coupled with the traditional equations of nonlinear elasticity. The coupled system of equations governs the configurational change in crystal structure, and elastic deformation driven by elastic, Eshelbian, and configurational stresses. As a result, numerical examples are presented to demonstrate interface motion as well as evolving microstructures of crystal structures.

36 MATERIALS SCIENCE↗
NASA NTRS1994

New description of charged particle propagation in random magnetic fields

When charged particles spiral along a large constant magnetic field, their trajectories are scattered by random components that are superposed on the guiding field. In the simplest analysis of this situation, scattering causes the particles to diffuse parallel to the guiding field. At the next level of approximation, moving pulses that correspond to a coherent mode of propagation are present, but they are represented by delta-functions whose infinitely narrow width makes no sense physically and is inconsistent with the finite duration of coherent pulses observed in solar energetic particle events. To derive a more realistic description, the transport problem is formulated in terms of 4 x 4 matrices, which derive from a representation of the particle distribution function in terms of eigenfunctions of the scattering operator, and which lead to useful approximations that give explicit predictions of the detailed evolution not only of the coherent pulses, but also of the diffusive wake. More specifically, the new description embodies a simple convolution of a narrow Gaussian with the solutions above that involve delta-functions, but with a slightly reduced coherent velocity. The validity of these approximations, which can easily be calculated on a desktop computer, has been exhaustively confirmed by comparison with results of Monte Carlo simulations which kept track of 50 million particles and which were carried out on the Maspar computer at Goddard Space Flight Center.

Earl, James A.↗
DOE OSTI2019

A multi-dimensional, moment-accelerated deterministic particle method for time-dependent, multi-frequency thermal radiative transfer problems

Thermal Radiative Transfer (TRT) is the dominant energy transfer mechanism in high-energy density physics with applications in inertial confinement fusion and astrophysics. The stiff interactions between the material and radiation fields make TRT problems challenging to model. In this study, we propose a multi-dimensional extension of the deterministic particle (DP) method. The DP method combines aspects from both particle and deterministic methods. If the emission source is known a priori, and no physical scattering is present, the intensity of a particle can be integrated analytically. This introduces no statistical noise compared to Monte-Carlo methods, while maintaining the flexibility of particle methods. The method is closely related to the popular method of long characteristics. The combination of the DP-method with a discretely-consistent, nonlinear, gray low-order system enables an efficient solution algorithm for multi-frequency TRT problems. Here, we demonstrate with numerical examples that the use of a linear-source approximation based on spatial moments improves the behavior of our method in the thick diffusion limit significantly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Blister resistant targets for nuclear reaction experiments with α-particle beams

Solid targets for nuclear measurements that use α-particle beams commonly experience a form of degradation known as blistering. The effect can prevent the use of solid targets for high intensity α-particle experiments, often necessitating complex gas target systems. To combat this problem, three different blister resistant target backings were designed for use in direct reaction measurements with high intensity α-particle beams. The blister resistant target designs utilize gas diffusive properties of fused silica, sintered metal, and porous evaporated metal. Each target was implanted with 22 Ne ions and bombarded with α-particle beam to test blister resistance. Targets were characterized and monitored using the 22 Ne(p, γ) 23 Na reaction to determine the degradation of implanted material, and compare them to typical implanted noble gas targets. We find that all targets studied exhibit resistance to blistering, with the porous evaporated metal targets displaying the least amount of target material degradation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Interaction of Planar Brushes with End-Functionalized Carboxylic or Amine Groups: Water Structure, Hydrogen Bonds, and Electrostatic Correlations

In this work we present a detailed analysis of two parallel plates grafted with polyethylene glycol (PEG) terminated with carboxylic or amine groups, and discuss the cases of opposite and same charge in the presence of electrolytes. We analyze the screening of charged groups by electrolytes, the role of hydrogen bonds in enhancing attraction or repulsion, hydration forces, competitive binding with counterions, and the hydrogen bond network, crucial to stabilizing the system. Our simulations reveal that, for large plate separation, the system is well described as a “Stern layer”, characterized by strong counterion binding and a subsequent weak diffuse layer. For short plate separation, the system is characterized by strong correlations among end groups and counterions. We also discuss the relevance of charged carboxylic/amine groups in different problems with a particular emphasis on nanoparticle assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1987

Analysis of non-premixed turbulent reacting flows

Studies of chemical reactions occurring in turbulent flows are important in the understanding of combustion and other applications. Current numerical methods are limited in their applications due to the numerical resolution required to completely capture all length scales, but, despite the fact that realistic combustion cannot be solved completely, numerical simulations can be used to give insight into the interaction between the processes of turbulence and chemical reaction. The objective was to investigate the effects of turbulent motion on the effects of chemical reaction to gain some insight on the interaction of turbulence, molecular diffusion, and chemical reaction to support modeling efforts. A direct turbulence simulation spectral code was modified to include the effects of chemical reaction and applied to an initial value problem of chemical reaction between non-premixed species. The influence of hydrodynamics on the instantaneous structure of the reaction was investigated. The local scalar dissipation rates and the local reaction rates were examined to determine the influence of vorticity or rate of strain on the reaction and the structure of the scalar field.

Leonard, Andy D.↗
NASA NTRS1988

Dynamics of two fluids under periodic acceleration

The evolution of the interface between two fluids confined in a rectangular cavity is investigated numerically to predict its transient behavior. These computations address mixing characteristics of fluids under microgravity conditions, particularly g-jitter conditions resulting from aircraft vibration or crew motion, and have applications in solution crystal growth. The two-dimensional formulation employs the Boussinesq approximation and treats the mixing of two fluids as an initial value problem with a prescribed concentration field. For fluid mixing of practical applications inside a cavity, it is shown that nonlinear convective transport can dominate over both viscous and molecular diffusion. However, viscous diffusion can become important for low Reynolds number or certain cavity sizes. A stacking phenomenon is shown to occur for aspect ratios (Ar) approximately between 5 to 10. For aspect ratios in the neighborhood of 0.1 to 0.2, inner cells evolve at the interface with increasing time. For a square cavity of Ar = 1, chaotic mixing of the fluid appears to occur for a Stokes-Reynolds number greater than about 5.

Duval, Walter M. B.↗
NASA NTRS1991

The Implications of Experimentally Controlled Gravitational Accelerations for Combustion Science

An overview of basic combustion problems which have been investigated under the condition of reduced gravity is presented to identify promising research directions. Attention is given to the broad categories of gas-jet diffusion flames, droplet combustions, particle clouds, flame spreading over liquid pools, smoldering, and flame spreading over solid fuels. Fire safety in spacecraft is the primary application that is addressed by the studies of combustion under microgravity. The need for more complete testing of the issues discussed in orbiting spacecraft is identified in the light of limited earth-based testing. Attention is also directed toward the need for advanced diagnostic methods for in-flight and other combustion investigations.

Sacksteder, Kurt R.↗
NASA NTRS1993

Influence of test configuration on the combustion characteristics of polymers as ignition sources

The experimental evaluation of polymers as ignition sources for metals was accomplished at the NASA White Sands Test Facility (WSTF) using a standard promoted combustion test. These tests involve the transient burning of materials in high-pressure oxygen environments. They have provided data from which design decisions can be made; data include video recordings of ignition and non-ignition for specific combinations of metals and polymers. Other tests provide the measured compositions of combustion products for polymers at select burn times and an empirical basis for estimating burn rates. With the current test configuration, the detailed analysis of test results requires modeling a three-dimensional, transient convection process involving fluid motion, thermal conduction and convection, the diffusion of chemical species, and the erosion of sample surface. At the high pressure extremes, it even requires the analysis of turbulent, transient convection where the physics of the problem are not well known and the computation requirements are not practical at this time. An alternative test configuration that can be analyzed with a relatively-simple convection model was developed during the summer period. The principal change constitutes replacing a large-diameter polymer disk at the end of the metal test rod with coaxial polymer cylinders that have a diameter nearer to that of the metal rod. The experimental objective is to assess the importance of test geometries on the promotion of metal ignition by testing with different lengths of the polymer and, with an extended effort, to analyze the surface combustion in the redesigned promoted combustion tests through analytical modeling of the process. The analysis shall use the results of cone-calorimeter tests of the polymer material to model primary chemical reactions and, with proper design of the promoted combustion test, modeling of the convection process could be conveniently limited to a quasi-steady boundary layer analysis where the economical solution of parabolic equations is involved. The products for the summer period are: (1) a conceptual-level redesign of the test apparatus, and (2) the development and use of an approximate integral boundary layer analysis to demonstrate the influence of geometry changes prior to testing. A computer code STAN5, an accurate numerical boundary layer model whose earlier versions were developed for the NASA Lewis Research Center by the Fellow, also was installed and validated on the WSTF and New Mexico State University computer systems as a starting point in the development of a more detailed fluid mechanics and combustion model.

Julien, Howard L.↗
NASA NTRS1999

Advanced Electromagnetic Methods for Aerospace Vehicles

The imminent destructive threats of Lightning on helicopters and other airborne systems has always been a topic of great interest to this research grant. Previously, the lightning induced currents on the surface of the fuselage and its interior were predicted using the finite-difference time-domain (FDTD) method as well as the NEC code. The limitations of both methods, as applied to lightning, were identified and extensively discussed in the last meeting. After a thorough investigation of the capabilities of the FDTD, it was decided to incorporate into the numerical method a subcell model to accurately represent current diffusion through conducting materials of high conductivity and finite thickness. Because of the complexity of the model, its validity will be first tested for a one-dimensional FDTD problem. Although results are not available yet, the theory and formulation of the subcell model are presented and discussed here to a certain degree. Besides lightning induced currents in the interior of an aircraft, penetration of electromagnetic fields through apertures (e.g., windows and cracks) could also be devastating for the navigation equipment, electronics, and communications systems in general. The main focus of this study is understanding and quantifying field penetration through apertures. The simulation is done using the FDTD method and the predictions are compared with measurements and moment method solutions obtained from the NASA Langley Research Center. Cavity-backed slot (CBS) antennas or slot antennas in general have many applications in aircraft-satellite type of communications. These can be flushmounted on the surface of the fuselage and, therefore, they retain the aerodynamic shape of the aircraft. In the past, input impedance and radiation patterns of CBS antennas were computed using a hybrid FEM/MoM code. The analysis is now extended to coupling between two identical slot antennas mounted on the same structure. The predictions are performed using both the hybrid FEM/MoM and the FDTD NEWS code. The results are compared with each other as well as with measurements performed in the ElectroMagnetic Anechoic Chamber (EMAC) of ASU. In addition to self and mutual impedances versus frequency, the comparisons include mutual coupling S(sub 12) as a function of distance for various slot orientations.

Balanis, Constantine A.↗
DOE OSTI2021

Towards robust autonomous impedance spectroscopy analysis: A calibrated hierarchical Bayesian approach for electrochemical impedance spectroscopy (EIS) inversion

Distribution-based analyses, such as the distribution of relaxation times (DRT) and the distribution of diffusion times (DDT), present model-free alternatives to equivalent circuit modeling for analysis of electrochemical impedance spectroscopy (EIS) data. However, reconstructing such distributions from noisy impedance data is an ill-posed problem that must be solved with specialized inversion algorithms, requiring careful control and tuning. Furthermore, most inversion algorithms developed to date can only solve problems of limited complexity. Herein, we present a new hierarchical Bayesian method for EIS inversion, leveraging efficient algorithms for optimization and Hamiltonian Monte Carlo (HMC) sampling to solve models of arbitrary complexity. We overcome the challenge of ad-hoc parameter tuning by encoding intrinsic characteristics of the DRT and DDT into flexible prior distributions and “pre-calibrating” the model to simulated data. This approach is versatile, highly robust to noise, and provides quantitative estimates of both the error structure of the data and the uncertainty in the recovered distributions. The model is validated with simulated data to demonstrate accurate recovery of the DRT and the DDT. The method also shows promise for simultaneous recovery of multiple distributions, raising the intriguing possibility of semi-autonomous EIS analysis and ad-hoc model construction. Finally, the practical utility of the method is illustrated with experimental data. Throughout, we draw comparisons to several recently published EIS inversion methodologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Physics-Informed Generative Adversarial Networks for Stochastic Differential Equations

We developed a new class of physics-informed generative adversarial networks (PI-GANs) to solve forward, inverse, and mixed stochastic problems in a unified manner based on a limited number of scattered measurements. Unlike standard GANs relying solely on data for training, here we encoded into the architecture of GANs the governing physical laws in the form of stochastic differential equations (SDEs) using automatic differentiation. In particular, we applied Wasserstein GANs with gradient penalty (WGAN-GP) for its enhanced stability compared to vanilla GANs. We first tested WGAN-GP in approximating Gaussian processes of different correlation lengths based on data realizations collected from simultaneous reads at sparsely placed sensors. We obtained good approximation of the generated stochastic processes to the target ones even if there is a mismatch between the input noise dimensionality and the effective dimensionality of the target stochastic processes. We also studied the overfitting issue for both the discriminator and the generator, and we found that overfitting occurs also in the generator in addition to the discriminator as previously reported. Subsequently, we considered the solution of elliptic SDEs requiring approximations of three stochastic processes, namely the solution, the forcing, and the diffusion coefficient. Here again, we assumed data collected from simultaneous reads at a limited number of sensors for the multiple stochastic processes. Furthermore, three generators were used for the PI-GANs: two of them were feed forward deep neural networks (DNNs), while the other one was the neural network induced by the SDE. For the case where we have one group of data, we employed one feed forward DNN as the discriminator, while for the case of multiple groups of data we employed multiple discriminators in PI-GANs. We solved forward, inverse, and mixed problems without changing the framework of PI-GANs, obtaining both the means and the standard deviations of the stochastic solution and the diffusion coefficient in good agreement with benchmarks. In this work, we have demonstrated the effectiveness of PI-GANs in solving SDEs for about 120 dimensions. In principle, PI-GANs could tackle very high dimensional problems given more sensor data with low-polynomial growth in computational cost.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Analysis of Integrated Fuselage Composite Seat for Small Aircraft in Crashworthiness Applications

The research paper is addressed to develop and validate an aircraft seat cushion to be used for the passenger seat directly integrated in the fuselage. For very light aircraft, these systems result to be a good choice in terms of weight reduction and space optimization, but can represent a challenging problem for crashworthiness design and a key energy-absorbing function for human protection from injuries so to improve the level of the safety. The first issue was dedicated to investigate the mechanical properties of four commercial foams, diffused in different industrial applications, which present nonlinear and time-dependent viscoelastic behavior, and finally to determine whether typical seat cushion materials are rate sensitive or not. Then, numerical simulations were validated with measurements to simulate the mechanical properties of foams and to evaluate the contact pressures, stress distribution, and the dependence of stress on the strain rate in order to provide adequate levels of both comfort and energy absorption.

36 MATERIALS SCIENCE↗