SEARCH · Engineering Papers

Results for “Mixing”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Influence of mixing on the generation of nanoparticles in cement systems

In this work, we study the effect of mixing intensity on the size and concentration of nanoparticles in quartz and cement systems. We vary the local shear rate in the studied suspensions induced by mixing by either changing solid volume fraction or mixing speed. Our results obtained on quartz, cement and quartz-cement suspensions provide a master curve relating the concentration of nanoparticles in the pore solution as a function of the local shearing conditions computed at a micromechanical level. Our results suggest, in accordance with literature, that high mixing intensity is at the origin of a detachment of nano-hydrates from the surface of cement particles. These nanoparticles can be stabilized and therefore detected in the pore-solution only in the presence of superplasticizers, which inhibit their agglomeration to other particles.

36 MATERIALS SCIENCE↗

DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effect of the numerical discretization scheme in Shock-Driven turbulent mixing simulations

In this work, we evaluate the effects of distinct numerical strategies in simulations of shock-driven turbulent mixing. An air-SF 6 -air gas curtain subjected to Mach 1.2 shock-waves is computed with Implicit Large-Eddy Simulation based on three numerical schemes: directional-split with Harten-Lax-van Leer (HLL) solver; directional-unsplit using HLL-Contact (HLLC) solver; and directional-unsplit with HLLC solver and a Low Mach number Correction (LMC). The results illustrate the importance of the numerical strategy to the accuracy of the predictions. Whereas both split and unsplit schemes result in a similar spatial development of the initial shock-driven instability, only the unsplit schemes can predict the turbulent mixing transition after reshock observed by the laboratory experiments. Such feature increases the mixing rate of the two fluids, this being particularly pronounced when the LMC is active due to i) the reduced flow Mach number, and ii) the larger effective Reynolds number. Since the selected mixing problem is driven by the deposition of vorticity at the fluids’ interface, the resultant flow physics is analyzed by investigating the contribution of distinct inviscid mechanisms to the production of vorticity: baroclinicity, stretching, and dilatation. As expected, it is observed that the production of vorticity is initially dominated by the baroclinicity mechanism. Yet, the relevance of the remainder mechanisms is enhanced after the reshock and may even surpass the baroclinicity term.

42 ENGINEERING↗

DOE OSTI2024

Machine learning assisted bayesian inference of mix and hot-spot conditions in NIF implosions

Experiments on the National Ignition Facility (NIF) have provided clear evidence of ablator material mixing into the Hot-Spot, leading to degraded performance. However, inferring the amount of mix and Hot-Spot conditions from typical experimental observations (e.g. x-ray spectra and images) is highly challenging. Here, we have developed an analysis method that utilizes machine learning assisted Bayesian inference to find the probability distributions of the Hot-Spot and mix conditions. This approach uses a neural network, trained on an idealized 2-dimensional representation of the Hot-Spot and mix distribution, and Bayesian inference to find the statistical distributions of Hot-Spot conditions that provide a match with observations. We have tested this method with synthetic data from simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Improved representation of black carbon mixing structures suggests stronger direct radiative heating

Black carbon significantly influences the Earth system because of its strong solar radiation absorption. However, its direct radiative effect remains poorly understood in current climate models, partly because current climate models oversimplify the diverse structures formed when black carbon mixes with other atmospheric components. Here we show that incorporating more realistic, multi-mixing-structure representations of black carbon increases the direct radiative effect. We find that aged black carbon particles, with thicker coatings and higher embedded fractions, enhance the direct radiative effect more efficiently. Using machine learning alongside the Community Earth System Model, we show that the direct radiative effect at the top of the atmosphere in regions with heavy black carbon pollution is 31.6% greater when multi-mixing structures are considered. These findings highlight the importance of modeling complex mixing structures of particle-resolved black carbon to accurately capture their warming impacts on global atmosphere, particularly in highly polluted regions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Pure- and mixed-gas transport properties of a microporous Tröger's Base polymer (PIM-EA-TB)

Polymers of intrinsic microporosity (PIMs) offer tantalizing combinations of high selectivity and permeability in initial gas permeation measurements. Here, we report characterization of pure- and mixed-gas permeation properties of a thick (~80 μm) PIM consisting of Tröger's Base (TB) and Ethanoanthracene (EA) films. The effects of feed pressure and temperature on pure-gas permeabilities of CH 4 , N 2 , O 2 , H 2 , and CO 2 were investigated. The physical aging behavior of the thick film was tracked via pure-gas O 2 , N 2 , and CH 4 permeability at 35 °C. Gas permeability decreased noticeably and selectivity increased as aging time increased. Particular attention was given to mixed-gas measurements of CO 2 and CH 4 (50/50) permeabilities at 35 °C and fugacities ranging from 2 to 18 atm to explore whether the rigid, bridged, bicyclic TB and EA units could resist CO 2 -induced plasticization. These results are presented along with pure-gas CO 2 and CH 4 results for membrane samples aged at different times. PIM-EA-TB aged for ~24 h did not show signs of plasticization across the fugacity range considered. Furthermore, dual-mode competitive sorption presumably caused the CO 2 /CH 4 mixed-gas selectivity to be slightly higher than its corresponding pure-gas selectivity. However, as aging time increased, aged films underwent progressively more rapid and extensive CO 2 -induced plasticization with increasing fugacity, suggesting a systematic relationship between physical aging and plasticization in PIMs. Consequently, physical aging caused less improvement in mixed-gas CO 2 /CH 4 selectivity than it did in pure-gas selectivity, due mainly to plasticization effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Real-Time Monitoring of Gas-Phase and Dissolved CO 2 Using a Mixed-Matrix Composite Integrated Fiber Optic Sensor for Carbon Storage Application

Novel chemical sensors that improve detection and quantification of CO 2 are critical to ensuring safe and cost-effective monitoring of carbon storage sites. Fiber optic (FO) based chemical sensor systems are promising field-deployable systems for real-time monitoring of CO 2 in geological formations for long-range distributed sensing. Here, a mixed-matrix composite integrated FO sensor system was developed with a purely optical readout that reliably operates as a detector for gas-phase and dissolved CO 2 . A mixed-matrix composite sensor coating consisting of plasmonic nanocrystals and hydrophobic zeolite embedded in a polymer matrix was integrated on the FO sensor. The mixed-matrix composite FO sensor showed excellent reversibility/stability in a high humidity environment and sensitivity to gas-phase CO 2 over a large concentration range. This remarkable sensing performance was enabled by using plasmonic nanocrystals to significantly enhance the sensitivity and a hydrophobic zeolite to effectively mitigate interference from water vapor. The sensor exhibited the ability to sense CO 2 in the presence of other geologically relevant gases, which is of importance for applications in geological formations. A prototype FO sensor configuration which possesses a robust sensing capability for monitoring dissolved CO 2 in natural water was demonstrated. Reproducibility was confirmed over many cycles, both in a laboratory setting and in the field. More importantly, we demonstrated on-line monitoring capabilities with a wireless telemetry system, which transferred the data from the field to a website. The combination of outstanding CO 2 sensing properties and facile coating processability makes this mixed-matrix composite FO sensor a good candidate suitable for practical carbon storage applications.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Silver(I) Supported Liquid Membranes for Selective Ethylene Recovery from Mixed-Gas Streams of Tandem CO 2 Electrolysis

Large-scale olefin separations from unreacted paraffins and other byproduct gases are primarily done by energy-intensive cryogenic distillation processes at refineries. Silver(I) supported liquid membranes (Ag SLMs) can be implemented at smaller production scales of ethylene (C 2 H 4 ), a critical industrial chemical, such as its electrocatalytic (EC) production from CO 2 . Challenges of EC C 2 H 4 production mainly stem from reducing gases like hydrogen (H 2 ), where the redox reactions pertinent to Ag(I) facilitators diminish olefin transport. Herein we report that aqueous Ag(I) solution in a composite Ag SLM can operate in mixed-gas conditions containing H 2 gas utilizing reduced titania compounds, such as titanium(III) oxide (Ti 2 O 3 ). Embedding Ti 2 O 3 in the polydimethylsiloxane layer of Ag SLM assisted in selectively separating C 2 H 4 from mixed-gas feed streams related to CO 2 electrolysis containing H 2 . The direct exposure of a mixed-gas stream containing C 2 H 4 , CO 2 , CO, N 2 , and CH 4 with as high as 50 vol % H 2 maintained excellent C 2 H 4 separations for 7 days of continuous operation. The Ag SLM provided effective separations of C 2 H 4 from CO (at detection limits), CH 4 (selectivity ratio (α) = 20–30), and H 2 (α = ∼20), but C 2 H 4 from CO 2 (α = 2–4) revealed a slightly lower separation. These data show that aqueous Ag(I) solution(s) used in the SLMs can separate C 2 H 4 for extended periods, even under highly reducing gas conditions. Also, we report C 2 H 4 recovery from the gas mixture produced in the EC CO 2 reduction process. In conclusion, the Ag SLM gave C 2 H 4 selective separation from a five-component complex mixed-gas stream, relevant to tandem CO 2 electrolysis.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

DOE OSTI2024

Surface and Bulk Stabilization of Silicon Anodes with Mixed-Multivalent Additives: Ca(TFSI) 2 and Mg(TFSI) 2

Here, silicon is drawing attention as the upcoming anode material for the next generation of lithium-ion batteries due to its higher capacity compared to commercial graphite. However, silicon anions formed during lithiation are highly reactive with binder and electrolyte components creating an unstable SEI layer and limiting the calendar life of silicon anodes. The reactivity of lithium silicide and the formation of an unstable SEI layer is mitigated by utilizing the use of a mixture of Ca and Mg multivalent cations as an electrolyte additive for Si anodes to improve their calendar life. The effect of mixed salts on the bulk and surface of silicon anodes was studied by multiple structural characterization techniques. Ca and Mg ions in the electrolyte formed relatively thermodynamically stable quaternary Li-Ca-Mg-Si Zintl phases in an in-situ fashion and more stable and denser SEI layer on the Si particles. These in turn protect silicon particles against side reactions with electrolytes in a coin cell. The full cell with the mixed cation electrolyte demonstrates enhanced calendar life performance with lower measured current and current leakage than that of the baseline electrolyte due to reduced side reactions. Electron Microscopy, HRXRD, and solid-state NMR results showed that electrodes with mixed cations tended to have less cracking on the electrode surface compared to Si electrodes with Gen2 + FEC and the presence of mixed cations enhances cation migration and formation of quaternary Zintl phases stabilizing bulk and forming a more stable SEI.

25 ENERGY STORAGE↗

DOE OSTI2021

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Differentiable Neural Architecture, Mixed Precision and Accelerator Co-Search

Quantization, effective Neural Network architecture, and efficient accelerator hardware are three important design paradigms to maximize accuracy and efficiency. Mixed Precision Quantization is a process of assigning different precision to different Neural Network layers for optimized inference. Neural Architecture Search (NAS) is a process of automatically designing the neural network for a task and can also be extended to search for the precision of each weight and activation matrix. In this paper, we develop the following three methods: (i) Fast Differentiable Hardware-aware Mixed Precision Quantization Search method to find optimal precision, (ii) Joint Differentiable hardware-aware Architecture and Mixed Precision Quantization Co-search, (iii) Joint Accelerator, Architecture, and Precision triple co-search to find best possibilities in all the three worlds. We demonstrate the effectiveness of our proposed methods targeting Bitfusion accelerator by searching mixed precision models on MobilenetV2. We achieve better accuracy-latency trade-off models than the manually designed and previously proposed search methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗

DOE OSTI2020

Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part I: Developing Clouds

Evolution of nonprecipitating cumulus clouds (Cu) at the developing stage under the influence of lateral entrainment and mixing is studied analytically using a minimalistic analytical model. We present a model of an ascending cloud volume (a model of developing Cu) whose structure is determined by the processes of droplet diffusion growth/evaporation and entrainment mixing in the horizontal direction. Spatial and time changes of liquid water content, the adiabatic fraction, droplet concentration, and the mean volume droplet radius are calculated. It is shown that the existence of a nondiluted core in a growing cumulus cloud significantly depends on the cloud width and vertical velocity. While at the updraft velocity of 2 m s −1 the core of a 400-m-wide cloud becomes diluted at distances of a few hundred meters above cloud base, the core of a cloud of 1000-m width remains nondiluted at distances up to 1500 m above cloud base. The explanation of this result is simple: the increase in cloud width and the decrease in the updraft velocity increase the time during which the cloud is diluted due to mixing. Since lateral mixing synchronously decreases both the cloud water content and droplet concentration, the variation of the mean volume droplet radius is low inside the cloud. The approximate quantitative condition for cloud formation in updraft is derived. It is shown that a cloud can arise when its vertical velocity exceeds a critical value. To produce clouds, narrow turbulent plumes should ascend at higher velocity as compared to wider plumes. High humidity of the environment air is favorable for formation of clouds from plumes. The comparison of the obtained results with previously published observational data indicates a reasonable agreement. The results can be useful for parameterization purposes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

M-Star ® Modeling of SME Mixing with Three Impeller Blades

Before a sample of the Slurry Mix Evaporator (SME) can be taken, the SME product sampling procedure requires that the agitator power be stable between 20 and 30 kW for at least one hour. The SME transfer to Melter Feed Tank (MFT) procedure also requires the SME agitator power be stabilized between 20 and 30 kW prior to transfer. These requirements are specified to ensure samples are homogeneous, as discussed in the Waste Form Qualification Report. During SME Batch 804, the agitator power dropped to 19 kW and struggled to achieve and maintain 20 kW. It was later discovered that the cause of the power drop was due to one of the bottom blades of the agitator breaking off. Both the sample and the transfer occurred without the power stabilizing between 20 and 30 kW. Therefore, the quality of SME Batch 804 is indeterminate, and homogeneity was questionable. To evaluate mixing in the SME with a broken agitator containing only three blades (90 degrees apart with space for a missing blade) on the bottom impeller, the Savannah River National Laboratory (SRNL) was requested to perform computer simulations of impeller mixing in the SME using M-Star ® software. The simulations will be used to determine if SME Batch 804 was well mixed and satisfies homogeneity requirements.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Consumer Preferences and Willingness to Pay for Potting Mix with Biochar

Biochar is a co-product of advanced biofuels production from feedstocks including food, agricultural, wood wastes, or dedicated energy crops. Markets for soil amendments using biochar are emerging, but little is known about consumer preferences and willingness to pay (WTP) for these products or the depth of the products’ market potential for this product. This research provides WTP estimates for potting mix amended with 25% biochar, conditioned on consumer demographics and attitudes about product information labeling. Data were collected with an online survey of 577 Tennessee home gardeners. WTP was elicited through a referendum contingent valuation. Consumer WTP for an 8.81 L bag of 25% biochar potting mix is $8.52; a premium of $3.53 over conventional potting mix. Demographics and attitudes toward biofuels and the environment influence WTP. Biochar amounts demanded are projected for the study area’s potential market. Optimal prices, profits, and market shares are estimated across different marginal costs of producing biochar potting mix.

09 BIOMASS FUELS↗

DOE OSTI2021

General Quantum Field Theory of Flavor Mixing and Oscillations

We review the canonical transformation in quantum physics known as the Bogoliubov transformation and present its application to the general theory of quantum field mixing and oscillations with an arbitrary number of mixed particles with either boson or fermion statistics. The mixing relations for quantum states are derived directly from the definition of mixing for quantum fields and the unitary inequivalence of the Fock space of energy and flavor eigenstates is shown by a straightforward algebraic method. The time dynamics of the interacting fields is then explicitly solved and the flavor oscillation formulas are derived in a unified general formulation with emphasis on antiparticle content and effect introduced by nontrivial flavor vacuum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗