Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DIRICHLET PROBLEM”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Numerical algorithms for finite element computations on arrays of microprocessors

The development of a multicolored successive over relaxation (SOR) program for the finite element machine is discussed. The multicolored SOR method uses a generalization of the classical Red/Black grid point ordering for the SOR method. These multicolored orderings have the advantage of allowing the SOR method to be implemented as a Jacobi method, which is ideal for arrays of processors, but still enjoy the greater rate of convergence of the SOR method. The program solves a general second order self adjoint elliptic problem on a square region with Dirichlet boundary conditions, discretized by quadratic elements on triangular regions. For this general problem and discretization, six colors are necessary for the multicolored method to operate efficiently. The specific problem that was solved using the six color program was Poisson's equation; for Poisson's equation, three colors are necessary but six may be used. In general, the number of colors needed is a function of the differential equation, the region and boundary conditions, and the particular finite element used for the discretization.

Ortega, J. M.↗
DOE OSTI2015

A numerical approach for simulating fluid structure interaction of flexible thin shells undergoing arbitrarily large deformations in complex domains

We present a new numerical methodology for simulating fluid–structure interaction (FSI) problems involving thin flexible bodies in an incompressible fluid. The FSI algorithm uses the Dirichlet–Neumann partitioning technique. The curvilinear immersed boundary method (CURVIB) is coupled with a rotation-free finite element (FE) model for thin shells enabling the efficient simulation of FSI problems with arbitrarily large deformation. Turbulent flow problems are handled using large-eddy simulation with the dynamic Smagorinsky model in conjunction with a wall model to reconstruct boundary conditions near immersed boundaries. The CURVIB and FE solvers are coupled together on the flexible solid–fluid interfaces where the structural nodal positions, displacements, velocities and loads are calculated and exchanged between the two solvers. Loose and strong coupling FSI schemes are employed enhanced by the Aitken acceleration technique to ensure robust coupling and fast convergence especially for low mass ratio problems. The coupled CURVIB-FE-FSI method is validated by applying it to simulate two FSI problems involving thin flexible structures: 1) vortex-induced vibrations of a cantilever mounted in the wake of a square cylinder at different mass ratios and at low Reynolds number; and 2) the more challenging high Reynolds number problem involving the oscillation of an inverted elastic flag. For both cases the computed results are in excellent agreement with previous numerical simulations and/or experiential measurements. Grid convergence tests/studies are carried out for both the cantilever and inverted flag problems, which show that the CURVIB-FE-FSI method provides their convergence. Finally, the capability of the new methodology in simulations of complex cardiovascular flows is demonstrated by applying it to simulate the FSI of a tri-leaflet, prosthetic heart valve in an anatomic aorta and under physiologic pulsatile conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Domain Optimization for an Acoustic Waveguide Scattering Problem

We consider a domain optimization problem of an unbounded domain, which models scattering of a time-harmonic acoustic wave at the junction of two closed waveguides. Solutions of our problem fulfill the Helmholtz equation with a real wavenumber, a modal radiation condition and homogeneous Dirichlet boundary conditions. We derive an a-priori bound for the solution on a certain class of domains (which is compact in the Hausdorff metric topology) and show that within this class the solution depends H{sup 1}-continuously on the domain. Furthermore, we show some numerical examples to illustrate our results, which were calculated using the domain (or shape) derivative of our problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Explicit block encodings of boundary value problems for many-body elliptic operators

Simulation of physical systems is one of the most promising use cases of future digital quantum computers. In this work we systematically analyze the quantum circuit complexities of block encoding the discretized elliptic operators that arise extensively in numerical simulations for partial differential equations, including high-dimensional instances for many-body simulations. When restricted to rectangular domains with separable boundary conditions, we provide explicit circuits to block encode the many-body Laplacian with separable periodic, Dirichlet, Neumann, and Robin boundary conditions, using standard discretization techniques from low-order finite difference methods. To obtain high-precision, we introduce a scheme based on periodic extensions to solve Dirichlet and Neumann boundary value problems using a high-order finite difference method, with only a constant increase in total circuit depth and subnormalization factor. We then present a scheme to implement block encodings of differential operators acting on more arbitrary domains, inspired by Cartesian immersed boundary methods. We then block encode the many-body convective operator, which describes interacting particles experiencing a force generated by a pair-wise potential given as an inverse power law of the interparticle distance. This work provides concrete recipes that are readily translated into quantum circuits, with depth logarithmic in the total Hilbert space dimension, that block encode operators arising broadly in applications involving the quantum simulation of quantum and classical many-body mechanics.

Kharazi, Tyler [University of California, Berkeley↗
DOE OSTI2013

Lower bounds for sums of eigenvalues of elliptic operators and systems

Two-term lower bounds of Berzin-Li-Yau type are obtained for the sums of eigenvalues of elliptic operators and systems with constant coefficients and Dirichlet boundary conditions. The polyharmonic operator, the Stokes system and its generalizations, the two-dimensional buckling problem, and also the Klein-Gordon operator are considered. Bibliography: 32 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Theoretical analysis and numerical simulation for a hyperbolic equation with Dirichlet and acoustic boundary conditions

This paper is concerned with a theoretical and numerical study for the initial-boundary value problem for a linear hyperbolic equation with variable coefficient and acoustic boundary conditions. On the theoretical results, we prove the existence and uniqueness of global solutions, and the uniform stability of the total energy. Numerical simulations using the finite element method associated with the finite difference method are employed, for one-dimensional and two-dimensional cases, to validate the theoretical results. In addition, numerically the uniform decay rate for energy and the order of convergence of the approximate solution are also shown.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A generalized Poisson solver for first-principles device simulations

Electronic structure calculations of atomistic systems based on density functional theory involve solving the Poisson equation. In this paper, we present a plane-wave based algorithm for solving the generalized Poisson equation subject to periodic or homogeneous Neumann conditions on the boundaries of the simulation cell and Dirichlet type conditions imposed at arbitrary subdomains. In this way, source, drain, and gate voltages can be imposed across atomistic models of electronic devices. Dirichlet conditions are enforced as constraints in a variational framework giving rise to a saddle point problem. The resulting system of equations is then solved using a stationary iterative method in which the generalized Poisson operator is preconditioned with the standard Laplace operator. The solver can make use of any sufficiently smooth function modelling the dielectric constant, including density dependent dielectric continuum models. For all the boundary conditions, consistent derivatives are available and molecular dynamics simulations can be performed. The convergence behaviour of the scheme is investigated and its capabilities are demonstrated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Numerical electromagnetic frequency domain analysis with discrete exterior calculus

In this paper, we perform a numerical analysis in frequency domain for various electromagnetic problems based on discrete exterior calculus (DEC) with an arbitrary 2-D triangular or 3-D tetrahedral mesh. We formulate the governing equations in terms of DEC for 3-D and 2-D inhomogeneous structures, and also show that the charge continuity relation is naturally satisfied. Then we introduce a general construction for signed dual volume to incorporate material information and take into account the case when circumcenters fall outside triangles or tetrahedrons, which may lead to negative dual volume without Delaunay triangulation. Then we examine the boundary terms induced by the dual mesh and provide a systematical treatment of various boundary conditions, including perfect magnetic conductor (PMC), perfect electric conductor (PEC), Dirichlet, periodic, and absorbing boundary conditions (ABC) within this method. An excellent agreement is achieved through the numerical calculation of several problems, including homogeneous waveguides, microstructured fibers, photonic crystals, scattering by a 2-D PEC, and resonant cavities. - Highlights: • The formulation of 2D and 3D electromagnetic equations with DEC are presented. • Characteristic Hodge stars are constructed without well-centered triangulation or Delaunay triangulation. • Systematical treatment of various boundary conditions, including PMC, PEC, periodic, and absorbing boundary conditions. • Switching between E to H field can change PEC from essential boundary condition to natural boundary condition. • Quality factor of an open dielectric sphere is considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Elastic waves trapped by a homogeneous anisotropic semicylinder

It is established that the problem of elastic oscillations of a homogeneous anisotropic semicylinder (console) with traction-free lateral surface (Neumann boundary condition) has no eigenvalues when the console is clamped at one end (Dirichlet boundary condition). If the end is free, under additional requirements of elastic and geometric symmetry, simple sufficient conditions are found for the existence of an eigenvalue embedded in the continuous spectrum and generating a trapped elastic wave, that is, one which decays at infinity at an exponential rate. The results are obtained by generalizing the methods developed for scalar problems, which however require substantial modification for the vector problem in elasticity theory. Examples are given and open questions are stated. Bibliography: 53 titles.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Calculation of vortex flows on complex configurations

The calculation of aerodynamic characteristics of complex configurations having strongly coupled vortex flows is a non-linear problem requiring iterative solution techniques. This paper discusses the use of a low-order panel method as a means of obtaining practical solutions to such problems. The panel method is based on piecewise constant source and doublet quadrilateral panels and uses the internal Dirichlet boundary condition of zero perturbation potential. The problems of predicting vortex/surface interaction and vortex separation are discussed. Some example calculations are included but further test cases have yet to be carried out, in particular for comparisons with experimental data. The problem of convergence on the iterative calculation for the shape of the free vortex sheet is addressed and a preprocessor routine, based on an unsteady, two-dimensional version of the panel method, is put forward as a cost-effective way of generating an initial vortex structure for use as a starting solution for general configurations.

Maskew, B.↗
DOE OSTI2018

Layer potential approach for fast eigenvalue characterization of the Helmholtz equation with mixed boundary conditions

Our goal is to propose an efficient approach to characterize the eigenvalues and eigenfunctions of the Helmholtz equation with mixed (Dirichlet and Neumann) boundary conditions. Our approach is based on layer potentials. We extend the eigenvalue characterization known for Neumann boundary conditions to the case of mixed boundary conditions. The problem is motivated by the need of such methods for real-time wave-field shaping by electronically tunable surfaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A stable partitioned FSI algorithm for incompressible flow and deforming beams

An added-mass partitioned (AMP) algorithm is described for solving fluid–structure interaction (FSI) problems coupling incompressible flows with thin elastic structures undergoing finite deformations. The new AMP scheme is fully second-order accurate and stable, without sub-time-step iterations, even for very light structures when added-mass effects are strong. The fluid, governed by the incompressible Navier–Stokes equations, is solved in velocity-pressure form using a fractional-step method; large deformations are treated with a mixed Eulerian-Lagrangian approach on deforming composite grids. The motion of the thin structure is governed by a generalized Euler–Bernoulli beam model, and these equations are solved in a Lagrangian frame using two approaches, one based on finite differences and the other on finite elements. The key AMP interface condition is a generalized Robin (mixed) condition on the fluid pressure. This condition, which is derived at a continuous level, has no adjustable parameters and is applied at the discrete level to couple the partitioned domain solvers. Special treatment of the AMP condition is required to couple the finite-element beam solver with the finite-difference-based fluid solver, and two coupling approaches are described. A normal-mode stability analysis is performed for a linearized model problem involving a beam separating two fluid domains, and it is shown that the AMP scheme is stable independent of the ratio of the mass of the fluid to that of the structure. A traditional partitioned (TP) scheme using a Dirichlet–Neumann coupling for the same model problem is shown to be unconditionally unstable if the added mass of the fluid is too large. A series of benchmark problems of increasing complexity are considered to illustrate the behavior of the AMP algorithm, and to compare the behavior with that of the TP scheme. The results of all these benchmark problems verify the stability and accuracy of the AMP scheme. Results for one benchmark problem modeling blood flow in a deforming artery are also compared with corresponding results available in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Numerical boundary condition procedure for the transonic axisymmetric inverse problem

Two types of boundary condition procedures for the axisymmetric inverse problem are described. One is a Neumann type boundary condition (analogous to the analysis problem) and the other is a Dirichlet type boundary conditon, both requiring special treatments to make the inverse scheme numerically stable. The dummy point concept is utilized in implementing both. Results indicate the Dirichlet type inverse boundary condition is more robust and conceptually simpler to implement than the Neumann type procedure. A few results demonstrating the powerful capability of the newly developed inverse method that can handle both shocked as well as shockless body design are included.

Shankar, V.↗
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.↗
DOE OSTI2016

Stability estimate for the aligned magnetic field in a periodic quantum waveguide from Dirichlet-to-Neumann map

In this article, we study the boundary inverse problem of determining the aligned magnetic field appearing in the magnetic Schrödinger equation in a periodic quantum cylindrical waveguide, by knowledge of the Dirichlet-to-Neumann map. We prove a Hölder stability estimate with respect to the Dirichlet-to-Neumann map, by means of the geometrical optics solutions of the magnetic Schrödinger equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Casimir effect for a scalar field via Krein quantization

In this work, we present a rather simple method to study the Casimir effect on a spherical shell for a massless scalar field with Dirichlet boundary condition by applying the indefinite metric field (Krein) quantization technique. In this technique, the field operators are constructed from both negative and positive norm states. Having understood that negative norm states are un-physical, they are only used as a mathematical tool for renormalizing the theory and then one can get rid of them by imposing some proper physical conditions. -- Highlights: • A modification of QFT is considered to address the vacuum energy divergence problem. • Casimir energy of a spherical shell is calculated, through this approach. • In this technique, it is shown, the theory is automatically regularized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Study of a mixed dispersal population dynamics model

In this study, we consider a mixed dispersal model with periodic and Dirichlet boundary conditions and its corresponding linear eigenvalue problem. This model describes the time evolution of a population which disperses both locally and non-locally. We investigate how long time dynamics depend on the parameter values. Furthermore, we study the minimization of the principal eigenvalue under the constraints that the resource function is bounded from above and below, and with a fixed total integral. Biologically, this minimization problem is motivated by the question of determining the optimal spatial arrangement of favorable and unfavorable regions for the species to die out more slowly or survive more easily. Our numerical simulations indicate that the optimal favorable region tends to be a simply-connected domain. Numerous results are shown to demonstrate various scenarios of optimal favorable regions for periodic and Dirichlet boundary conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Low- and high-order accurate boundary conditions: From Stokes to Darcy porous flow modeled with standard and improved Brinkman lattice Boltzmann schemes

The present contribution focuses on the accuracy of reflection-type boundary conditions in the Stokes–Brinkman–Darcy modeling of porous flows solved with the lattice Boltzmann method (LBM), which we operate with the two-relaxation-time (TRT) collision and the Brinkman-force based scheme (BF), called BF-TRT scheme. In parallel, we compare it with the Stokes–Brinkman–Darcy linear finite element method (FEM) where the Dirichlet boundary conditions are enforced on grid vertices. In bulk, both BF-TRT and FEM share the same defect: in their discretization a correction to the modeled Brinkman equation appears, given by the discrete Laplacian of the velocity-proportional resistance force. This correction modifies the effective Brinkman viscosity, playing a crucial role in the triggering of spurious oscillations in the bulk solution. While the exact form of this defect is available in lattice-aligned, straight or diagonal, flows; in arbitrary flow/lattice orientations its approximation is constructed. At boundaries, we verify that such a Brinkman viscosity correction has an even more harmful impact. Already at the first order, it shifts the location of the no-slip wall condition supported by traditional LBM boundary schemes, such as the bounce-back rule. For that reason, this work develops a new class of boundary schemes to prescribe the Dirichlet velocity condition at an arbitrary wall/boundary-node distance and that supports a higher order accuracy in the accommodation of the TRT-Brinkman solutions. For their modeling, we consider the standard BF scheme and its improved version, called IBF; this latter is generalized in this work to suppress or to reduce the viscosity correction in arbitrarily oriented flows. Our framework extends the one- and two-point families of linear and parabolic link-wise boundary schemes, respectively called B-LI and B-MLI, which avoid the interference of the Brinkman viscosity correction in their closure relations. The performance of LBM and FEM is thoroughly evaluated in three benchmark tests, which are run throughout three distinctive permeability regimes. The first configuration is a horizontal porous channel, studied with a symbolic approach, where we construct the exact solutions of FEM and BF/IBF with different boundary schemes. The second problem refers to an inclined porous channel flow, which brings in as new challenge the formation of spurious boundary layers in LBM; that is, numerical artefacts that arise due to a deficient accommodation of the bulk solution by the low-accurate boundary scheme. The third problem considers a porous flow past a periodic square array of solid cylinders, which intensifies the previous two tests with the simulation of a more complex flow pattern. The ensemble of numerical tests provides guidelines on the effect of grid resolution and the TRT free collision parameter over the accuracy and the quality of the velocity field, spanning from Stokes to Darcy permeability regimes. It is shown that, with the use of the high-order accurate boundary schemes, the simple, uniform-mesh-based TRT-LBM formulation can even surpass the accuracy of FEM employing hardworking body-fitted meshes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗