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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1981

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗

NASA NTRS1993

The NASA/GISS Mars general circulation model: Preliminary experiments

The NASA/GISS Mars General Circulation Model (GCM) is an adapted version of the GISS Global Climate/Middle Atmosphere Model, specifically developed for the diagnostic validation and objective analysis of measured atmospheric temperatures from the Mars Observer Pressure Modulator Infrared Radiometer (PMIRR) experiment. The GISS Mars GCM has 23 vertical layers extending from the surface to approximately 80 km altitude, representing a vertical resolution of about 0.3 scale heights. The primitive (vertically hydrostatic) equations are solved in finite difference form on the Krakawa B grid, with a horizontal resolution of 8 deg x 10 deg (latitude-longitude). The model includes a diurnal solar cycle, heat transport within a two-layer ground, and a high-order 'slopes-scheme' for the advection of heat in the upper atmosphere. The radiative transfer scheme is based on the correlated k distribution method for the treatment of nongray gaseous absorption thermal emission, and multiple scattering, including options for suspended dust. A special feature of the model of particular importance for Mars is a parameterization of gravity-wave-induced drag incorporating orographic forcing, wind shear, convection, and radiative damping. The implementation of the GISS Mars model includes global maps of topography, roughness, and albedo.

Allison, Michael↗

NASA NTRS2000

Model Study of the Sources of Upper Tropospheric Ozone Over The Atlantic, with Emphasis on Interhemispheric Asymmetries

Measurements of ozone and aerosols from ships transecting the Atlantic Ocean reveal considerable interhemispheric asymmetries in the vertical structure of the ozone distribution. There is much more upper tropospheric ozone to the South of the Equator than in the northern hemisphere. These profiles suggest that the air in the southern hemispheric tropics and subtropics has a source in the stratosphere, while that in the northern hemisphere originates in the boundary layer. A model experiment designed to investigate this problem will be discussed. They are based on a new general circulation model, based on a Finite-Volume dynamical core and the physical parametrizations from the Community Climate Model at the National Center for Atmospheric Research (NCAR). The model simulates a climatological state in reasonable agreement with observations, comparable to other contemporary climate models, but has a much more accurate transport capability than most other present-day models. The experiment involves using a suite of tracers, designed to highlight air from different source regions. Some of these (water vapor and ozone) represent real atmospheric quantities, while others are markers which indicate whether the air originates in the stratosphere, the boundary layer, or has traversed regions of convection (for example). The experiments were run for periods of up to one month, initialized from snapshots from a long climate integration of the model. A sample of results will be presented, chosen to illustrate the different sources of air in the upper troposphere and the physical mechanisms by which these air parcels arrived there.

Pawson, Steven↗

DOE OSTI2016

Thermo-mechanical simulations of CO{sub 2} laser–fused silica interactions

CO{sub 2} laser heating of silica glass is used in many scientific and industrial applications. Particularly, localized CO{sub 2} laser heating of silica glass has demonstrated its ability to mitigate surface damage on optics used for high power laser applications. To develop such applications, the control of temperature, heat affected area, and resulting mechanical stresses are critical. Therefore, it is necessary to understand the silica transformation, the material ejection, and the thermo-mechanical stresses induced by the laser heating and subsequent cooling. In this paper, we detail the development of comprehensive thermo-mechanical numerical simulations of these physical processes, based on finite-element method. The approach is developed for 2D or 3D cases to tackle the case of a moving beam at the surface of the sample, and we particularly discuss the choice of the different parameters based on bibliographic inputs. The thermal and mechanical numerical results have been compared to different dedicated experimental studies: infrared thermography measurements at the surface of the irradiated area, optical profilometry measurements of the laser-processed sites, and photo-elastic measurements. Very consistent results are obtained between numerical and experimental results for the description of the temperature gradients, the material ejection, and the residual stresses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

A concurrent irradiation-mechanics multiscale coupling model

The vast majority of our current knowledge regarding the basic mechanisms controlling irradiation effect on mechanical properties is based almost entirely on the results of post-irradiation experiments or theoretical models. However, the concurrent effects of irradiation, mechanical stress, and thermal damage on the failure phenomena of materials and components remain largely unexplored due to its internal multiscale-multiphysics coupling nature. Here, we present here a concurrent irradiation-mechanics multiscale coupling model. The concurrent evolutions of nanoscale irradiation defect clusters, microscale dislocation configurations, and mechanical responses are well captured based on coupling cluster dynamics, discrete dislocation dynamics, and the finite element methods using an effective time marching scheme. Model predictions of defect densities and size are in general agreement with experimental observations. Irradiation hardening is shown to take place also in samples undergoing concurrent irradiation-mechanical loading, similar to samples tested post-irradiation. However, the occurrence of plastic flow localization and dislocation channel formation is not accompanied with apparent yield drop (softening) under concurrent irradiation-mechanical loading conditions, which is different from the post-irradiation case.

42 ENGINEERING↗

DOE OSTI2024

Extracting Topological Orders of Generalized Pauli Stabilizer Codes in Two Dimensions

In this paper, we introduce an algorithm for extracting topological data from translation invariant generalized Pauli stabilizer codes in two-dimensional systems, focusing on the analysis of anyon excitations and string operators. The algorithm applies to Z d qudits, including instances where d is a nonprime number. This capability allows the identification of topological orders that differ from the Z d toric codes. It extends our understanding beyond the established theorem that Pauli stabilizer codes for Z p qudits (with p being a prime) are equivalent to finite copies of Z p toric codes and trivial stabilizers. The algorithm is designed to determine all anyons and their string operators, enabling the computation of their fusion rules, topological spins, and braiding statistics. The method converts the identification of topological orders into computational tasks, including Gaussian elimination, the Hermite normal form, and the Smith normal form of truncated Laurent polynomials. Furthermore, the algorithm provides a systematic approach for studying quantum error-correcting codes. We apply it to various codes, such as self-dual CSS quantum codes modified from the two-dimensional honeycomb color code and non-CSS quantum codes that contain the double semion topological order or the six-semion topological order. Published by the American Physical Society 2024

Physics↗

NASA NTRS1984

Variational objective analyses for cyclone studies

The basic analysis equations, i.e., the two horizontal momentum equations, the hydrostatic equation, and the integrated continuity equation were derived for the nonlinear vertical coordinate, nondimensionalized, and expressed in finite differences on a staggered grid. Special care was taken to transform the hydrostatic equation and the pressure gradient terms of the horizontal momentum equations to nearly eliminate truncation error over steeply sloping terrain. This formulation also eliminated explicit reference to orographically induced variations in the thermodynamic variables so that the variational adjustments are on the scale of the meteorological perturbations. The analysis equations were subjected to the Euler-Lagrange operations as expressed for finite differences and an additional set of five partial differential equations was derived, bringing to nine the number of equations in Model I. Higher order terms, terms containing observed quantities, and terms containing none of the variables to be adjusted were grouped into forcing functions and the equations were solved for the zero order terms. Zero order variables were eliminated between these equations and there resulted two diagnostic equations which take the form of general linear second order partial differential equations with nonconstant coefficients.

Achtemeier, G. L.↗

DOE OSTI2016

Energy/dissipation-preserving Birkhoffian multi-symplectic methods for Maxwell's equations with dissipation terms

In this study, we propose two new energy/dissipation-preserving Birkhoffian multi-symplectic methods (Birkhoffian and Birkhoffian box) for Maxwell's equations with dissipation terms. After investigating the non-autonomous and autonomous Birkhoffian formalism for Maxwell's equations with dissipation terms, we first apply a novel generating functional theory to the non-autonomous Birkhoffian formalism to propose our Birkhoffian scheme, and then implement a central box method to the autonomous Birkhoffian formalism to derive the Birkhoffian box scheme. We have obtained four formal local conservation laws and three formal energy global conservation laws. We have also proved that both of our derived schemes preserve the discrete version of the global/local conservation laws. Furthermore, the stability, dissipation and dispersion relations are also investigated for the schemes. Theoretical analysis shows that the schemes are unconditionally stable, dissipation-preserving for Maxwell's equations in a perfectly matched layer (PML) medium and have second order accuracy in both time and space. Numerical experiments for problems with exact theoretical results are given to demonstrate that the Birkhoffian multi-symplectic schemes are much more accurate in preserving energy than both the exponential finite-difference time-domain (FDTD) method and traditional Hamiltonian scheme. Finally, we also solve the electromagnetic pulse (EMP) propagation problem and the numerical results show that the Birkhoffian scheme recovers the magnitude of the current source and reaction history very well even after long time propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Impact of built-in fields and contact configuration on the characteristics of ultra-thin GaAs solar cells

We discuss the effects of built-in fields and contact configuration on the photovoltaic characteristics of ultra-thin GaAs solar cells. The investigation is based on advanced quantum-kinetic simulations reaching beyond the standard semi-classical bulk picture concerning the consideration of charge carrier states and dynamics in complex potential profiles. The thickness dependence of dark and photocurrent in the ultra-scaled regime is related to the corresponding variation of both, the built-in electric fields and associated modification of the density of states, and the optical intensity in the films. Losses in open-circuit voltage and short-circuit current due to the leakage of electronically and optically injected carriers at minority carrier contacts are investigated for different contact configurations including electron and hole blocking barrier layers. The microscopic picture of leakage currents is connected to the effect of finite surface recombination velocities in the semi-classical description, and the impact of these non-classical contact regions on carrier generation and extraction is analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

A finite element program for postbuckling calculations (PSTBKL)

The object of the research reported herein was to develop a general mathematical model and solution methodologies for analyzing the structural response of thin, metallic shell structures under large transient, cyclic, or static thermochemical loads. This report describes the computer program resulting from the research. Among the system responses associated with these loads and conditions are thermal buckling, creep buckling, and ratcheting. Thus geometric and material nonlinearities (of high order) have been anticipated and are considered in developing the mathematical model. The methodology is demonstrated through different problems of extension, shear, and of planar curved beams. Moreover, importance of the inclusion of large strains is clearly demonstrated, through the chosen applications.

Simitses, G. T.↗

DOE OSTI2023

Shedding light on shadow generalized parton distributions

The feasibility of extracting generalized parton distributions (GPDs) from deeply virtual Compton scattering (DVCS) data has recently been questioned because of the existence of an infinite set of so-called “shadow GPDs” (SGPDs). These SGPDs depend on the process and manifest as multiple solutions (at a fixed scale 𝑄 2 ) to the inverse problem that needs to be solved to infer GPDs from DVCS data. SGPDs therefore pose a significant challenge for extracting GPDs from DVCS data. With this motivation we study the extent to which QCD evolution can provide constraints on SGPDs. This is possible because the known classes of SGPDs begin to contribute to observables after evolution, and can then be constrained (at the input scale 𝑄$^{2}_{0}$) by data that has a finite 𝑄 2 range. The impact that SGPDs could have on determining the total angular momentum, pressure and sheer force distributions, and tomography is also discussed. We find that 𝑄 2 evolution, given empirical data over a wide range of skewness (𝜉 ) and 𝑄 2 scales, constrains the SGPDs studied at low 𝑥 and low 𝜉. If this finding generalizes to all SGPDs—and we have not been able to find any SGPDs with qualitatively different behavior—then this potentially makes extraction of GPDs from DVCS possible over a limited range of 𝑥 and 𝜉.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1996

Multigrid Methods for EHL Problems

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of the application of such methods is the recent Effective Influence Method of Dowson and Wang. Multigrid methods have also been used with great success by Venner and Venner and Lubrecht with a good summary being given by Venner. As both these finite difference discretization based approaches appear to provide an efficient way of solving EHL problems, it is important to understand their relative merits. This paper is a first attempt at providing such an understanding in the context of EHL point contact problem, (contact of two spheres), in which the contact zone is a point and an ellipse or circle for unloaded and loaded dry contacts respectively. Since the film thickness and the contact width are generally small compared to the local radius of curvature of the two surfaces, the reduced geometry of the surfaces in the contact area can be accurately approximated to the contact between a paraboloid and a flat surface. The layout of the remainder of this paper is as follows. In section 2 we introduce the form of the equations to be solved. The Effective Influence Newton Method is described in Section 3 while Section 4 describes the Multigrid method to be used. Sections 5 and 6 describe the test problems to be used in the comparison between the two methods and compare the performance of the two methods. Section 7 concludes the paper with an argument of the two methods and suggests some future research directions.

Nurgat, Elyas↗

NASA NTRS2001

Processing and Testing of Thermoplastic Composite Cylindrical Shells Fabricated by Automated Fiber Placement

Two 61-cm-diameter eight-ply quasi-isotropic IM7/PEEK cylindrical shells were fabricated by automated fiber placement the NASA Langley Research Center using only infrared radiant heat to preheat the substrate and incoming composite uni-tape. The shells were characterized by ultrasonic c-scans for overall consolidation quality, and by optical microscopy and acid digestion for void content. Compression tests were also performed. Although the material used in the study was of generally poor quality due to numerous splits and dry fiber regions, the process was able to achieve a net reduction in void content in the as-placed component. Microscopy of the composite shells revealed well-consolidated, void-free interfaces. The two cylinders were then tested in uni-axial compression in a 1334 kN-capacity hydraulic test machine until buckling occurred. A geometrically nonlinear finite element analysis was conducted, and the differences between the predicted and measured values were 18.0 and 25.8%, respectively. Inclusion of measured imperfections of the cylinder into the analysis is expected to reduce these differences.

Hulcher, Anthony Bruce↗

NASA NTRS1979

Generalization of the Manley-Rowe relations for non-periodic signals constrained in bandwidth with finite energy

General relations for the energy spectral densities of nonperiodic signals constrained in bandwidth in nonlinear reactances are derived. It is assumed that the charge-voltage characteristics of the reactances are polynomials. In addition certain conditions for the frequency bands of the signals must be met. By contrast with the Manley-Rowe relations integrals with respect to frequency appear instead of the ratios power over frequency. The integrands are energy spectral densities divided by frequency. For parametric devices inequalities can be derived for the ratios of the energy levels in the different circuits and the limiting frequencies of the energy spectral densities. With these inequalities it is possible to determine limits for the energy levels.

Proksch, P.↗

DOE OSTI2019

Role of structural anisotropy in geometric birefringence of high-birefringence index-guiding PCFs: analysis by moment of inertia

We have studied the role of structural anisotropy on the geometric birefringence of four typical high-birefringence (HB) index-guiding photonic crystal fibers (PCFs) with the same air-silica structure and circular, ellipse, and rectangle air holes, respectively. The normalized difference between the moment of inertia ΔI in x and y direction of the PCF structure is utilized to quantificationally describe the structural anisotropy. We demonstrate that the phase birefringence B increases monotonously when the normalized ΔI rises in different cases of PCFs structures. In order to obtain high birefringence B, it can be an effective method to optimize the parameters of the structure to increase the normalized ΔI. This work can be not only a reference of founding the essence of geometric birefringence and proposing the optimal design of HB-IG-PCFs, but also an enlightenment for bringing the different method (compared to the finite element) into the investigation of PCFs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1979

Failure analysis of composite laminates with free edges

Linear elastic stress distributions obtained from a refined finite element mesh are used in conjunction with the tensor polynomial failure criterion to predict the initiation of failure in symmetric, finite-width graphite-epoxy laminates under tensile loading. Results are presented for a wide variety of laminates including: (+ and - theta)s angle-ply; cross-ply (0/90)s and (90/0)s; and quasi-isotropic (90/0/+ and - 45)s and (+ and -45/0/90)s. It is shown that the elastic stress distributions generally compare favorably with other published results, but also indicate improved satisfaction of the stress-free boundary conditions and indicate some differences in the singular behavior of selected stress components at the free edge. The tensor polynomial failure criterion is used to predict the location and mode of first failure. Examination of the individual terms of the polynomial indicates different modes of failure depending upon the laminate configuration.

Herakovich, C. T.↗

DOE OSTI2016

Study on Modeling of Nuclear Fuel End Drop during a Hypothetical Accident

It is well known that shell elements have been used to model the cladding of single nuclear fuel pins. These have been used, for example, to analyze the response of a fuel assembly by assuming that all fuel pins in a fuel assembly deform in the same pattern. It is increasingly significant to model the entire fuel assembly in order to capture more realistically the behavior of each individual nuclear fuel during a hypothetical accident in a manner that can account for interactions among single pins. However, modeling each individual fuel pin in an entire fuel assembly using shell elements will result in models with large amount of elements and degrees of freedom that take correspondingly larger amount of computational times to execute. Therefore, modeling the nuclear fuel assembly with beam elements may become an option. This paper presents the results of a set of analysis with the LS-DYNA finite element software that investigated the differences and similarities of modeling nuclear fuel assemblies using either shell elements or beam elements. This study provides additional insights on the appropriateness of using beam element models of the fuel pins in the simulation of a hypothetical accident. In this study, two single pin models, a single pin model using shell elements and a single beam model using beam elements, are constructed to analyze end drop scenarios with the LS-DYNA code. Each model is then analyzed for a 19 ft end drop scenario (linear behavior) and 30 ft end drop scenario (non-linear behavior). Comparisons are made between the results obtained with the single pin shell model and with the single pin beam model to study the differences and similarities resulting from the two modeling approaches. In this study, a B and W 15x15 PWR fuel is analyzed for both 30 ft and 19 ft end hypothetical accident drops using both single pin finite element model with shell elements and beam elements for the fuel cladding. It is found that for the 19 ft end drop, the difference in the Max/Min total axial strains is about 12% between the two models. For the 30 ft end drop, the maximum total axial strains calculated with both models is similar, and a difference of about 12% is observed for the minimum total axial strains calculated by the two models. Advantages of the single pin shell model are that it permits consideration of internal pressures and may capture local buckling. Therefore, use of beam elements to model the fuel cladding needs to take into account the possible significance of local buckling and internal pressurization for the problem being analyzed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2018

Coupled complex Ginzburg–Landau systems with saturable nonlinearity and asymmetric cross-phase modulation

We formulate and study dynamics from a complex Ginzburg–Landau system with saturable nonlinearity, including asymmetric cross-phase modulation (XPM) parameters. Such equations can model phenomena described by complex Ginzburg–Landau systems under the added assumption of saturable media. When the saturation parameter is set to zero, we recover a general complex cubic Ginzburg–Landau system with XPM. We first derive conditions for the existence of bounded dynamics, approximating the absorbing set for solutions. We use this to then determine conditions for amplitude death of a single wavefunction. We also construct exact plane wave solutions, and determine conditions for their modulational instability. In a degenerate limit where dispersion and nonlinearity balance, we reduce our system to a saturable nonlinear Schrödinger system with XPM parameters, and we demonstrate the existence and behavior of spatially heterogeneous stationary solutions in this limit. Using numerical simulations we verify the aforementioned analytical results, while also demonstrating other interesting emergent features of the dynamics, such as spatiotemporal chaos in the presence of modulational instability. In other regimes, coherent patterns including uniform states or banded structures arise, corresponding to certain stable stationary states. For sufficiently large yet equal XPM parameters, we observe a segregation of wavefunctions into different regions of the spatial domain, while when XPM parameters are large and take different values, one wavefunction may decay to zero in finite time over the spatial domain (in agreement with the amplitude death predicted analytically). We also find a collection of transient features, including transient defects and what appear to be rogue waves, while in two spatial dimensions we observe highly localized pattern formation. While saturation will often regularize the dynamics, such transient dynamics can still be observed – and in some cases even prolonged – as the saturability of the media is increased, as the saturation may act to slow the timescale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗