Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference discretization”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical investigation of the thermal stratification in cryogenic tanks subjected to wall heat flux

The flow pattern and thermal stratification of a cryogenic cylindrical tank are numerically studied. The tank sidewall is subjected to either a uniform heat-flux or two discrete levels of uniform heat-flux at the upper and lower halves of the tank wall. The tank bottom is kept at a constant temperature controlled by the heat exchanger of a thermodynamic vent system. The tank pressure is also assumed constant resulting in a constant saturation temperature at the interface which is higher than the tank bottom temperature. The effects of vapor motion and vapor superheat on the mass and heat transfer processes at the interface are assumed negligible such that the calculations of liquid region can be decoupled from the vapor region. Dimensionless steady-state conservation equations are solved by a finite-difference method. The effects of modified Rayleigh number, Prandtl number, tank aspect ratio, wall heat-flux parameter, and wall heat-flux distribution on the liquid velocity and temperature fields are investigated. Also, their effects on the rate of heat transfer through the interface and the tank bottom are examined.

Lin, Chin-Shun↗
NASA NTRS1995

Slot spiral antenna modeling using hybrid/mixed FE-BI technique

This report is concerned with the numerical simulation of the printed slot spiral antenna (SSA) using the finite element - boundary integral (FE-BI) analysis. It has been reported that the FE-BI technique is suited for modeling microstrip patch antennas of any shape, printed on layered planar structure or cylindrical platform and fed with a coaxial cable or a microstripline network underneath the radiating elements. However, direct application of the hybrid FE-BI technique to thin slot spirals requires excessive sampling rates to accurately simulate the geometry. To alleviate the meshing/modeling difficulties encountered with cavity-backed slot antennas, we describe a mixed finite element-boundary integral formulation. As in the past, the boundary integral is used to describe the radiation of the slot above the cavity. However the cavity is now modeled using a suitable mix of edge and node elements. The latter are used only at the aperture of the thin slot so that the nodes follow the center line of the slot. In this manner, regular size elements can be used regardless of the slot's width and any meshing restrictions are substantially relaxed. The proposed mixed element FE-BI formulation introduces three different computational regions and as expected this complicates the generation of the discrete system.

Gong, Jian↗
NASA NTRS1998

Numerical Modeling of Physical Vapor Transport in Contactless Crystal Growth Geometry

Growth from the vapor under conditions of limited contact with the walls of the growth ampoule is beneficial for the quality of the growing crystal due to reduced stress and contamination which may be caused by interactions with the growth container. The technique may be of a particular interest for studies on crystal growth under microgravity conditions: elimination of some factors affecting the crystal quality may make interpretation of space-conducted processes more conclusive and meaningful. For that reason, and as a part of our continuing studies on 'contactless' growth technique, we have developed a computational model of crystal growth process in such system. The theoretical model was built, and simulations were performed using the commercial computational fluid dynamics code, (CFD) ACE. The code uses an implicit finite volume formulation with a gray discrete ordinate method radiation model which accounts for the diffuse absorption and reflection of radiation throughout the furnace. The three-dimensional model computes the heat transfer through the crystal, quartz, and gas both inside and outside the ampoule, and mass transport from the source to the crystal and the sink. The heat transport mechanisms by conduction, natural convection, and radiation, and mass transport by diffusion and convection are modeled simultaneously and include the heat of the phase transition at the solid-vapor interfaces. As the thermal boundary condition, temperature profile along the walls of the furnace is used. For different thermal profiles and furnace and ampoule dimensions, the crystal growth rate and development of the crystal-vapor and source-vapor interfaces (change of the interface shape and location with time) are obtained. Super/under-saturation in the ampoule is determined and critical factors determining the 'contactless' growth conditions are identified and discussed. The relative importance of the ampoule dimensions and geometry, the furnace dimensions and its temperature, and the properties of the grown material are analyzed. The results of the simulations are compared with related experimental results on growth of CdTe, CdZnTe, ZnTe, PbTe, and PbSnTe crystals by this technique.

Palosz, W.↗
DOE OSTI2022

Impact of bulk-edge coupling on observation of anyonic braiding statistics in quantum Hall interferometers

Abstract Quantum Hall interferometers have been used to probe fractional charge and statistics of quasiparticles. We present measurements of a small Fabry–Perot interferometer in which the electrostatic coupling constants which affect interferometer behavior can be determined experimentally. Near the center of the ν = 1/3 state this device exhibits Aharonov–Bohm interference interrupted by a few discrete phase jumps, and Φ 0 oscillations at higher and lower magnetic fields, consistent with theoretical predictions for detection of anyonic statistics. We estimate the electrostatic parameters K I and K I L by two methods: using the ratio of oscillation periods in compressible versus incompressible regions, and from finite-bias conductance measurements. We find that the extracted K I and K I L can account for the deviation of the phase jumps from the theoretical anyonic phase θ a = 2 π /3. At integer states, we find that K I and K I L can account for the Aharonov–Bohm and Coulomb-dominated behavior of different edge states.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1996

A Full Multi-Grid Method for the Solution of the Cell Vertex Finite Volume Cauchy-Riemann Equations

The system of inhomogeneous Cauchy-Riemann equations defined on a square domain and subject to Dirichlet boundary conditions is considered. This problem is discretised by using the cell vertex finite volume method on quadrilateral meshes. The resulting algebraic problem is overdetermined and the solution is defined in a least squares sense. By this approach a consistent algebraic problem is obtained which differs from the original one by O(h(exp 2)) perturbations of the right-hand side. A suitable cell-based convergent smoothing iteration is presented which is naturally linked to the least squares formulation. Hence, a standard multi-grid algorithm is reported which combines the given smoother and cell-based transfer operators. Some remarkable reduction properties of these operators are shown. A full multi-grid method is discussed which solves the discrete problem to the level of truncation error by employing one multi-grid cycle at each current level of discretisation. Experiments and applications of the full multi-grid scheme are presented.

Borzi, A.↗
NASA NTRS1990

Low-frequency ionization-driven instability of a discrete auroral arc model

The low-frequency (time scales of tens of seconds) dynamics of the auroral ionosphere differs from that of the nonauroral ionosphere by the presence of strong, unstable space- and time-dependent ionospheric ionization produced by precipitating auroral electrons. If recombination is relatively unimportant (as at high ionospheric heights), it is shown that, in general, transport processes cannot remove this ionization as fast as it is created, and no equilibrium is possible. These nonequilibrium phenomena are investigated in the context of a nonlinear adiabatic auroral model, which has previously been studied in static situations. A linearized local perturbation analysis is given of what amounts to a current-driven E x B gradient-drift instability with an ionization source, as well as some exact nonlinear solutions valid in a finite but limited spatial region. These solutions show continuing motion of auroral potential and plasma density, as the aurora tries to shift its ionization problems from one place to another. The analysis gives clues to the possibility of generation of chaos and of fine-scale spatial structure.

Cornwall, John M.↗
NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗
NASA NTRS1995

Flow simulations about steady-complex and unsteady moving configurations using structured-overlapped and unstructured grids

The limiting factor in simulating flows past realistic configurations of interest has been the discretization of the physical domain on which the governing equations of fluid flow may be solved. In an attempt to circumvent this problem, many Computational Fluid Dynamic (CFD) methodologies that are based on different grid generation and domain decomposition techniques have been developed. However, due to the costs involved and expertise required, very few comparative studies between these methods have been performed. In the present work, the two CFD methodologies which show the most promise for treating complex three-dimensional configurations as well as unsteady moving boundary problems are evaluated. These are namely the structured-overlapped and the unstructured grid schemes. Both methods use a cell centered, finite volume, upwind approach. The structured-overlapped algorithm uses an approximately factored, alternating direction implicit scheme to perform the time integration, whereas, the unstructured algorithm uses an explicit Runge-Kutta method. To examine the accuracy, efficiency, and limitations of each scheme, they are applied to the same steady complex multicomponent configurations and unsteady moving boundary problems. The steady complex cases consist of computing the subsonic flow about a two-dimensional high-lift multielement airfoil and the transonic flow about a three-dimensional wing/pylon/finned store assembly. The unsteady moving boundary problems are a forced pitching oscillation of an airfoil in a transonic freestream and a two-dimensional, subsonic airfoil/store separation sequence. Accuracy was accessed through the comparison of computed and experimentally measured pressure coefficient data on several of the wing/pylon/finned store assembly's components and at numerous angles-of-attack for the pitching airfoil. From this study, it was found that both the structured-overlapped and the unstructured grid schemes yielded flow solutions of comparable accuracy for these simulations. This study also indicated that, overall, the structured-overlapped scheme was slightly more CPU efficient than the unstructured approach.

Newman, James C., III↗
NASA NTRS1995

Shape optimization of single- and two-element airfoils on multiblock grids

A multiblock, discrete sensitivity analysis method is used to couple a direct optimization method and a flow analysis method. The domain is divided into smaller subdomains for which the sensitivities are obtained separately. Then, an effective sensitivity equation is solved to complete the coupling of all the sensitivity information. The flow analysis is based on the thin-layer Navier-Stokes equations solved by an implicit, upwind-biased, finite-volume method. The method of feasible directions is used for the present gradient-based optimization approach. First, a transonic airfoil is optimized to investigate the behavior of the method in highly nonlinear flows as well as the effect of different blocking strategies on the procedure. A supercritical airfoil is produced from an initially symmetric airfoil with multiblocking affecting the path but not the final shape. Secondly, a two-element airfoil is shape optimized in subsonic flow to demonstrate the present method's capability of shaping aerodynamically interfering elements simultaneously. For a very low and a very high Reynolds number cases, the shape of the main airfoil and the flap are optimized to yield improved lift-to-drag ratios.

Lacasse, James M.↗
DOE OSTI2021

Multi-scale Coupled Processes Modeling of Fractures as Porous, Interfacial and Granular Systems from Rock Images with the Numerical Manifold Method

Abstract The greatest challenges of rigorously modeling coupled hydro-mechanical processes in fractured rocks at different scales are associated with computational geometry. In addition, selections of continuous or discontinuous models, physical laws, and coupling priorities at different scales based on different geometric features determine the applicability of a numerical model for a certain type of problem. In this study, we present our multi-scale modeling capabilities that have been developed based on the numerical manifold method for analyzing coupled hydro-mechanical processes in fractured rocks. Based on their geometric features, the fractures are modeled as continua—finite-thickness porous zones, and discontinua—discontinuous interfaces and microscale asperities and granular systems. Different governing equations, physical laws, coupling priorities, and approaches for addressing fracture intersections and shearing are then applied to describe these. We applied these models to simulate coupled processes in fractured rocks using realistic geometry obtained from rock images at different scales. We first calculated shearing of a single fracture with different models and demonstrated the impacts of asperities on shearing. We then applied the continuous and discontinuous models to simulate a network of rough fractures, demonstrating that contact dynamics contribute significantly to the geometric, multi-physical evolution of systems where rough fractures are not mineral filled. For a discrete fracture network, our coupled processes modeling demonstrates that shearing of the discrete fractures can have a major impact on stress and pore pressure distribution. Lastly, we applied the discontinuous granular model to simulate evolution of a complex granular system with a deformation band, demonstrating that the deformation band can dominate contact dynamics, the structural and the stress evolution of the granular system.

58 GEOSCIENCES↗
NASA NTRS1974

Flutter of panels on discrete flexible supports

The supersonic flutter of wide panels on discrete flexible supports is investigated for three different panel-support configurations. The first study examines the effect of support stiffness on the flutter behavior of a specific five-support configuration with leading-and trailing-edge overhangs; this configuration was recently under consideration for a heat shield on an entry body. The second study investigates the effect of support stiffness on flutter of panels with various numbers of equally spaced supports. The third study examines the effect of center-support location on the flutter of panels with three supports. Results are presented in nondimensional form. The analysis is based on wide-plate structural theory and Ackeret aerodynamics. Finite differences are employed to obtain solutions for flutter pressure. A computer program based on this analysis and including a direct solution technique is presented. The program can be used to find the flutter pressure of wide panels of variable thickness supported by any number of flexible supports.

Stein, M.↗
NASA NTRS1992

User's manual for three dimensional FDTD version D code for scattering from frequency-dependent dielectric and magnetic materials

The Penn State Finite Difference Time Domain Electromagnetic Scattering Code version D is a 3-D numerical electromagnetic scattering code based upon the finite difference time domain technique (FDTD). The manual provides a description of the code and corresponding results for several scattering problems. The manual is organized into 14 sections: introduction; description of the FDTD method; operation; resource requirements; version D code capabilities; a brief description of the default scattering geometry; a brief description of each subroutine; a description of the include file; a section briefly discussing Radar Cross Section computations; a section discussing some scattering results; a sample problem setup section; a new problem checklist; references and figure titles. The FDTD technique models transient electromagnetic scattering and interactions with objects of arbitrary shape and/or material composition. In the FDTD method, Maxwell's curl equations are discretized in time-space and all derivatives (temporal and spatial) are approximated by central differences.

Beggs, John H.↗
NASA NTRS1993

Application of multigrid methods to the solution of liquid crystal equations on a SIMD computer

We will describe a finite difference code for computing the equilibrium configurations of the order-parameter tensor field for nematic liquid crystals in rectangular regions by minimization of the Landau-de Gennes Free Energy functional. The implementation of the free energy functional described here includes magnetic fields, quadratic gradient terms, and scalar bulk terms through the fourth order. Boundary conditions include the effects of strong surface anchoring. The target architectures for our implementation are SIMD machines, with interconnection networks which can be configured as 2 or 3 dimensional grids, such as the Wavetracer DTC. We also discuss the relative efficiency of a number of iterative methods for the solution of the linear systems arising from this discretization on such architectures.

Farrell, Paul A.↗
NASA NTRS1985

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗
NASA NTRS1994

Accurate finite difference methods for time-harmonic wave propagation

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

Harari, Isaac↗
NASA NTRS1990

A computational procedure for the dynamics of flexible beams within multibody systems

The dynamic analysis of three dimensional elastic beams which experience large rotational and large deformational motions are examined. The beam motion is modeled using an inertial reference for the translational displacements and a body-fixed reference for the rotational quantities. Finite strain rod theories are then defined in conjunction with the beam kinematic description which accounts for the effects of stretching, bending, torsion, and transverse shear deformations. A convected coordinate representation of the Cauchy stress tensor and a conjugate strain definition is introduced to model the beam deformation. To treat the beam dynamics, a two-stage modification of the central difference algorithm is presented to integrate the translational coordinates and the angular velocity vector. The angular orientation is then obtained from the application of an implicit integration algorithm to the Euler parameter/angular velocity kinematical relation. The combined developments of the objective internal force computation with the dynamic solution procedures result in the computational preservation of total energy for undamped systems. The present methodology is also extended to model the dynamics of deployment/retrieval of the flexible members. A moving spatial grid corresponding to the configuration of a deployed rigid beam is employed as a reference for the dynamic variables. A transient integration scheme which accurately accounts for the deforming spatial grid is derived from a space-time finite element discretization of a Hamiltonian variational statement. The computational results of this general deforming finite element beam formulation are compared to reported results for a planar inverse-spaghetti problem.

Downer, Janice Diane↗
DOE OSTI2025

A Simple, Scalable Large Deformation Solid Mechanics Implementation in the MOOSE Framework

This article describes a large deformation solid mechanics solver implemented as part of the freely available and open source MOOSE finite element simulation framework. The article documents the choices made in developing the solid mechanics framework and describes novel formulations for the gradient operator and constitutive modeling framework made to simplify implementations of different coordinate systems, stabilized gradient operators, and different constitutive model inputs and outputs. In the process, the article describes a new formulation that casts objective integration of the Cauchy stress as a linear transformation of the small stress rate. Finally, the article presents key implementation details and examines the parallel efficiency of the solid mechanics solver implemented in MOOSE. The implementation retains a good weak scaling efficiency beyond 1,000 parallel processes. The article includes a discussion of the factors limiting the parallel efficiency of implicit, large deformation solid mechanics codes on current high-performance computers, with the main current limitation being the scalability of the algebraic multigrid methods used to solve the linearized equilibrium equations.

Applied computing → Computer-aided design↗