Results for “conservation laws” — Page 46

NASA NTRS2019

Universal Expansion with Spatially Varying G

We calculate the expansion of the Universe under the assumptions that G varies in space and the radial size r of the Universe is very large (we call this the MOND regime of varying-G gravity). The inferred asymptotic behaviour turns out to be different from that found by McCrea & Milne in 1934 and our equations bear no resemblance to those of the relativistic case. In this cosmology, the scale factor R(t) increases linearly with time t, the radial velocity is driven by inertia, and gravity is incapable of hindering the expansion. Yet, Hubble’s law is borne out without any additional assumptions. When we include a repulsive acceleration a(sub de) due to dark energy, the resulting universal expansion is then driven totally by this new term and the solutions for a(sub de) → 0 do not reduce to those of the a(sub de) ≡ 0 case. This is a realization of a new Thom catastrophe: The inclusion of the new term alters the conservation of energy and the dark energy solutions are not reducible to those in the case without dark energy.

Christodoulou, Dimitris M.↗

DOE OSTI2023

Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution

Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.

58 GEOSCIENCES↗

DOE OSTI2021

Thermodynamically consistent semi-compressible fluids: a variational perspective

This paper presents (Lagrangian) variational formulations for single and multicomponent semi-compressible fluids with both reversible (entropy-conserving) and irreversible (entropy-generating) processes. Semi-compressible fluids are useful in describing low-Mach dynamics, since they are soundproof. These models find wide use in many areas of fluid dynamics, including both geophysical and astrophysical fluid dynamics. Specifically, the Boussinesq, anelastic and pseudoincompressible equations are developed through a unified treatment valid for arbitrary Riemannian manifolds, thermodynamic potentials and geopotentials. By design, these formulations obey the 1st and 2nd laws of thermodynamics, ensuring their thermodynamic consistency. This general approach extends and unifies existing work, and helps clarify the thermodynamics of semi-compressible fluids. To further this goal, evolution equations are presented for a wide range of thermodynamic variables: entropy density s, specific entropy η, buoyancy b, temperature T, potential temperature θ and a generic entropic variable χ; along with a general definition of buoyancy valid for all three semicompressible models and arbitrary geopotentials. Finally, the elliptic equation for the pressure perturbation (the Lagrange multiplier that enforces semi-compressibility) is developed for all three equation sets in the case of reversible dynamics, and for the Boussinesq/anelastic equations in the case of irreversible dynamics; and some discussion is given of the difficulty in formulating it for the pseudoincompressible equations with irreversible dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

REopt: Energy Decision Analysis

REopt transforms complex energy project decisions into actionable results for building owners, utilities, and industry. Based on decades of NREL decision-support expertise, REopt guides investment in economic, resilient energy technologies. This presentation provides an overview of how REopt can help building owners, utilities, developers, and industry optimize the economic and resilience benefits of distributed energy resources.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS2022

Enhancement of a global lake and reservoir database to aid climate studies and resource monitoring utilizing satellite radar altimetry

Here, we describe the expansion and enhancement of a large (surface area ≥100 sq. km) lake and reservoir database (1.GREALD). These efforts have also resulted in the spin-off of two additional databases, one containing lakes and reservoirs (2.GREALD, area 10–99 sq. km), and the other containing ephemeral lakes (3.GELD, area ≥100 sq. km). These databases are unique sources for projects that utilize satellite radar altimeter data to monitor surface water levels. While 1.GREALD aims to be a complete catalog, 2.GREALD focuses on reservoirs in response to applied sciences programs that monitor water and energy resources. The creation of 3.GELD has climate change objectives as well as water resources and ecosystem conservation applications. The recording of information pertaining to the potential overpasses (waterbody crossings) of the current and archive satellite altimeters is a primary objective as is the need to highlight any form of controlled water level variation. The permanent water databases now contain 6282 entries, half experience some form of anthropogenic influence and ∼430 have been identified as potential climatically sensitive terminal lakes. The revised integral surface area distribution is a power law with exponent −1.016. Statistics reveal that with altimetric repeat visit times of 10-day to monthly, at least 80% of the permanent water bodies (≥10 sq. km) have been overflown at some period since the 1990s. Current information on water use and reservoir formation date show that the primary use of the reservoir class is hydroelectric power, and that China, Brazil, India, Turkey, and Vietnam dominate the dam building in recent decades.

Charon M. Birkett↗

DOE OSTI2024

Particle Acceleration in Relativistic Alfvénic Turbulence

Abstract Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0 . If δ B 0 ≪ B 0 , the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0 , interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

$|\Delta \mathcal{B}| =2$: A State of the Field, and Looking Forward--A brief status report of theoretical and experimental physics opportunities

The origin of the matter-antimatter asymmetry apparently obligates the laws of physics to include some mechanism of baryon number ($\mathcal{B}$) violation. Searches for interactions violating $\mathcal{B}$ and baryon-minus-lepton number $\mathcal{(B-L)}$ represent a rich and underutilized opportunity. These are complementary to the existing, broad program of searches for $\mathcal{L}$-violating modes such as neutrinoless double $\beta$-decay which could provide deeper understandings of the plausibility of leptogenesis, or $\mathcal{B}$-violating, $\mathcal{(B-L)}$-conserving processes such as proton decay. In particular, a low-scale, post-sphaleron violation mechanism of $\mathcal{(B-L)}$ could provide a \textit{testable} form of baryogenesis. Though theoretically compelling, searches for such $\mathcal{(B-L)}$-violating processes like $\Delta\mathcal{B}=2$ dinucleon decay and $n\rightarrow\bar{n}$ remain relatively underexplored experimentally compared to other rare processes. By taking advantage of upcoming facilities such as the Deep Underground Neutrino Experiment and the European Spallation Source, this gap can be addressed with new intranuclear and free searches for neutron transformations with very high sensitivity, perhaps greater than three orders of magnitude higher than previous experimental searches. This proceedings reports on recent theoretical and experimental advances and sensitivities of next-generation searches for neutron transformations were detailed as part of the Amherst Center for Fundamental Interactions Workshop, "Theoretical Innovations for Future Experiments Regarding Baryon Number Violation," directly coordinated with the Rare Processes and Precision Measurements Frontier.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

The conserved bridging domain on HCV E1E2 glycoprotein complex is targeted by neutralizing antibodies from diverse lineages

Explore the source record for details and available documents.

Chen, Fang (ORCID:0000000234813363)↗

NASA NTRS1978

Effect of wall temperature on a supersonic turbulent boundary layer

Measurement of mean flow profiles in a fully developed Mach 3 turbulent boundary layer with negligible pressure gradient is reported. Data were acquired at several streamwise locations for wall-to-total temperature ratios of 0.94, 0.71, and 0.54. The results demonstrate that the velocity defect formulation of the law-of-the-wake, which successfully correlates compressible, adiabatic boundary layers, is also valid for nonadiabatic flows. It is also shown that for adiabatic walls, the linear Crocco relation between total temperature and velocity does not provide a valid test of the nature of the boundary-layer flow for practical cases where the Prandtl number departs from unity. Finally, the turbulent shear stress, mixing length, and eddy viscosity were extracted from the 'time-averaged' conservation equations using the measured mean flow profiles and found to be insensitive to wall temperature. In particular, the latter properties are in good agreement with earlier compressible, adiabatic correlations of turbulent transport properties.

Laderman, A. J.↗

DOE OSTI2024

Error field predictability and consequences for ITER

Abstract ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER’s coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of n = 1 error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scalign laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2024

Equipment Qualification Report Environmental Qualification of GNB Absolyte Valve Regulated Lead Acid (VRLA) 1600 Ah100G33 Battery Rack Assembly (24590-QL-POA-EDB0-00001-11-00002_00A)

Greenberry Environmental Qualification Report 550001.001-35.0.5 provides basis for assignment of qualified life for: GNB Absolyte Valve Regulated Lead Acid (VRLA) 1600 Ah 100G33 Battery Rack Assembly in accordance with the requirements specified 24590-WTP-3PS-G000-T0015 (Rev 2) and Environmental Qualification Plan 550001.001-35.0.1 (Rev. 3). The equipment qualification basis represents the most conservative capability of the equipment. The analysis performed for the qualification is not less conservative than the bounding environmental conditions detailed in contract documents issued to Greenberry in contract 24509-QL-POA-EDB0-00001 Rev.0. The qualified life of 10 years has been established based upon an end-condition objective of the equipment condition indicators that correlate to the ability of equipment to perform its safety function. The VRLA Battery Cell Assembly was aged by 10 years (minimum) in accordance with conditions specified by Bechtel Equipment Qualification Datasheet 24590-LAW-EUQ-UPE 00003 Rev. 3 and the process conditions specified by the Instrument Data Sheet 24590-LAW EUD-UPE-00009 Rev. 2. Greenberry Industrial has contracted with GNB Industrial Battery Co located at 4115 S Zero St, Fort Smith, AR 72908 to perform age conditioning, monitoring, and capacity testing in accordance with Environmental Qualification Plan 550001.001-35.0.1 Rev. 1. The required process at the stated conditions set by the parameters established by the plan were completed satisfactorily. The details of the of the test process observed by Greenberry is detailed in the attached Seismic Test Log 550001.001-7.0.3, including examples of the objective evidence collected during the qualification process.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

West Valley Demonstration Project (WVDP) Annual Site Environmental Report (ASER) for Calendar Year 2020

West Valley Demonstration Project (WVDP) Annual Site Environmental Report (ASER) for Calendar Year 2020. The report, prepared for the U.S. Department of Energy West Valley Demonstration Project office (DOE-WVDP), summarizes the environmental protection program at the WVDP for calendar year (CY) 2020. Monitoring and surveillance of the facilities used by the DOE are conducted to verify protection of public health and safety and the environment. The report is a key component of DOE’s effort to keep the public informed of environmental conditions at the WVDP. The quality assurance protocols applied to the environmental monitoring program ensure the validity and accuracy of the monitoring data. In addition to demonstrating compliance with environmental laws, regulations, and directives, evaluation of data collected in 2020 continued to indicate that WVDP activities pose no threat to public health or safety, or to the environment.

Resource Conservation and Recovery Act↗

NASA NTRS1995

Active galaxies observed during the Extreme Ultraviolet Explorer all-sky survey

We present observations of active galactic nuclei (AGNs) obtained with the Extreme Ultraviolet Explorer (EUVE) during the all-sky survey. A total of 13 sources were detected at a significance of 2.5 sigma or better: seven Seyfert galaxies, five BL Lac objects, and one quasar. The fraction of BL Lac objects is higher in our sample than in hard X-ray surveys but is consistent with the soft X-ray Einstein Slew Survey, indicating that the main reason for the large number of BL Lac objects in the extreme ulktraviolet (EUV) and soft X-ray bands is their steeper X-ray spectra. We show that the number of AGNs observed in both the EUVE and ROSAT Wide Field Camera surveys can readily be explained by modelling the EUV spectra with a simple power law in the case of BL Lac objects and with an additional EUV excess in the case of Seyferts and quasars. Allowing for cold matter absorption in Seyfert galaxy hosts drive up the inferred average continuum slope to 2.0 +/- 0.5 (at 90% confidence), compared to a slope of 1.0 usually found from soft X-ray data. If Seyfert galaxies without EUV excesses form a significant fraction of the population, then the average spectrum of those with bumps should be even steeper. We place a conservative limit on neutral gas in BL Lac objects: N(sub H) less than 10(exp 20)/sq cm.

Marshall, H. L.↗

NASA NTRS1996

Non-Parabolic Hydrodynamic Formulations for the Simulation of Inhomogeneous Semiconductor Devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models cannot fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations or the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship ((hk)(exp 2)/2m = W(1 + alphaW). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(exp y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships; parabolic, Kane dispersion and power law dispersion.

Smith, A. W.↗

NASA NTRS1995

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Smith, Arlynn W.↗