Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Geometry”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗
NASA NTRS1975

Stress analysis and stress intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses and the stress intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress intensity factors are also obtained for a through-thickness, center cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress intensity factors for commonly encountered surface crack geometries in finite solids.

Gyekenyesi, J. P.↗
DOE OSTI2020

Using Differential Forms to Find Symmetries in the Noh Problem for an Ideal Gas in a Spherical System

The Euler equations describe inviscid, compressible gas flow and are derived using the laws of conservation of mass, momentum, and energy. From these equations a common test problem, the Noh problem, can be solved for invariant solutions in planar, cylindrical, and spherical geometries. Using Lie group analysis methods, we will solve the Noh Problem for an ideal gas in a spherical system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1977

Stress analysis and stress-intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses, and the stress-intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress-intensity factors are also obtained for a through-thickness, center-cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress-intensity factors for commonly encountered surface crack geometries in finite solids

Gyekenyesi, J. P.↗
DOE OSTI2015

Highly temperature insensitive, low threshold-current density (λ = 8.7–8.8 μm) quantum cascade lasers

By stepwise tapering, both the barrier heights and quantum-well depths in the active regions of 8.7–8.8 μm-emitting quantum-cascade-laser (QCL) structures, virtually complete carrier-leakage suppression is achieved. Such step-taper active-region-type QCLs possess, for 3 mm-long devices with high-reflectivity-coated back facets, threshold-current characteristic temperature coefficients, T{sub 0}, as high as 283 K and slope-efficiency characteristic temperature coefficients, T{sub 1}, as high as 561 K, over the 20–60 °C heatsink-temperature range. These high T{sub 0} and T{sub 1} values reflect at least a factor of four reduction in carrier-leakage current compared to conventional 8–9 μm-emitting QCLs. Room temperature, pulsed, threshold-current densities are 1.58 kA/cm{sup 2}; values comparable to those for 35-period conventional QCLs of similar injector-region doping level. Superlinear behavior of the light-current curves is shown to be the result of the onset of resonant extraction from the lower laser level at a drive level of ∼1.3× threshold. Maximum room-temperature slope efficiencies are 1.23 W/A; that is, slope efficiency per period values of 35 mW/A, which are 37%–40% higher than for same-geometry conventional 8–9 μm-emitting QCLs. Since the waveguide-loss coefficients are very similar, we estimate that the internal differential efficiency is at least 30% higher than in conventional QCLs. Such high internal differential efficiency values reflect the combined effect of nearly complete carrier-leakage suppression and high differential efficiency of the laser transition (∼90%), due to resonant extraction from the lower laser level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

SU-F-19A-08: Optimal Time Release Schedule of In-Situ Drug Release During Permanent Prostate Brachytherapy

Purpose: Permanent prostate brachytherapy spacers can be used to deliver sustained doses of radiosentitizing drug directly to the target, in order to enhance the radiation effect. Implantable nanoplatforms for chemo-radiation therapy (INCeRTs) have a maximum drug capacity and can be engineered to control the drug release schedule. The optimal schedule for sensitization during continuous low dose rate irradiation is unknown. This work studies the optimal release schedule of drug for both traditional sensitizers, and those that work by suppressing DNA repair processes. Methods: Six brachytherapy treatment plans were used to model the anatomy, implant geometry and calculate the spatial distribution of radiation dose and drug concentrations for a range of drug diffusion parameters. Three state partial differential equations (cells healthy, damaged or dead) modeled the effect of continuous radiation (radiosensitivities α,β) and cellular repair (time tr) on a cell population. Radiosensitization was modeled as concentration dependent change in α,β or tr which with variable duration under the constraint of fixed total drug release. Average cell kill was used to measure effectiveness. Sensitization by means of both enhanced damage and reduced repair were studied. Results: Optimal release duration is dependent on the concentration of radiosensitizer compared to the saturation concentration (csat) above which additional sensitization does not occur. Long duration drug release when enhancing α or β maximizes cell death when drug concentrations are generally over csat. Short term release is optimal for concentrations below saturation. Sensitization by suppressing repair has a similar though less distinct trend that is more affected by the radiation dose distribution. Conclusion: Models of sustained local radiosensitization show potential to increase the effectiveness of radiation in permanent prostate brachytherapy. INCeRTs with high drug capacity produce the greatest benefit with drug release over weeks. If in-vivo drug concentrations are not able to approach saturation concentration, durations of days is optimal. DOD 1R21CA16977501; A. David Mazzone Awards Program 2012PD164.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-I-51: CT/MR Image Deformation: The Clinical Assessment QA in Target Delineation

Purpose: To study the deformation effects in CT/MR image registration of head and neck (HN) cancers. We present a clinical indication in guiding and simplifying registration procedures of this process while CT images possessed artifacts. Methods: CT/MR image fusion provides better soft tissue contrast in intracranial GTV definition with artifacts. However, whether the fusion process should include the deformation process is questionable and not recommended. We performed CT/MR image registration of a HN patient with tonsil GTV and nodes delineation on Varian Velocity™ system. Both rigid transformation and deformable registration of the same CT/MR imaging data were processed separately. Physician’s selection of target delineation was implemented to identify the variations. Transformation matrix was shown with visual identification, as well as the deformation QA numbers and figures were assessed. Results: The deformable CT/MR images were traced with the calculated matrix, both translation and rotational parameters were summarized. In deformable quality QA, the calculated Jacobian matrix was analyzed, which the min/mean/max of 0.73/0/99/1.37, respectively. Jacobian matrix of right neck node was 0.84/1.13/1.41, which present dis-similarity of the nodal area. If Jacobian = 1, the deformation is at the optimum situation. In this case, the deformation results have shown better target delineation for CT/MR deformation than rigid transformation. Though the root-mean-square vector difference is 1.48 mm, with similar rotational components, the cord and vertebrae position were aligned much better in the deformable MR images than the rigid transformation. Conclusion: CT/MR with/without image deformation presents similar image registration matrix; there were significant differentiate the anatomical structures in the region of interest by deformable process. Though vendor suggested only rigid transformation between CT/MR assuming the geometry remain similar, our findings indicated with patient positional variations, deformation registration is needed to generate proper GTV coverage, which will be irradiated more accurately in the following boost phase.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗
NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗
NASA NTRS1980

Gigahertz analog modulation and differential delay of GaAlAs lasers - Temperature and current behavior

Gigahertz analog modulation characteristics of broad-area commercially available GaAlAs lasers have been investigated as a function of temperature and current in the vicinity of the upper frequency limit, where the resonance phenomena occur. The optimum temperature for small-signal amplitude modulation was found to be around -15 C for our particular broad-stripe geometry double-heterostructure laser. The Q was found to increase by a factor of 2 and the bandwidth by about 2%; the external quantum efficiency was maximized in this range. The optimum dc current bias was about 2% above the threshold current. Differential delays have also been measured down to a few picosecond accuracy by a unique phase-angle measurement method using a vector voltmeter. Some of the temperature effects observed may be related to mode changes and multimode and superradiance behavior.

Eng, S. T.↗
DOE OSTI2018

Interference in electron–molecule elastic scattering

General formulas describing the multiple scattering of electron by polyatomic molecules have been derived within the framework of the model of non-overlapping atomic potentials. These formulas are applied to different carbon molecules, both for fixed-in-space and randomly oriented molecules. Here, the molecular continuum wave function is represented as a plane wave plus a linear combination of the Green’s functions for free motion and the derivatives of these functions. Far from the target the electron spherical waves interfere as in the case of the Young’s double slits experiment. This interference manifests itself as diffraction oscillations in the differential and total cross sections for elastic electron scattering. The amplitude of electron scattering by a molecule is defined by the phase shifts for each of the atoms forming the target and its geometry. The numerical calculation of the scattering amplitude in closed form (rather than in the form of S-matrix expansion) is reduced to solving a system of algebraic equations. The number of atoms in a molecule and the atomic phase shifts for the orbital angular momenta included in the calculation define the number of equations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

A multiscale restriction-smoothed basis method for high contrast porous media represented on unstructured grids

A wide variety of multiscale methods have been proposed in the literature to reduce runtime and provide better scaling for the solution of Poisson-type equations modeling flow in porous media. We present a new multiscale restricted-smoothed basis (MsRSB) method that is designed to be applicable to both rectilinear grids and unstructured grids. Like many other multiscale methods, MsRSB relies on a coarse partition of the underlying fine grid and a set of local prolongation operators (multiscale basis functions) that map unknowns associated with the fine grid cells to unknowns associated with blocks in the coarse partition. These mappings are constructed by restricted smoothing: Starting from a constant, a localized iterative scheme is applied directly to the fine-scale discretization to compute prolongation operators that are consistent with the local properties of the differential operators. The resulting method has three main advantages: First of all, both the coarse and the fine grid can have general polyhedral geometry and unstructured topology. This means that partitions and good prolongation operators can easily be constructed for complex models involving high media contrasts and unstructured cell connections introduced by faults, pinch-outs, erosion, local grid refinement, etc. In particular, the coarse partition can be adapted to geological or flow-field properties represented on cells or faces to improve accuracy. Secondly, the method is accurate and robust when compared to existing multiscale methods and does not need expensive recomputation of local basis functions to account for transient behavior: Dynamic mobility changes are incorporated by continuing to iterate a few extra steps on existing basis functions. This way, the cost of updating the prolongation operators becomes proportional to the amount of change in fluid mobility and one reduces the need for expensive, tolerance-based updates. Finally, since the MsRSB method is formulated on top of a cell-centered, conservative, finite-volume method, it is applicable to any flow model in which one can isolate a pressure equation. Herein, we only discuss single and two-phase incompressible models. Compressible flow, e.g., as modeled by the black-oil equations, is discussed in a separate paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A hybrid full-wave Markov chain approach to calculating radio-frequency wave scattering from scrape-off layer filaments

The interaction of radio-frequency (RF) waves with edge turbulence modifies the incident wave spectrum, and can significantly affect RF heating and current drive in tokamaks. Previous lower hybrid (LH) scattering models have either used the weak-turbulence approximation, or treated more realistic, filamentary turbulence in the ray tracing limit. In this study, a new model is introduced which retains full-wave effects of RF scattering in filamentary turbulence. First, a Mie-scattering technique models the interaction of an incident wave with a single Gaussian filament. Next, an effective differential scattering width is derived for a statistical ensemble of filaments. Lastly, a Markov chain solves for the transmitted wave spectrum in slab geometry. This model is applied to LH launching for current drive. The resulting wave spectrum is asymmetrically broadened in angular wavenumber space. This asymmetry is not accounted for in previous LH scattering models. The modified wave spectrum is coupled to a ray tracing/Fokker–Planck solver (GENRAY/CQL3D) to study its impact on current drive. The resulting current profile is greatly altered, and there is significant increase in the on-axis current and decrease in the off-axis peaks. This is attributed to a portion of the modified wave spectrum that is strongly dampened on-axis during the first pass.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Fault geometry and slip distribution of the 2008 M w 7.9 Wenchuan, China earthquake, inferred from GPS and InSAR measurements

We revisit the problem of coseismic rupture of the 2008 M w 7.9 Wenchuan earthquake. Precise determination of the fault structure and slip distribution provides critical information about the mechanical behaviour of the fault system and earthquake rupture. We use all the geodetic data available, craft a more realistic Earth structure and fault model compared to previous studies, and employ a nonlinear inversion scheme to optimally solve for the fault geometry and slip distribution. Compared to a homogeneous elastic half-space model and laterally uniform layered models, adopting separate layered elastic structure models on both sides of the Beichuan fault significantly improved data fitting. Our results reveal that: (1) The Beichuan fault is listric in shape, with near surface fault dip angles increasing from ~36° at the southwest end to ~83° at the northeast end of the rupture. (2) The fault rupture style changes from predominantly thrust at the southwest end to dextral at the northeast end of the fault rupture. (3) Fault slip peaks near the surface for most parts of the fault, with~8.4mthrust and ~5 m dextral slip near Hongkou and ~6 m thrust and ~8.4 m dextral slip near Beichuan, respectively. (4) The peak slips are located around fault geometric complexities, suggesting that earthquake style and rupture propagation were determined by fault zone geometric barriers. Such barriers exist primarily along restraining left stepping discontinuities of the dextralcompressional fault system. (5) The seismic moment released on the fault above 20 km depth is 8.2 × 10 21 N m, corresponding to an M w 7.9 event. The seismic moments released on the local slip concentrations are equivalent to events of M w 7.5 at Yingxiu-Hongkou, M w 7.3 at Beichuan-Pingtong, M w 7.2 near Qingping, M w 7.1 near Qingchuan, and M w 6.7 near Nanba, respectively. (6) The fault geometry and kinematics are consistent with a model in which crustal deformation at the eastern margin of the Tibetan plateau is decoupled by differential motion across a decollement in the mid crust, above which deformation is dominated by brittle reverse faulting and below which deformation occurs by viscous horizontal shortening and vertical thickening.

58 GEOSCIENCES↗
DOE OSTI2016

Atomic layer deposition of Al{sub 2}O{sub 3} for single electron transistors utilizing Pt oxidation and reduction

This work describes the fabrication of single electron transistors using electron beam lithography and atomic layer deposition to form nanoscale tunnel transparent junctions of alumina (Al{sub 2}O{sub 3}) on platinum nanowires using either water or ozone as the oxygen precursor and trimethylaluminum as the aluminum precursor. Using room temperature, low frequency conductance measurements between the source and drain, it was found that devices fabricated using water had higher conductance than devices fabricated with ozone. Subsequent annealing caused both water- and ozone-based devices to increase in conductance by more than 2 orders of magnitude. Furthermore, comparison of devices at low temperatures (∼4 K) showed that annealed devices displayed much closer to the ideal behavior (i.e., constant differential conductance) outside of the Coulomb blockade region and that untreated devices showed nonlinear behavior outside of the Coulomb blockade region (i.e., an increase in differential conductance with source-drain voltage bias). Transmission electron microscopy cross-sectional images showed that annealing did not significantly change device geometry, but energy dispersive x-ray spectroscopy showed an unusually large amount of oxygen in the bottom platinum layer. This suggests that the atomic layer deposition process results in the formation of a thin platinum surface oxide, which either decomposes or is reduced during the anneal step, resulting in a tunnel barrier without the in-series native oxide contribution. Furthermore, the difference between ozone- and water-based devices suggests that ozone promotes atomic layer deposition nucleation by oxidizing the surface but that water relies on physisorption of the precursors. To test this theory, devices were exposed to forming gas at room temperature, which also reduces platinum oxide, and a decrease in resistance was observed, as expected.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Preparation for a Measurement of Charge Asymmetry in the Bethe-Heitler Process

We have prepared a measurement of the energy asymmetry in wide- and medium-angle electron/positron pair production off protons and heavy targets. This asymmetry is caused by the interference between the first- and second-order Born diagrams and the Compton scattering diagram. It directly probes aspects of QED, as well as providing a direct measurement of the real part of the Compton amplitude. It will be conducted at the HI??S facility at Duke University, using a 60 MeV photon beam. This dissertation serves as documentation of the preparation stage of the Bethe-Heitler experiment. The major was the recommissioning of the vertical drift chambers previously used in the Q-weak experiment at the Jefferson Lab. Cosmic test runs were conducted, drift time data were collected and efficiency plateaus were measured. We made modifications to the JLAB Hall A analyzer to suit the geometry and drift characteristics of these wire chambers. The analyzer was used for the reconstruction of the trajectories of cosmic ray test runs with the results confirmed by direct measurement of trigger geometry. Spatial and angular resolution is estimated to ~300?? and 0.17° respectively. Geant 4 simulations with generated Bethe-Heitler pairs satisfying theoretical differential cross sections. It was used to check detector acceptance, optimize apparatus layout, and estimate measurable energy asymmetry. The measurable asymmetries from electron/positron pairs with polar angles around between approximately 5° and 8°, azimuthal angles differing by 180°, and energy differing by approximately 9 MeV to 15 MeV are predicted to be above 10%. The kinematics of primary vertices are reconstructed using the data from wire chambers in the simulation. The energy resolution is determined to be better than 1MeV.

Chen, Haoyu↗
DOE OSTI2025

Mechanical Properties and Interfacial Bonding of Overmolded Polymer Composite Lattice Structures

This study investigates the mechanical properties and interfacial bonding of Polyamide 12 filled with Glass Fiber (PA12/GF) lattice structures overmolded with Nylon 66 (PA66) and Thermoplastic Polyurethane (TPU). The PA12/GF lattices, designed in Gyroid, Isotruss, and Octahedral geometries, were produced using the Selective Laser Sintering (SLS) technique. The thermal characteristics of both the lattice and overmolding materials were analyzed using Differential Scanning Calorimetry (DSC) and Thermogravimetric Analysis (TGA). The mechanical performance and interfacial bonding were evaluated through mechanical testing and microstructural analysis. The results indicate that overmolding improved flexural strength and impact resistance of the PA12/GF lattices. For instance, the plain PA12/GF Isotruss lattice exhibited a flexural strength of 15.5 MPa. By overmolding it with TPU the flexural strength increased by up to 137%, while PA66-overmolding resulted in an increase of 371%, achieving a flexural strength of 73 MPa. In terms of impact resistance, TPU-overmolding improved performance significantly, with an increase of 1800% (Gyroid) compared to the plain lattice. Microstructural analysis revealed good adhesion at the interface, especially when both the lattice and overmolding materials were thoroughly dried prior to the overmolding process to minimize interfacial porosity. This study highlights the potential of using overmolding to enhance lattice structures performance, offering lightweight and cost-efficient solutions for automotive applications.

Talabi, Isaac [ORNL] (ORCID:0000000340215594)↗
NASA NTRS1995

Zero side force volute development

Collector scrolls on high performance centrifugal pumps are currently designed with methods which are based on very approximate flowfield models. Such design practices result in some volute configurations causing excessive side loads even at design flowrates. The purpose of this study was to develop and verify computational design tools which may be used to optimize volute configurations with respect to avoiding excessive loads on the bearings. The new design methodology consisted of a volute grid generation module and a computational fluid dynamics (CFD) module to describe the volute geometry and predict the radial forces for a given flow condition, respectively. Initially, the CFD module was used to predict the impeller and the volute flowfields simultaneously; however, the required computation time was found to be excessive for parametric design studies. A second computational procedure was developed which utilized an analytical impeller flowfield model and an ordinary differential equation to describe the impeller/volute coupling obtained from the literature, Adkins & Brennen (1988). The second procedure resulted in 20 to 30 fold increase in computational speed for an analysis. The volute design analysis was validated by postulating a volute geometry, constructing a volute to this configuration, and measuring the steady radial forces over a range of flow coefficients. Excellent agreement between model predictions and observed pump operation prove the computational impeller/volute pump model to be a valuable design tool. Further applications are recommended to fully establish the benefits of this new methodology.

Anderson, P. G.↗