Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current conservation”

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 541 records · Page 30

Numerical modeling of an arcjet thruster

A numerical model of an arcjet thruster was developed and the results of calculations using this model are compared to the experimental data from a 30 kW arcjet using ammonia as a propellant. The model contains equations for the conservation of mass, radial, axial, and azimuthal momentum, energy, and the radial component of the magnetic field and provides a solution for three components of velocity, temperature, and the axial and radial current as a function of position in an axisymmetric flowfield. The model predicts 12 percent more specific impulse than the experiment when the mass flow and power input are constrained to the experimental values. A major conclusion of this study is that improved procedures for calculating transport properties are necessary if the accuracy of the model is to be improved.

Rhodes, Robert P.↗

A Prognostics Framework Development for Swarm Satellite Formations

Prognostics is the science of predicting the failure(s) of a component or a system and understanding how the performance will change in the event of a failure or degradation mechanism. With accurate predictions of possible failures, autonomous mitigative actions can be taken to correct/repair any issues or alert human operators of a failure threshold exceedance requiring condition-based maintenance. Although there is extensive research on failure predictions for a component or a system, there are significantly more opportunities to foray into failure predictions and prognostics for a system of systems such as an airspace consisting of multiple aircraft, a fleet of unmanned aerial vehicles, and a swarm of intelligent satellite systems. Failure prediction and mitigation are particularly important in autonomous systems such as satellite swarm systems that need effective resource management and minimal human interactions. Based on NASA's decadal survey, there is a clear need to prioritize the development of satellite swarm technology for studies of space physics and Earth science. The science community will propose future missions that return in-situ measurements from a 3-D (three-dimensional) volume of space, with relative spacecraft motion and inter-satellite baselines controlled according to the mission objectives. For such multi-spacecraft missions, it is required that ground operations resources do not scale with the number of satellites, thus compromising the swarm or leading to inefficiencies in resource allocation. Swarms of tens or hundreds of small satellites will require autonomy in attitude control, navigation and failure. Although significant research has been conducted in the areas of autonomous formation flying algorithms, less attention has been given to the development of resilient systems robust to failures.The focus of this research paper is the integration of model-based prognostics into the swarm dynamics control and decision-making algorithms. We simulate swarm management strategies for a subsystem failure to demonstrate the importance of failure predictions by comparing two cases: (i) no health information is provided to the system and utilized in the decision-making process and (2) system health information is obtained using prognostics and employed by the control system. One example scenario presented is for the GPS (Global Positioning System) system of an individual satellite to perform off-nominally due to increasing estimated error. In this scenario, the keep-out zone for that satellite would become more conservative, thereby decreasing the risk of collision. This is achieved via tuning the individual artificial repulsive functions assigned to each satellite.This paper is structured as follows. First we provide an overview of current swarm technology development, where we specifically use the term swarm to define multiple satellites flying in formation in similar orbits, with cross-link communication and station-keeping capabilities. Second, we give an introduction to the Swarm Orbital Dynamics Advisor (SODA), a tool that accepts high-level configuration commands and provides the orbital maneuvers required to achieve the prescribed formation configuration. Third, we provide the details of the model-based prognostics algorithm implementation in SODA. Finally, we present different case studies for potential component/subsystem failures and the swarm responses based with and without failure prediction information.

prognostics↗

The Atmosphere-Space Transition Region Explorer (ASTRE) – A Low Perigee Satellite to Investigate the Coupling of the Earth’s Upper Atmosphere and Magnetosphere

The Atmosphere-Space Transition Region Explorer (ASTRE) is a mission concept designed to carry out an unprecedented study of the interaction between the Earth’s atmosphere and the ionized gases of space within the atmosphere-space transition region. By gathering direct measurements of the coupling of ion and neutral gases in this region, ASTRE provides the critical missing link in our knowledge of the transfer, dissipation, and regulation of energy and momentum between the sun and the upper atmosphere. ASTRE provides the first detailed, systematic investigation of this important unexplored region, vastly improves and constrains models of the upper atmosphere, and fills a critical gap in our understanding of how the coupled lower ionosphere/upper atmosphere “works” as a system. To achieve its science objectives, ASTRE gathers accurate measurements of plasma and neutral gases, electric and magnetic fields, and energetic particles using well-proven, in situ instruments with excellent flight heritage. Furthermore, as described herein, all the instruments have been designed to perform well in the low perigee environment, including altitudes of 150 km and lower. ASTRE utilizes a three-axis stabilized satellite that uses on-board propulsion to carry out systematic, low perigee measurements at high latitudes. The satellite design incorporates a “form follows function” approach with a cylindrical shape and conductive body-mounted solar arrays to minimize drag and perturbations to the space environment. Atomic-oxygen resistant materials are utilized and a passive thermal design with heat pipes and radiator panels minimizes the impact of aero-heating. ASTRE launches into a 250 km × 1500 km elliptical insertion orbit with an inclination of 83°. Perigee precesses from its highest northern latitude to its highest southern latitude every ~60 days. Hydrazine propulsion provides over 2000 orbits with perigee below 200 km, with a significant fraction as low as 150 km, during high latitude, two-week campaigns when perigee precesses to either the northern or southern high latitude region. At mid and low latitudes, the perigee is near 225 km. Because very conservative assumptions were made with respect to the orbital analysis and drag, including continuous 3-sigma “worst-case” solar flux and atmospheric density, when additional propellant is included (accommodated in the current design) and some of the stringent assumptions are relaxed, the ASTRE mission described herein may be expected to include repeated excursions to altitudes of ~130 km or even lower. This paper presents an overview of the ASTRE mission, its science motivation, and objectives. It includes a discussion of the science-driven requirements and traceability, followed by a “proof-of concept” implementation that includes notional instruments and a straightforward spacecraft design. Three key points are demonstrated: 1) There is a critical knowledge gap in the high latitude, atmosphere-space transition region below 250 km; 2) The instrument and measurement techniques needed to obtain the ASTRE measurements are well-proven and function well in the low-altitude environment; and 3) A mature spacecraft design, flight dynamics analysis, and concept of operations have been developed that demonstrate that the ASTRE mission can be achieved in a straightforward manner using current technologies.

Ionosphere↗

The evolution of rotating stars. III - Predicted surface rotation velocities for stars which conserve total angular momentum

Predicted surface rotation velocities for Population I stars at 10, 7, 5, 3, and 1.5 solar masses are presented. The surface velocities were computed for angular momentum with no radial redistribution, complete redistribution, and partial redistribution as predicted by consideration of circulation currents in rotating stars. Near the main sequence, rotational effects can reduce the moment of inertia of a star, so nonrotating models underestimate the expected velocities for evolving stars. On the red giant branch, angular momentum redistribution reduces the surface velocity by a factor of 2 or more, relative to the velocity expected for no radial redistribution. This removes the discrepancy between predicted and observed rotation rates for the K giants and makes it unlikely that these stars lose significant amounts of angular momentum by stellar winds. Calculations indicate that improved observations of the red giants in the Hyades cluster can be used to determine how angular momentum is redistributed by convection

Endal, A. S.↗

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗

Overview of the SPARC tokamak

The SPARC tokamak is a critical next step towards commercial fusion energy. SPARC is designed as a high-field ( $B_0 = 12.2$ T), compact ( $R_0 = 1.85$ m, $a = 0.57$ m), superconducting, D-T tokamak with the goal of producing fusion gain $Q>2$ from a magnetically confined fusion plasma for the first time. Currently under design, SPARC will continue the high-field path of the Alcator series of tokamaks, utilizing new magnets based on rare earth barium copper oxide high-temperature superconductors to achieve high performance in a compact device. The goal of $Q>2$ is achievable with conservative physics assumptions ( $H_{98,y2} = 0.7$ ) and, with the nominal assumption of $H_{98,y2} = 1$ , SPARC is projected to attain $Q \approx 11$ and $P_{\textrm {fusion}} \approx 140$ MW. SPARC will therefore constitute a unique platform for burning plasma physics research with high density ( $\langle n_{e} \rangle \approx 3 \times 10^{20}\ \textrm {m}^{-3}$ ), high temperature ( $\langle T_e \rangle \approx 7$ keV) and high power density ( $P_{\textrm {fusion}}/V_{\textrm {plasma}} \approx 7\ \textrm {MW}\,\textrm {m}^{-3}$ ) relevant to fusion power plants. SPARC's place in the path to commercial fusion energy, its parameters and the current status of SPARC design work are presented. This work also describes the basis for global performance projections and summarizes some of the physics analysis that is presented in greater detail in the companion articles of this collection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Intermittent Turbulence and SOC Dynamics in a 2-D Driven Current-Sheet Model

Borovsky et al. have shown that Earth's magnetotail plasma sheet is strongly turbulent. More recently, Borovsky and Funsten have shown that eddy turbulence dominates and have suggested that the eddy turbulence is driven by fast flows that act as jets in the plasma. Through basic considerations of energy and magnetic flux conservation, these fast flows are thought to be localized to small portions of the total plasma sheet and to be generated by magnetic flux reconnection that is similarly localized. Angelopoulos et al., using single spacecraft Geotail data, have shown that the plasma sheet turbulence exhibits signs of intermittence and Weygand et al., using four spacecraft Cluster data, have confirmed and expanded on this conclusion. Uritsky et al., using Polar UVI image data, have shown that the evolution of bright, nightside, UV auroral emission regions is consistent with many of the properties of systems in self-organized criticality (SOC). Klimas et al. have suggested that the auroral dynamics is a reflection of the dynamics of the fast flows in the plasma. sheet. Their hypothesis is that the transport of magnetic fludenergy through the magnetotail is enabled by scale-free avalanches of localized reconnection whose SOC dynamics are reflected in the auroral UV emission dynamics. A corollary of this hypothesis is that the strong, intermittent, eddy turbulence of the plasma sheet is closely related to its critical dynamics. The question then arises: Can in situ evidence for the SOC dynamics be found in the properties of the plasma sheet turbulence? A 2-dimensional numerical driven current-sheet model of the central plasma sheet has been developed that incorporates an idealized current-driven instability with a resistive MHD system. It has been shown that the model can evolve into SOC in a physically relevant parameter regime. Initial results from a study of intermittent turbulence in this model and the relationship of this turbulence to the model's known SOC dynamics will be discussed.

Klimas, A. J.↗

Accurately simulating core-collapse self-interacting dark matter halos

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations have motivated current efforts to explain surprisingly compact DM halos. If DM is not collisionless, but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we sought to better understand the numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially, our aim was to better understand the evolution of satellite halos. To do so, we ran SIDM N-body simulations of a low-mass halo in isolation and within an external gravitational potential. The simulation set-up was motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a SIDM kernel size that is too large can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our most highly resolved simulation (N = 5 × 10 7 ) we have made the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, for example that the energy conservation error is better kept well below 1%. This allows simulations to be run of halos that become concentrated enough to explain observations of GD-1-like stellar streams or strong gravitational lensing systems.

dark matter↗

General aviation energy-conservation research programs

A review is presented of nonturbine general aviation engine programs underway at the NASA-Lewis Research Center. The program encompasses conventional, lightweight diesel, and rotary engines. Its three major thrusts are: (1) reduced SFC's; (2) improved fuels tolerance; and (3) reducing emissions. Current and planned future programs in such areas as lean operation, improved fuel management, advanced cooling techniques, and advanced engine concepts, are described. These are expected to lay the technology base, by the mid to latter 1980's, for engines whose total fuel costs are as much as 30% lower than today's conventional engines.

Willis, E. A.↗

Whole-genome demography of COVID-19 virus during its pandemic period and on “panvalent” vaccine design

With over 16 million submitted genomic sequences, the SARS-CoV-2 (SC2) virus, the cause of the most recent worldwide COVID-19 pandemic, has become the most sequenced genome of all known viruses, revealing, for example, a vast number of expanding viral lineages. Since the pandemic phase appears to be over, we performed a retrospective re-examination of the demographic grouping pattern and their genomic characteristics during the entire pandemic period up to the peak of the last pandemic wave. For our study, we extracted from the NCBI only unique viral sequences and converted each sequence data to a relational vector, indicating the presence/absence of each variational event compared to a “reference” sequence. Our study revealed several genomic features that are unexpected or different from those of previous studies. For example, approximately 44,000 variants with unique sequences emerged during the pandemic period; they group into only four major viral-genomic groups and each has a set of mostly unique highly-conserved variant-genotypes (HCVGs); and a small set from the first (“ancestral”) group was inherited by the three (“descendant”) groups, suggesting that HCVGs in the next group may be predictable from the current group(s). Such a concept may be potentially important in designing “panvalent” vaccines against the current and future waves of viral infections.

60 APPLIED LIFE SCIENCES↗

Xyce(™) Parallel Electronic Simulator v.7.5

The Xyce Parallel Electronic Simulator simulates electronic circuit behavior in DC, AC, HB, MPDE and transient mode using standard analog (DAE) and/or device (PDE) device models including several age and radiation aware devices. It supports a variety of computing platforms (both serial and parallel) computers. Lastly, it uses a variety of modern solution algorithms dynamic parallel load-balancing and iterative solvers.! ! Xyce is primarily used to simulate the voltage and current behavior of a circuit network (a network of electronic devices connected via a conductive network). As a tool, it is mainly used for the design and analysis of electronic circuits.! ! Kirchoff's conservation laws are enforced over a network using modified nodal analysis. This results in a set of differential algebraic equations (DAEs). The resulting nonlinear problem is solved iteratively using a fully coupled Newton method, which in turn results in a linear system that is solved by either a standard sparse-direct solver or iteratively using Trilinos linear solver packages, also developed at Sandia National Laboratories.

Source record↗

Dialogue Between Lighting and HVAC Systems: Improving building system integration

Lighting systems have long been capable of sensing when someone enters or exits a room and using that knowledge to turn lights on or off. More recently, connected lighting systems with sensors integrated into every luminaire have become broadly available, facilitating highly granular occupancy detection. Similarly, HVAC systems have long been able to use an understanding of building occupancy to adjust temperature setpoints and reduce energy use without significant impacts to occupant comfort. Energy codes (e.g., ASHRAE/IES Standard 90.1, IECC, Title 24) now require “occupied standby HVAC control,” whereby systems adjust both temperature and ventilation setpoints in zones that are determined to be unoccupied during normal occupancy hours. Here, this article discusses current issues that stymie the integration of Lighting and HVAC systems, and DOE activities focused on addressing them.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reactor Pressure Vessel Fracture Mechanics Development and Concrete Application Testing for Grizzly

The Grizzly code is being developed to address degradation issues in nuclear reactor structures and components. For light-water reactors, Grizzly currently has capabilities to simulate degradation processes and their effects on structural integrity in two key areas: reactor pressure vessels (RPVs) and reinforced concrete structures. This report documents improvements made to Grizzly’s ability to address both of these structural systems. For RPVs, the reduced-order models (ROMs) used in fracture mechanics calculations have been expanded to allow their application over a broader range of the parameter space than was permitted by the previous models. The ROMs currently used in Grizzly for the evaluation of flaws that are fully embedded within the RPV (as opposed to surface-breaking flaws) are based on a model that is known to be conservative, indicating higher stress intensity factors than would be obtained from direct simulations. A more accurate model that eliminates these excess conservatisms has been recently included in the American Society of Mechanical Engineers Boiler and Pressure Vessel Code but was not applicable for flaws near the RPV surface, which is where the most critical flaws are usually located. That model has recently been extended for increased applicability in this near-surface region. The ROMs for embedded flaws in the Grizzly code have been expanded to include these recent extensions, which permit their use in a much broader set of cases than previously possible. Direct 3D simulations have been used to check these ROMs and have shown good agreement in most cases, although there are still some cases that need further investigation. There are considerable benefits to using these these more accurate and less conservative ROMs for embedded flaws. On a benchmark probabilistic fracture mechanics problem tested here, the conditional probability of fracture initiation computed for a population of flaws in a single plate in an RPV decreased by over a factor of 3. To address aging in reinforced concrete structures, a capability to simulate multiple degradation mechanisms, including alkali-silica reaction and radiation-induced volumetric expansion has been developed in Grizzly over the past several years. This had previously been demonstrated on laboratory-scale specimens but not on full-scale nuclear concrete structures with reinforcement. To demonstrate the applicability of Grizzly to the analysis of large-scale structures of interest, a full 3D model of a representative reinforced concrete structure, including a complex arrangement of reinforcing bars, was developed and demonstrated in Grizzly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computer Science Research Needs for Parallel Discrete Event Simulation (PDES)

Historically, scientific computing efforts have demonstrated the clear need for, and effective use of, supercomputing with traditional time-stepped simulations. Nevertheless, there are several areas in the mission spaces of the U.S. Department of Energy and other agencies waiting to tap advanced computing research using a different, discrete event style of modeling, simulation, and analysis. These span a wide spectrum of applications including energy grid resilience, urban planning and policy, transportation science, building technologies, emergency response and planning, environmental impact analysis, computational epidemiology, Internet communications, cyber security, and cyber-physical systems, to name only a few. Even within traditional scientific applications, the role of discrete event modes of execution is increasing in the form of new event-based mathematical solvers such as quantized state integration methods and discrete-continuous hybrid system solvers. Co-design of advanced supercomputing hardware systems is another area that exploits discrete event simulation at its core for effective analyses. Complex systems, entity behaviors and interconnections play a significant role in all these applications, which are mapped to large-scale models with discrete event formulations. To make advancements in all the aforementioned scientific areas, many technical aspects need to be more thoroughly studied and deeply understood in parallel discrete event simulation (PDES). The unique dynamics inherent in a discrete event modeling approach, by their very nature, intersect and influence the entire stack of the computing system, including (a) the unique nature of the instruction sets exercised in PDES workloads without a predominance of high-precision floating point operations, (b) virtual time-constrained multi-threaded execution of many logical processes per processor, (c) extremely variable and difficult to predict network traffic characteristics, (d) interfaces and inter-dependencies with machine learning and artificial intelligence codes at higher software layers, and (e) highly challenging load balancing needs, especially in effectively accounting for accelerated/extremely heterogeneous computing in current and future high-performance computing systems. Efficient and accurate parallel execution of PDES workloads is also dominated by challenges in dealing with their asynchronous concurrency fundamentally present at the model level. Conservative synchronization, optimistic/speculative synchronization, and their hybrid schemes open new questions in fundamental computer science with respect to reversibility of computation and prediction (lookahead) of behaviors inherent within model codes. On the implementation front, there are relatively few scalable, general-purpose parallel discrete event simulators in the world, and even fewer have been studied on emerging hardware platforms. To enable scientific advances using PDES, the research needs in computer science must also be pursued and met in the intersection of the algorithmic and hardware-aware aspects of scalable PDES engines. This report is aimed at capturing a computer science-oriented view of this important area of research in PDES, presenting a sample of important applications with their inherent discrete event technology elements. Needs are outlined in core areas of parallel discrete event research as well as cross-cutting directions in computer science research that positively impact scientific advancements across several important application areas. A selection of priority research opportunities in advanced computing for PDES is identified to serve as reference for key research topics and their order of importance for scientific advancements.

97 MATHEMATICS AND COMPUTING↗

Effects of Borate and Organics on U(VI) Solubility in WIPP Brine

The solubility of uranium (VI) in Waste Isolation Pilot Plant (WIPP)-relevant brine was determined to support ongoing WIPP recertification activities. This research was performed by the Los Alamos National Laboratory Carlsbad Operations (LANL-CO) Actinide Chemistry and Repository Science Program (ACRSP). The WIPP Actinide Source Term Program (ASTP) did not develop a model for the solubility of actinides in the VI oxidation state. The solubility of UO 2 2+ , in the absence of WIPP specific data, is presently set to be equal to a conservatively high 1 mM within the WIPP Performance Assessment (PA) for all expected WIPP conditions (SOTERM, 2019) as selected at the recommendation of the Environment Protection Agency (EPA) (EPA, 2005). According to the current WIPP chemistry model assumptions and conditions, the expected pC H+ is about 9.5 and controlled by MgO buffering CO 3 2- . The goal of this study is to perform screening experiments that account for the contributions of organics and borate on uranium solubility. In this report, the solubility of U(VI) was determined at pC H+ 9 WIPP brine in the absence or presence of borate and organics at under-saturation approach. Experiments were equilibrated for about 135 days. Organic compounds present in WIPP waste can form strong complexes with actinides and can affect the oxidation states of actinides. The organic compounds addressed in WIPP performance assessment include EDTA (Ethylenediaminetetraacetic Acid), oxalate, citrate, and acetate (SOTERM, 2019). These data quantify the effects of WIPP-relevant concentrations of borate and organics effects on the solubility of U(VI) to challenge the predictions of the WIPP actinide model and inform decisions and recommendations made in the upcoming recertification of the WIPP (CRA-2024). The experiments performed were done according to the U.S. Department of Energy (DOE) approved Test Plan entitled “Experimental Strategy to Challenge Actinide Solubility Predictions” and designated LCO-ACP-26. All data reported were obtained under the LANL-CO Quality Assurance Program, which is compliant with the DOE Carlsbad Field Office, Quality Assurance Program Document (CBFO/QAPD) (QAPD, 2017).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Buckling of thin-walled circular cylinders

This monograph indicates current practices for predicting buckling of uniform stiffened and unstiffened circular cylindrical shells under various types of static loading, and suggests the procedures that yield estimates of static buckling loads considered to be conservative. The buckling of truncated conical shells and shells of double curvature will be treated in separate monographs.

V. I. Weingarten↗

Wave propagation in a quasi-chemical equilibrium plasma

Wave propagation in a quasi-chemical equilibrium plasma is studied. The plasma is infinite and without external fields. The chemical reactions are assumed to result from the ionization and recombination processes. When the gas is near equilibrium, the dominant role describing the evolution of a reacting plasma is played by the global conservation equations. These equations are first derived and then used to study the small amplitude wave motion for a near-equilibrium situation. Nontrivial damping effects have been obtained by including the conduction current terms.

Fang, T.-M.↗