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 649 records · Page 36

Directional finite difference method for directly solving 3D gyrokinetic field equations with enhanced accuracy

The gyrokinetic (GK) field equation is a three-dimensional (3D) elliptic equation, but it is often simplified to a set of two-dimensional (2D) equations by assuming that the field does not vary along a specific direction. However, this simplification can introduce inevitable 0th-order numerical errors, as nonlinear mode coupling in toroidal geometry can produce undesirable harmonic modes that violate the assumption. In this work, we propose a novel directional finite difference method (FDM) with a local coordinate transformation to better resolve the target field of interest. The directional FDM can accurately solve 3D GK field equations without simplifications, which can overcome the limitations of conventional methods. The accuracy and efficiency of different FDMs are analyzed in great detail for a variety of geometries, from simple 2D Cartesian coordinates to realistic 3D curvilinear coordinates. The 0th-order numerical errors of simplified 2D GK equations were found to be more problematic for low-harmonic modes and low aspect ratio geometries such as spherical tokamaks. On the other hand, the directional 3D FDM can accurately resolve a much wider range of harmonic modes aligned to the direction of interest, including the low-harmonic modes. In conclusion, we demonstrate that the directional 3D FDM is a highly effective algorithm for solving the 3D GK field equations, achieving accuracy improvements of 10 to 100 times or more, particularly for low-harmonic modes in spherical tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Coupled fcGCM-GCE Modeling System: A 3D Cloud Resolving Model and a Regional Scale Model

Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional single-column models in simulating various types of clouds and cloud systems from different geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a super-parameterization or multi-scale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and ore sophisticated physical parameterization. NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. The Goddard MMF is based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM), and it has started production runs with two years results (1998 and 1999). Also, at Goddard, we have implemented several Goddard microphysical schemes (21CE, several 31CE), Goddard radiation (including explicity calculated cloud optical properties), and Goddard Land Information (LIS, that includes the CLM and NOAH land surface models) into a next generation regional scale model, WRF. In this talk, I will present: (1) A Brief review on GCE model and its applications on precipitation processes (microphysical and land processes), (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), and preliminary results (the comparison with traditional GCMs), (3) A discussion on the Goddard WRF version (its developments and applications), and (4) The characteristics of the four-dimensional cloud data sets (or cloud library) stored at Goddard.

Tao, Wei-Kuo↗

A Coupled GCM-Cloud Resolving Modeling System, and A Regional Scale Model to Study Precipitation Processes

Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional single-column models in simulating various types of clouds and cloud systems from different geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a super-parameterization or multi-scale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. The Goddard MMF is based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM), and it has started production runs with two years results (1998 and 1999). Also, at Goddard, we have implemented several Goddard microphysical schemes (21CE, several 31CE), Goddard radiation (including explicitly calculated cloud optical properties), and Goddard Land Information (LIS, that includes the CLM and NOAH land surface models) into a next generation regional scale model, WRF. In this talk, I will present: (1) A brief review on GCE model and its applications on precipitation processes (microphysical and land processes), (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), and preliminary results (the comparison with traditional GCMs), and (3) A discussion on the Goddard WRF version (its developments and applications).

Tao, Wei-Kuo↗

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

The purpose of this paper is to achieve more versatile, convenient stability criteria for a wide class of finite-difference approximations to initial boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter-plane x greater than or equal to 0, t greater than or equal to 0. With these criteria, stability is easily established for a large number of examples, thus incorporating and generalizing many of the cases studied in recent literature.

Goldberg, M.↗

Numerical Analysis of Temperature Gradients and Interface Shape During Directional Solidification of Al and Al-Cu Alloy Under Microgravity Conditions

Numerical modeling was undertaken to analyze the influence of radial thermal gradient on solid/liquid (s/1) interface shape and convection patterns during solidification of pure Al and Al-4 wt% Cu alloy. The objective of the numerical task was to predict the influence of convective velocity on an insoluble particle near a s/l interface. These predictions would then be used to define the minimum gravity level (g) required to investigate the fundamental physics of interaction between a particle and a s/I interface. To satisfy this objective, steady state calculations were performed for different gravity levels and orientations with the gravity vector. ne furnace configuration used in this analysis is the proposed International Space Station Furnace, Quench Module Insert (QMI) 1. Results from a thermal model of the furnace core were used as initial boundary conditions for solidification modeling. General model of binary alloy solidification was based on the finite element code FIDAP. It was found that for the worst case orientation of 90 degrees with the gravity vector and a g level of 10(exp -4)g(sub o) (g(sub o) = 9.8 m/s(exp 2)) the dominant forces acting on the particle would be the fundamental drag and interfacial forces.

Bune, Andris V.↗

Building Thermal Models

This presentation is meant to be an overview of the model building process It is based on typical techniques (Monte Carlo Ray Tracing for radiation exchange, Lumped Parameter, Finite Difference for thermal solution) used by the aerospace industry This is not intended to be a "How to Use ThermalDesktop" course. It is intended to be a "How to Build Thermal Models" course and the techniques will be demonstrated using the capabilities of ThermalDesktop (TD). Other codes may or may not have similar capabilities. The General Model Building Process can be broken into four top level steps: 1. Build Model; 2. Check Model; 3. Execute Model; 4. Verify Results.

Modeling↗

Influences of Quench Cooling Rate on Microstructure and Corrosion Resistance of Al-Cu-Mg Alloy Based on the End-Quenching Test

To investigate the effects of the quench cooling rate on corrosion resistance of Al-Cu-Mg alloy, an end-quenching test was conducted and the microstructures at different cooling rates were observed by SEM and TEM. Additionally, the corrosion resistance was characterized by an intergranular corrosion test and electrochemical test. Moreover, the finite element method was applied to simulate the end quenching process. The results indicate that the actual end quenching process can be approximated as one-dimensional heat transfer, and the cooling rate varies at different cooling distances. By affecting the microstructures, decreasing the cooling rate leads to a decline in the corrosion properties. Low cooling rates coarsen the constituent particles and grain boundary particles, resulting in a wide precipitation-free zone and an increase in the intensity of corrosion reactions. A high cooling rate concentrates on the intragranular precipitant, which can reduce the pitting depth and represents a conversion from localized corrosion to general corrosion.

36 MATERIALS SCIENCE↗

Scalable algorithms for three-field mixed finite element coupled poromechanics

We introduce a class of block preconditioners for accelerating the iterative solution of coupled poromechanics equations based on a three-field formulation. The use of a displacement/velocity/pressure mixed finite-element method combined with a first order backward difference formula for the approximation of time derivatives produces a sequence of linear systems with a 3 x 3 unsymmetric and indefinite block matrix. The preconditioners are obtained by approximating the two-level Schur complement with the aid of physically-based arguments that can be also generalized in a purely algebraic approach. A theoretical and experimental analysis is presented that provides evidence of the robustness, efficiency and scalability of the proposed algorithm. In conclusion, the performance is also assessed for a real-world challenging consolidation experiment of a shallow formation.

42 ENGINEERING↗

Accelerated iterative calculation of transonic nacelle flowfields

A method is presented for the calculation of inviscid, supercritical flowfields about axisymmetric inlet cowls. A finite-difference calculation is performed in a simple, rectangular domain obtained from the nacelle geometry by a nearly-conformal mapping procedure. Type-dependent finite-differences are constructed using a coordinate-independent, 'rotated' differencing scheme. Methods of accelerating convergence of the iterative solution are demonstrated including a hybrid fast-Poisson-solver/relaxation scheme and an extrapolated relaxation procedure. Calculated pressure distributions are compared with experimental data for a variety of Mach numbers and mass-flow ratios, and show generally good agreement.

Caughey, D. A.↗

Computational design of the basic dynamical processes of the UCLA general circulation model

The 12-layer UCLA general circulation model encompassing troposphere and stratosphere (and superjacent 'sponge layer') is described. Prognostic variables are: surface pressure, horizontal velocity, temperature, water vapor and ozone in each layer, planetary boundary layer (PBL) depth, temperature, moisture and momentum discontinuities at PBL top, ground temperature and water storage, and mass of snow on ground. Selection of space finite-difference schemes for homogeneous incompressible flow, with/without a free surface, nonlinear two-dimensional nondivergent flow, enstrophy conserving schemes, momentum advection schemes, vertical and horizontal difference schemes, and time differencing schemes are discussed.

Arakawa, A.↗

An algorithm for maximum likelihood estimation using an efficient method for approximating sensitivities

An algorithm for maximum likelihood (ML) estimation is developed primarily for multivariable dynamic systems. The algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). The method determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort compared with integrating the analytically determined sensitivity equations or using a finite-difference method. Different surface-fitting methods are discussed and demonstrated. Aircraft estimation problems are solved by using both simulated and real-flight data to compare MNRES with commonly used methods; in these solutions MNRES is found to be equally accurate and substantially faster. MNRES eliminates the need to derive sensitivity equations, thus producing a more generally applicable algorithm.

Murphy, P. C.↗

Modeling and analysis of the space shuttle nose-gear tire with semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The Space Shuttle Orbiter nose gear tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynominals in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell. Numerical results of the Space Shuttle Orbiter nose gear tire model are compared with experimental measurements of the tire subjected to inflation loading.

Kim, Kyun O.↗

Analysis of aircraft tires via semianalytic finite elements

A computational procedure is presented for the geometrically nonlinear analysis of aircraft tires. The tire was modeled by using a two-dimensional laminated anisotropic shell theory with the effects of variation in material and geometric parameters included. The four key elements of the procedure are: (1) semianalytic finite elements in which the shell variables are represented by Fourier series in the circumferential direction and piecewise polynomials in the meridional direction; (2) a mixed formulation with the fundamental unknowns consisting of strain parameters, stress-resultant parameters, and generalized displacements; (3) multilevel operator splitting to effect successive simplifications, and to uncouple the equations associated with different Fourier harmonics; and (4) multilevel iterative procedures and reduction techniques to generate the response of the shell.

Noor, Ahmed K.↗

Impact of Corrections to the Spallings Volume Calculation on Waste Isolation Pilot Plant Performance Assessment - 16135

The numerical code DRSPALL (from direct release spallings) is written to calculate the volume of Waste Isolation Pilot Plant solid waste subject to material failure and transport to the surface (i.e., spallings) as a result of a hypothetical future inadvertent drilling intrusion into the repository. An error in the implementation of the DRSPALL finite difference equations was discovered and documented in a software problem report in accordance with the quality assurance procedure for software requirements. This paper describes the corrections to DRSPALL and documents the impact of the new spallings data from the modified DRSPALL on previous performance assessment calculations. Updated performance assessments result in more simulations with spallings, which generally translates to an increase in spallings releases to the accessible environment. Total normalized radionuclide releases using the modified DRSPALL data were determined by forming the summation of releases across each potential release pathway, namely borehole cuttings and cavings releases, spallings releases, direct brine releases, and transport releases. Because spallings releases are not a major contributor to the total releases, the updated performance assessment calculations of overall mean complementary cumulative distribution functions for total releases are virtually unchanged. Therefore, the corrections to the spallings volume calculation did not impact Waste Isolation Pilot Plant performance assessment calculation results. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Atomistically derived cohesive zone model of intergranular fracture in polycrystalline graphene

Pristine single crystal graphene is the strongest known two-dimensional material, and its nonlinear anisotropic mechanical properties are well understood from the atomic length scale up to a continuum description. However, experiments indicate that grain boundaries in the polycrystalline form reduce the mechanical behavior of polycrystalline graphene. Herein, we perform atomistic-scale molecular dynamics simulations of the deformation and fracture of graphene grain boundaries and express the results as continuum cohesive zone models (CZMs) that embed notions of the grain boundary ultimate strength and fracture toughness. To facilitate energy balance, we employ a new methodology that simulates a quasi-static controlled crack propagation which renders the kinetic energy contribution to the total energy negligible. We verify good agreement between Griffith's critical energy release rate and the work of separation of the CZM, and we note that the energy of crack edges and fracture toughness differs by about 35%, which is attributed to the phenomenon of bond trapping. This justifies the implementation of the CZM within the context of the finite element method (FEM). To enhance computational efficiency in the FEM implementation, we discuss the use of scaled traction-separation laws (TSLs) for larger element sizes. As a final result, we have established that the failure characteristics of pristine graphene and high tilt angle bicrystals differ by less than 10%. This result suggests that one could use a unique or a few typical TSLs as a good approximation for the CZMs associated with the mechanical simulations of the polycrystalline graphene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Collaboration with NSTX in Calculations of Radiofrequency and Neutral Beam Heating and Current Drive

The overall objective of this project was to advance the understanding of heating and current drive necessary for toroidal magnetic confinement fusion energy reactors, and to provide accurate models for these sources to aid in understanding magneto-hydrodynamic, transport and boundary physics issues. Increased understanding is provided by computer simulation using the underlying physics laws, and will lead to greater reliability and increased efficiency of a power producing fusion reactor. The proposed work entailed applications of existing CompX Fokker-Planck, ray tracing, full-wave microwave, and particle diffusion calculator codes to the National Spherical Torus Experiment at Princeton Plasma Physics Laboratory. Special features of NSTX include that large ion orbits are emphasized in its geometry and thus offer a major opportunity within fusion research for focused exploration of their ubiquitous effects, including comparison with simulation codes. The CompX CQL3D Fokker-Planck code was be further adapted to the finite-orbit-width physics. Comparison with experiment ensured that the computer codes have sufficient accuracy to accurately model the experiments. Nonthermal distributions of energetic ions are produced in experiments-- and modeled with the codes-- exhibiting the signatures of auxiliary heating, transport, and radiation processes within the plasma. Using physics-based modeling, in conjunction with experimental data, nonthermal distributions can provide keys to validation and quantification of auxiliary heating physics. Nonthermal distributions also hold particular information on radial transport processes, not obscured by the smoothing Maxwellization effects of Coulomb collisions. Generally, radial transport is not well understood but is a key to efficient fusion energy production. By comparison with Fokker-Planck modeling, the velocity dependence of radial transport coefficients becomes experimentally accessible, particularly in cases of well-modeled sources such as neutral beams, fast wave power, or the toroidal electric field. It is clear from modeling of fast-wave interactions with ions that finite-orbit-width effects play a significant role. For example, a different set of ions is computed to intersect a given cyclotron resonance depending on whether or not radial orbit drifts are taken into account. Also, radial transport, associated bootstrap current, and alpha channeling depend on radial drifts. Ion cyclotron heating above the first harmonic interacts preferentially with higher velocity ions, giving large radial drifts. The CQL3D Fokker-Planck code, maintained by the proposers, is a physics based model used extensively within the US fusion community to compare with and interpret experimental observations. At the beginning of this grant, bounce-averages of the rf quasilinear and other operators have been performed neglecting the radial drifts which give finite orbit widths to the particle orbits. The proposed work included a major upgrade to CQL3D to include finite-orbit-width effects be accounted for by performing the bounce-averages of rf and neutral beam fast ions along guiding center orbits. This CQL3D enhancement provided microwave and neutral beam contributions to neoclassical radial transport coefficients for nonthermal distributions, and substantially increased the comprehensiveness of the ion model. Additional work validated Fokker-Planck and ray tracing modeling of neutral beam and high harmonic fast wave absorption and current drive against the NSTX experiment, as had been very successful for electron cyclotron current drive in the DIII-D tokamak at General Atomics and lower hybrid waves on C-Mod at MIT. The validation work was in conjunction with staff and graduate students at Princeton Plasma Physics Laboratory. This DOE grant supported fusion energy research, a potential long-term solution to the world's energy needs. The research grant under consideration, DE-SC0006614, provided research support at the 0.3 Full-Time-Equivalent level for application of computer models to aid in understanding and projecting efficacy of heating and current drive sources in the National Spherical Torus Experiment, a tokamak variant, at the Princeton Plasma Physics Laboratory. The NSTX experiment explores the physics of very tight aspect ratio, almost spherical tokamaks, aiming at producing steady-state fusion plasmas. The current drive is an integral part of the steady-state concept, maintaining the magnetic geometry in the steady-state tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PolyPole-1: An accurate numerical algorithm for intra-granular fission gas release

This paper describes the development of a new numerical algorithm (called PolyPole-1) to efficiently solve the equation for intra-granular fission gas release in nuclear fuel. The work was carried out in collaboration with Politecnico di Milano and Institute for Transuranium Elements. The PolyPole-1 algorithms is being implemented in INL's fuels code BISON code as part of BISON's fission gas release model. The transport of fission gas from within the fuel grains to the grain boundaries (intra-granular fission gas release) is a fundamental controlling mechanism of fission gas release and gaseous swelling in nuclear fuel. Hence, accurate numerical solution of the corresponding mathematical problem needs to be included in fission gas behaviour models used in fuel performance codes. Under the assumption of equilibrium between trapping and resolution, the process can be described mathematically by a single diffusion equation for the gas atom concentration in a grain. In this work, we propose a new numerical algorithm (PolyPole-1) to efficiently solve the fission gas diffusion equation in time-varying conditions. The PolyPole-1 algorithm is based on the analytic modal solution of the diffusion equation for constant conditions, with the addition of polynomial corrective terms that embody the information on the deviation from constant conditions. The new algorithm is verified by comparing the results to a finite difference solution over a large number of randomly generated operation histories. Furthermore, comparison to state-of-the-art algorithms used in fuel performance codes demonstrates that the accuracy of the PolyPole-1 solution is superior to other algorithms, with similar computational effort. Finally, the concept of PolyPole-1 may be extended to the solution of the general problem of intra-granular fission gas diffusion during non-equilibrium trapping and resolution, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Three-dimensional rotational plasma flows near solid surfaces in an axial magnetic field

A rotational flow of a conducting viscous medium near an extended dielectric disk in a uniform axial magnetic field is analyzed in the magnetohydrodynamic (MHD) approach. An analytical solution to the system of nonlinear differential MHD equations of motion in the boundary layer for the general case of different rotation velocities of the disk and medium is obtained using a modified Slezkin–Targ method. A particular case of a medium rotating near a stationary disk imitating the end surface of a laboratory device is considered. The characteristics of a hydrodynamic flow near the disk surface are calculated within the model of a finite-thickness boundary layer. The influence of the magnetic field on the intensity of the secondary flow is studied. Calculations are performed for a weakly ionized dense plasma flow without allowance for the Hall effect and plasma compressibility. An MHD flow in a rotating cylinder bounded from above by a retarding cap is considered. The results obtained can be used to estimate the influence of the end surfaces on the main azimuthal flow, as well as the intensities of circulating flows in various devices with rotating plasmas, in particular, in plasma centrifuges and laboratory devices designed to study instabilities of rotating plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗