SEARCH · Engineering Papers

Results for “Elliptic interface problems”

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.

18 records

All sources NASA NTRS DOE OSTI

NASA NTRS1993

On the shear stabilization of capillary break-up of finite liquid bridges

In this paper we consider an isothermal finite liquid bridge under zero-gravity. A sinusoidal interfacial shear stress drives a nonparallel flow in the bridge. The linear stability of this flow to three-dimensional disturbances (which may deform the cylindrical gas-liquid interface) is determined numerically by solving an elliptic eigenvalue problem. Previous results on shear stabilization of capillary break-up of axially unbounded cylindrical interfaces containing a parallel flow are hereby extended to a nonparallel flow. The turning flow regions influence the stability significantly. However, for the particular cases considered, a small area in parameter space remains where the capillary instability is suppressed through interfacial shear. Second, non-axisymmetric oscillatory instabilities are found which originate from an interaction of the interface deformation and the basic flow. These instabilities may be the isothermal limit of the oscillatory instabilities observed in float-zone crystal growth.

Dijkstra, Henk A.↗

NASA NTRS1987

Isoparametric finite element methods to estimate discontinuous coefficients in two-dimensional elliptic equations

A summary of findings is presented on the problem of estimating discontinuous coefficients, including the locations of interfaces at which discontinuities occur, for elliptic equations on a two-dimensional domain. Finite-element approximations are constructed that are based on an underlying grid that adaptively moves to follow locations of interfaces as they are iteratively updated in the numerical estimation process.

Lamm, Patricia K.↗

NASA NTRS1996

A Multigrid Algorithm for Immersed Interface Problems

Many physical problems involve interior interfaces across which the coefficients in the problem, the solution, its derivatives, the flux, or the source term may have jumps. These interior interfaces may or may not align with a underlying Cartesian grid. Zhilin Li, in his dissertation, showed how to discretize such elliptic problems using only a Cartesian grid and the known jump conditions to second order accuracy. In this paper, we describe how to apply the full multigrid algorithm in this context. In particular, the restriction, interpolation, and coarse grid problem will be described. Numerical results for several model problems are given to demonstrate that good rates can be obtained even when jumps in the coefficients are large and do not align with the grid.

Adams, Loyce↗

NASA NTRS1992

Boundary-conforming mapping applied to computations of highly deformed solidification interfaces

A new boundary-conforming mapping is developed for the calculation of highly deformed cellular solidification interfaces in a model of directional solidification of a binary alloy. The mapping is derived through a variational formulation that is designed so that the grid penetrates the grooves between cells along the interface without causing a loss of ellipticity of the mapping equations. A finite element/Newton method is presented for simultaneous solution of the free boundary problem described by the solutal model of directional solidification and the mapping equations. Results are compared to previous calculations and demonstrate the importance of accurate representation of the interface shape for understanding the solution structure.

Tsiveriotis, K.↗

NASA NTRS1993

A comparison of locally adaptive multigrid methods: LDC, FAC and FIC

This study is devoted to a comparative analysis of three 'Adaptive ZOOM' (ZOom Overlapping Multi-level) methods based on similar concepts of hierarchical multigrid local refinement: LDC (Local Defect Correction), FAC (Fast Adaptive Composite), and FIC (Flux Interface Correction)--which we proposed recently. These methods are tested on two examples of a bidimensional elliptic problem. We compare, for V-cycle procedures, the asymptotic evolution of the global error evaluated by discrete norms, the corresponding local errors, and the convergence rates of these algorithms.

Khadra, Khodor↗

NASA NTRS1996

First-Order System Least-Squares for Second-Order Elliptic Problems with Discontinuous Coefficients

The first-order system least-squares methodology represents an alternative to standard mixed finite element methods. Among its advantages is the fact that the finite element spaces approximating the pressure and flux variables are not restricted by the inf-sup condition and that the least-squares functional itself serves as an appropriate error measure. This paper studies the first-order system least-squares approach for scalar second-order elliptic boundary value problems with discontinuous coefficients. Ellipticity of an appropriately scaled least-squares bilinear form of the size of the jumps in the coefficients leading to adequate finite element approximation results. The occurrence of singularities at interface corners and cross-points is discussed. and a weighted least-squares functional is introduced to handle such cases. Numerical experiments are presented for two test problems to illustrate the performance of this approach.

Manteuffel, Thomas A.↗

NASA NTRS2000

Discontinuous dual-primal mixed finite elements for elliptic problems

We propose a novel discontinuous mixed finite element formulation for the solution of second-order elliptic problems. Fully discontinuous piecewise polynomial finite element spaces are used for the trial and test functions. The discontinuous nature of the test functions at the element interfaces allows to introduce new boundary unknowns that, on the one hand enforce the weak continuity of the trial functions, and on the other avoid the need to define a priori algorithmic fluxes as in standard discontinuous Galerkin methods. Static condensation is performed at the element level, leading to a solution procedure based on the sole interface unknowns. The resulting family of discontinuous dual-primal mixed finite element methods is presented in the one and two-dimensional cases. In the one-dimensional case, we show the equivalence of the method with implicit Runge-Kutta schemes of the collocation type exhibiting optimal behavior. Numerical experiments in one and two dimensions demonstrate the order accuracy of the new method, confirming the results of the analysis.

Bottasso, Carlo L.↗

NASA NTRS1999

A Gas-Kinetic Method for Hyperbolic-Elliptic Equations and Its Application in Two-Phase Fluid Flow

A gas-kinetic method for the hyperbolic-elliptic equations is presented in this paper. In the mixed type system, the co-existence and the phase transition between liquid and gas are described by the van der Waals-type equation of state (EOS). Due to the unstable mechanism for a fluid in the elliptic region, interface between the liquid and gas can be kept sharp through the condensation and evaporation process to remove the "averaged" numerical fluid away from the elliptic region, and the interface thickness depends on the numerical diffusion and stiffness of the phase change. A few examples are presented in this paper for both phase transition and multifluid interface problems.

Xu, Kun↗

NASA NTRS1998

A Parallel Non-Overlapping Domain-Decomposition Algorithm for Compressible Fluid Flow Problems on Triangulated Domains

This paper considers an algebraic preconditioning algorithm for hyperbolic-elliptic fluid flow problems. The algorithm is based on a parallel non-overlapping Schur complement domain-decomposition technique for triangulated domains. In the Schur complement technique, the triangulation is first partitioned into a number of non-overlapping subdomains and interfaces. This suggests a reordering of triangulation vertices which separates subdomain and interface solution unknowns. The reordering induces a natural 2 x 2 block partitioning of the discretization matrix. Exact LU factorization of this block system yields a Schur complement matrix which couples subdomains and the interface together. The remaining sections of this paper present a family of approximate techniques for both constructing and applying the Schur complement as a domain-decomposition preconditioner. The approximate Schur complement serves as an algebraic coarse space operator, thus avoiding the known difficulties associated with the direct formation of a coarse space discretization. In developing Schur complement approximations, particular attention has been given to improving sequential and parallel efficiency of implementations without significantly degrading the quality of the preconditioner. A computer code based on these developments has been tested on the IBM SP2 using MPI message passing protocol. A number of 2-D calculations are presented for both scalar advection-diffusion equations as well as the Euler equations governing compressible fluid flow to demonstrate performance of the preconditioning algorithm.

Barth, Timothy J.↗

NASA NTRS1988

Domain decomposition algorithms and computation fluid dynamics

In the past several years, domain decomposition was a very popular topic, partly motivated by the potential of parallelization. While a large body of theory and algorithms were developed for model elliptic problems, they are only recently starting to be tested on realistic applications. The application of some of these methods to two model problems in computational fluid dynamics are investigated. Some examples are two dimensional convection-diffusion problems and the incompressible driven cavity flow problem. The construction and analysis of efficient preconditioners for the interface operator to be used in the iterative solution of the interface solution is described. For the convection-diffusion problems, the effect of the convection term and its discretization on the performance of some of the preconditioners is discussed. For the driven cavity problem, the effectiveness of a class of boundary probe preconditioners is discussed.

Chan, Tony F.↗

NASA NTRS1974

On the flame spreading over a polymer surface

The characteristics of flame propagation along the interface between a solid fuel and a gaseous oxidizer are considered. It is found that difficulties in a study of this problem are related to the eigenvalue character of the relative velocity of the flame with respect to the oxidizer flow. In the proposed model of flame spreading the difficulties are overcome by retaining the ellipticity of the governing equations and by taking the surface pyrolysis into account. The solution procedure makes use of a sink-source approach and the Fourier-transform method. Flame propagation velocity data computed on the basis of the described approach are compared with experimental results.

Ohki, Y.↗

NASA NTRS1996

MGLab: An Interactive Multigrid Environment

MGLab is a set of Matlab functions that defines an interactive environment for experimenting with multigrid algorithms. The package solves two-dimensional elliptic partial differential equations discretized using either finite differences or finite volumes, depending on the problem. Built-in problems include the Poisson equation, the Helmholtz equation, a convection-diffusion problem, and a discontinuous coefficient problem. A number of parameters controlling the multigrid V-cycle can be set using a point-and-click mechanism. The menu-based user interface also allows a choice of several Krylov subspace methods, including CG, GMRES(k), and Bi-CGSTAB, which can be used either as stand-alone solvers or as multigrid acceleration schemes. The package exploits Matlab's visualization and sparse matrix features and has been structured to be easily extensible.

Bordner, James↗

NASA NTRS1988

Boundary element method for 3-D cracks in a plate

Fundamental solutions which automatically satisfy boundary conditions at the interfaces of an elastic plate perfectly bonded to two elastic halfspaces are implemented in a three-dimensional BEM for crack problems. The BEM features a new integration scheme for highly singular kernels. The capability is achieved through a part analytic and part numerical integration procedure, such that the analytic part of the integration is similar for all slip/opening variations. Part-through elliptic cracks in an elastic plate with traction-free surfaces are analyzed and the SIF values along the crack front are found to compare favorably with the numerical SIF results of Raju and Newman (1979).

Fares, N.↗

NASA NTRS1985

Two-dimensional modelling of compressive failure in delaminated laminates

An analytical model is developed to assess the compressive strength criticality of near-surface interlaminar defects in laminated composites. The delaminated region is elliptic in shape, separating a thick isotropic plate from a thin orthotropic layer whose material axes coincide with the ellipse axes. The growth conditions and growth behavior of this defect are studied by breaking the overall problem into an elastic stability problem and a fracture problem. Post-buckling solution for the elliptic section is obtained using the Rayleigh-Ritz method while an energy balance criterion based on a self-similar disbond growth governs the fracture. The parameters controlling the growth or arrest of the delamination damage are identified as the fracture energy, disbond depth and elastic properties of the materials from both sides of the delaminating interface. By varying the degree of material anisotropy relative to the loading axis a range in growth behavior was found including stable or unstable crack growth parallel to or normal to the loading axis.

Chai, H.↗

NASA NTRS1995

Key design issues for user terminals operating with HEO satellites

Coverage of high-latitude areas in the northern hemisphere using geostationary satellites can not be achieved with good elevation angles, and highly inclined elliptical orbits (HEO) present one possible way to overcome this problem. The high elevation angles they achieve and the relative motion of the satellites across the sky, including satellite hand-over, present a specific set of requirements for the design of user terminals, be they for communication or broadcast-reception purposed. This paper reviews key design issues for user terminals operating with satellites in highly inclined elliptical orbits. While all aspects of user terminal design are broadly covered, special emphasis is given to the antenna and RF front-end design. The types of terminals discussed include communications terminals: the hand-held, the notebook, and the vehicle-mounted types, as well as the vehicle-mounted digital audio broadcasting (DAB) receiver terminal. A number of solutions are outlined, based on tradeoffs with respect to the most significant system interface parameters such as G/T and equivalent isotropically radiated power (EIRP). Candidate antenna and front-end configurations were studied, in particular patch and helical antennas with the principal aim of deriving low cost, yet effective solutions. Issues related to achieving best EIRP for the communications terminals are presented, taking into account radiation safety limits. The differences in the design of receive only versus receive/transmit RF front-ends for HEO system terminals are highlighted and solutions for both proposed. This paper summarizes the results of internal ESTEC studies and external ESA-sponsored work concerning major user terminal parameters and design aspects performed for Archimedes -- a HEO satellite constellation for multimedia broadcasting and communications.

Stojkovic, I.↗

NASA NTRS2003

Flow Field Measurement of Mixing Driven by Buoyancy

Mixing driven by buoyancy-induced flows inside a cavity consists of stretching and folding of an interface. Measurement of the flow field using particle imaging velocimetry shows that during stretching the flow field has a single elliptic point, thus dominated by a single vortex. However, global bifurcation that results in folding introduces a hyperbolic point whereby the flow field degenerates to multiple vortex interactions. The short-lived coherent structure observed during mixing which results in the Rayleigh- Taylor morphology is attributed to vortex interactions. The mixing characteristics of non-homogeneous fluids driven by buoyancy are important towards understanding transport phenomenon in a microgravity environment. Mixing consists of stretching and folding of an interface due to a flow field whose intensity depends on the body force. For miscible liquids, the characteristic of the flow field determines whether mass transport is governed by diffusion or bulk stirring which induces mixing. For technologically important processes, transport of mass is governed by the coupling of the body force to scalar gradients such as concentration and or temperature' 2 3 . In order to lend insight into these classes of problems we consider a model experimental system to study mixing driven by buoyancy-induced flows. The characteristics of mixing is addressed from detail measurements of the flow field using particle imaging velocimetry (PIV), and its corresponding interface dynamics using image processing techniques.

Batur, C.↗

NASA NTRS2010

GOATS - Orbitology Component

The GOATS Orbitology Component software was developed to specifically address the concerns presented by orbit analysis tools that are often written as stand-alone applications. These applications do not easily interface with standard JPL first-principles analysis tools, and have a steep learning curve due to their complicated nature. This toolset is written as a series of MATLAB functions, allowing seamless integration into existing JPL optical systems engineering modeling and analysis modules. The functions are completely open, and allow for advanced users to delve into and modify the underlying physics being modeled. Additionally, this software module fills an analysis gap, allowing for quick, high-level mission analysis trades without the need for detailed and complicated orbit analysis using commercial stand-alone tools. This software consists of a series of MATLAB functions to provide for geometric orbit-related analysis. This includes propagation of orbits to varying levels of generalization. In the simplest case, geosynchronous orbits can be modeled by specifying a subset of three orbit elements. The next case is a circular orbit, which can be specified by a subset of four orbit elements. The most general case is an arbitrary elliptical orbit specified by all six orbit elements. These orbits are all solved geometrically, under the basic problem of an object in circular (or elliptical) orbit around a rotating spheroid. The orbit functions output time series ground tracks, which serve as the basis for more detailed orbit analysis. This software module also includes functions to track the positions of the Sun, Moon, and arbitrary celestial bodies specified by right ascension and declination. Also included are functions to calculate line-of-sight geometries to ground-based targets, angular rotations and decompositions, and other line-of-site calculations. The toolset allows for the rapid execution of orbit trade studies at the level of detail required for the early stage of mission concept development.

Haber, Benjamin M.↗

NASA NTRS1992

The kinematics of buoyancy induced mixing

A model problem to study mixing driven by buoyancy flow fields due to a time dependent body force which consists of a steady and an oscillatory component is considered. Flow fields generated by a time dependent body force, which consists of a steady and an oscillatory component from vibrations or g-jitter, mix two fluids inside a cavity by stretching and folding the interface into various morphological patterns. These patterns show two basic structures depending on the dominant components of the body force. These structures are whirls and tendrils, corresponding to the dominance of either the steady or oscillatory component, respectively. A combination of these structures also occurs in the interaction region. In this region the effect of the steady component opposes the oscillatory component and the result is to smooth the interface, and prevent breakup of the interface. For the parametric range considered, two basic flow regimes occur: convective and chaotic. In the convective regime, the morphological patterns are topological. In the chaotic regime, mixing occurs by a repetitive sequence of bubble formation, necking, and breakup which spreads the concentration field throughout the cavity. The length stretch of the interface increases exponentially with time. The dynamical system characteristics show that in the convective regime, at a point in space, the flow field is oscillatory with elliptical phase space trajectories, and its power spectrum indicates that the flow field responds to the corresponding input frequency. Whereas, in the chaotic regime, the flow field is aperiodic, the phase space trajectories show irregular patterns, and the power spectrum shows a broadband distribution.

Duval, Walter M. B.↗