Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Numerical calculation of the transonic flow past a swept wing

A numerical method is presented for analyzing the transonic potential flow past a lifting, swept wing. A finite difference approximation to the full potential equation is solved in a coordinate system which is nearly conformally mapped from the physical space in planes parallel to the symmetry plane, and reduces the wing surface to a portion of one boundary of the computational grid. A coordinate invariant, rotated difference scheme is used, and the difference equations are solved by relaxation. The method is capable of treating wings of arbitrary planform and dihedral, although approximations in treating the tips and vortex sheet make its accuracy suspect for wings of small aspect ratio. Comparisons of calculated results with experimental data are shown for examples of both conventional and supercritical transport wings. Agreement is good for both types, but it was found necessary to account for the displacement effect of the boundary layer for the supercritical wing, presumably because of its greater sensitivity to changes in effective geometry.

Jameson, A.↗
DOE OSTI2019

Novel refractive index sensor based on fiber bragg grating in nano-bore optical fiber

In this paper, we propose a novel refractive index sensor with large dynamic range based on fiber Bragg grating in the nano-bore optical fiber. The spectral responses of grating on the filling index in the bore with different bore diameters and index ranges are numerically studied by the finite element method and coupled mode theory based transmission matrix method. It is found that the Bragg wavelength of grating is sensitive to the filling index in the bore. The influences of bore diameter on the wavelength sensitivity, resolution and dynamic range of the index sensor are analyzed. There exists a trade off between the linear dynamic range and sensitivity of index when the bore diameter varies. The optimized wavelength sensitivity and resolution of filling index are 53.6923 nm/index unit and 9.3123 × 10{sup −6} index under a linear index range of 1.0–1.4 using 2 µm bore diameter. The optimized linear dynamic range of index is 1.0–1.48 with wavelength sensitivity of 13.3056 nm/index unit and resolution of 3.7578 × 10{sup −5} index using 1 µm bore diameter. Our proposed index sensor has the technical advantages of simple all-fiber structure, large index range, and potential capability to synchronously detect index, temperature and strain.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Environmental influence on cracking and debonding of electrically conductive adhesives

Electrically conductive adhesives (ECAs) are starting to replace metallic solders in recent designs of photovoltaic (PV) modules. This transition represents a significant material change, and a proper understanding of the durability and reliability of the new interconnect needs to be established. This paper presents our continued work on developing a degradation model for ECA interconnects in PV modules. Here, we characterize the fracture mechanics properties of an epoxy-based ECA, for both critical and subcritical, mode I and mode II loading conditions. Emphasis is on the influence of different environmental conditions such as temperature and humidity. We use the Finite Element Method to account for residual stresses, induced by temperature changes and moisture absorption, and correct the apparent fracture toughness. We found that high moisture levels not only can weaken the fracture resistance of the ECA interconnect, but can also promote subcritical debonding at significantly lower driving forces than in dry environments.

42 ENGINEERING↗
DOE OSTI2017

Metal nano-grids for transparent conduction in solar cells

A general procedure for predicting metal grid performance in solar cells was developed. Unlike transparent conducting oxides (TCOs) or other homogeneous films, metal grids induce more resistance in the neighbor layer. The resulting balance of transmittance, neighbor and grid resistance was explored in light of cheap lithography advances that have enabled metal nano-grid (MNG) fabrication. The patterned MNGs have junction resistances and degradation rates that are more favorable than solution-synthesized metal nanowires. Neighbor series resistance was simulated by the finite element method, although a simpler analytical model was sufficient in most cases. Finite-difference frequency-domain transmittance simulations were performed for MNGs with minimum wire width (w) of 50 nm, but deviations from aperture transmittance were small in magnitude. Depending on the process, MNGs can exhibit increased series resistance as w is decreased. However, numerous experimental reports have already achieved transmittance-MNG sheet resistance trade-offs comparable to TCOs. The transmittance, neighbor and MNG series resistances were used to parameterize a grid fill factor for a solar cell. In conclusion, this new figure of merit was used to demonstrate that although MNGs have only been employed in low efficiency solar cells, substantial gains in performance are predicted for decreased w in all high efficiency absorber technologies.

14 SOLAR ENERGY↗
DOE OSTI2016

Compact spherical neutron polarimeter using high-T{sub c} YBCO films

We describe a simple, compact device for spherical neutron polarimetry measurements at small neutron scattering angles. The device consists of a sample chamber with very low (<0.01 G) magnetic field flanked by regions within which the neutron polarization can be manipulated in a controlled manner. This allows any selected initial and final polarization direction of the neutrons to be obtained. We have constructed a prototype device using high-T{sub c} superconducting films and mu-metal to isolate regions with different magnetic fields and tested device performance in transmission geometry. Finite-element methods were used to simulate the device’s field profile and these have been verified by experiment using a small solenoid as a test sample. Measurements are reported using both monochromatic and polychromatic neutron sources. The results show that the device is capable of extracting sample information and distinguishing small angular variations of the sample magnetic field. As a more realistic test, we present results on the characterization of a 10 μm thick Permalloy film in zero magnetic field, as well as its response to an external magnetic field.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Inverse Design of Plasma Metamaterial Devices for Optical Computing

We apply inverse-design methods to produce two-dimensional plasma metamaterial (PMM) devices. Backpropagated finite-difference frequency-domain (FDFD) simulations are used to design waveguides and demultiplexers operating under both TE and TM modes. Demultiplexing and waveguiding are demonstrated for devices composed of plasma elements with reasonable plasma densities of approximately 7 GHz, allowing for future in situ training and experimental realization of these designs. We also explore the possible applicability of PMMs to nonlinear boolean operations for use in optical computing. Functionally complete logical connectives (or and and) are achieved in the TM mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

A computer program for the calculation of laminar and turbulent boundary layer flows

The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.

Dwyer, H. A.↗
NASA NTRS1974

Finite-amplitude thermal convection in a spherical shell

The properties of finite-amplitude thermal convection for a Boussinesq fluid contained in a spherical shell are investigated. All nonlinear terms are retained in the equations, and both axisymmetric and nonaxisymmetric solutions are studied. The velocity is expanded in terms of poloidal and toroidal vectors. Spherical surface harmonics resolve the horizontal structure of the flow, but finite differences are used in the vertical. With a few modifications, the transform method developed by Orszag (1970) is used to calculate the nonlinear terms, while Green's function techniques are applied to the poloidal equation and diffusion terms.

Young, R. E.↗
NASA NTRS1974

Subsonic and supersonic indicial aerodynamics and aerodynamic transfer function for complex configurations

A general theory for indicial-potential-compressible aerodynamics around complex configurations is presented. The motion is assumed to consist of constant subsonic or supersonic speed (steady state) and small perturbations around the steady state. Using the finite-element method to discretize the space problem, a set of differential-difference equations in time relating the potential to its normal derivative on the surface of the body was obtained. The aerodynamics transfer function was derived by using standard method of operational calculus.

Morino, L.↗
NASA NTRS1977

An efficient numerical technique for calculating thermal spreading resistance

An efficient numerical technique for solving the equations resulting from finite difference analyses of fields governed by Poisson's equation is presented. The method is direct (noniterative)and the computer work required varies with the square of the order of the coefficient matrix. The computational work required varies with the cube of this order for standard inversion techniques, e.g., Gaussian elimination, Jordan, Doolittle, etc.

Gale, E. H., Jr.↗
NASA NTRS1981

A remark on the theory of measuring thermal diffusivity by the modified Angstrom's method

A theory of the measurement of the thermal diffusivity of a sample by the modified Angstrom method is developed for the case in which radiative heat loss from the end surface of the sample is not negligible, and applied to measurements performed on lunar samples. Formulas allowing sample thermal diffusivity to be determined from the amplitude decay and phase lag of a temperature wave traveling through the sample are derived for a flat disk sample for which only heat loss from the end surface is important, and a sample of finite diameter and length for which heat loss through the end and side surfaces must be considered. It is noted that in the case of a flat disk, measurements at a single angular frequency of the temperature wave are sufficient, while the sample of finite diameter and length requires measurements at two discrete angular frequencies. Comparison of the values of the thermal diffusivities of two lunar samples of dimensions approximately 1 x 1 x 2 cm derived by the present methods and by the Angstrom theory for a finite bar reveals them to differ by not more than 5%, and indicates that more refined data are required as the measurement theory becomes more complicated.

Horai, K.-I.↗
NASA NTRS1982

Radiative transfer with partial frequency redistribution in inhomogeneous atmospheres - Application to the Jovian aurora

A direct finite difference numerical solution for the equation of radiative transfer by the Feautrier method is developed for use in planetary atmospheres. The procedure described here uses a plane-parallel atmosphere, and can treat partial frequency redistribution, inhomogeneity, external or internal sources, and various boundary conditions. Isotropic scattering is assumed, but in the case of no frequency redistribution, Rayleigh scattering can also be handled. A program utilizing this method is tested in a variety of situations against more powerful and elaborate methods. The case of the Lyman alpha aurora on Jupiter is then considered, where the effects of partial frequency redistribution are shown to be of great importance. New results for the detailed line profiles for Lyman alpha in the Jovian aurora are presented. The method is quite versatile, and should be especially useful in studying a wide range of problems related to auroral or dayglow emissions in planetary atmospheres.

Gladstone, G. R.↗
NASA NTRS1984

An efficient coordinate transformation technique for unsteady, transonic aerodynamic analysis of low aspect-ratio wings

An efficient coordinate transformation technique is presented for constructing grids for unsteady, transonic aerodynamic computations for delta-type wings. The original shearing transformation yielded computations that were numerically unstable and this paper discusses the sources of those instabilities. The new shearing transformation yields computations that are stable, fast, and accurate. Comparisons of those two methods are shown for the flow over the F5 wing that demonstrate the new stability. Also, comparisons are made with experimental data that demonstrate the accuracy of the new method. The computations were made by using a time-accurate, finite-difference, alternating-direction-implicit (ADI) algorithm for the transonic small-disturbance potential equation.

Guruswamy, G. P.↗
NASA NTRS1987

A computational method for the determination of the response of a linear system

A number of issues are discussed relating to the analysis of a linear damped oscillator equation with a forcing term where the right side of the equation, g(t), is known only at fixed, equal time intervals. (A particular class of such equations has application to civil earthquake engineering.) It is noted that the equation is not a differential equation because the condition of uniqueness does not obtain. A general computational method is presented for calculating x(t), based on an unconditionally stable finite difference technique. The work generalizes the results given in the recent paper of Ly (1984).

Mickens, R. E.↗
NASA NTRS1989

Integral equation solution of the full potential equation for transonic flows

An integral equation method for solving the full potential equation has been developed for arbitrary configurations in two- or three-dimensional transonic flows. This method is capable of capturing shocks using Murman-Cole type of finite difference scheme and is capable of predicting accurate and force-free wake shape as well. A rectangular grid combined with a technique of local grid refinement greatly improved the computational efficiency.

Chu, Li-Chuan↗
NASA NTRS1989

Solution of electromagnetic scattering problems using time domain techniques

New methods are developed to calculate the electromagnetic diffraction or scattering characteristics of objects of arbitrary material and shape. The methods extend the efforts of previous researchers in the use of finite-difference and pulse response techniques. Examples are given of the scattering from infinite conducting and nonconducting cylinders, open channel, sphere, cone, cone sphere, coated disk, open boxes, and open and closed finite cylinders with axially incident waves.

Britt, Charles L.↗