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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Finite element buckling and postbuckling solutions for multilayered composite panels

A study is made of the buckling and postbuckling responses of flat, unstiffened composite panels subjected to various combinations of mechanical and thermal loads. The analysis is based on a first-order shear deformation von Karman-type plate theory. A mixed formulation is used with the fundamental unknowns consisting of the strain components, stress resultants and the generalized displacements of the plate. The stability boundary, postbuckling response and the sensitivity coefficients are evaluated. The sensitivity coefficients measure the sensitivity of the buckling and postbuckling responses to variations in the different lamination and material parameters of the panel. Numerical results are presented for both solid panels and panels with central circular cutouts. The results show the effects of the variations in the fiber orientation angels, aspect ratio of the panel, and the hole diameter (for panels with cutouts) on the stability boundary, postbuckling response and sensitivity coefficients.

Noor, Ahmed K.↗
DOE OSTI2018

Light Scattering by Small Multilayer Particles: A Generalized Separation of Variables Method

A Rayleigh approximation is constructed for light scattering by small multilayer axisymmetric particles, in which their polarizability is determined by the generalized separation of variables method (SVM). In this method, scalar potentials, the gradients of which yield the electric-field strengths, are represented as expansions in spherical harmonics of the Laplace equation. Unknown coefficients of expansions are determined from the boundary conditions, which are reduced to infinite systems of linear algebraic equations (ISLAEs), since the separation of variables is incomplete. The T matrix of the electrostatic problem, principal element T{sub 11} of which is proportional to the particle polarizability, is determined. The necessary condition for the ISLAEs solvability for the SVM coincides with the condition of the correct application of the extended boundary conditions method (EBCM). However, numerical calculations in which finite-dimensional (i.e., reduced) systems are solved, yield different results in ranges of variation of parameters that are close to the boundary of the range of applicability. An analysis of the numerical calculations of the scattering and absorption cross sections for two-layer confocal spheroids, an exact solution for which can be obtained using spheroidal harmonics, shows that the SVM is preferable to the EBCM. It turned out that the proposed method yields workable results in a wider range of variation of parameters. Even outside the range of applicability, in which it should be regarded as a certain approximate solution, its use in a number of cases is quite acceptable. Additional calculations for three-layer nonconfocal spheroids, as well as for three-layer similar pseudospheroids and Pascal’s snails, which can be obtained from spheroids as a result of the inversion with respect to the coordinate origin and one of the foci, respectively, confirm these inferences. We note that, for certain values of the parameters, the shapes of the latter particles are nonconvex.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

Here, we perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient D tr (x), at various positions x within the pore. D tr (x) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.↗
NASA NTRS1985

Prediction of HiMAT RPRV component interference effects at transonic speeds

Detailed surface flow field analyses of the NASA Highly Maneuverable Aircraft Technology Remotely Piloted Research Vehicle (HiMAT RPRV) are described. Computational predictions are provided by a nonlinear finite difference transonic wing-body code. Component interference effects are obtained through modifications in lifting surface twist and camber boundary conditions. Complex aircraft interactions created by wing-tip fins, landing gear, and lower surface wing actuator fairings are included. Lifting surface viscous flow effects are modeled using a two-dimensional 'strip' boundary layer method. Correlations with flight test data presented demonstrate the high degree of accuracy achieved for this realistic aircraft configuration. This methodology and approach appears to be generally applicable to other multisurface configurations.

Laiosa, J.↗
NASA NTRS1991

The P1-RKDG method for two-dimensional Euler equations of gas dynamics

A class of nonlinearly stable Runge-Kutta local projection discontinuous Galerkin (RKDG) finite element methods for conservation laws is investigated. Two dimensional Euler equations for gas dynamics are solved using P1 elements. The generalization of the local projections, which for scalar nonlinear conservation laws was designed to satisfy a local maximum principle, to systems of conservation laws such as the Euler equations of gas dynamics using local characteristic decompositions is discussed. Numerical examples include the standard regular shock reflection problem, the forward facing step problem, and the double Mach reflection problem. These preliminary numerical examples are chosen to show the capacity of the approach to obtain nonlinearly stable results comparable with the modern nonoscillatory finite difference methods.

Cockburn, Bernardo↗
NASA NTRS1992

Multidisciplinary tailoring of hot composite structures

A computational simulation procedure is described for multidisciplinary analysis and tailoring of multilayered multimaterial hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loadings including the effect of aggressive environments. The simulation is based on a 3D finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., the laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermomechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loadings lead to different tailored designs, even those opposite of each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.

Singhal, Surendra N.↗
NASA NTRS1993

Multidisciplinary tailoring of hot composite structures

A computational simulation procedure is described for multidisciplinary analysis and tailoring of layered multi-material hot composite engine structural components subjected to simultaneous multiple discipline-specific thermal, structural, vibration, and acoustic loads. The effect of aggressive environments is also simulated. The simulation is based on a three-dimensional finite element analysis technique in conjunction with structural mechanics codes, thermal/acoustic analysis methods, and tailoring procedures. The integrated multidisciplinary simulation procedure is general-purpose including the coupled effects of nonlinearities in structure geometry, material, loading, and environmental complexities. The composite material behavior is assessed at all composite scales, i.e., laminate/ply/constituents (fiber/matrix), via a nonlinear material characterization hygro-thermo-mechanical model. Sample tailoring cases exhibiting nonlinear material/loading/environmental behavior of aircraft engine fan blades, are presented. The various multidisciplinary loads lead to different tailored designs, even those competing with each other, as in the case of minimum material cost versus minimum structure weight and in the case of minimum vibration frequency versus minimum acoustic noise.

Singhal, Surendra N.↗
NASA NTRS2002

Application of the generalized vertical coordinate ocean model for better representing satellite data

It is found that two adaptive parametric functions can be introduced into the basic ocean equations for utilizing the optimal or hybrid features of commonly used z-level, terrain- following, isopycnal, and pressure coordinates in numerical ocean models. The two parametric functions are formulated by combining three techniques: the arbitrary vertical coordinate system of Kasahara (1 974), the Jacobian pressure gradient formulation of Song (1 998), and a newly developed metric factor that permits both compressible (non-Boussinesq) and incompressible (Boussinesq) approximations. Based on the new formulation, an adaptive modeling strategy is proposed and a staggered finite volume method is designed to ensure conservation of important physical properties and numerical accuracy. Implementation of the combined techniques to SCRUM (Song and Haidvogel1994) shows that the adaptive modeling strategy can be applied to any existing ocean model without incurring computational expense or altering the original numerical schemes. Such a generalized coordinate model is expected to benefit diverse ocean modelers for easily choosing optimal vertical structures and sharing modeling resources based on a common model platform. Several representing oceanographic problems with different scales and characteristics, such as coastal canyons, basin-scale circulation, and global ocean circulation, are used to demonstrate the model's capability for multiple applications. New results show that the model is capable of simultaneously resolving both Boussinesq and non-Boussinesq, and both small- and large-scale processes well. This talk will focus on its applications of multiple satellite sensing data in eddy-resolving simulations of Asian Marginal Sea and Kurosio. Attention will be given to how Topex/Poseidon SSH, TRMM SST; and GRACE ocean bottom pressure can be correctly represented in a non- Boussinesq model.

numerical ocean model generalized coordinate syste↗
NASA NTRS2011

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Schmerr, Nicholas C.↗
NASA NTRS1983

A selection principle in Benard-type convection

In a Benard-type convection problem, the stationary flows of an infinite layer of fluid lying between two rigid horizontal walls and heated uniformly from below are determined. As the temperature difference across the layer increases beyond a certain value, other convective motions appear. These motions areoften cellular in character in that their streamlines are confined to certain well-defined cells having, for example, the shape of rolls or hexagons. A selection principle that explains why hexagonal cells seem to be preferred for certain ranges of the parameters is formulated. An operator-theoretical formulation of one generalized Bernard problem is given. The infinite dimensional problem is reduced to one of solving a finite dimensional system of equations, namely, the selection equations. These equations are solved and a linearized stability analysis of the resultant stationary flows is presented.

Knightly, G. H.↗
NASA NTRS1993

Nonreflective boundary conditions for high-order methods

A different approach to nonreflective boundary conditions for the Euler equations is presented. This work is motivated by a need for in and outflow boundary conditions that do not limit the useful accuracy of high-order accurate methods. The primary interest is in the propagation and convection of continuous acoustic and convective waves. This new approach employs the exact solution to finite waves to relate interior values and ambient conditions to boundary values. The method is first presented in one dimension and then generalized to multidimensions. Grid refinement studies are used to demonstrate high-order convergence for both one-dimensional and two-dimensional flows.

Atkins, H. L.↗
NASA NTRS1994

Nonreflective boundary conditions for high-order methods

A different approach to nonreflective boundary conditions for the Euler equations is presented. This work is motivated by a need for inflow and outflow boundary conditions that do not limit the useful accuracy of high-order accurate methods. The primary interest is in the propagation and convection of continuous acoustic and convective waves. This new approach employs the exact solution to finite waves to relate interior values and ambient conditions to boundary values. The method is first presented in one dimension and then generalized to multidimensions. Grid refinement studies are used to demonstrate high-order convergence for both one-dimensional and two-dimensional flows.

Atkins, H.↗
NASA NTRS2007

Automatic Conversion of Conceptual Geometry to CFD Geometry for Aircraft Design

Conceptual aircraft design is usually based on simple analysis codes. Its objective is to provide an overall system performance of the developed concept, while preliminary aircraft design uses high-fidelity analysis tools such as computational fluid dynamics (CFD) analysis codes or finite element structural analysis codes. In some applications, such as low-boom supersonic concept development, it is important to be able to explore a variety of drastically different configurations while using CFD analysis to check whether a given configuration can be tailored to have a low-boom ground signature. It poses an extremely challenging problem of integrating CFD analysis in conceptual design. This presentation will discuss a computer code, called iPatch, for automatic conversion of conceptual geometry to CFD geometry. In general, conceptual aircraft geometry is not as well-defined as a CAD geometry model. In particular, a conceptual aircraft geometry model usually does not define the intersection curves for the connecting surfaces. The computer code iPatch eliminates the gap between conceptual geometry and CFD geometry by accomplishing the following three tasks automatically: (1) use bicubic B-splines to extrapolate (if necessary) each surface in a conceptual geometry so that all the independently defined geometry components (such as wing and fuselage) can be intersected to form a watertight CFD geometry, (2) compute the intersection curves of surface patches at any resolution (up to 10-7 accuracy) specified by users, and (3) write the B-spline surface patches and the corresponding boundary points for the watertight CFD geometry in the format that can be directly exported to the meshing tool VGRID in the CFD software TetrUSS. As a result, conceptual designers can get quick feedback on the aerodynamic characteristics of their concepts, which will allow them to understand some subtlety in their concepts and to be able to assess their concepts with a higher degree of confidence. This integration of CFD analysis in conceptual aircraft design will greatly eliminate some uncertainty due to simple analysis codes used to develop the concepts and improve the feasibility/credibility of the final concept. The presentation will highlight the mathematical challenges of accomplishing the aforementioned three tasks and the computational algorithms used by iPatch.

Li, Wu↗
NASA NTRSDate not supplied

Modeling and Sensitivity Analysis of Sandwich Composite Cylinders with Geometric Imperfections

It is well known that real shell structures can have significantly lower buckling loads and even different mode shapes than the theoretical predictions for perfect structures. Much of this difference can be attributed to geometric and loading imperfections, and geometrically nonlinear effects. The realistic buckling response can be investigated using geometrically nonlinear finite element analyses and including radial imperfections. Such analyses are used in the NASA Engineering and Safety Center Shell Buckling Knockdown Factor Project, which has the goal of developing new analysis-based buckling design recommendations for select classes of cylindrical shell structures under uniaxial compressive load. The approach for modeling several sandwich composite cylinders with two-dimensional general-purpose shell elements and the influence of the element type selection and element size on the buckling load is discussed. The influence of geometric imperfections of various magnitudes on buckling behavior of a sandwich composite cylinder was also investigated.

Structural Modeling↗
NASA NTRS1997

An Application of the Difference Potentials Method to Solving External Problems in CFD

Numerical solution of infinite-domain boundary-value problems requires some special techniques that would make the problem available for treatment on the computer. Indeed, the problem must be discretized in a way that the computer operates with only finite amount of information. Therefore, the original infinite-domain formulation must be altered and/or augmented so that on one hand the solution is not changed (or changed slightly) and on the other hand the finite discrete formulation becomes available. One widely used approach to constructing such discretizations consists of truncating the unbounded original domain and then setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The role of the ABC's is to close the truncated problem and at the same time to ensure that the solution found inside the finite computational domain would be maximally close to (in the ideal case, exactly the same as) the corresponding fragment of the original infinite-domain solution. Let us emphasize that the proper treatment of artificial boundaries may have a profound impact on the overall quality and performance of numerical algorithms. The latter statement is corroborated by the numerous computational experiments and especially concerns the area of CFD, in which external problems present a wide class of practically important formulations. In this paper, we review some work that has been done over the recent years on constructing highly accurate nonlocal ABC's for calculation of compressible external flows. The approach is based on implementation of the generalized potentials and pseudodifferential boundary projection operators analogous to those proposed first by Calderon. The difference potentials method (DPM) by Ryaben'kii is used for the effective computation of the generalized potentials and projections. The resulting ABC's clearly outperform the existing methods from the standpoints of accuracy and robustness, in many cases noticeably speed up the multigrid convergence, and at the same time are quite comparable to other methods from the standpoints of geometric universality and simplicity of implementation.

Ryaben 'Kii, Victor S.↗
DOE OSTI2016

Enhancing the light extraction efficiency of AlGaN deep ultraviolet light emitting diodes by using nanowire structures

The performance of conventional AlGaN deep ultraviolet light emitting diodes has been limited by the extremely low light extraction efficiency (<10%), due to the unique transverse magnetic (TM) polarized light emission. Here, we show that, by exploiting the lateral side emission, the extraction efficiency of TM polarized light can be significantly enhanced in AlGaN nanowire structures. Using the three-dimensional finite-difference time domain simulation, we demonstrate that the nanowire structures can be designed to inhibit the emission of guided modes and redirect trapped light into radiated modes. A light extraction efficiency of more than 70% can, in principle, be achieved by carefully optimizing the nanowire size, nanowire spacing, and p-GaN thickness.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗