Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution multiplicity”

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 1,405 records · Page 78

Pilot-Scale Testing of an Integrated Circuit for the Extraction of Rare Earth Minerals and Elements from Coal and Coal Byproducts Using Advanced Separation Technologies

The primary objective of this project was to develop and demonstrate an integrated pilot-scale circuitry for recovering high-value rare earth elements (REEs) from coal and coal byproducts. The target performance was to produce a mixed REE product with content of at least two percent by weight on a dry mass basis in a cost-effective and environmentally benign manner. During the first nine months of the project period (Phase 2 Budget Period 2), pilot plant construction was completed including all field site startup activities such as permitting, engineering design, procurement/bidding, unit fabrication, site construction, equipment installation, module assembly, safety training, and circuit shakedown. During the remaining 21 months of project period (Phase 2 Budget Period 3), detailed field-testing activities were performed including feedstock sample collection and preparation, exploratory testing, circuit modification, detailed parametric study, and performance optimization. A detailed techno-economic analysis was performed based on the pilot plant testing findings which provided various scenarios for REE production. The project successfully accomplished the proposed target performance by producing mixed rare earth oxide (REO) with greater than 90% purity by weight in a continuous pilot scale operation from two distinctly different coarse refuse materials (i.e., West Kentucky No. 13 and Fire Clay coal seams), and at least three secondary sources (i.e., heap leach process and naturally formed acid mine drainage system). Project partners included the University of Kentucky, Virginia Tech, West Virginia University, Alliance Coal, Blackhawk Mining, Mineral Refining Company, and Mineral Separation Technologies. The pilot scale test facility was constructed at a former mining complex owned by Alliance Natural Resource Partners (Alliance Coal). The site was rehabilitated to accommodate the equipment installation, construction and fabrication, electrical power requirement, water line management and containment. The process units constructed and installed included X-ray sorting unit, crushing and grinding unit, physical separation unit, acid leaching unit, solvent extraction unit, and wastewater management unit. A rare earth mineral concentration unit was constructed as a standalone unit for flexible operation. A detailed environmental assessment and control plan was carried out to identify and quantify any potential impacts of the pilot-scale processing circuitry on the human and eco-system health and well-being. Corresponding mitigation strategies and control measures were provided. A conceptual flowsheet was developed to effectively remove thorium and uranium from high purity rare earth oxide mix or any potential radionuclide enriched stream. The two distinct feedstock materials were secured from the Blackhawk Mining Complex in eastern Kentucky where the Fire Clay (Hazard No. 4) seam is processed. The West Kentucky No. 13 (Baker) coarse refuse material was collected from an active process stream at an Alliance coal preparation plant located in western Kentucky. Characterization analysis indicated that both of feed materials generated from the two sources contained >300 ppm of TREEs on a dry whole mass basis which met the requirements for a qualified feed stock. The two feedstocks were further upgraded using a dual x-ray sorter to prepare the feed material for hydrometallurgical circuit. Thermal treatment on feed material prior to leaching was found to: 1) improve the leaching recovery of REEs, 2) increase the leaching kinetics, and 3) allow the leaching reaction to occur at lower acidity. Roasting at 600°C was selected as the pre-treatment condition for both West Kentucky No. 13 and Fire Clay coarse refuse material. Over 40% of leaching recovery was achieved by roasting West Kentucky No. 13 material having a top particle size of 3 mm in the pilot scale operation using 1.2M sulfuric acid leaching at 75OC. Initial pilot scale testing involved continuous operation of the pilot plant for 94 hours. The leaching unit was operated at solid-to-liquid ratio of 1 to 10 (w/v) using 0.5M sulfuric acid solution at a temperature of 75°C. The continuous solvent extraction circuit utilized rougher and cleaner units with DEHPA and TBP as the extractants. An innovative stripping circuit was developed to accumulate the REE concentration in the stripping solution to a level above 600 ppm. A bleed stream from the recycled strip solution was treated using oxalic acid precipitation which produced a high grade rare earth oxalate. The oxalate product was roasted to remove the oxalate which produced a rare earth oxide product having a purity greater 90%. Due to high concentrations of contaminant ions in the pregnant leach solution (PLS), a modified flowsheet was developed that involved pre-concentration of the REEs using multiple stages of precipitation and redissolution. The advantage of this process was improved removal of contamination before the downstream purification process and a significant cost reduction relative to the circuit that utilized the solvent extraction process. The modified circuitry included processes involving leaching, multistage precipitation, redissolution, and oxalate precipitation followed by roasting of the oxalate product. The circuit produced a mixed REO that was 92.96% pure from the initial test. A detailed parametric test plan was carried out which involved varying key parameters including solids feed rate, acid flowrate, acid concentration, multistage precipitation pH, redissolution pH, oxalate precipitation dosage and pH. The response variables included REE recovery, contaminant recovery, REO product grade and overall chemical consumption. Test results indicated that the acid-to-solid ratio is the key parameter to leaching efficiency as performance deteriorated with an increase in solids concentration. The optimal pH determined for REE precipitation and redissolution was 6.5 and 2.5, respectively. Additional tests were conducted to further improve the flowsheet. Recirculating a portion of the PLS to the feed of the leach tanks improved the leaching performance by lowering the pH of the leaching system and reducing the contamination recovery by shortening the residence time. Moreover, the removal of Al prior to REE precipitation significantly reduced the oxalic acid consumption in the oxalate precipitation circuit. The modified circuit produced over 90% grade REO by weight from both West Kentucky No. 13 and Fire Clay coarse refuse material in pilot scale continuous test programs. A case analysis model was developed to project the REE and major contaminants concentration in each PLS stream based on the leaching condition, pH cut point, and oxalic acid dosage. A correlation was established using empirical and semi-empirical models. Using the models, chemical consumption required for each stage was predicted based on the projected performance of the hydrometallurgy circuit. After identifying the optimum conditions, validation tests were carried out for the treatment of both West Kentucky No. 13 and Fire Clay coarse refuse materials in the pilot plant. The actual circuit performance and chemical consumptions were very close to the model predictions. Other than the two coarse refuse sources, several secondary feedstocks were also tested in the pilot plant facility. A “heap leach” system was constructed using the coarse refuse material generated from cleaning the West Kentucky No. 13 seam coal. Using the two stage SX rougher and cleaner circuit, a concentrate with a grade >90% REO was produced while recovering >97% of the REEs from the heap leach PLS. Naturally generated acid mine drainage (AMD) from West Kentucky No.13 mine was processed using the multistage precipitation circuit in the pilot plant in a test conducted for a period of 32 hours. The final grade of the mix RE oxide produced from the AMD was 90.84% with an overall circuit recovery of 64%. The primary source of REE was the selective precipitation steps involving iron and aluminum rejection. The hydrophobic-hydrophilic separation (HHS) process was proven to effectively recover coal from fine waste materials. For REM recovery, the HHS process was able to produce concentrates at grades of approximately 1.8% REE on an ash basis; however, recovery values were typically low, <10%, under the optimal conditions determined in the laboratory-scale testing. Staged testing of the pilot-scale HHS process for coal recovery and semi-continuous laboratory testing for REM testing showed that a total concentration ratio of more than 15x was observed for the REM recovery process. A circuit simulation package was developed for REE extraction and purification using a spreadsheet-based platform (Microsoft Excel). The REESim circuit simulation package is configured to track the mass and volume flows of components passing through a series of unit operations specified and configured by the user. The mass rates can then be utilized by the user to determine important performance indicators such as product mass yields, concentrate purity levels, element-by-element recoveries, and so forth. The techno-economic analysis showed that the roasting and leaching operations were the most expensive capital items, each contributing approximately 30% to the total capital cost. One notable contributor to the high production costs was the low REE recovery observed in the pilot scale trials. The product basket price was shown to have a strong influence on the economic viability of the scenarios, with the scandium price being the most significant influencer. Operating cost was shown to be extremely sensitive to REE recovery, REE feed grade, and leaching acid consumption. An analysis of ten different scenarios for a 500 t/h commercial operation revealed that three were economically favorable, producing internal rates of return varying from 27.7% to 33.1% and payback periods of 4 to 5 years. The project successfully developed and demonstrated a process to recover REEs from coal and coal byproducts in a pilot-plant operation which consistently produced over 90% grade REO mix from varies types of feedstocks. Commercialization analysis showed that the technology readiness level successfully achieved TRL 6 at the end of the project and demonstrated the need and the potential for scaling the process to further advance the technologies toward the goal of providing a domestic supply of REEs at a commercial scale.

01 COAL, LIGNITE, AND PEAT↗

DOC-WMF - A Generic Model for Assessment of Waste Disposal Facilities - 19278

Numerical models and computer codes play an important role in the safety assessment of nuclear waste disposal facilities, as well as in the regulatory review of proponents' safety assessments. An important part of a safety assessment is the radiation dose calculation. This paper presents a computer code called DOC-WMF for calculating radiation dose from nuclear waste disposal facilities. DOC-WMF is a generic code that could be used for dose calculations for deep geological or near surface disposal facilities of nuclear wastes. DOC-WMF is much faster than complicated 3-D numerical models in terms of computational speed since the model uses the commercial code GoldSim as the computation engine which is based on analytical solutions for contaminant transport. As a result, it enables a user to quickly assess the dose consequences of multiple options or scenarios for a nuclear waste disposal facility. DOC-WMF was adapted from, and an improvement to, the Scoping of Options and Analyzing Risk (SOAR) model developed by the United States Nuclear Regulatory Commission (US NRC). A key improvement over the SOAR model is that DOC-WMF allows the calculation of dose contributions from other pathways such as food ingestion and air inhalation in addition to the drinking water pathway which is the sole dose contribution considered in the SOAR model. DOC-WMF also allows for: the consideration of 55 radionuclides with option to add more if required; the choice of source release method; various pathways to the biosphere; and dose calculations including an extended array of exposure pathways. In addition, DOC-WMF can simulate more complex hydrogeological settings and waste forms than SOAR does. DOC-WMF was compared against known analytical solutions for verification. It was then applied to two test cases and the results were compared. The first was a hypothetical deep geological disposal repository. The second test case was a near-surface nuclear waste disposal facility. DOC-WMF independently modeled these two test cases and achieved similar results as the more complicated numerical models did. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.4 is a program that obtains a steady-state flow and temperature solution, in geometrical extent, for a nuclear reactor core, or a single fuel assembly, including the effect of manufacturing and modelling uncertainties by the means of hot channel factors. It is based on an evolutionary sequence of codes originally used for calculating plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over the past 30 years. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset of nucleate boiling ratio (ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability for select design and safety parameters is available to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiple fields in stochastic inflation

Stochastic effects in multi-field inflationary scenarios are investigated. A hierarchy of diffusion equations is derived, the solutions of which yield moments of the numbers of inflationary e-folds. Solving the resulting partial differential equations in multi-dimensional field space is more challenging than the single-field case. A few tractable examples are discussed, which show that the number of fields is, in general, a critical parameter. When more than two fields are present for instance, the probability to explore arbitrarily large-field regions of the potential, otherwise inaccessible to single-field dynamics, becomes non-zero. In some configurations, this gives rise to an infinite mean number of e-folds, regardless of the initial conditions. Another difference with respect to single-field scenarios is that multi-field stochastic effects can be large even at sub-Planckian energy. This opens interesting new possibilities for probing quantum effects in inflationary dynamics, since the moments of the numbers of e-folds can be used to calculate the distribution of primordial density perturbations in the stochastic-δN formalism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multiple origins of extra electron diffractions in fcc metals

Diffuse intensities in the electron diffraction patterns of concentrated face-centered cubic solid solutions have been widely attributed to chemical short-range order, although this connection has been recently questioned. This article explores the many nonordering origins of commonly reported features using a combination of experimental electron microscopy and multislice diffraction simulations, which suggest that diffuse intensities largely represent thermal and static displacement scattering. A number of observations may reflect additional contributions from planar defects, surface terminations incommensurate with bulk periodicity, or weaker dynamical effects.

36 MATERIALS SCIENCE↗

Project H - A Complete Spaceport Hydrogen Solution

This slide presentation reviews Project H, and its importance in the development of Kennedy Space Center (KSC) as a Spaceport capable of multiple launches. It is known that current KSC cryogenic technology results in only approximately 55 % of purchased hydrogen being used. The rest is lost at various points in the process: transfer from transporting vehicle to tank, storage tank boil off, and from the tank to the intended propulsion tanks. Project H's goals would be to have local hydrogen production and liquifaction capability, and to increase the efficiency of hydrogen operations to greater than 80 %. The project envisions two phases: Phase 1 will build a smaller scale demonstration system, and phase 2 will build a full scale spaceport system. This initial project has proposed ideas for local hydrogen production, gaseous distribution, integrated refrigeration and storage, and high efficiency transfer lines that merit further investigation.

Notardonato, William↗

Linearity-Preserving Limiters on Irregular Grids

This paper examines the behavior of flux and slope limiters on non-uniform grids in multiple dimensions. We note that on non-uniform grids the scalar formulation in standard use today sacrifices k-exactness, even for linear solutions, impacting both accuracy and convergence. We rewrite some well-known limiters in a n way to highlight their underlying symmetry, and use this to examine both traditional and novel limiter formulations. A consistent method of handling stretched meshes is developed, as is a new directional formulation in multiple dimensions for irregular grids. Results are presented demonstrating improved accuracy and convergence using a combination of model problems and complex three-dimensional examples.

Berger, Marsha↗

Future Low-Carbon Technology Options for a 100 Percent Decarbonized Power Sector

Recent studies have modeled 100 percent decarbonized or renewable power sectors and proposed technical solutions to address costs, reliability, and other challenges, specifically related to meeting the last 10 percent of energy demand and eliminating carbon dioxide (CO 2 ) and other greenhouse gas (GHG) emissions. Here, this article summarizes multiple studies and provides references to the detailed and nuanced scenarios of reliably and cost-effectively operating 100 percent clean and/or renewable power systems. It also outlines future technology solutions that might be needed to meet this simply stated but deeply complex goal.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Dynamic properties of different liquid states in systems with competing interactions studied with lysozyme solutions

Recent studies of colloidal systems with a short-range attraction and long-range repulsion (SALR) have been demonstrated to have a generalized phase diagram with multiple liquid states defined by their structures. In this paper, we identify the different liquid states of previous experimentally studied lysozyme samples within this proposed generalized state diagram and explore the dynamic properties of each liquid state. We show that most lysozyme samples studied here and previously at low and intermediate concentrations are dispersed fluids while a few high concentration samples are randomly percolated liquids. In the dispersed fluid region, the short-time diffusion coefficient measured by neutron spin echo agrees well with the long time diffusion coefficient estimated with the solution viscosity. This dynamic feature is maintained even for some samples in the random percolated region. However, the short-time and long-time diffusion coefficients of random percolated fluids deviate at larger concentration and attraction strength. At high enough concentrations, the mean square displacement can be as slow as those of many glassy colloidal systems at time scales near the characteristic diffusion time even though these lysozyme samples remain in liquid states at the long-time limit. Here, we thus identify the region in the generalized phase diagram where these equilibrium states with extremely slow local dynamics exist relative to bulk percolation and kinetic arrest (gel and glassy) transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Workflow to Rapidly Interrogate Multiscale Model Simulation Results Across Multiple Length Scales

Many tools can be used to visualize field and state variables for a single scale analysis so that the influence of relevant mechanisms can be evaluated. Finite element software is often utilized to simulate a unit cell of a material and visualize results at that scale. Material properties can be homogenized from individual constituents and local deformation, damage, and failure mechanisms can be evaluated within the unit cell due to globally applied boundary conditions. Such solutions can produce satisfactory results if a user is only interested in analyzing a single scale. But materials in general contain features across multiple length scales, and assumptions must be made when attempting to account for lower length scale phenomena within a higher length scale model. Multiscale modeling is an attractive means to model materials because detailed material responses can be tracked across multiple disparate length scales while reducing the amount of required assumptions. However, as the complexity of these models increases, a large amount of data can be produced, and data traceability can become increasingly more difficult. Field and state variables, which are naturally dependent on spatial position, may themselves be calculated from one or more lower length scale unit cell models each with their own appropriate field and state variables. The NASA Multiscale Analysis Tool (NASMAT) is one software that can be used to perform a multiscale analysis efficiently and output requested data at all length scales in the analysis. A companion open-source Python software, NASMAT PrePost, can be used to visualize NASMAT model results and rapidly interrogate multiscale data across multiple length scales. This presentation will demonstrate some of the key features of NASMAT PrePost on two multiscale problems by quickly displaying and demonstrating connectivity among multiscale results from large datasets.

Python↗

How structural differences influence cross-model consistency: An electric sector case study

Multiple models are often employed to describe a range of possible outcomes for one or more scenarios, yielding insights into causal relationships and their uncertainties. Electric sector capacity expansion scenarios are a common topic of such efforts due to the economic influence of the electric sector, but model results typically span a broad solution space despite efforts to harmonize input assumptions, making decision implications difficult to discern. This study investigates the relationship between input harmonization and cross-model scenario consistency under disparate electric sector scenarios. We compare cross-model consistency between two state-of-the-art electric sector capacity expansion models (GCAM-USA and ReEDS) for six electric sector scenarios comprising alternate assumptions about fossil fuel resource availability, technology innovation, and long-term economy-wide transitions under four harmonization configurations varying model representations of electricity demand, fuel prices, renewable resources, and capacity retirements. These comparisons reveal that cross-model consistency can vary across scenarios under a given harmonization configuration, suggesting that harmonization efforts must often be scenario-specific if comparable cross-model consistency is desired. Model structural differences can hinder consistency, and the impact of these differences can depend on the scenario. Ultimately, thorough harmonization can reveal insights into cross-model consistency, which can be used to tighten uncertainty bounds and improve the decision-making implications of multi-model activities.

Cohen, Stuart↗

Maximal charge injection of consecutive electron pulses with uniform temporal pulse separation

A charge sheet model is proposed for the study of the space-charge limited density of consecutive electron pulses injected to in a diode with uniform temporal pulse separation. Based on the model, an analytical formula is derived for expressing the dependency of the charge density limit on the gap spacing, gap voltage, and pulse separation. The theoretical results are verified by numerical solutions up to electron energy of a few MeV, including relativistic effects. The model can be applied to the design of multiple-pulse electron beams for time resolved electron microscopy and free electron lasers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Crystallization by particle attachment in synthetic, biogenic, and geologic environments

Field and laboratory observations show that crystals commonly form by the addition and attachment of particles that range from multi-ion complexes to fully formed nanoparticles. The particles involved in these nonclassical pathways to crystallization are diverse, in contrast to classical models that consider only the addition of monomeric chemical species. We review progress toward understanding crystal growth by particle-attachment processes and show that multiple pathways result from the interplay of free-energy landscapes and reaction dynamics. Much remains unknown about the fundamental aspects, particularly the relationships between solution structure, interfacial forces, and particle motion. Developing a predictive description that connects molecular details to ensemble behavior will require revisiting long-standing interpretations of crystal formation in synthetic systems, biominerals, and patterns of mineralization in natural environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Multilevel Solution Algorithm for Markov Chains

We discuss the recently introduced multilevel algorithm for the steady-state solution of Markov chains. The method is based on an aggregation principle which is well established in the literature and features a multiplicative coarse-level correction. Recursive application of the aggregation principle, which uses an operator-dependent coarsening, yields a multi-level method which has been shown experimentally to give results significantly faster than the typical methods currently in use. When cast as a multigrid-like method, the algorithm is seen to be a Galerkin-Full Approximation Scheme with a solution-dependent prolongation operator. Special properties of this prolongation lead to the cancellation of the computationally intensive terms of the coarse-level equations.

Horton, Graham↗

Determination of critical experiment correlations for experiments involving highly enriched uranium solutions

Quantifying the correlations among critical experiments has become a topic of interest over the last 10-15 years. The correlations arise from the shared use of fissile material, non-fissile materials, fixtures, and procedures in multiple experiments or multiple experiment series. A Monte Carlo sampling technique proposed as a viable technical method to evaluate these correlations has been implemented within the SCALE 6.2 Sampler sequence. This paper presents the investigation of the use of this new tool and technique for highly enriched uranium solution experiments. One of the first attempts to quantify the correlations among critical experiments examined a series of highly enriched solution experiments, including the HST-001 series. These correlations are included in the International Criticality Safety Benchmark Evaluation Project Handbook, and they form a basis for methodology comparisons. Results generated using the Sampler implementation are compared with these previously published results. Additional studies involving the HST-001 experiment series are also presented. A range of different assumptions of the independence of the geometry components, enrichments, and other material compositions are implemented to investigate the impact of these assumptions on the correlation coefficients. These effects are investigated with these experiments because the calculation time is shorter for the solution experiments than for the low-enriched lattice experiments discussed in a companion paper. The magnitude of the effects of different independence assumptions on the correlation coefficients is important to establish which parameters must be clearly identified as shared or unique for accurate correlation determinations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A perpendicular ion beam instability - Solutions to the linear dispersion relation

A 200-eV Xe(+) ion beam directed perpendicular to the terrestrial magnetic field in the F region ionosphere produced very narrow band electrostatic emissions just above multiples of the hydrogen cyclotron frequency. Although the plasma conditions associated with the ion beam were undoubtedly very complex, a simple ion beam in a background ionosphere is considered first. The dispersion relation for flute mode waves and an unmagnetized perpendicular ion beam is solved for a diffuse H(+) plasma and then for a combination of dense O(+) and diffuse H(+). These solutions account for most of the wave properties, including the observation of narrow spectral peaks separated by the hydrogen cyclotron frequency and the observation of no spectral peaks below 2000 Hz. We cannot dismiss field-aligned currents associated with the Xe(+) beam as an alternate source of free energy for the narrow band emissions. However, our intention here is to examine closely the Xe(+) beam as a source for directly exciting the plasma waves.

Kintner, P. M.↗

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗

A moving mesh finite difference method for equilibrium radiation diffusion equations

An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor–corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗