Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “control volume methods”

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 181 records · Page 10

The value of the 4-day headdown bedrest model for screening countermeasures

In order to evaluate the benefits of periodic exposure to the +G(z) vector as a countermeasure to the physiological responses to minus 6 degree head down bedrest (HDT), we considered a two-tiered approach: (a) to use 4 days HDT as a quick and inexpensive means of screening countermeasures, (b) to use a 60 day HDT to validate the most promising candidates. The approach and results of a 4 day study are described here. Methods: Nine males were admitted to our Human Research Facility for one ambulatory control day followed by 4 days HDT and were released on the next day after completion of a peak oxygen consumption test (VO(sub 2 peak)). A battery of tests was selected and standardized to evaluate the known early effects of HDT on plasma volume, early bone markers, orthostatic tolerance, physical performance, and fluid and electrolytes and their hormone regulation. Fluid sodium (Na) and potassium (K) intake and output in the urine were monitored throughout. Plasma volume was determined with a modified Evans Blue method and orthostatic tolerance with a 60 degree head-up tilt test for 30 minutes - both of which were determined on the ambulatory control day and on day 4 of HDT. Immediately after completion of the tilt test subjects were returned to the minus 6 degree position until the next morning when a VO(sub 2 peak) (horizontal ergometer) was done. This was compared to a similar control test determined on 2 separate occasions before subject admission. Results: Four hours after going HDT produced significant decreases (p less than 0.05) in the circulating concentration of fluid and electrolyte regulating hormones. Plasma volume, orthostatic tolerance and VO(sub 2 peak) changed significantly after 4 days HDT. There was also the expected natriuresis on day 1 of HDT but no significant diuresis. The consistency of the pre-bedrest VO(sub 2 peak) tilt tests and plasma volumes was remarkable. Conclusions: The 4 day HDT model seems highly promising for screening a variety of countermeasures alone and in combination before validating their benefits in extended bedrest or flight experiments.

Vernikos, J.↗

Computationally-efficient stochastic cluster dynamics method for modeling damage accumulation in irradiated materials

An improved version of a recently developed stochastic cluster dynamics (SCD) method (Marian and Bulatov, 2012) [6] is introduced as an alternative to rate theory (RT) methods for solving coupled ordinary differential equation (ODE) systems for irradiation damage simulations. SCD circumvents by design the curse of dimensionality of the variable space that renders traditional ODE-based RT approaches inefficient when handling complex defect population comprised of multiple (more than two) defect species. Several improvements introduced here enable efficient and accurate simulations of irradiated materials up to realistic (high) damage doses characteristic of next-generation nuclear systems. The first improvement is a procedure for efficiently updating the defect reaction-network and event selection in the context of a dynamically expanding reaction-network. Next is a novel implementation of the τ-leaping method that speeds up SCD simulations by advancing the state of the reaction network in large time increments when appropriate. Lastly, a volume rescaling procedure is introduced to control the computational complexity of the expanding reaction-network through occasional reductions of the defect population while maintaining accurate statistics. The enhanced SCD method is then applied to model defect cluster accumulation in iron thin films subjected to triple ion-beam (Fe{sup 3+}, He{sup +} and H{sup +}) irradiations, for which standard RT or spatially-resolved kinetic Monte Carlo simulations are prohibitively expensive.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The MR appearance of volume overload in the lower extremities

PURPOSE: Our goal was to describe the MR findings of volume overload (VO) in the lower extremities. METHOD: Fifteen individuals were studied, including eight healthy controls and seven patients with VO (four cardiac, three renal). MR evaluation included various SE techniques. Edema detection, localization, and symmetry were assessed subjectively. Relaxation time estimates were also made of the subcutaneous tissue, marrow, and three muscles. RESULTS: Subcutaneous tissue was markedly edematous in seven of seven patients and asymmetric in four of seven, whereas marrow was normal in all patients. Muscle edema was mild and asymmetric in six and two of seven patients, respectively. Perifascial fluid collections were identified in six of seven patients. CONCLUSION: Subcutaneous tissue edema is the dominant feature of VO in the lower extremities. Perifascial fluid is common but does not necessarily distribute symmetrically. Muscle edema is relatively mild. These findings should aid in identifying VO as the potential cause of swelling in patients with swollen legs.

Non-NASA Center↗

System for detecting operating errors in a variable valve timing engine using pressure sensors

A method and control module includes a pressure sensor data comparison module that compares measured pressure volume signal segments to ideal pressure volume segments. A valve actuation hardware remedy module performs a hardware remedy in response to comparing the measured pressure volume signal segments to the ideal pressure volume segments when a valve actuation hardware failure is detected.

47 OTHER INSTRUMENTATION↗

GSRP/David Marshall: Fully Automated Cartesian Grid CFD Application for MDO in High Speed Flows

With the renewed interest in Cartesian gridding methodologies for the ease and speed of gridding complex geometries in addition to the simplicity of the control volumes used in the computations, it has become important to investigate ways of extending the existing Cartesian grid solver functionalities. This includes developing methods of modeling the viscous effects in order to utilize Cartesian grids solvers for accurate drag predictions and addressing the issues related to the distributed memory parallelization of Cartesian solvers. This research presents advances in two areas of interest in Cartesian grid solvers, viscous effects modeling and MPI parallelization. The development of viscous effects modeling using solely Cartesian grids has been hampered by the widely varying control volume sizes associated with the mesh refinement and the cut cells associated with the solid surface. This problem is being addressed by using physically based modeling techniques to update the state vectors of the cut cells and removing them from the finite volume integration scheme. This work is performed on a new Cartesian grid solver, NASCART-GT, with modifications to its cut cell functionality. The development of MPI parallelization addresses issues associated with utilizing Cartesian solvers on distributed memory parallel environments. This work is performed on an existing Cartesian grid solver, CART3D, with modifications to its parallelization methodology.

Source record↗

Cuprous Oxide Scale up: Gram Production via Bulk Synthesis using Classic Solvents at Low Temperatures

Cuprous oxide is a p-type semiconducting material that has been highly researched for its interesting properties. Many small-scale syntheses have exhibited excellent control over size and morphology. As the demand for cuprous oxide grows, the synthesis method need to evolve to facilitate large-scale production. This paper supplies a facile bulk synthesis method for Cu₂O on average, 1-liter reaction volume can produce 1 gram of particles. In order to study the shape and size control mechanisms on such a scale, the reaction volume was diminished to 250 mL producing on average 0.3 grams of nanoparticles per batch. Well-shaped nanoparticles have been synthesized using an aqueous solution of CuCl₂, NaOH, SDS surfactant, and NH₂OH-HCl at mild temperatures. The time allotted between the addition of NaOH and NH₂OH-HCl was determined to be critical for Cu(OH)2 production, an important precursor to the final produce The effects of stirring rates on a large scale was also analyzed during reagent addition and post reagent addition. A morphological change from rhombic dodecahedra to spheres occurred as the stirring speed was increased. The effects of NH₂OH-HCl concentration were also studied to control the etching effects of the final product.

36 MATERIALS SCIENCE↗

A compilation and analysis of helicopter handling qualities data. Volume 2: Data analysis

A compilation and an analysis of helicopter handling qualities data are presented. Multiloop manual control methods are used to analyze the descriptive data, stability derivatives, and transfer functions for a six degrees of freedom, quasi static model. A compensatory loop structure is applied to coupled longitudinal, lateral and directional equations in such a way that key handling qualities features are examined directly.

Heffley, R. K.↗

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗

SU-E-T-199: How Number of Control Points Influences the Dynamic IMRT Plan Quality and Deliverability

Purpose: To study the influence of number of control points on plan quality and deliverability. Methods: Five previously treated patients of carcinoma of rectum were selected. Planning target volume (PTV) and organs at risk (OARs) i.e. bladder and bowel were contoured. Dynamic IMRT plans (6MV, 7-fields, 45Gy/25 fractions and prescribed at 95% isodose) were created in Eclipse (Varian medical system, Palo Alto, CA) treatment planning system (TPS) for Varian CL2300C/D linear-accelerator. Base plan was calculated with 166 control points, variable mode (Eclipse Default). For generating other plans, all parameters were kept constant, only number of control points (Fixed mode) was varied as follows: 100, 166 and 200. Then, plan quality was analyzed in terms of maximum and mean dose received by the PTV and OARs. For plan deliverability, TPS calculated fluence was verified with I’matriXX (IBA Dosimetry, Germany) array and compared with TPS dose-plane using gamma index criteria of 3% dose difference and 3mm distance to agreement (DTA). Total number of monitor units (MU) required to deliver a plan was also noted. Results: The maximum variation for the PTV maximum with respect to eclipse default control point (166) was 0.28% (0.14Gy). Similarly, PTV mean varied only up to 0.22 %( 0.11Gy). Bladder maximum and bladder mean varied up to 0.51% (0.24Gy) and 0.16% (0.06Gy). The variation for the bowel maximum and bowel mean was also only 0.39% (0.19Gy) and 0.33% (0.04Gy). Total MU was within 0.32 % (4MU). Average gamma pass rate using different control points for five patients are 98.75±0.33%, 99.37±0.09%, 99.29±0.12%, 98.14±0.13% and 99.25±0.14% respectively. Conclusion: Slight variation (<1%) in PTV and OARs maximum and mean doses was observed with varying number of control points. Monitor unit was also not varied much. Reducing number of control points did not showed any comprise in plan deliverability in terms of gamma index pass rate.

60 APPLIED LIFE SCIENCES↗

Active Control of Fan Noise-Feasibility Study. Volume 2: Canceling Noise Source-Design of an Acoustic Plate Radiator Using Piezoceramic Actuators

The feasibility of using acoustic plate radiators powered by piezoceramic thin sheets as canceling sources for active control of aircraft engine fan noise is demonstrated. Analytical and numerical models of actuated beams and plates are developed and validated. An optimization study is performed to identify the optimum combination of design parameters that maximizes the plate volume velocity for a given resonance frequency. Fifteen plates with various plate and actuator sizes, thicknesses, and bonding layers were fabricated and tested using results from the optimization study. A maximum equivalent piston displacement of 0.39 mm was achieved with the optimized plate samples tested with only one actuator powered, corresponding to a plate deflection at the center of over 1 millimeter. This is very close to the deflection required for a full size engine application and represents a 160-fold improvement over previous work. Experimental results further show that performance is limited by the critical stress of the piezoceramic actuator and bonding layer rather than by the maximum moment available from the actuator. Design enhancements are described in detail that will lead to a flight-worthy acoustic plate radiator by minimizing actuator tensile stresses and reducing nonlinear effects. Finally, several adaptive tuning methods designed to increase the bandwidth of acoustic plate radiators are analyzed including passive, active, and semi-active approaches. The back chamber pressurization and volume variation methods are investigated experimentally and shown to be simple and effective ways to obtain substantial control over the resonance frequency of a plate radiator. This study shows that piezoceramic-based plate radiators can be a viable acoustic source for active control of aircraft engine fan noise.

Pla, F. G.↗

The method of planes pressure tensor for a spherical subvolume

Various formulas for the local pressure tensor based on a spherical subvolume of radius, R, are considered. An extension of the Method of Planes (MOP) formula of Todd et al. [Phys. Rev. E 52, 1627 (1995)] for a spherical geometry is derived using the recently proposed Control Volume formulation [E. R. Smith, D. M. Heyes, D. Dini, and T. A. Zaki, Phys. Rev. E 85, 056705 (2012)]. The MOP formula for the purely radial component of the pressure tensor is shown to be mathematically identical to the Radial Irving-Kirkwood formula. Novel offdiagonal elements which are important for momentum conservation emerge naturally from this treatment. The local pressure tensor formulas for a plane are shown to be the large radius limits of those for spherical surfaces. The radial-dependence of the pressure tensor computed by Molecular Dynamics simulation is reported for virtual spheres in a model bulk liquid where the sphere is positioned randomly or whose center is also that of a molecule in the liquid. The probability distributions of angles relating to pairs of atoms which cross the surface of the sphere, and the center of the sphere, are presented as a function of R. The variance in the shear stress calculated from the spherical Volume Averaging method is shown to converge slowly to the limiting values with increasing radius, and to be a strong function of the number of molecules in the simulation cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Navier-Stokes analysis of three-dimensional flow and heat transfer inside turbine blade rows

A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.

Hah, C.↗

Effect of natural ageing on volume stability of MSW and wood waste incineration residues

Highlights: ► Natural weathering on BA from MSW and wood waste incineration was evaluated. ► Type of mineral phases, pH and volume stability were considered. ► Weathering reactions effect in improved stability of the materials. - Abstract: This paper presents the results of a study on the effect of natural weathering on volume stability of bottom ash (BA) from municipal solid waste (MSW) and wood waste incineration. BA samples were taken at different steps of treatment (fresh, 4 weeks and 12 weeks aged) and then characterised for their chemical and mineralogical composition and for volume stability by means of the mineralogical test method (M HMVA-StB), which is part of the German quality control system for using aggregates in road construction (TL Gestein-StB 04). Changes of mineralogical composition with the proceeding of the weathering treatment were also monitored by leaching tests. At the end of the 12 weeks of treatment, almost all the considered samples resulted to be usable without restrictions in road construction with reference to the test parameter volume stability.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-dimensional thermal modeling of power monolithic microwave integrated circuits (MMIC's)

Numerical simulations of the two-dimensional temperature distributions for a typical GaAs MMIC circuit are conducted, aiming at understanding the heat conduction process of the circuit chip and providing temperature information for device reliability analysis. The method used is to solve the two-dimensional heat conduction equation with a control-volume-based finite difference scheme. In particular, the effects of the power dissipation and the ambient temperature are examined, and the criterion for the worst operating environment is discussed in terms of the allowed highest device junction temperature.

Fan, Mark S.↗

Validation of deformable image registration algorithms on CT images of ex vivo porcine bladders with fiducial markers

Purpose: The spatial accuracy of deformable image registration (DIR) is important in the implementation of image guided adaptive radiotherapy techniques for cancer in the pelvic region. Validation of algorithms is best performed on phantoms with fiducial markers undergoing controlled large deformations. Excised porcine bladders, exhibiting similar filling and voiding behavior as human bladders, provide such an environment. The aim of this study was to determine the spatial accuracy of different DIR algorithms on CT images ofex vivo porcine bladders with radiopaque fiducial markers applied to the outer surface, for a range of bladder volumes, using various accuracy metrics. Methods: Five excised porcine bladders with a grid of 30–40 radiopaque fiducial markers attached to the outer wall were suspended inside a water-filled phantom. The bladder was filled with a controlled amount of water with added contrast medium for a range of filling volumes (100–400 ml in steps of 50 ml) using a luer lock syringe, and CT scans were acquired at each filling volume. DIR was performed for each data set, with the 100 ml bladder as the reference image. Six intensity-based algorithms (optical flow or demons-based) implemented in theMATLAB platform DIRART, a b-spline algorithm implemented in the commercial software package VelocityAI, and a structure-based algorithm (Symmetric Thin Plate Spline Robust Point Matching) were validated, using adequate parameter settings according to values previously published. The resulting deformation vector field from each registration was applied to the contoured bladder structures and to the marker coordinates for spatial error calculation. The quality of the algorithms was assessed by comparing the different error metrics across the different algorithms, and by comparing the effect of deformation magnitude (bladder volume difference) per algorithm, using the Independent Samples Kruskal-Wallis test. Results: The authors found good structure accuracy without dependency on bladder volume difference for all but one algorithm, and with the best result for the structure-based algorithm. Spatial accuracy as assessed from marker errors was disappointing for all algorithms, especially for large volume differences, implying that the deformations described by the registration did not represent anatomically correct deformations. The structure-based algorithm performed the best in terms of marker error for the large volume difference (100–400 ml). In general, for the small volume difference (100–150 ml) the algorithms performed relatively similarly. The structure-based algorithm exhibited the best balance in performance between small and large volume differences, and among the intensity-based algorithms, the algorithm implemented in VelocityAI exhibited the best balance. Conclusions: Validation of multiple DIR algorithms on a novel physiological bladder phantom revealed that the structure accuracy was good for most algorithms, but that the spatial accuracy as assessed from markers was low for all algorithms, especially for large deformations. Hence, many of the available algorithms exhibit sufficient accuracy for contour propagation purposes, but possibly not for accurate dose accumulation.

60 APPLIED LIFE SCIENCES↗

Three-dimensional flow analysis inside turbomachinery stages with steady and unsteady Navier-Stokes method

This study presents a numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed in the present study. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow field of a subsonic turbine stage and the unsteady flow field inside a transonic, high-through-flow, axial compressor stage. The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid.

Copenhaver, W. W.↗

SU-C-BRE-02: BED Vs. Local Control: Radiobiological Effect of Tumor Volume in Monte Carlo (MC) Lung SBRT Planning

Purpose: SBRT with hypofractionated dose schemata has emerged a compelling treatment modality for medically inoperable early stage lung cancer patients. It requires more accurate dose calculation and treatment delivery technique. This report presents the relationship between tumor control probability(TCP) and size-adjusted biological effective dose(sBED) of tumor volume for MC lung SBRT patients. Methods: Fifteen patients who were treated with MC-based lung SBRT to 50Gy in 5 fractions to PTVV100%=95% were studied. ITVs were delineated on MIP images of 4DCT-scans. PTVs diameter(ITV+5mm margins) ranged from 2.7–4.9cm (mean 3.7cm). Plans were generated using non-coplanar conformal arcs/beams using iPlan XVMC algorithm (BrainLABiPlan ver.4.1.2) for Novalis-TX with HD-MLCs and 6MVSRS(1000MU/min) mode, following RTOG-0813 dosimetric guidelines. To understand the known uncertainties of conventional heterogeneities-corrected/uncorrected pencil beam (PBhete/ PB-homo) algorithms, dose distributions were re-calculated with PBhete/ PB-homo using same beam configurations, MLCs and monitor units. Biologically effective dose(BED10) was computed using LQ-model with α/β=10Gy for meanPTV and meanITV. BED10-c*L, gave size-adjusted BED(sBED), where c=10Gy/cm and L=PTV diameter in centimeter. The TCP model was adopted from Ohri et al.(IJROBP, 2012): TCP = exp[sBEDTCD50]/ k /(1.0 + exp[sBED-TCD50]/k), where k=31Gy corresponding to TCD50=0Gy; and more realistic MC-based TCP was computed for PTV(V99%). Results: Mean PTV PB-hete TCP value was 6% higher, but, mean PTV PB-homo TCP value was 4% lower compared to mean PTV MC TCP. Mean ITV PB-hete/PB-homo TCP values were comparable (within ±3.0%) to mean ITV MC TCP. The mean PTV(V99%)had BED10=90.9±3.7%(median=92.2%),sBED=54.1±8.2%(median=53.5%) corresponding to mean MC TCP value of 84.8±3.3%(median=84.9%) at 2- year local control. Conclusion: The TCP model which incorporates BED10 and tumor diameter indicates that radiobiological effect of target volume and dose calculation algorithm significantly affects TCP for lung SBRT patients. Dose calculation using MC-based algorithm is more realistic with tissue heterogeneities and is routinely performed in our clinic. Patients will be followed up to determine whether TCP prediction correlate clinical outcomes.

60 APPLIED LIFE SCIENCES↗

AERACEPT (AErosol Rapid Analysis Combined Entry Probe/sonde Technology) for the Nephele Venus Cloud Mission Concept

Nephele is a small atmospheric probe mission concept to determine whether the Venus clouds contain organic matter. Nephele complements larger missions that seek to conduct Venus atmospheric analyses, such as DAVINCI, Venera-D, and the Morning Star missions, by directly sampling and analyzing cloud and haze particles. AERACEPT is an aerosol sampling technology designed to minimize the mass, volume, and complexity required for in situ planetary atmosphere characterization. A single aeroshell functions as a combined entry vehicle, descent probe, and instrumented sonde, combining recent advances in thermal protection materials with heritage technology from planetary and airborne science. The velocity of descent drives aerosol through inlets embedded in the aeroshell’s thermal protection material; inertial separation isolates the particles from the gas stream; and the particle capture surface doubles as the analysis substrate. This method can work without heat shield separation, descent control, or active sample handling; these advantages trade against a smaller sample volume, higher risk of sample heating, and shorter time for analysis and data return. AERACEPT further benefits from modeling of the internal and external flow to understand the particle capture efficiency and bias for size and concentration distributions. AERACEPT is well suited for a Venus mission such as Nephele, where a passive descent trajectory is both subsonic over the altitudes of interest and provides sufficient velocity and pressure differential for the particle sizes of interest. A thermal and flow toolchain developed in-house (Traj→ DPLR→ PATO→ openFOAM) was used to model an 80-cm aeroshell with science operations 63–39 km. These models indicate that AERACEPT will provide ~24 µL accumulated sample of particles ≥ 0.2 µm, substantially above the 2 µL threshold requirement. The worst-case particle heating will be ~14 K above ambient, for the smallest particles at the lowest altitude; most particles will experience ≤ 1 K rise, at which almost all analytes of interest will be unaffected. Plasma jet testing of a small-scale aeroshell-material inlet was recently conducted to validate the thermal model parameters; wind tunnel testing to validate the flow model parameters is underway.

AERACEPT↗