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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Numerical simulation using vorticity-vector potential formulation

An accurate and efficient computational method is needed for three-dimensional incompressible viscous flows in engineering applications. On solving the turbulent shear flows directly or using the subgrid scale model, it is indispensable to resolve the small scale fluid motions as well as the large scale motions. From this point of view, the pseudo-spectral method is used so far as the computational method. However, the finite difference or the finite element methods are widely applied for computing the flow with practical importance since these methods are easily applied to the flows with complex geometric configurations. However, there exist several problems in applying the finite difference method to direct and large eddy simulations. Accuracy is one of most important problems. This point was already addressed by the present author on the direct simulations on the instability of the plane Poiseuille flow and also on the transition to turbulence. In order to obtain high efficiency, the multi-grid Poisson solver is combined with the higher-order, accurate finite difference method. The formulation method is also one of the most important problems in applying the finite difference method to the incompressible turbulent flows. The three-dimensional Navier-Stokes equations have been solved so far in the primitive variables formulation. One of the major difficulties of this method is the rigorous satisfaction of the equation of continuity. In general, the staggered grid is used for the satisfaction of the solenoidal condition for the velocity field at the wall boundary. However, the velocity field satisfies the equation of continuity automatically in the vorticity-vector potential formulation. From this point of view, the vorticity-vector potential method was extended to the generalized coordinate system. In the present article, we adopt the vorticity-vector potential formulation, the generalized coordinate system, and the 4th-order accurate difference method as the computational method. We present the computational method and apply the present method to computations of flows in a square cavity at large Reynolds number in order to investigate its effectiveness.

Tokunaga, Hiroshi↗
NASA NTRSDate not supplied

The Effect of Fiber-Angle Fidelity on the Linear Response of Tow-steered Composite Plates

Composite laminate tailoring is traditionally performed by uniformly changing the in-plane ply orientation to obtain the desired mechanical performance. The emergence of tow-steered plies, where the fibers follow a prescribed curvilinear path, have increased the tailorability of composite laminates. However, characterizing the behavior of tow-steered laminates using finite element analysis is challenging because additional, and often numerous, orientation definitions may be required. The additional orientation definitions detrimentally increase the computational cost and hinder the use of advanced analysis techniques, such as Monte Carlo or uncertainty quantification, in the design process. To reduce computational cost without adversely affecting mechanics-based performance, the results of a parametric study that was used to investigate the effects of fiber-angle fidelity, element size, and tow-steered radius-of-curvature on various loading scenarios for tow-steered composite plates are presented. The three loading scenarios that were analyzed using finite element analysis included an axial tension load, an applied constant through-thickness-direction pressure load, and an axial compression load. Results for mechanics-based metrics of interest are presented and discussed for each loading scenario. Little sensitivity (less than one percent difference) to the effect of fiber-angle fidelity is observed in the mechanics-based metrics until the coarse-end of the considered range. Sensitivity to element size generally dominates the observed results for the mechanics-based metrics. Notable reductions in the preprocessing time are observed for increasingly coarse element size and fiber-rounding parameters. The preprocessing times decreased up to three orders of magnitude from a few thousand seconds to a few seconds without loss of accuracy in the mechanics-based metrics. Such increased computational performance is of particular interest to the structural design and analysis communities that may be conducting large counts of finite-element analyses, such as in other parametric studies, Monte-Carlo analyses, uncertainty quantification, or tow-steered optimization.

tow-steered composites↗
DOE OSTI2018

Exact sum rules for quantum billiards of arbitrary shape

Highlights: Sum rules of order one for domains of arbitrary shape are considered. The sum rules are proved to be finite. Exact expressions for number of examples are obtained. These expressions are verified with accurate numerical tests. -- Abstract: We have derived explicit expressions for the sum rules of order one of the eigenvalues of the negative Laplacian on two dimensional domains of arbitrary shape. Taking into account the leading asymptotic behavior of these eigenvalues, as given from Weyl’s law, we show that it is possible to define sum rules that are finite, using different prescriptions. We provide the explicit expressions and test them on a number of non trivial examples, comparing the exact results with precise numerical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Quasinormal modes, echoes and the causal structure of the Green's function

Quasinormal modes describe the return to equilibrium of a perturbed system, in particular the ringdown phase of a black hole merger. But as globally-defined quantities, the quasinormal spectrum can be highly sensitive to global structure, including distant small perturbations to the potential. In what sense are quasinormal modes a property of the resulting black hole? We explore this question for the linearized perturbation equation with two potentials having disjoint bounded support. We give a composition law for the Wronskian that determines the quasinormal frequencies of the combined system. We show that over short time scales the evolution is governed by the quasinormal frequencies of the individual potentials, while the sensitivity to global structure can be understood in terms of echoes. We introduce an echo expansion of the Green’s function and show that, as expected on general grounds, at any finite time causality limits the number of echoes that can contribute. We illustrate our results with the soluble example of a pair of -function potentials. We explicate the causal structure of the Green’s function, demonstrating under what conditions two very different quasinormal spectra give rise to very similar ringdown waveforms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Design of monolithic piezoelectric bimorph mirrors made from lithium niobate

Recent advances in deformable mirrors based on monolithic piezoelectric substrates, such as lithium niobate, have the potential to improve the image quality of X-ray optical systems on synchrotron and free-electron laser beamlines and microscopes. However, the quantitative relationship between design parameters and the deformed shape has not been readily available in the literature. We present an analytical model, validated through finite element analysis, enabling calculation of tangential and sagittal curvatures based on mirror dimensions, crystallographic orientation, and applied voltage. We demonstrate that through the selection of material orientation, it is possible to achieve different deformed shapes (sphere, cylinder, or hyperbolic paraboloid). This methodology can be generalized to other piezoelectric materials and substrate-integrated actuator systems.

Marzari, Francesco [University of Trento (Italy); ↗
DOE OSTI2022

An Exact and Near-Exact Distribution Approach to the Behrens–Fisher Problem

The Behrens–Fisher problem occurs when testing the equality of means of two normal distributions without the assumption that the two variances are equal. This paper presents approaches based on the exact and near-exact distributions for the test statistic of the Behrens–Fisher problem, depending on different combinations of even or odd sample sizes. We present the exact distribution when both sample sizes are odd and the near-exact distribution when one or both sample sizes are even. The near-exact distributions are based on a finite mixture of generalized integer gamma (GIG) distributions, used as an approximation to the exact distribution, which consists of an infinite series. The proposed tests, based on the exact and the near-exact distributions, are compared with Welch’s t-test through Monte Carlo simulations, in particular for small and unbalanced sample sizes. The results show that the proposed approaches are competent solutions to the Behrens–Fisher problem, exhibiting precise sizes and better powers than Welch’s approach for those cases. Numerical studies show that the Welch’s t-test tends to be a bit more conservative than the test statistics based on the exact or near-exact distribution, in particular when sample sizes are small and unbalanced, situations in which the proposed exact or near-exact distributions obtain higher powers than Welch’s t-test.

Hong, Serim↗
NASA NTRS1973

Application of the generalized Galerkin method to the computation of fluid flows.

The purpose of this paper is to show that most existing methods for the calculation of fluid flows can be interpreted as special applications of a single mathematical formalism. This formalism, called the generalized Galerkin method, then provides a single conceptual framework for comparing various methods in terms of convergence and accuracy. The method is presented in sufficient mathematical detail to permit its interpretation as a projection in function space. In order to demonstrate the basic thesis, finite difference, finite element, strip integral, and classical integral methods of boundary layer theory are developed by application of the method.

Murphy, J. D.↗
NASA NTRS1978

Transonic flow over the NACA 64A006 with an oscillating flap-calculations based on the Euler equations

Exploratory calculations of transonic flows over the airfoil with a quarter-chord oscillating flap were made using a program which obtains approximate solutions to the Euler equations with an explicit, shock-capturing, finite-difference scheme. The calculations, essentially inviscid and for the airfoil at zero angle-of-attack in a free-stream, are at Mach numbers and reduced frequencies which were tested in experiments by Tijdeman. The oscillatory lifts from analogus calculations by various investigators generally agree with one another better than they agree with Tijdeman's data. Inclusion in the calculations of an approximate modelling of boundary conditions expected at slotted wind tunnel walls tends to shift some of the results closer to the experimental values.

Magnus, R. J.↗
NASA NTRS1980

On the solution of the unsteady Navier-Stokes equations for hypersonic flow about axially-symmetric blunt bodies

A formulation of the complete Navier-Stokes problem for a viscous hypersonic flow in general curvilinear coordinates is presented. This formulation is applicable to both the axially symmetric and three dimensional flows past bodies of revolution. The equations for the case of zero angle of attack were solved past a circular cylinder with hemispherical caps by point SOR finite difference approximation. The free stream Mach number and the Reynolds number for the test case are respectively 22.04 and 168883. The whole algorithm is presented in detail along with the preliminary results for pressure, temperature, density and velocity distributions along the stagnation line.

Warsi, Z. U. A.↗
NASA NTRS1985

Structural optimization by generalized, multilevel decomposition

The developments toward a general multilevel optimization capability and results for a three-level structural optimization are described. The method partitions a structure into a number of substructuring levels where each substructure corresponds to a subsystem in the general case of an engineering system. The method is illustrated by a portal framework that decomposes into individual beams. Each beam is a box that can be further decomposed into stiffened plates. Substructuring for this example spans three different levels: (1) the bottom level of finite elements representing the plates; (2) an intermediate level of beams treated as substructures; and (3) the top level for the assembled structure. The three-level case is now considered to be qualitatively complete.

Sobieszczanski-Sobieski, J.↗
NASA NTRS1993

The uncertainty principle and quantum chaos

The conception of quantum chaos is described in some detail. The most striking feature of this novel phenomenon is that all the properties of classical dynamical chaos persist here but, typically, on the finite and different time scales only. The ultimate origin of such a universal quantum stability is in the fundamental uncertainty principle which makes discrete the phase space and, hence, the spectrum of bounded quantum motion. Reformulation of the ergodic theory, as a part of the general theory of dynamical systems, is briefly discussed.

Chirikov, Boris V.↗
NASA NTRS1993

The annual pressure cycle on Mars: Results from the LMD Martian atmospheric general circulation model

We have been developing a General Circulation Model (GCM) of the martian atmosphere since 1989. The model has been described rather extensively elsewhere and only the main characteristics are given here. The dynamical part of the model, adapted from the LMD terrestrial climate model, is based on a finite-difference formulation of the classical 'primitive equations of meteorology.' The radiative transfer code includes absorption and emission by CO2 (carefully validated by comparison to line-by-line calculations) and dust in the thermal range and absorption and scattering by dust in the visible range. Other physical parameterizations are included: modeling of vertical turbulent mixing, dry convective adjustment (in order to prevent vertical unstable temperature profiles), and a multilayer model of the thermal conduction in the soil. Finally, the condensation-sublimation of CO2 is introduced through specification of a pressure-dependent condensation temperature. The atmospheric and surface temperatures are prevented from falling below this critical temperature by condensation and direct precipitation onto the surface of atmospheric CO2. The only prespecified spatial fields are the surface thermal inertia, albedo, and topography.

Hourdin, Frederic↗
NASA NTRS2001

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗
DOE OSTI2016

Non-singular Brans–Dicke collapse in deformed phase space

We study the collapse process of a homogeneous perfect fluid (in FLRW background) with a barotropic equation of state in Brans–Dicke (BD) theory in the presence of phase space deformation effects. Such a deformation is introduced as a particular type of non-commutativity between phase space coordinates. For the commutative case, it has been shown in the literature (Scheel, 1995), that the dust collapse in BD theory leads to the formation of a spacetime singularity which is covered by an event horizon. In comparison to general relativity (GR), the authors concluded that the final state of black holes in BD theory is identical to the GR case but differs from GR during the dynamical evolution of the collapse process. However, the presence of non-commutative effects influences the dynamics of the collapse scenario and consequently a non-singular evolution is developed in the sense that a bounce emerges at a minimum radius, after which an expanding phase begins. Such a behavior is observed for positive values of the BD coupling parameter. For large positive values of the BD coupling parameter, when non-commutative effects are present, the dynamics of collapse process differs from the GR case. Finally, we show that for negative values of the BD coupling parameter, the singularity is replaced by an oscillatory bounce occurring at a finite time, with the frequency of oscillation and amplitude being damped at late times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

A mixed pseudospectral/finite difference method for a thermally driven fluid in a nonuniform gravitational field

A numerical study of the steady, axisymmetric flow in a heated, rotating spherical shell is conducted to model the Atmospheric General Circulation Experiment (AGCE) proposed to run aboard a later shuttle mission. The AGCE will consist of concentric rotating spheres confining a dielectric fluid. By imposing a dielectric field across the fluid a radial body force will be created. The numerical solution technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is based on the incompressible Navier-Stokes equations. In the method a pseudospectral technique is used in the latitudinal direction, and a second-order accurate finite difference scheme discretizes time and radial derivatives. This paper discusses the development and performance of this numerical scheme for the AGCE which has been modelled in the past only by pure FD formulations. In addition, previous models have not investigated the effect of using a dielectric force to simulate terrestrial gravity. The effect of this dielectric force on the flow field is investigated as well as a parameter study of varying rotation rates and boundary temperatures. Among the effects noted are the production of larger velocities and enhanced reversals of radial temperature gradients for a body force generated by the electric field.

Macaraeg, M. G.↗
DOE OSTI2020

Optical analogues to the equatorial Kerr–Newman black hole

Optical analogues to black holes allow the investigation of general relativity in a laboratory setting. Previous works have considered analogues to Schwarzschild black holes in an isotropic coordinate system; the major drawback is that required material properties diverge at the horizon. We present the dielectric permittivity and permeability tensors that exactly reproduce the equatorial Kerr–Newman metric, as well as the gradient-index material that reproduces equatorial Kerr–Newman null geodesics. Importantly, the radial profile of the scalar refractive index is finite along all trajectories except at the point of rotation reversal for counter-rotating geodesics. Construction of these analogues is feasible with available ordinary materials. A finite-difference frequency-domain solver of Maxwell’s equations is used to simulate light trajectories around a variety of Kerr–Newman black holes. For reasonably sized experimental systems, ray tracing confirms that null geodesics can be well-approximated in the lab, even when allowing for imperfect construction and experimental error.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Boundary Corrections for Kernel Approximation to Differential Operators

The kernel-based approach to operator approximation for partial differential equations has been shown to be unconditionally stable for linear PDEs and numerically exhibit unconditional stability for non-linear PDEs. These methods have the same computational cost as an explicit finite difference scheme but can exhibit order reduction at boundaries. In previous work on periodic domains, order reduction was addressed, yielding high-order accuracy. The issue addressed in this work is the elimination of order reduction of the kernel-based approach for a more general set of boundary conditions. Further, we consider the case of both first and second order operators. To demonstrate the theory, we provide not only the mathematical proofs but also experimental results by applying various boundary conditions to different types of equations. The results agree with the theory, demonstrating a systematic path to high order for kernel-based methods on bounded domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

On the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers: Selective mode-type excitation and specific mode enhancement

Here, we investigate the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers. Non-uniform polarization is represented as a non-uniform electromechanical coupling coefficient using a polarization function. Governing equations are derived for the forced in-plane vibrations of a thin wafer under the assumption of generalized plane stress. The effect of non-uniform polarization is explicitly obtained in the forcing terms of the governing equations. These equations are then recast into a variational weak form that is then solved using the finite element method to obtain the displacement fields for different modes. The electromechanical response of a non-uniformly polarized piezoelectric wafer is derived in terms of the out-of-plane displacement profile on the surface of the wafer. Using the derived analytical expression, a necessary and sufficient condition for the presence/absence of a vibrational mode in the electromechanical impedance spectrum is obtained. Based on this condition, criteria for selective mode-type excitation and specific mode enhancement of vibrational modes in the electromechanical impedance spectrum are postulated. Selective mode-type excitation of in-plane extensional and shear modes is demonstrated for a square wafer and that of in-plane bending modes is demonstrated for a rectangular wafer. Specific mode enhancement is demonstrated for both square and rectangular wafers. In addition, it is also demonstrated how the criteria can be used to suppress specific vibrational modes in the electromechanical impedance spectrum. The proposed methodology of using non-uniformly polarized piezoelectric wafers finds application in the design of single element transducers with multi-frequency operation, frequency-tuned receivers/sensors, acoustic holograms, designing acoustic beams of prescribed shape/lobes, and other non-traditional applications such as information storage.

36 MATERIALS SCIENCE↗