Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flows”

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 901 records · Page 50

Development of a Multiphase Adjoint Capability in OpenFOAM

Thermal hydraulic predictions of flow regimes within nuclear reactor cores require significant computational resources and accuracy to ensure that the core design does not violate the thermal limits of materials. Multiphase computational fluid dynamics (CFD) has become a valuable computational capability for predicting flow characteristics as they relate to applications such as nuclear reactor safety and mixing phenomena around spacer grids within fuel assemblies. Multiphase CFD requires closure relationships for both the effects of turbulence as well as the transfer of mass, momentum, and energy between phases. Understanding the sensitivity of the parameters used in closure relationships is important for improving accuracy, but performing rigorous sensitivity analysis is computationally expensive. Adjoint mathematics have been utilized for reactor physics and fluid dynamics applications with respect to verification, sensitivity analysis, adaptive mesh refinement, and optimization. The usage of adjoint methods to predict the sensitivities of a variety of closure relationship parameters within a given problem has been shown to return sensitivity responses with significantly fewer computations than normal 'large-scale' methods. As the complexity of multiphase CFD methods increases, the demand for accurate verification tools and improvements in computational efficiency also increases, making the development of adjoint capabilities an attractive option. Closure relationships that describe interphase momentum physics in multiphase CFD modelling rely heavily on prescribed constants determined from experimental data. Due to the computational efficiency of an adjoint approach, the time to derive these closure relationships could be dramatically reduced, and data assimilation techniques could be used to return probability density functions of closure models derived using experimental data. This work develops a multiphase adjoint capability for the Eulerian-Eulerian two field, two phase equations. Specifically, the adjoint capability has been developed for the sub-cooled nucleate boiling solver boilEulerFoam developed by Dr. Nam Dinh of NC State University within the open-source CFD library OpenFOAM.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Numerical studies of dilute and dense spray characteristics

Several issues involving the improvement of physical submodels and the computational efficiency in modeling dilute and dense spray combustion are discussed. First, the implementations of a dispersion width approach accounting for turbulent dispersion within each computational parcel is discussed. This is essentially a statistical transport model and the testings of this model confirm the capability of accurately representing dispersion in nearly-homogeneous and inhomogeneous turbulent flows with improved efficiency over the delta function stochastic separated flow model. To account for the dense spray effects, an existing drop collision and coalescence model and a Taylor analogy breakup (TAB) model were employed. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Several examples including nonevaporating, evaporating, and burning dense spray cases were studied. The numerical results show reasonably good comparisons with available experimental data in terms of spray penetration, drop sizes, and overall configuration of a spray flame.

Kim, Y. M.↗

Studies on fuel spray characteristics in high-pressure environment

The present study deals with several issues involving the improvement of physical submodels and the computational efficiency in modeling dense fuel sprays. To improve the computational efficiency, a parcel PDF approach is implemented which can account for turbulent dispersion within each computational parcel. The advantage of a parcel PDF tracking method is to reduce the number of computational parcels representing the spray dynamics as well as to obtain grid-independent solutions for two-phase flows. To account for the dense spray effects, an existing drop collision and coalescence model, two breakup models, and a Reitz's wave instability model were used. These models were incorporated into a state-of-the-art multiphase all-speed transient flow solution procedure. Comparative performance for two breakup models as well as the turbulence modulation effects are also studied. Validation cases include the nonevaporating and evaporating solid-cone dense sprays. The predictions show a reasonably good agreement with available experimental results in terms of spray penetration, drop sizes, gas and drop mean velocities, and gas and drop rms velocities. The numerical results indicate that the present parcel PDF model has the capability of accurately representing drop dispersion in dense sprays with manageable number of computational parcels.

Shang, H. M.↗

FEARCE: Fast, Easy, Accurate, and Robust Continuum Engineering Improving fuel efficiency and reducing emissions in combustion engines 2019 R&D 100 Awards

FEARCE investigates and simulates turbulent chemically reactive flows in fuel combustion, such as fuel injection and fuel-air mixing, thermodynamic combustion losses, and combustion–emission formation processes in engines. Our objective is to understand turbulent flows and how they react and interact with other complicating variables, including highly nonlinear heat and mass transfer, small-scale velocity, and chemical kinetics or reactions. Our emphasis is on the fluid dynamics and chemical reactions that take place in all kinds of internal combustion engines, including the reactions between multiple species (such as air, gasoline, carbon dioxide, water, and hundreds of other species) and multiphase (liquid and gas) fluid flow.

42 ENGINEERING↗

Modeling Three-Dimensional Fluid-Driven Propagation of Multiple Fractures using TOUGH-FEMM

Here, a novel numerical simulation tool named TOUGH-FEMM is presented and applied to model three-dimensional (3D) hydraulic fracturing in porous rock. The fluid flow in both fractures and porous rock is modeled using TOUGH2, which is a well-established code for analysis of multiphase and multi-component fluid flow. Rock deformations associated with fracture propagation are modeled using finite element-meshfree method (FEMM). FEMM is an approach to simulate fracture propagation without remeshing, in which the fracture path does not need to be predetermined. Fracture mechanics with mixed-mode stress intensity factors are employed to detect fracture instability and determine the direction of fracture propagation. TOUGH-FEMM is confirmed for modeling fluid-driven fracture propagation in 3D through a number of simulation examples, including modeling of hydraulic fracturing laboratory experiments and by comparison to independent numerical simulation results for multiple interacting hydraulic fractures at ten to hundred meter scale.

58 GEOSCIENCES↗

Robust hierarchical dispatch for residential distribution network management considering home thermal flexibility and model predictive control

In the transactive energy (TE) paradigm, the devices of participative consumers, or prosumers, may be aggregated and employed to drive operational objectives at the network level. Home heating, ventilation and air-conditioning (HVAC) systems in particular are well-suited to modulate their behaviours based on both home thermal flexibility and requests from the utility grid. This paper develops a robust, hierarchical power dispatch scheme in the context of a residential distribution network. The formulation couples a unique, multiphase linear distribution optimal power flow (OPF) at the upper level with model predictive control (MPC)-based HVAC fleet controllers at the lower level. The proposed approach is tested on nearly 2000 homes with a three-phase distribution network in an intraday market setting, where two major applications are explored and analysed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Novel Membrane-Based Corrosion Sensors For Operating In Natural Gas Pipelines

Membrane-based electrochemical sensors (MBES) provide an attractive option for real-time continuous humidity and corrosion rate measurements inside of natural gas pipelines. The MBES, can use electrochemical techniques such as linear polarization resistance and electrochemical impedance spectroscopy even when a conductive aqueous phase is intermittent or unavailable. We now present field-test performance data of our MBES operated inside a multiphase water and natural gas flow facility. The MBES measured internal corrosion and conductivity at the top, bottom, and side of the pipe spool. Tests were carried out with flow rates of 4.6, 9.1, and 18.3 m/s, temperatures of 21, 29, and 35 °C, and a pressure of 1000 psia. All sensors followed similar trends in open circuit potential, corrosion rate, and conductivity over time. Neither flow rate, temperature, nor build-up of metal-based solids on the sensor surfaces impacted the quality of the electrochemical measurements. The field-test results confirm that the membrane-based electrochemical sensors are suitable for monitoring internal corrosion of operational natural-gas pipelines.

Duffy, Timothy↗

DNS of turbulent flow with hemispherical wall roughness

The present study of the effect of roughness density on the mean flow turbulence parameters is motivated by the need for new generation of boundary conditions for multiphase computational multiphase fluid dynamics (CMFD) models applied to boiling flows. Effect of roughness element density on the turbulent flow in a channel is quantified through direct numerical simulations (DNSs). The Navier--Stokes equations are solved using finite element method and bubbles are approximated as rigid near-hemispherical obstacles at the wall. Six different cases were analysed including channel flow with smooth wall and channel flow with rough wall for five different bubble nucleation site densities. Friction factor and the law of the wall was calculated and compared with the previously published results. Existing correlations for nucleating bubble site density dependency on a wall heat flux were used to obtain a relation between the heat flux and the friction factor, leading to the law of the wall dependency on the heat flux. This separate effect study provides new guidelines on how the heat flux in subcooled boiling regime affects the turbulence behaviour near the wall and guides the computational fluid dynamics model development for boiling two-phase flows.

Mechanics↗

Entropic multiple-relaxation-time multirange pseudopotential lattice Boltzmann model for two-phase flow

Here, an entropic multiple-relaxation-time lattice Boltzmann approach is coupled to a multirange Shan-Chen pseudopotential model to study the two-phase flow. Compared with previous multiple-relaxation-time multiphase models, this model is stable and accurate for the simulation of a two-phase flow in a much wider range of viscosity and surface tension at a high liquid-vapor density ratio. A stationary droplet surrounded by equilibrium vapor is first simulated to validate this model using the coexistence curve and Laplace’s law. Then, two series of droplet impact behavior, on a liquid film and a flat surface, are simulated in comparison with theoretical or experimental results. Droplet impact on a liquid film is simulated for different Reynolds numbers at high Weber numbers. With the increase of the Sommerfeld parameter, onset of splashing is observed and multiple secondary droplets occur. The droplet spreading ratio agrees well with the square root of time law and is found to be independent of Reynolds number. Moreover, shapes of simulated droplets impacting hydrophilic and superhydrophobic flat surfaces show good agreement with experimental observations through the entire dynamic process. The maximum spreading ratio of a droplet impacting the superhydrophobic flat surface is studied for a large range of Weber numbers. Results show that the rescaled maximum spreading ratios are in good agreement with a universal scaling law. This series of simulations demonstrates that the proposed model accurately captures the complex fluid-fluid and fluid-solid interfacial physical processes for a wide range of Reynolds and Weber numbers at high density ratios.

42 ENGINEERING↗

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗

Temperature Homogeneity under Selective and Localized Microwave Heating in Structured Flow Reactors

Selective heating of different phases of multiphase systems via microwaves can result in energy savings and suppression of side reactions. However, materials properties and operating conditions that maximize temperature gradients are poorly understood. Here we utilize computational fluid dynamics (CFD) computations and temperature measurements in structured flow reactors (monoliths) in a monomodal microwave cavity to assess the temperature difference between the walls and the fluid and develop a simple lumped model to estimate when temperature gradients exist. We also explore the material’s thermal and electrical properties of structured reactors for isothermal catalyst conditions. We propose that CFD simulations can be used as a nonintrusive, predictive tool of temperature homogeneity. Importantly, we demonstrate that localized heating in the bed under several conditions rather than selective heating is responsible for the selectivity enhancement. Furthermore, our results indicate that structured beds made of high thermal conductivity materials avoid arcing and enable temperature homogeneity and low electrical conductivity materials allow microwaves to penetrate the domain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A smoothed particle hydrodynamics model for droplet and film flow on smooth and rough fracture surfaces

Flow on fracture surfaces has been identified by many authors as an important flow process in unsaturated fractured rock formations. Given the complexity of flow dynamics on such small scales, robust numerical methods have to be employed in order to capture the highly dynamic interfaces and flow intermittency. In this work we present microscale free-surface flow simulations using a three-dimensional multiphase Smoothed Particle Hydrodynamics (SPH) code. Pairwise solid-fluid and fluid-fluid interaction forces are used to control the wetting behavior and cover a wide range of static and transient contact angles as well as Reynolds numbers encountered in droplet flow on rock surfaces. We validate our model via comparison with existing empirical and semi-analyical solutions for droplet flow. We use the model to investigate the occurence of adsorbed trailing films of droplets under various flow conditions and its importance for the flow dynamics when films and droplets coexist. We show that flow velocities are higher on prewetted surfaces covered by a thin film which is qualitatively attributed to the enhanced dynamic wetting and dewetting at the trailing and advancing contact line.

unsaturated flow, fractured rocks, smoothed partic↗

Preface to the Special Issue on TOUGH Symposium 2015

The TOUGH Symposium 2015 was held in Berkeley, California, September 28-30, 2015. The TOUGH family of codes, developed at the Energy Geosciences Division of Lawrence Berkeley National Laboratory (LBNL), is a suite of computer programs for the simulation of multiphase and multicomponent fluid and heat flows in porous and fractured media with applications in many geosciences fields, such as geothermal reservoir engineering, nuclear waste disposal, geological carbon sequestration, oil and gas reservoirs, gas hydrate research, vadose zone hydrology and environmental remediation. Since the first release in the 1980s, many modifications and enhancements have been continuously made to TOUGH and its various descendants (iTOUGH2, TOUGH+, TOUGH-MP, TOUGHREACT, TOUGH+HYDRATE, TMVOC...), at LBNL and elsewhere. Today, these codes are used worldwide in academia, government organizations and private companies in problems involving coupled hydrological, thermal, biogeochemical and geomechanical processes. The Symposia, organized every 2-3 years, bring together developers and users for an open exchange on recent code enhancements and applications. In 2015, the Symposium was attended by one hundred participants, representing thirty-four nationalities. This Special Issue in Computers & Geosciences gathers extended versions of selected Symposium proceedings related to (i) recent enhancements to the TOUGH family of codes and (ii) coupled flow and geomechanics processes modeling. In this special issue, a paper by Jung et al. provides a new base version of TOUGH, TOUGH3. This latest version addresses the computational robustness and efficiency needed for the increasing complexity of the simulated processes and the growing size of model domains. Moreover, it consolidates the serial and parallel versions into a single code, hence, inconsistencies between different versions are less likely to occur, and code maintenance and development become more efficient. Finsterle et al. present a review of the capabilities of iTOUGH2 and illustrate the workflow of an integrated simulation-optimization framework. Liu et al. present an approach to reduce the computational cost of performing Bayesian inverse modeling and demonstrate their approach using a synthetic experiment simulated with TOUGH2. Regarding the mesh generation process, the paper by Doughty presents a generalization of the grid-generation algorithm embedded in TOUGH2 (MESHMAKER) to create one-column grids with fractal dimension, useful for representing flow through fractures or highly heterogeneous media with fractal geometry, and also for inverse methods. Sentís and Gable present a module to interface the open source code Lagrit with TOUGH2, that assists in the processes of mesh generation and model setup, and that can be particularly useful for complex geometries and large models. Bonduà et al. present a pre-processing tool (VORO2MESH) that features a full 3D Voronoi approach for structured and unstructured grids, and an enhanced post-processor (TOUGH2Viewer) for the visualization of 3D Voronoi grids. The literature review shows that coupled fluid flow and geomechanics has gained increasing attention in the last years in many geosciences fields. In this context, several codes of the TOUGH family have been successfully coupled with geomechanics software. The paper by Rutqvist presents a review of some of these TOUGH-based geomechanics models, of which the TOUGH-FLAC simulator – presented for the first time in the TOUGH Symposium 2003 – is well-known worldwide and has been used in a large variety of applications, encompassing the application domains of the TOUGH codes. Blanco-Martín et al. present an extension of TOUGH-FLAC to the finite strain framework, motivated by the existence of large strains associated with rock salt deformation in the context of heat-generating nuclear waste disposal. Following a paper published in the Special Issue of the TOUGH Symposium 2012 in this journal, Kim et al. present the current status of the TOUGH-RSBN simulator for coupled hydro-mechanical modeling of hydraulic fracturing using a discrete fracture network approach, and present some validation and demonstration cases. Rinaldi and Nespoli present TOUGH-SEED, an injection-induced seismicity modeling approach using TOUGH2 and a geomechanics stochastic model, and apply it to the Basel EGS project. Finally, Rinaldi et al. use the PEST protocol implemented in iTOUGH2, combined with TOUGH-FLAC as the forward model, to perform inverse modeling for the In Salah CO2 demonstration site in Algeria. Using pore pressures and InSAR displacements as observations, this approach allows for a multi-physics joint inversion and takes advantage of the inverse capabilities of iTOUGH2 (parameter estimation, sensitivity analysis, uncertainty estimation …), applied to coupled fluid flow and geomechanics. On behalf of the members of the TOUGH Symposium 2015 Organizing Committee, I hope that this collection of papers will provide the readers of Computers & Geosciences with a recent view of the capabilities of the TOUGH family of codes, as well as with current modeling research domains in the geosciences. Last but not least, I would like to express my gratitude to all the reviewers who kindly contributed to this Special Issue, as well as the Journal staff who have made this publication possible.

99 GENERAL AND MISCELLANEOUS↗

Research requirements for a real-time flight measurements and data analysis system for subsonic transport high-lift research

A multiphased research program to obtain detailed flow characteristics on a multielement high-lift flap system is being conducted on the Transport Systems Research Vehicle (B737-100 aircraft) at NASA Langley Research Center. Upcoming flight tests have required the development of a highly capable and flexible flight measurement and data analysis instrumentation system. This instrumentation system will be more comprehensive than any of the systems used on previous high-lift flight experiment at NASA Langley. The system will provide the researcher near-real-time information for decision making needed to modify a flight test in order to further examine unexpected flow conditions. This paper presents the research requirements and instrumentation design concept for an upcoming flight experiment for the subsonic transport high-lift research program. The flight experiment objectives, the measurement requirements, the data acquisition system, and the onboard data analysis and display capabilities are described.

Whitehead, Julia H.↗

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗

Fault-Related CO 2 Degassing, Geothermics, and Fluid Flow in Southern California Basins--Physiochemical Evidence and Modeling

Our studies have had an important impact on societal issues. Experimental and field observations show that CO 2 degassing, such as might occur from stored CO 2 reservoir gas, can result in significant stable isotopic disequilibrium. In the offshore South Ellwood field of the Santa Barbara channel, we show how oil production has reduced natural seep rates in the area, thereby reducing greenhouse gases. Permeability is calculated to be ~20-30 millidarcys for km-scale fault-focused fluid flow, using changes in natural gas seepage rates from well production, and poroelastic changes in formation pore-water pressure. In the Los Angeles (LA) basin, our characterization of formation water chemistry, including stable isotopic studies, allows the distinction between deep and shallow formations waters. Our multiphase computational-based modeling of petroleum migration demonstrates the important role of major faults on geological-scale fluid migration in the LA basin, and show how petroleum was dammed up against the Newport-Inglewood fault zone in a “geologically fast” interval of time (less than 0.5 million years). Furthermore, these fluid studies also will allow evaluation of potential cross-formational mixing of formation fluids. Lastly, our new study of helium isotopes in the LA basin shows a significant leakage of mantle helium along the Newport Inglewood fault zone (NIFZ), at flow rates up to 2 cm/yr. Crustal-scale fault permeability (~60 microdarcys) and advective versus conductive heat transport rates have been estimated using the observed helium isotopic data. The NIFZ is an important deep-seated fault that may crosscut a proposed basin decollement fault in this heavily populated area, and appears to allow seepage of helium from the mantle sources about 30 km beneath Los Angeles. The helium study has been widely cited in recent weeks by the news media, both in radio and on numerous web sites.

58 GEOSCIENCES↗

Computational study of the shock driven instability of a multiphase particle-gas system

This paper considers the interaction of a shock wave with a multiphase particle-gas system which creates an instability somewhat similar to the Richtmyer-Meshkov instability but with a larger parameter space. Because this parameter space is large, we only present an introductory survey of the effects of many of these parameters. We highlight the effects of particle-gas coupling, incident shock strength, particle size, effective system density differences, and multiple particle relaxation time effects. We focus on dilute flows with mass loading up to 40% and do not attempt to cover all parametric combinations. Instead, we vary one parameter at a time leaving additional parametric combinations for future work. In this work, the simulations are run with the Ares code, developed at Lawrence Livermore National Laboratory, which uses a multiphase particulate transport method to model two-way momentum and energy coupling. A brief validation of these models is presented and coupling effects are explored. It is shown that even for small particles, on the order of 1μm, multi-phase coupling effects are important and diminish the circulation deposition on the interface by up to 25%. These coupling effects are shown to create large temperature deviations from the dusty gas approximation, up to 20% greater, especially at higher shock strengths. It is also found that for a multiphase instability, the vortex sheet deposited at the interface separates into two sheets. In conclusion, depending on the particle and particle-gas Atwood numbers, the instability may be suppressed or enhanced by the interactions of these two vortex sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗