Results for “flow conservation” — Page 18

DOE OSTI2021

A Circularity Assessment for Silicon Solar Panels Based on Dynamic Material Flow Analysis

Solar photovoltaics (PV) are the fastest growing renewable energy technology for clean, inexpensive, and sustainable electricity generation. Along with numerous technical roadmaps to improve system cost, performance and reliability, the PV industry should also plan to handle large volumes of silicon panel waste, which is initially estimated to be ~13 million metric tons (MT) by 2050 in the U.S. alone. Understanding the magnitude of material needs and how material flows throughout the PV panel life cycle could respond to design, operational and different end-of-life (EOL) circular pathways will help transition into a circular, resource-conserving economy. Herein, we introduce a dynamic material flow analysis (DMFA) framework based on electricity generation to quantify time-series stocks and flows of bulk PV materials (e.g., solar glass and aluminum frames) throughout the life cycles of utility-scale silicon PV systems in the U.S. in the period 2000-2100. We apply the model to a range of scenarios to understand how material demands depend on selected PV-related parameters, different material circularity strategies, and recent module design trends (e.g., bifacial, frameless). We found that float glass and aluminum in PV installations would likely reach 100 million MT and 12 million MT by 2100, respectively, in the baseline scenario. The most influential parameters for PV installation and subsequent waste reduction are found to be module lifetime, module efficiency, annual degradation, and material reduction. Module recycling and component remanufacturing were found to be the most effective material circularity strategies for waste minimization. Panel reuse has negligible savings on waste under current module efficiencies compared to replacements with newer generations with higher efficiency. Ongoing trends to produce larger power frameless modules could save 10 million MT of glass and ~9 million MT of aluminum. Our results enable advanced planning for future materials needs and provide insight into potential opportunities to minimize waste.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

A Circularity Assessment for Silicon Solar Panels Based on Dynamic Material Flow Analysis: Preprint

Solar photovoltaics (PV) are the fastest growing renewable energy technology for clean, inexpensive, and sustainable electricity generation. Along with numerous technical roadmaps to improve system cost, performance and reliability, the PV industry should also plan to handle large volumes of silicon panel waste, which is initially estimated to be ~13 million metric tons (MT) by 2050 in the U.S. alone. Understanding the magnitude of material needs and how material flows throughout the PV panel life cycle could respond to design, operational and different end-of-life (EOL) circular pathways will help transition into a circular, resource-conserving economy. Herein, we introduce a dynamic material flow analysis (DMFA) framework based on electricity generation to quantify time-series stocks and flows of bulk PV materials (e.g., solar glass and aluminum frames) throughout the life cycles of utility-scale silicon PV systems in the U.S. in the period 2000-2100. We apply the model to a range of scenarios to understand how material demands depend on selected PV-related parameters, different material circularity strategies, and recent module design trends (e.g., bifacial, frameless). We found that float glass and aluminum in PV installations would likely reach 100 million MT and 12 million MT by 2100, respectively, in the baseline scenario. The most influential parameters for PV installation and subsequent waste reduction are found to be module lifetime, module efficiency, annual degradation, and material reduction. Module recycling and component remanufacturing were found to be the most effective material circularity strategies for waste minimization. Panel reuse has negligible savings on waste under current module efficiencies compared to replacements with newer generations with higher efficiency. Ongoing trends to produce larger power frameless modules could save 10 million MT of glass and ~9 million MT of aluminum. Our results enable advanced planning for future materials needs and provide insight into potential opportunities to minimize waste.

circular economy↗

DOE OSTI2021

Preliminary Implementation of Two-Dimensional Cartesian Solver in CTF-R

Sub-channel codes are one of the the modeling and simulation tools used for thermal-hydraulic analysis of nuclear reactors. A few examples of such sub-channel codes are the COolant Boiling in Rod Arrays (COBRA) family of codes. The approximations that are used to simplify the fluid conservation equations into sub-channel form, mainly that of axially-dominated flow, lead to noticeable limitations on sub-channels solvers for problems with significant flow in lateral directions. In this report, a two-dimensional Cartesian solver is developed and implemented within CTF-R, which is the residual solver in the North Carolina State University version of COBRA-TF (CTF). The new solver will enable CTF to simulate flow that is not axially-dominated. The appropriate Cartesian forms of the conservation equations are derived and implemented in the solver. Once the conservation equations are established, the process of constructing the matrix system was altered to solve a two-dimensional staggered grid system. A simple case was used to test that the two-dimensional Cartesian solver is accurate. The test problem does not include any source terms or flow in the lateral direction. The results show that the solver was able to run the simple case and converge to a steady-state solution. Future work will focus on testing existing capabilities by using test cases that include transients and equation cross-terms. Future work will also include adding additional capabilities such as enabling the solver to include cases with source terms and three dimensional cases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1990

Upward ion flow in ionospheric holes on Venus

The nature of ion flow within two ionospheric holes on the nightside of Venus was investigated using ion composition measurements made by the ion mass spectrometer on the Pioneer Venus Orbiter. A comparison of the altitude profiles of the observed ion densities with those expected under diffusive equilibrium conditions indicates that the major ions O(+), NO(+), and O2(+), as well as the minor ions H(+) and He(+) flow upward and away from Venus along the axes of the holes. This result agrees with a quantitative evaluation of the ion flow speeds appearing in expressions derived from the equations for conservation of mass and momentum of the ions and electrons. The analysis shows that all ion species flow upward in the holes because the upward force produced by the plasma pressure gradient exceeds all downward forces. However, the nature of the ion source required to maintain such flow is not known.

Hartle, R. E.↗

NASA NTRS1983

Numerical study of scramjet and ramjet flow fields

Two computer programs have been developed to numerically calculate complex, two-dimensional flow fields in scramjets. The first program is written for inlet analysis whereas the second program is written primarily for combustor analysis. Both programs solve the full two-dimensional Navier-Stokes equations by a well-known explicit, predictor-corrector technique. Turbulence is modeled by an algebraic eddy-viscosity model. The combustor program also includes one or more species conservation equations to calculate mixing and reacting flows. The hydrogen/air chemistry in this program is modeled by a complete reaction model. The combustor program has been recently modified to analyze axisymmetric ramjet dump combustor flow field. Results from these computer programs are presented that predict the flow in several scramjet inlet configurations, two model scramjet engine configurations, and in a dump combustor simulator. Computed results are also compared with available experimental data to allow assessment of the programs.

Kumar, A.↗

NASA NTRS1986

The construction and use of divergence free vector expansions for incompressible fluid flow calculations

For incompressible fluids the law of mass conservation reduces to a constraint on the velocity vector, namely that it be divergence free. This constraint has long been a source of great difficulty to the numericist seeking to discretize the Navier-Stokes and Euler equations. A spectral method is discussed which overcomes this difficulty. Its efficacy is demonstrated on some simple problems. The velocity is approximated by a finite sum of divergence free vectors, each of which satisfies the same boundary conditions as the velocity. Projecting the governing equation onto the space of inviscid vector fields eliminates the pressure term and produces a set of ordinary differential equations that must be solved for the coefficents in the velocity. The pressure can then be recovered if it is needed.

Mhuiris, N. M. G.↗

NASA NTRS2000

A Unified Wall Boundary Treatment for Viscous and Inviscid Flows in the CE/SE Method

In the setting of the conservation element-solution element (CE/SE) method, a new and unified wall boundary treatment for the Navier-Stokes and Euler Equations is proposed. In essence, the shear stress exerted on the fluid by a wall is modeled as a source term as a part of local spacetime flux conservation in the vicinity of a wall boundary. When the fluid is inviscid, the source term vanishes and the boundary condition reduces to the usual 'slip' condition. On the other hand, when the fluid is viscous, the source-term effect is consistent with the traditional no-slip condition. Numerical results show that the new treatment is robust, efficient, and accurate for viscous and inviscid flows.

Chang, Sin-Chung↗

NASA NTRS1992

Resonant behaviour of MHD waves on magnetic flux tubes. III - Effect of equilibrium flow

The Hollweg et al. (1990) analysis of MHD surface waves in a stationary equilibrium is extended. The conservation laws and jump conditions at Alfven and slow resonance points obtained by Sakurai et al. (1990) are generalized to include an equilibrium flow, and the assumption that the Eulerian perturbation of total pressure is constant is recovered as the special case of the conservation law for an equilibrium with straight magnetic field lines and flow along the magnetic field lines. It is shown that the conclusions formulated by Hollweg et al. are still valid for the straight cylindrical case. The effect of curvature is examined.

Goossens, Marcel↗

NASA NTRS2001

A Space-Time Conservation Element and Solution Element Method for Solving the Two- and Three-Dimensional Unsteady Euler Equations Using Quadrilateral and Hexahedral Meshes

In this paper, we report a version of the Space-Time Conservation Element and Solution Element (CE/SE) Method in which the 2D and 3D unsteady Euler equations are simulated using structured or unstructured quadrilateral and hexahedral meshes, respectively. In the present method, mesh values of flow variables and their spatial derivatives are treated as independent unknowns to be solved for. At each mesh point, the value of a flow variable is obtained by imposing a flux conservation condition. On the other hand, the spatial derivatives are evaluated using a finite-difference/weighted-average procedure. Note that the present extension retains many key advantages of the original CE/SE method which uses triangular and tetrahedral meshes, respectively, for its 2D and 3D applications. These advantages include efficient parallel computing ease of implementing non-reflecting boundary conditions, high-fidelity resolution of shocks and waves, and a genuinely multidimensional formulation without using a dimensional-splitting approach. In particular, because Riemann solvers, the cornerstones of the Godunov-type upwind schemes, are not needed to capture shocks, the computational logic of the present method is considerably simpler. To demonstrate the capability of the present method, numerical results are presented for several benchmark problems including oblique shock reflection, supersonic flow over a wedge, and a 3D detonation flow.

Zhang, Zeng-Chan↗

NASA NTRS1994

Green's formula and variational principles for cosmic-ray transport with application to rotating and shearing flows

Green's theorem and Green's formula for the diffusive cosmic-ray transport equation in relativistic flows are derived. Green's formula gives the solution of the transport equation in terms of the Green's function of the adjoint transport equation, and in terms of distributed sources throughout the region R of interest, plus terms involving the particle intensity and streaming on the boundary. The adjoint transport equation describes the time-reversed particle transport. An Euler-Lagrange variational principle is then obtained for both the mean scattering frame distribution function f, and its adjoint f(dagger). Variations of the variational functional with respect to f(dagger) yield the transport equation, whereas variations of f yield the adjoint transport equation. The variational principle, when combined with Noether's theorem, yields the conservation law associated with Green's theorem. An investigation of the transport equation for steady, azimuthal, rotating flows suggests the introduction of a new independent variable H to replace the comoving frame momentum variable p'. For the case of rigid rotating flows, H is conserved and is shown to be analogous to the Hamiltonian for a bead on a rigidly rotating wire. The variable H corresponds to a balance between the centrifugal force and the particle inertia in the rotating frame. The physical interpretation of H includes a discussion of nonrelativistic and special relativistic rotating flows as well as the cases of aziuthal, differentially rotating flows about Schwarzs-child and Kerr black holes. Green's formula is then applied to the problem of the acceleration of ultra-high-energy cosmic rays by galactic rotation. The model for galactic rotation assumes an angular velocity law Omega = Omega(sub 0)(omega(sub 0)/omega), where omega denotes radial distance from the axis of rotation. Green's functions for the galactic rotation problem are used to investigate the spectrum of accelerated particles arising from monoenergetic and truncated power-law sources. We conclude that it is possible to accelerate particles beyond the knee by galactic rotation, but not in sufficient number to adequately explain the observed spectrum.

Webb, G. M.↗

NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗

NASA NTRS1981

Computation of simple three-dimensional wing-vortex interaction in transonic flow

An implicit delta form finite-difference algorithm for Euler equations in conservation law form has been used in preliminary calculations of three-dimensional wing-vortex interactions. Both steady and unsteady transonic flow wing-vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two- and three-dimensions.

Srinivasan, G. R.↗

NASA NTRS1980

An alternating-direction-implicit algorithm for the unsteady potential equation in conservation form

An implicit finite difference scheme for an efficient computation of unsteady potential flow about airfoils is presented. The formulation uses density and velocity potential as dependent variables, and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time varying sheared rectilinear coordinate transformation is employed. Calculated time history solutions on a pulsating airfoil are compared with the results of another unsteady transonic code. It is concluded that the method has excellent numerical stability and gives accurate solutions with sharply resolved shocks.

Chipman, R. R.↗

NASA NTRS1991

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.↗

NASA NTRS1976

Implicit finite-difference procedures for the computation of vortex wakes

Implicit finite-difference procedures for the primitive form of the incompressible Navier-Stokes and the compressible Euler equations are used to compute vortex wake flows. The partial differential equations in strong conservation-law form are transformed to cluster grid points in regions with large changes in vorticity. In addition to clustering, fourth-order accurate, spatial difference operators are used to help resolve the flow-field gradients. The use of implicit time-differencing permits large time steps to be taken since temporal variations are typically small. Computational efficiency is achieved by approximate factorization. Both two-dimensional and preliminary three-dimensional calculations are described and qualitatively compared with existing experimental data.

Steger, J. L.↗

NASA NTRS1986

Calculation of external-internal flow fields for mixed-compression inlets

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Chyu, W. J.↗

NASA NTRS1987

Calculation of external-internal flow fields for mixed-compression inlets

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Chyu, W. J.↗

DOE OSTI2020

A reconstructed discontinuous Galerkin method for compressible flows in Lagrangian formulation

In this work, a high-order accurate reconstructed discontinuous Galerkin (rDG) method is developed for solving two-dimensional hydrodynamic problems in cell-centered updated Lagrangian formulation. This method is the Lagrangian limit of the unsplit rDG-ALE formulation, and is obtained by assuming the equality of the grid velocity to the fluid velocity only at cell boundaries. The conservative variables and the Taylor basis defined on the time-dependent moving mesh, provide the piece-wise polynomial expansion in the updated Lagrangian formulation. A multi-directional nodal Riemann solver is implemented for computing the grid velocity at the vertices and the numerical flux at the cell boundaries. A characteristic limiting procedure is extended from the primitive variable version to the conservative variable version, and its performance is compared with the limiter on physical variables. A number of benchmark test cases are conducted to assess the accuracy, robustness, and non-oscillatory property of the DG(P0), DG(P1) and rDG(P1P2) methods. The numerical experiments demonstrate that the developed rDG method is able to attain the designed order of accuracy and the characteristic limiting procedure outperforms the limiter on physical variables in terms of the monotonicity and symmetry preservation for shock problems.

42 ENGINEERING↗