Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical analysis”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Reduced quantum dynamics with arbitrary bath spectral densities: Hierarchical equations of motion based on several different bath decomposition schemes

We investigated applications of the hierarchical equation of motion (HEOM) method to perform high order perturbation calculations of reduced quantum dynamics for a harmonic bath with arbitrary spectral densities. Three different schemes are used to decompose the bath spectral density into analytical forms that are suitable to the HEOM treatment: (1) The multiple Lorentzian mode model that can be obtained by numerically fitting the model spectral density. (2) The combined Debye and oscillatory Debye modes model that can be constructed by fitting the corresponding classical bath correlation function. (3) A new method that uses undamped harmonic oscillator modes explicitly in the HEOM formalism. Methods to extract system-bath correlations were investigated for the above bath decomposition schemes. We also show that HEOM in the undamped harmonic oscillator modes can give detailed information on the partial Wigner transform of the total density operator. Theoretical analysis and numerical simulations of the spin-Boson dynamics and the absorption line shape of molecular dimers show that the HEOM formalism for high order perturbations can serve as an important tool in studying the quantum dissipative dynamics in the intermediate coupling regime.

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

Analysis of Complex Valve and Feed Systems

A numerical framework for analysis of complex valve systems supports testing of propulsive systems by simulating key valve and control system components in the test loop. In particular, it is designed to enhance the analysis capability in terms of identifying system transients and quantifying the valve response to these transients. This system has analysis capability for simulating valve motion in complex systems operating in diverse flow regimes ranging from compressible gases to cryogenic liquids. A key feature is the hybrid, unstructured framework with sub-models for grid movement and phase change including cryogenic cavitations. The multi-element unstructured framework offers improved predictions of valve performance characteristics under steady conditions for structurally complex valves such as pressure regulator valve. Unsteady simulations of valve motion using this computational approach have been carried out for various valves in operation at Stennis Space Center such as the split-body valve and the 10-in. (approx.25.4-cm) LOX (liquid oxygen) valve and the 4-in. (approx.10 cm) Y-pattern valve (liquid nitrogen). Such simulations make use of variable grid topologies, thereby permitting solution accuracy and resolving important flow physics in the seat region of the moving valve. An advantage to this software includes possible reduction in testing costs incurred due to disruptions relating to unexpected flow transients or functioning of valve/flow control systems. Prediction of the flow anomalies leading to system vibrations, flow resonance, and valve stall can help in valve scheduling and significantly reduce the need for activation tests. This framework has been evaluated for its ability to predict performance metrics like flow coefficient for cavitating venturis and valve coefficient curves, and could be a valuable tool in predicting and understanding anomalous behavior of system components at rocket propulsion testing and design sites.

Ahuja, Vineet↗
NASA NTRS1989

SINFAC simulation of a high-power hybrid CPL

The systems improved numerical fluids analysis code (SINFAC) was used to model a high-power hybrid capillary-pumped loop. A comparison was made between the simulation results and actual experimental test data for the mechanical-pump-assist mode of operation. The hybrid capillary-pumped loop consisted of a capillary-pumped loop with a mechanical pump in the liquid return line with a bypass value so that the system could be operated either as a capillary-only loop or a series hybrid loop. Excellent agreement was obtained between numerical and experimental data.

Benner, Steve M.↗
DOE OSTI2020

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Bayesian analysis for kaon photoproduction

We have investigated contribution of the nucleon resonances in the kaon photoproduction process by using an established statistical decision making method, i.e. the Bayesian method. This method does not only evaluate the model over its entire parameter space, but also takes the prior information and experimental data into account. The result indicates that certain resonances have larger probabilities to contribute to the process.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Implications of theoretical analysis to explore the functional core dimension in one-rod core microstructured optical fibers

A new development in optical fiber technology that has attracted interest for multitude of novel and potential applications is an all-in-silica microstructured optical fibre (MOF) with one-rod core. We aim to estimate the functional core dimension for one-rod core triangular MOF (T-MOF) from its geometrical characteristics by using an alternative analytical field model. Moreover, to scrutinize the validity of the choice of effective core dimension (by adapting V-parameter, i.e., the normalized frequency) for T-MOF, we have explored the rudimentary propagation characteristics of the T-MOF in single mode guidance realm. We made comparisons with those which are based on the experimental and numerical studies. Relative errors are also reported.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗
DOE OSTI2021

A Butterfly-Accelerated Volume Integral Equation Solver for Broad Permittivity and Large-Scale Electromagnetic Analysis

In this work, a butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly (HOD-BF) scheme to construct the system matrix and obtain its approximate inverse, used as a preconditioner. Complexity analysis and numerical experiments validate the O(N log 2 N) construction cost of the HOD-BF-compressed system matrix and O(N log 1.5 N) inversion cost for the preconditioner, where N is the number of unknowns in the high-frequency EM scattering problem. For many practical scenarios, the proposed VIE solver requires less memory and computational time to construct the system matrix and obtain its approximate inverse compared to a H matrix-accelerated VIE solver. The accuracy and efficiency of the proposed solver have been demonstrated via its application to the EM analysis of large-scale canonical and real-world structures comprising of broad permittivity values and involving millions of unknowns.

42 ENGINEERING↗
NASA NTRS1976

Research on MLS airborne antenna

Numerical solutions for the radiation patterns of antennas mounted on aircraft are developed. The airborne antenna problems associated with the Microwave Landing System (MLS) are emphasized. Based on the requirements of the MLS, volumetric pattern solutions are essential. Previous attempts at solving for the volumetric patterns were found to be far too complex and very inefficient. However as a result of previous efforts, it is possible to combine the elevation and roll plane pattern solutions to give the complete volumetric pattern. This combination is described as well as the aircraft simulation models used in the analysis. A numerical technique is presented to aid in the simulation of the aircraft studied. Finally, a description of the input data used in the computer code is given.

Yu, C. L.↗
DOE OSTI2021

Lithium source for focused ion beam implantation and analysis

In this work, we present a new Li source for focused ion beam applications. Based on an AuSi eutectic alloy, Li is added as an impurity to minimize effects from degradation when exposed to air. We show the source is stable over the course of an hour and spot sizes ≲10 nm can be achieved. The Li beam can achieve hundreds of nanometer ranges in semiconductors with minimal damage being generated along the path length. The source performance is evaluated through a high-resolution ion beam induced charge collection experiment on an Si-based detector. Further application of the source for ion beam analysis is numerically explored; the example investigated is based on probing a semiconductor heterostructure through a Rutherford backscattering experiment, where the Li beam can reveal information that is inaccessible with either low energy or high energy He projectiles used as probes.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Effect of sputtered lanthanum hexaboride film thickness on field emission from metallic knife edge cathodes

We report experiments and analysis of field emission from metallic knife-edge cathodes, which are sputter-coated with thin films of lanthanum hexaboride (LaB{sub 6}), a low-work function material. The emission current is found to depend sensitively on the thickness of the LaB{sub 6} layer. We find that films thinner than 10 nm greatly enhance the emitted current. However, cathodes coated with a thicker layer of LaB{sub 6} are observed to emit less current than the uncoated metallic cathode. This result is unexpected due to the higher work function of the bare metal cathode. We show, based on numerical calculation of the electrostatic potential throughout the structure, that the external (LaB{sub 6}/vacuum) barrier is reduced with respect to uncoated samples for both thin and thick coatings. However, this behavior is not exhibited at the internal (metal/LaB{sub 6}) barrier. In thinly coated samples, electrons tunnel efficiently through both the internal and external barrier, resulting in current enhancement with respect to the uncoated case. In contrast, the thick internal barrier in thickly coated samples suppresses current below the value for uncoated samples in spite of the lowered external barrier. We argue that this coating thickness variation stems from a relatively low (no higher than 10{sup 18} cm{sup -3}) free carrier density in the sputtered polycrystalline LaB{sub 6}.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Impacts of environmental exposure on thermal and mycological characteristics of insulation wools

Highlights: • Insulation materials (IM) suffer expositions to weather on energy efficient buildings. • Induced exposure to weather of IM reveals differences in ATR-FTIR spectra. • Mycological analysis shows a significant effect of exposure condition on IM. • Mold species detected may induce human health risk and degradation of IM. • Thermal conductivity of IM according humidity of air may impact energy performance. - Abstract: Fieldwork was conducted on 138 energy efficient buildings, to check if building site practices induced exposure to weather of thermal insulation. In nearly 80% of cases, insulation was exposed to weather due to insufficient protection. Therefore insulating wools were exposed, divided into three parts (control, exposure under shelter, exposure to weather). After the exposure period these samples were submitted to thermal conductivity measurements and infrared spectroscopy. A numerical study was carried out to evaluate the influence of weather conditions on thermal performance. Identifications of molds were realized by direct cultures of each sample. Mycological analyses provided the most significant results, showing a significant effect of exposure condition, with possible consequences in human health (toxinogenic, allergenic or pathogenic species) and material properties (cellulolytic species). Infrared spectrometry also showed some changes in sensitivity to water. These results show the interest to study in more detail the variations in the sensitivity, to moisture and mold, of insulation material throughout its ageing on site, and the interactions with the practices of the building process.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Electromagnetic energy flow in confined regions of evanescent waves: wavelength-scale analysis by the method of single expression

An alternative expression for the Poynting vector for monitoring electromagnetic energy flow in confined regions of evanescent waves is presented. This expression is an intrinsic integrable function of the boundary problem solution in the method of single expression (MSE). The peculiarity of the MSE is the fact that it does not exploit the superposition principle by representing solutions of the Helmholtz equation as a sum of counter propagating waves in the medium, but rather considers resulting wave with certain amplitude and phase by solving special set of the MSE equations with backward calculation approach. This permits to have a universal expression for the Poynting vector appropriate both for media of a positive product of permittivity and permeability and of a negative product of these values that is relevant to the region of evanescent waves. For the boundary problem solution carried out numerically by the MSE it is possible to observe distributions of not only electric and magnetic field amplitudes, but also the Poynting vector in any confined media including regions of evanescent waves. The presented alternative expression for the Poynting vector in the MSE is in complete agreement with the traditional representation of the Poynting vector for the medium where it is defined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Conjugate heat and mass transfer in the lattice Boltzmann equation method

An interface treatment for conjugate heat and mass transfer in the lattice Boltzmann equation method is proposed based on our previously proposed second-order accurate Dirichlet and Neumann boundary schemes. The continuity of temperature (concentration) and its flux at the interface for heat (mass) transfer is intrinsically satisfied without iterative computations, and the interfacial temperature (concentration) and their fluxes are conveniently obtained from the microscopic distribution functions without finite-difference calculations. The present treatment takes into account the local geometry of the interface so that it can be directly applied to curved interface problems such as conjugate heat and mass transfer in porous media. For straight interfaces or curved interfaces with no tangential gradient, the coupling between the interfacial fluxes along the discrete lattice velocity directions is eliminated and thus the proposed interface schemes can be greatly simplified. Several numerical tests are conducted to verify the applicability and accuracy of the proposed conjugate interface treatment, including (i) steady convection-diffusion in a channel containing two different fluids, (ii) unsteady convection-diffusion in the channel, (iii) steady heat conduction inside a circular domain with two different solid materials, and (iv) unsteady mass transfer from a spherical droplet in an extensional creeping flow. The accuracy and order of convergence of the simulated interior temperature (concentration) field, the interfacial temperature (concentration), and heat (mass) flux are examined in detail and compared with those obtained from the "half-lattice division" treatment in the literature. The present analysis and numerical results show that the half-lattice division scheme is second-order accurate only when the interface is fixed at the center of the lattice links, while the present treatment preserves second-order accuracy for arbitrary link fractions. For curved interfaces, the present treatment yields second-order accurate interior and interfacial temperatures (concentrations) and first-order accurate interfacial heat (mass) flux. An increase of order of convergence by one degree is obtained for each of these three quantities compared with the half-lattice division scheme. The surface-averaged Sherwood numbers computed in test (iv) agree well with published results.

Li, LK↗
NASA NTRS1996

Variational Methods in Sensitivity Analysis and Optimization for Aerodynamic Applications

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equation is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

Ibrahim, A. H.↗
NASA NTRS2007

Core Physics and Kinetics Calculations for the Fissioning Plasma Core Reactor

Highly efficient, compact nuclear reactors would provide high specific impulse spacecraft propulsion. This analysis and numerical simulation effort has focused on the technical feasibility issues related to the nuclear design characteristics of a novel reactor design. The Fissioning Plasma Core Reactor (FPCR) is a shockwave-driven gaseous-core nuclear reactor, which uses Magneto Hydrodynamic effects to generate electric power to be used for propulsion. The nuclear design of the system depends on two major calculations: core physics calculations and kinetics calculations. Presently, core physics calculations have concentrated on the use of the MCNP4C code. However, initial results from other codes such as COMBINE/VENTURE and SCALE4a. are also shown. Several significant modifications were made to the ISR-developed QCALC1 kinetics analysis code. These modifications include testing the state of the core materials, an improvement to the calculation of the material properties of the core, the addition of an adiabatic core temperature model and improvement of the first order reactivity correction model. The accuracy of these modifications has been verified, and the accuracy of the point-core kinetics model used by the QCALC1 code has also been validated. Previously calculated kinetics results for the FPCR were described in the ISR report, "QCALC1: A code for FPCR Kinetics Model Feasibility Analysis" dated June 1, 2002.

Butler, C.↗
NASA NTRS2014

Instantaneous Wavenumber Estimation for Damage Quantification in Layered Plate Structures

This paper illustrates the application of instantaneous and local wavenumber damage quantification techniques for high frequency guided wave interrogation. The proposed methodologies can be considered as first steps towards a hybrid structural health monitoring/ nondestructive evaluation (SHM/NDE) approach for damage assessment in composites. The challenges and opportunities related to the considered type of interrogation and signal processing are explored through the analysis of numerical data obtained via EFIT simulations of damage in CRFP plates. Realistic damage configurations are modeled from x-ray CT scan data of plates subjected to actual impacts, in order to accurately predict wave-damage interactions in terms of scattering and mode conversions. Simulation data is utilized to enhance the information provided by instantaneous and local wavenumbers and mitigate the complexity related to the multi-modal content of the plate response. Signal processing strategies considered for this purpose include modal decoupling through filtering in the frequency/wavenumber domain, the combination of displacement components, and the exploitation of polarization information for the various modes as evaluated through the dispersion analysis of the considered laminate lay-up sequence. The results presented assess the effectiveness of the proposed wavefield processing techniques as a hybrid SHM/NDE technique for damage detection and quantification in composite, plate-like structures.

Mesnil, Olivier↗
DOE OSTI2014

Fourier analysis of conductive heat transfer for glazed roofing materials

For low-rise buildings, roof is the most exposed surface to solar radiation. The main mode of heat transfer from outdoor via the roof is conduction. The rate of heat transfer and the thermal impact is dependent on the thermophysical properties of roofing materials. Thus, it is important to analyze the heat distribution for the various types of roofing materials. The objectives of this paper are to obtain the Fourier series for the conductive heat transfer for two types of glazed roofing materials, namely polycarbonate and polyfilled, and also to determine the relationship between the ambient temperature and the conductive heat transfer for these materials. Ambient and surface temperature data were collected from an empirical field investigation in the campus of Universiti Teknologi MARA Shah Alam. The roofing materials were installed on free-standing structures in natural ventilation. Since the temperature data are generally periodic, Fourier series and numerical harmonic analysis are applied. Based on the 24-point harmonic analysis, the eleventh order harmonics is found to generate an adequate Fourier series expansion for both glazed roofing materials. In addition, there exists a linear relationship between the ambient temperature and the conductive heat transfer for both glazed roofing materials. Based on the gradient of the graphs, lower heat transfer is indicated through polyfilled. Thus polyfilled would have a lower thermal impact compared to polycarbonate.

42 ENGINEERING↗