Search NASA⌕ Search

SEARCH · Search NASA

Results for “method of cells”

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 37 records · Page 2

Simulation of three-dimensional incompressible flows with a vortex-in-cell method

A new method for the numerical simulation of three-dimensional incompressible flows is described. The vortex-in-cell (VIC) method presented traces the motion of the vortex filaments in the velocity field which these filaments create. The velocity field is not calculated directly by the Biot-Savart law of interaction but by creating a mesh-record of the vorticity field, then integrating a Poisson's equation via the fast Fourier transform to generate a mesh-record of the velocity field. The computed scales of motion are assumed to be essentially inviscid. Viscous of subgrid-scale effects are incorporated into a filtering procedure in wave vector space. Results of tracing a periodic array of single vortex rings are compared with a Green's function calculation. The agreement is very good.

Couet, B.↗

Calibration of solar cells by the reference cell method - The spectral mismatch problem

The calibration of solar cells by means of solar simulators and calibrated reference cells has always been faced with the problem of errors due to source spectral irradiance mismatch and cell spectral response mismatch. A general expression is derived in this paper to enable calculation of mismatch error. The information required to make this calculation is the relative spectra responses of the test and reference cells, and the relative spectral irradiance of the simulator (over the spectral range defined by cell response). The spectral irradiance of the solar standard is assumed given. Experimental results are presented.

Seaman, C. H.↗

Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow

We present preliminary development of an approach for simulating high Reynolds number steady compressible flow in two space dimensions using a Cartesian cut-cell finite volume method. We consider both laminar and turbulent flow with both low and high cell Reynolds numbers near the wall. The approach solves the full Navier-Stokes equations in all cells, and uses a wall model to address the resolution requirements near boundaries and to mitigate mesh irregularities in cut cells. We present a quadratic wall model for low cell Reynolds numbers. At high cell Reynolds numbers, the quadratic is replaced with a newly developed analytic wall model stemming from solution of a limiting form of the Spalart-Allmaras turbulence model which features a forward evaluation for flow velocity and exactly matches characteristics of the SA turbulence model in the field. We develop multigrid operators which attain convergence rates similar to inviscid multigrid. Investigations focus on preliminary verification and validation of the method. Flows over flat plates and compressible airfoils show good agreement with both theoretical results and experimental data. Mesh convergence studies on sub- and transonic airfoil flows show convergence of surface pressures with wall spacings as large as approx.0.1% chord. With the current analytic wall model, one or two additional refinements near the wall are required to obtain mesh converged values of skin friction.

Berger, Marsha↗

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Particle-in-cell Method for Plasmas with A Generalized Momentum Formulation, Part III: A family of Gauge Conserving Methods

In this paper, we introduce a new family of spatially co-located field solvers for particle-in-cell applications which evolve the potential formulation of Maxwell’s equations under the Lorenz gauge. Our recent work [2] introduced the concept of time-consistency, which connects charge conservation to the preservation of the gauge at the semi-discrete level. It will be shown that there exists a large family of time discretizations which satisfy this property. Additionally, it will be further shown that for large classes of time marching methods, the satisfaction of the gauge condition automatically implies the satisfaction of Gauss’s law for electricity, with the potential formulation ensuring that that Gauss’s law for magnetism is satisfied by definition. We focus on popular time marching methods including centered differences, backward differences, and diagonally-implicit Runge-Kutta methods, which are coupled to a spectral discretization in space. We demonstrate the theory by testing the methods on a relativistic Weibel instability and a drifting cloud of electrons.

97 MATHEMATICS AND COMPUTING↗

Accuracy of the explicit energy-conserving particle-in-cell method for under-resolved simulations of capacitively coupled plasma discharges

The traditional explicit electrostatic momentum-conserving particle-in-cell algorithm requires strict resolution of the electron Debye length to deliver numerical stability and accuracy. The explicit electrostatic energy-conserving particle-in-cell algorithm alleviates this constraint with minimal modification to the traditional algorithm, retaining its simplicity, ease of parallelization, and acceleration on modern supercomputing architectures. In this article, we apply the algorithm to model a one-dimensional radio frequency capacitively coupled plasma discharge relevant to industrial applications. The energy-conserving approach closely matches the results from the momentum-conserving algorithm and retains accuracy even for cell sizes up to 8 times the electron Debye length. For even larger cells, the algorithm loses accuracy due to poor resolution of steep gradients within the radio frequency sheath. Accuracy can be recovered by adopting a non-uniform grid, which resolves the sheath and allows for cell sizes up to 32 times the electron Debye length in the quasi-neutral bulk of the discharge. The effect is an up to 8 times reduction in the number of required simulation cells, an improvement that can compound in higher-dimensional simulations. We therefore consider the explicit energy-conserving algorithm as a promising approach to significantly reduce the computational cost of full-scale device simulations and a pathway to delivering kinetic simulation capabilities of use to industry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A simple centrifuge cell method for ex situ quantification of electrical conductivity of slurry electrode materials

We present the design, experimental procedure, and experimental evaluation of a system for fast, simple, and ex situ characterization of electrical conductivity of slurry electrode materials. The system uses a custom-designed electrochemical cell compatible with centrifugation in a swing-bucket centrifuge. The cell features cylindrical graphite electrodes that are partially sheathed so as to expose only 2 mm of the electrode surface to the bottom region of the packed particulate pellet. Also presented is a conduction model that provides a shape factor for estimating effective conductivity. We tested aqueous solutions of carbon black (CB), activated carbon (AC), and mixtures thereof. These particles were dispersed in 0.0 and 0.5 M NaCl solutions. Measurements show that the effective conductivity initially increases linearly with pellet mass and then saturates at higher masses. Notably, CB exhibited a fivefold increase in conductivity than AC at equal pellet masses. CB/AC mixtures at a fixed pellet mass were tested with CB mass fractions of 0 to 100%. Interestingly, the mixture conductivity was found to be a non-monotonic function of CB mass fraction, with a maximum conductivity at about 60 % CB mass fraction. At this maximum, the mixture conductivity is approximately 30 % higher than pure CB. NaCl concentration in the slurry solution had no effects on conductivity. These results highlight the interactions between slurry electrode composition and compaction, offering insights for optimizing slurry electrodes. Furthermore, the system and results may also be applicable to evaluation of particulate materials (including slurries) used for Li-ion batteries, capacitive deionization, fuel cells, and flow electrodes.

Capacitive deionization↗

Three-dimensional simulation of the free shear layer using the vortex-in-cell method

We present numerical simulations of the evolution of a mixing layer from an initial state of uniform vorticity with simple two- and three-dimensional small perturbations. A new method for tracing a large number of three-dimensional vortex filaments is used in the simulations. Vortex tracing by Biot-Savart interaction originally implied ideal (non-viscous) flow, but we use a 3-d mesh, Fourier transforms and filtering for vortex tracing, which implies 'modeling' of subgrid scale motion and hence some viscosity. Streamwise perturbations lead to the usual roll-up of vortex patterns with spanwise uniformity maintained. Remarkably, spanwise perturbations generate streamwise distortions of the vortex filaments and the combination of both perturbations leads to patterns with interesting features discernable in the movies and in the records of enstrophy and energy for the three components of the flow.

Couet, B.↗

FLIP MHD - A particle-in-cell method for magnetohydrodynamics

The fluid-implicit-particle, or 'FLIP' method presently extended to 2D and 3D MHD flow incorporates a Lagrangian field representation and yields a grid magnetic Reynolds number of up to 16 while preserving contact continuities that retain the Galilean invariance of the MHD flow equations. Analytical arguments and numerical examples demonstrate the conservation of mass, momentum, magnetic flux, and energy; 2D calculation results for the illustrative cases of contact discontinuity convection, Rayleigh-Taylor unstable flow.

Brackbill, J. U.↗

Estimation of Effective Elastic Properties of General Multifunctional Honeycomb Structures Using a Unit Cell Method

Sandwich composite structures are ideal configurations in which to incorporate additional functionality beyond load carrying capabilities. The inner core can be layered to facilitate other functions such as power storage for a battery. In this work we investigate an assemblage of analytical tools to compute effective properties that allow complex layered core architectures to be homogenized into a single continuum layer. This provides a great increase in computational efficiency to numerically simulate the structural response of multifunctional sandwich structures under applied loads.

Krishnamurthy, T.↗

Investigating Local Buckling and Plastic Collapse in Multifunctional Sandwich Composite Cores using a Unit Cell Method

The primary goal of developing analytical models of composite sandwich cores has been to obtain simple expressions for elastic and failure response to improve the rapid analysis of composite structures. An earlier investigation by the authors was made to verify the accuracy of widely-used analytical models that calculate equivalent elastic moduli of a uniform homogeneous single layer of material in place of the complex honeycomb core properties in sandwich composites. This transformation of elastic properties permits a coarsening of the finite element model (FEM) of the core to a fraction of the computational cost in simulating the original complex core architecture. In addition, the core was assumed constructed of multi-layered isotropic walls to support potential additional functionality beyond simple load-carrying capability.

Erik Saether↗