Energy Technology Data Exchange (ETDEWEB)

Fully coupled, Newton-Krylov algorithms are investigated for solving strongly coupled, nonlinear systems of partial differential equations arising in the field of computational fluid dynamics. Primitive variable forms of the steady incompressible and compressible Navier-Stokes and energy equations that describe the flow of a laminar Newtonian fluid in two-dimensions are specifically considered. Numerical solutions are obtained by first integrating over discrete finite volumes that compose the computational mesh. The resulting system of nonlinear algebraic equations are linearized using Newton`s method. Preconditioned Krylov subspace based iterative algorithms then solve these linear systems on each Newton iteration. Selected Krylov algorithms include the Arnoldi-based Generalized Minimal RESidual (GMRES) algorithm, and the Lanczos-based Conjugate Gradients Squared (CGS), Bi-CGSTAB, and Transpose-Free Quasi-Minimal Residual (TFQMR) algorithms. Both Incomplete ...

Energy Technology Data Exchange (ETDEWEB)

No abstract prepared.

Energy Technology Data Exchange (ETDEWEB)

The control volume finite element method (CVFEM) was developed to combine the local numerical conservation property of control volume methods with the unstructured grid and generality of finite element methods (FEMs). Most implementations of CVFEM include mass-lumping and upwinding techniques typical of control volume schemes. In this work we compare, via numerical error analysis, CVFEM and FEM utilizing consistent and lumped mass implementations, and stabilized Petrov-Galerkin streamline upwind schemes in the context of advection-diffusion processes. For this type of problem, we find no apparent advantage to the local numerical conservation aspect of CVFEM as compared to FEM. The stabilized schemes improve accuracy and degree of positivity on coarse grids, and also reduce iteration counts for advection-dominated problems.

British Library Electronic Table of Contents (United Kingdom)

In this paper, the unsteady MHD free convection heat and mass transfer of viscous fluid flowing through a Darcian porous regime adjacent to a moving vertical semi-infinite plate under Soret and Dufour effect have been examined. Viscous dissipation effects are included in the energy equation. A uniform magnetic field is applied transversely to the direction of the flow. The differential equations governing the problem have been transformed by a similarity transformation into a system of non-dimensional differential equations which are solved numerically by element free Galerkin method. The influence of Grashof number (Gr), magnetic parameter (M), heat absorption parameter (Q), permeability parameter (K), Schmidt number (Sc), Soret number (Sr), and Dufour number (Du) on the velocity, tempera...

CERN Document Server

We consider the numerical discretization of the time-domain Maxwell's equations with an energy-conserving discontinuous Galerkin finite element formulation. This particular formulation allows for higher order approximations of the electric and magnetic field. Special emphasis is placed on an efficient implementation which is achieved by taking advantage of recurrence properties and the tensor-product structure of the chosen shape functions. These recurrences have been derived symbolically with computer algebra methods reminiscent of the holonomic systems approach.

British Library Electronic Table of Contents (United Kingdom)

We present modeling, simulation, and characterization for the dynamic response of clamped-clamped microbeams under mechanical shock. A Galerkin-based reduced-order model is utilized and its results are verified by comparing to finite-element results. The results indicate that the response of a microbeam to mechanical shock is inherently non-linear because of the dominating effect of mid-plane stretching. The effect of the shock pulse shape is investigated. It is concluded that the shape of the shock pulse can result in significant dynamic amplification in the response of the microbeam even in cases where the shock load is considered quasi-static.The combined effect of the electrostatic force and mechanical shock is investigated. The results show that this combined effect can lead to early ...

Energy Technology Data Exchange (ETDEWEB)

The solution of the governing steady transport equations for momentum, heat and mass transfer in flowing fluids can be very difficult. These difficulties arise from the nonlinear, coupled, nonsymmetric nature of the system of algebraic equations that results from spatial discretization of the PDEs. In this manuscript the authors focus on evaluating a proposed nonlinear solution method based on an inexact Newton method with backtracking. In this context they use a particular spatial discretization based on a pressure stabilized Petrov-Galerkin finite element formulation of the low Mach number Navier-Stokes equations with heat and mass transport. The discussion considers computational efficiency, robustness and some implementation issues related to the proposed nonlinear solution scheme. Computational results are presented for several challenging CFD benchmark problems as well as two large scale 3D flow simulations.

International Nuclear Information System (INIS)

The study of natural convection flow and heat transfer within a cylindrical annulus has received considerable attention because of its numerous applications, such as in nuclear reactor design, electronic component cooling, thermal storage systems, energy conservation, energy storage, and energy transmission. Here, the effects of multiple geometric perturbations on the inner and outer cylinders of an annulus with impermeable end walls are investigated in this work. A three-dimensional study was done using a numerical scheme based on a Galerkin method of finite element formulation. The nature of the buoyancy-induced flow field has been analyzed in detail. The flow fields for the cases considered were found to be qualitatively similar, and the introduction of each additional perturbation altered the flow field in a regular and recurring manner. The introduction of each perturbation on the outer cylinder causes clockwise and counterclock-wise rotating patterns on ...

International Nuclear Information System (INIS)

Although B-spline techniques have been used to solve two-point boundary value problems with Dirac Hamiltonians for more than 20 years, the treatment of boundary conditions is still a matter of controversy. Spurious, non-physical, solutions are endemic when boundary conditions are not handled correctly. These pathological problems are absent when traditional finite difference methods are used as in computer packages such as GRASP. Accurate approximation using both finite differences and B-splines depends on controlling local approximation errors, and this common property suggests no a priori reason to suppose that B-spline algorithms should be more prone to generate spurious solutions. The relativistic Bloch operators of [24], when added to the Dirac differential operator, permit the construction of a self-adjoint differential operator for the two-point boundary value problem on a finite interval. Approximate solution of this problem exploiting the properties of B-splines in variational ...

Science.gov (United States)

In this study, adjoint sensitivity theory is developed for equations of two-dimensional steady-state flow in a confined aquifer. Both the primary flow equation and the adjoint sensitivity equation are solved using the Galerkin finite element method. The developed computer code is used to investigate the regional flow parameters of the Leadville Formation of the Paradox Basin in Utah and the Wolcamp carbonate/sandstone aquifer of the Palo Duro Basin in the Texas Panhandle. Two performance measures are evaluated, local heads and velocity in the vicinity of potential high-level nuclear waste repositories. The results illustrate the sensitivity of calculated local heads to the boundary conditions. Local velocity-related performance measures are more sensitive to hydraulic conductivities. The uncertainty in the performance measure is a function of the parameter sensitivity, parameter variance and the correlation between parameters. Given a parameter covariance matrix, ...

Science.gov (United States)

In the self-cooling blankets of the Tokamak fusion reactor, a liquid metal, namely liquid lithium, is pumped through a system of ducts to transfer heat and capture neutrons. One of the blanket designs proposed in Argonne National Laboratory's Blanket Comparison and Selection study uses a combination of poloidal and toroidal ducts in order to maximize heat transfer while minimizing net pressure drop. In the design, the poloidal and toroidal ducts meet at sharp, abrupt corners. They were modelled as two identical, straight, semi-infinite, thin-walled, rectangular ducts with 45{degree} miters and joined at a 90{degree} angle in the plane of a strong, uniform magnetic field. While in the toroidal containment vessel (i.e. the blanket), the liquid lithium is subjected to a large electromagnetic body force due to the presence of a strong magnetic field. This body force so dominates the flow as to make the inertial and viscous forces negligible everywhere, except in thin boundary or ...

Energy Technology Data Exchange (ETDEWEB)

Singlet oxygen generators are multiphase flow chemical reactors used to generate energetic oxygen to be used as a fuel for chemical oxygen iodine lasers. In this paper, a theoretical model of the generator is presented along with its solutions over ranges of parameter space and oxygen maximizing optimizations. The singlet oxygen generator (SOG) is a low-pressure, multiphase flow chemical reactor that is used to produce molecular oxygen in an electronically excited state, i.e. singlet delta oxygen. The primary product of the reactor, the energetic oxygen, is used in a stage immediately succeeding the SOG to dissociate and energize iodine. The gas mixture including the iodine is accelerated to a supersonic speed and lased. Thus the SOG is the fuel generator for the chemical oxygen iodine laser (COIL). The COIL has important application for both military purposes--it was developed by the US Air Force in the 1970s--and, as the infrared beam is readily absorbed by metals, industrial cutting ...