Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Modeling of noise-assisted quantum transfer between donor and acceptor with finite bandwidths

We present an analytical and numerical study of noise-assisted quantum exciton (electron) transfer (ET) in a dimer consisting of donor and acceptor, modeled by interacting continuous electron bands of finite widths. We start with discrete energy levels for both bands. The interaction with the environment is modeled by a stationary stochastic process (noise) acting on all the donor and acceptor energy levels. Then, by using a continuous limit for the electron spectra, we derive integro-differential equations for the ET dynamics between the two bands. Finally, we derive from these equations rate-type differential equations for the ET dynamics. We formulate the conditions of validity of the rate-type equations. We consider different regions of parameters characterizing the widths of the donor and acceptor bands and the strength of the dimer-noise interaction. For a simplified model with a single energy level donor and a continuous acceptor band, we derive a generalized simple analytic expression and provide numerical simulations for the ET rate. They are consistent with Wigner–Weisskopf, Förster-type, and Marcus-type expressions, in their corresponding regime of parameters. For a weak dimer-noise interaction, our approach leads to the Wigner–Weisskopf ET rate, or to the Förster-type ET rate, depending on other parameters. In the limit of strong dimer-noise interaction, our approach is non-perturbative in the dimer-noise interaction constant, and it recovers the Marcus-type ET rate. Our analytic results are confirmed by numerical simulations. We demonstrate how our theoretical results are modified when both the donor and the acceptor are described by finite bands. We also show that, for a relatively wide acceptor band, the efficiency of the ET from donor to acceptor can be close to 100% for a broad range of noise amplitudes, for both 'downhill' and 'uphill' ET, for sharp and flat redox potentials, and for reasonably short times. We discuss possible applications of our approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

A design procedure for fan-out improvement in all-optical photonic crystal logic design

In this paper, a general method is proposed to improve the fan-out parameter for all-optical logic gates and functions. Two different types of devices are designed to increase the fan-out of a logic operation and its inverted logic. Nonlinear cavities are used to design driving units. Silicon nanocrystal is used as the nonlinear material to create the required frequency shift for different values of input power. Plane wave expansion and finite difference time domain methods are used to simulate and analyze the proposed structures. The propagation delay of the proposed structure is less than 1.5 ps and the maximum required power for fan-out of three is 15 W.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

High resolution applications of the Osher upwind scheme for the Euler equations

The 'Osher' scheme was introduced by Osher (1981). It represents an upwind finite-difference method for hyperbolic systems of conservation laws, including the Euler equations. In studies conducted by Osher (1981) and Osher and Solomon (1982), the method was applied to the nonisentropic form of the Euler equations in one dimension and the isentropically restricted form in two spatial dimensions, both in Cartesian coordinates. Chakravarthy and Osher (1982) have shown an approach for extending the Osher scheme to the Euler equations written for general geometries, taking into account the use of mappings to arbitrary curvilinear coordinate systems. The present investigation is concerned with the high resolution extension of the Osher scheme to second-order accuracy. Results are presented for several example problems, giving attention to quasi-one-dimensional Laval nozzle flow, a one-dimensional shock tube problem, and supersonic flow over a cylinder.

Chakravarthy, S. R.↗
NASA NTRS1990

Numerical derivative techniques for trajectory optimization

The adoption of robust numerical optimization techniques in trajectory simulation programs has resulted in powerful design and analysis tools. These trajectory simulation/optimization programs are widely used, and a representative list includes the GTS system, the POST program, and newer collocation methods such as OTIS and FONPAC. All of these programs rely on optimization algorithms which require objective function and constraint gradient data during the iteration process. However, most trajectory optimization problems lack simple analytical expressions for these derivatives. In the general case a function evaluation involves integrating aerodynamic, propulsive, and gravity forces over multiple trajectory phases with complex control models. With the newer collocation methods, the integration is replaced by defect constraints and cubic approximations for the state. While analytic gradient expressions can sometimes be derived for trajectory optimization problems, the derivation is cumbersome, time consuming, and prone to mistakes. Fortunately, an alternate method exists for the gradient evaluation, namely finite difference approximations. In this paper some finite difference gradient techniques developed for use with the GTS system are presented. These techniques include methods for computing first and second partial derivatives of single and multiple sets of functions. A key feature of these methods is an error control mechanism which automatically adjusts the perturbation size to obtain accurate derivative values.

Hallman, Wayne P.↗
NASA NTRS1994

Accurate finite difference methods for time-harmonic wave propagation

Finite difference methods for solving problems of time-harmonic acoustics are developed and analyzed. Multidimensional inhomogeneous problems with variable, possibly discontinuous, coefficients are considered, accounting for the effects of employing nonuniform grids. A weighted-average representation is less sensitive to transition in wave resolution (due to variable wave numbers or nonuniform grids) than the standard pointwise representation. Further enhancement in method performance is obtained by basing the stencils on generalizations of Pade approximation, or generalized definitions of the derivative, reducing spurious dispersion, anisotropy and reflection, and by improving the representation of source terms. The resulting schemes have fourth-order accurate local truncation error on uniform grids and third order in the nonuniform case. Guidelines for discretization pertaining to grid orientation and resolution are presented.

Harari, Isaac↗
NASA NTRS2001

Evaluation of Stratospheric Transport in New 3D Models Using the Global Modeling Initiative Grading Criteria

The Global Modeling Initiative (GMI) Team developed objective criteria for model evaluation in order to identify the best representation of the stratosphere. This work created a method to quantitatively and objectively discriminate between different models. In the original GMI study, 3 different meteorological data sets were used to run an offline chemistry and transport model (CTM). Observationally-based grading criteria were derived and applied to these simulations and various aspects of stratospheric transport were evaluated; grades were assigned. Here we report on the application of the GMI evaluation criteria to CTM simulations integrated with a new assimilated wind data set and a new general circulation model (GCM) wind data set. The Finite Volume Community Climate Model (FV-CCM) is a new GCM developed at Goddard which uses the NCAR CCM physics and the Lin and Rood advection scheme. The FV-Data Assimilation System (FV-DAS) is a new data assimilation system which uses the FV-CCM as its core model. One year CTM simulations of 2.5 degrees longitude by 2 degrees latitude resolution were run for each wind data set. We present the evaluation of temperature and annual transport cycles in the lower and middle stratosphere in the two new CTM simulations. We include an evaluation of high latitude transport which was not part of the original GMI criteria. Grades for the new simulations will be compared with those assigned during the original GMT evaluations and areas of improvement will be identified.

Strahan, Susan E.↗
NASA NTRS2020

OGLE-2013-BLG-0911Lb: A Secondary on the Brown-dwarf Planet Boundary around an M Dwarf

We present the analysis of the binary-lens microlensing event OGLE-2013-BLG-0911. The best-fit solutions indicate the binary mass ratio of q 0.03, which differs from that reported in Shvartzvald et al. The event suffers from the well-known close/wide degeneracy, resulting in two groups of solutions for the projected separation normalized by the Einstein radius of s ~ 0.15 or s ~ 7. The finite source and the parallax observations allow us to measure the lens physical parameters. The lens system is an M dwarf orbited by a massive Jupiter companion at very close (M(host) = 0.30(sup +0.08, sub -0.06)Mʘ, M(comp) = 10.1 (sup +2.9, sub -2.2 M(jup), a(exp) = 0.40 (sup +0.05, sub -0.40)au) or wide (M(host) = 0.28 (sup +0.10, sub -0.80)Mʘ, M(comp) = 9.9 (sup +3.8, sub -3.5)M(jup), a(exp) = 18.0 (sup +3.2, sub -3.2)au) separation. Although the mass ratio is slightly above the planet-brown dwarf (BD) mass-ratio boundary of q = 0.03, which is generally used, the median physical mass of the companion is slightly below the planet-BD mass boundary of 13M(Jup). It is likely that the formation mechanisms for BDs and planets are different and the objects near the boundaries could have been formed by either mechanism. It is important to probe the distribution of such companions with masses of ~13M(Jup) in order to statistically constrain the formation theories for both BDs and massive planets. In particular, the microlensing method is able to probe the distribution around low-mass M dwarfs and even BDs, which is challenging for other exoplanet detection methods.

Shota Miyazaki↗
DOE OSTI2017

Preconditioned steepest descent methods for some nonlinear elliptic equations involving p-Laplacian terms

We describe and analyze preconditioned steepest descent (PSD) solvers for fourth and sixth-order nonlinear elliptic equations that include p-Laplacian terms on periodic domains in 2 and 3 dimensions. The highest and lowest order terms of the equations are constant-coefficient, positive linear operators, which suggests a natural preconditioning strategy. Such nonlinear elliptic equations often arise from time discretization of parabolic equations that model various biological and physical phenomena, in particular, liquid crystals, thin film epitaxial growth and phase transformations. The analyses of the schemes involve the characterization of the strictly convex energies associated with the equations. We first give a general framework for PSD in Hilbert spaces. Based on certain reasonable assumptions of the linear pre-conditioner, a geometric convergence rate is shown for the nonlinear PSD iteration. We then apply the general theory to the fourth and sixth-order problems of interest, making use of Sobolev embedding and regularity results to confirm the appropriateness of our pre-conditioners for the regularized p-Lapacian problems. Our results include a sharper theoretical convergence result for p-Laplacian systems compared to what may be found in existing works. We demonstrate rigorously how to apply the theory in the finite dimensional setting using finite difference discretization methods. Numerical simulations for some important physical application problems – including thin film epitaxy with slope selection and the square phase field crystal model – are carried out to verify the efficiency of the scheme.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

PLATO (Parallel Load Assignment Tool): A Parallel Workload Partitioner for Particle-based Methods

A key first step in many large-scale high-performance scientific computing applications is the distribution of the computational domain across all participating processors of a parallel computing platform such as a shared memory machine, cluster or a tightly-coupled supercomputer. The primary goal of such a partitioning step is to ensure that the computational work performed by each process remains load-balanced throughout the execution of the parallel application. In this project, a general-purpose parallel partitioning tool, code-named PLATO, was designed and implemented. PLATO assumes that the workload is a set of points in an n-dimensional metric space. These points are abstractions of computational units that may represent point-like computational units such as celestial objects in astrophysics, atoms in materials sciences and molecules in biology/chemistry, or may represent more abstract computational units such as grid points in finite difference methods, grid elements in finite element methods or n-dimensional points in the feature space of a host of other applications. PLATO supports several key and useful capabilities, namely: (a) multiple load partitioning strategies, (b) support for arbitrary user-defined processor topology, (c) pre-fetching data from neighboring processors based on user-defined parameters, and (d) support for arbitrary number and types of attributes associated with each point in the computational domain. PLATO is implemented in the C programming language using Message Passing Interface (MPI) for inter-processor communications and its parallel performance tested on up to 48 processors of a shared memory machine. This report includes a detailed performance study of PLATO using parallel speedup and iso-efficiency as the two main metrics of parallel performance. In addition, a user manual and DOXYGEN generated code documentation have been included for completeness.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Thermal buoyancy and Marangoni convection in a two fluid layered system - A numerical study

Thermal buoyancy and surface tension driven convection is numerically investigated in a system with two, immiscible fluids. The geometry investigated, has an open cavity configuration with the lighter fluid situated on top of the heavier fluid, forming a stable layered system. The upper fluid meniscus and the interface are assumed to be flat and undeformable in the calculations. The governing equations and boundary/interface conditions are solved by a control volume based finite difference scheme for two pairs of immiscible fluids; the water-hexadecane system and a so called generic system. The steady state calculations show, that dramatically different flows are predicted when the interfacial tension effects are included or excluded from the system model. These differences are particularly appreciable in surface tension dominated flows, that are typical of microgravity situations. Complex flow patterns, with induced secondary flows are noticed in both the fluids. In general, the overall system heat transfer is found to increase with increases in the thermal buoyancy and surface tension effects, but the behavior of the system flow and thermal fields is not easy to characterize, when different combinations of these forces are considered.

Ramachandran, N.↗
NASA NTRS1990

An easily implemented static condensation method for structural sensitivity analysis

A black-box approach to static condensation for sensitivity analysis is presented with illustrative examples of a cube and a car structure. The sensitivity of the structural response with respect to joint stiffness parameter is calculated using the direct method, forward-difference, and central-difference schemes. The efficiency of the various methods for identifying joint stiffness parameters from measured static deflections of these structures is compared. The results indicate that the use of static condensation can reduce computation times significantly and the black-box approach is only slightly less efficient than the standard implementation of static condensation. The ease of implementation of the black-box approach recommends it for use with general-purpose finite element codes that do not have a built-in facility for static condensation.

Gangadharan, S. N.↗
NASA NTRS1991

Design and analysis of grid stiffened concepts for aircraft composite primary structural applications

An approach to buckling resistant design of general grid stiffened flat plates based on smeared stiffener theory for combined inplane loading is discussed. Some results from parametric studies performed to assess the validity of smeared stiffener for practical stiffener configurations and to illustrate the benefits of different stiffening concepts are given. Details of a design study are discussed where the present analysis method is used to design a grid stiffened panel for a fuselage application and verified using a finite element analysis results.

Ambur, Damodar R.↗
NASA NTRS1998

Flow/Damage Surfaces for Fiber-Reinforced Metals Having Different Periodic Microstructures

Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/ titanium composite system having rectangular, hexagonal, and square diagonal fiber packing arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics, using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue; for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.

Lissenden, Cliff J.↗
NASA NTRS1998

Flow/Damage Surfaces for Fiber-Reinforced Metals having Different Periodic Microstructures

Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/ titanium composite system having rectangular, hexagonal, and square diagonal fiber packing, arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics. using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber-matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue, for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.

Lissenden, Cliff J.↗
DOE OSTI2016

Analysis of the static magnetic field-dependent optical transmission of Ni nanorod colloidal suspensions

The magnetic field-dependent optical transmission of dilute Ni nanorod aqueous suspensions was investigated. A series of four samples of nanorods were synthesized using the AAO template method and processed to stable colloids. The distributions of their length and diameter were characterized by analysis of TEM images and revealed average diameters of ∼25 nm and different lengths in the range of 60 nm–1100 nm. The collinear magnetic and optical anisotropy was studied by static field-dependent transmission measurements of linearly polarized light parallel and perpendicular to the magnetic field direction. The experimental results were modelled assuming the field-dependent orientation distribution function of a superparamagnetic ensemble for the uniaxial ferromagnetic nanorods in liquid dispersion and extinction cross sections for longitudinal and transversal optical polarization derived from different approaches, including the electrostatic approximation and the separation of variables method, both applied to spheroidal particles, as well as finite element method simulations of spheroids and capped cylindrical particles. The extinction cross sections were compared to reveal the differences associated with the approximations of homogeneous polarization and/or particle shape. The consequences of these approximations for the quantitative analysis of magnetic field-dependent optical transmission measurements were investigated and a reliable protocol derived. Furthermore, the changes in optical cross sections induced by electromagnetic interaction between two nanorods in parallel end-to-end and side-by-side configuration as a function of their separation were studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Superconductivity in multiorbital systems with repulsive interactions: Hund's pairing versus spin-fluctuation pairing

Hund's pairing refers to Cooper pairing generated by onsite interactions that become attractive due to large Hund's exchange J. This is possible in multiorbital systems even when all local bare interactions are repulsive, since attractions in specific channels are given by certain linear combinations of interaction parameters. On the other hand, pairing processes such as the exchange of spin fluctuations are also present. We compare mean-field Hund's pairing and spin-fluctuation-mediated pairing using electronic bands appropriate for different classes of multiorbital systems over a wide range of interaction parameters. We find that, for systems without clear nesting features, the superconducting state generated by the Hund's mechanism agrees well with that from the full fluctuation exchange vertex when Hund's exchange and spin-orbit coupling are sufficiently large. On the other hand, for systems characterized by a peaked finite-momentum particle-hole susceptibility, spin-fluctuation pairing generally dominates over Hund's pairing. From this perspective Hund's pairing states appear unlikely to be realized in systems like Sr 2 RuO 4 and generic iron-based superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Simulating coupled surface–subsurface flows with ParFlow v3.5.0: capabilities, applications, and ongoing development of an open-source, massively parallel, integrated hydrologic model

Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.

58 GEOSCIENCES↗