Optimization of Nonlinear, Coupled Fluid-Thermal Systems

Содержание

Слайд 2

Presentation Outline Overview Project Goals Microgravity Research MGFLO Optimization Theory Previous

Presentation Outline

Overview
Project Goals
Microgravity Research
MGFLO
Optimization Theory
Previous Work

Code Details
Overview
Validation
Applications

Conclusions
Recommendations
Слайд 3

Project Goals To Design and Implement an optimization algorithm for a

Project Goals
To Design and Implement an optimization algorithm for a fluid-thermal

simulator
MGFLO
Boundary Condition Manipulation
Слайд 4

Microgravity Fluid Research Surface Tension Smallest Surface Area Possible Dominated on

Microgravity Fluid Research

Surface Tension
Smallest Surface Area Possible
Dominated on Earth by Gravity,

which Makes Surfaces Flat

Liquid Bridges
ALEX: A Liquid Electrohydrodynamics eXperiment
Surface Tension Dominates with Decreased Electric Field

In a microgravity environment, surface tension and
thermocapillary effects can be dominant.

Слайд 5

Strong Field Weak Field

Strong
Field

Weak
Field

Слайд 6

Microgravity Test Facilities Drop Towers Evacuated tubes used to expose experiments

Microgravity Test Facilities

Drop Towers
Evacuated tubes used to expose experiments to several

seconds of microgravity
Only short durations of microgravity are achieved
Слайд 7

Test Facilities NASA’s KC-135 “Vomit Comet” Parabolic flight pattern can produce

Test Facilities

NASA’s KC-135 “Vomit Comet”
Parabolic flight pattern can produce up to

30 seconds of microgravity
Several periods of microgravity in one flight
Слайд 8

Test Facilities Sounding Rockets Also flown in a parabolic flight path

Test Facilities

Sounding Rockets
Also flown in a parabolic flight path to produce

microgravity
Can provide 6-7 minutes of microgravity
Слайд 9

Microgravity Simulation Computational Fluid Dynamics (CFD) allows cost-effective microgravity simulation Advances

Microgravity Simulation

Computational Fluid Dynamics (CFD) allows cost-effective microgravity simulation
Advances in parallel

supercomputing allow large problems to be solved
Слайд 10

Incompressible Navier-Stokes Equations: Energy Equation: Governing Equations

Incompressible Navier-Stokes Equations:
Energy Equation:

Governing Equations

Слайд 11

MGFLO Developed Under NASA-Grand Challenge Support Parallel, Finite Element Formulation of

MGFLO

Developed Under NASA-Grand Challenge Support
Parallel, Finite Element Formulation of Navier-Stokes and

Energy Equations
Allows for Coupled and Uncoupled Solution
Systems Optimized Through Matlab Using Existing Algorithms
Слайд 12

Optimization Theory Attempt to find “best value” of a merit function

Optimization Theory

Attempt to find “best value” of a merit function within

defined constraints
Gradient versus non-gradient methods
Gradient methods can be complex and require several merit function evaluations
Non-gradient methods optimize based on a sample set of merit function values
Nelder-Mead Simplex Search Algorithm
Слайд 13

Nelder and Mead’s Method Efficient search method for minimizing a merit

Nelder and Mead’s Method

Efficient search method for minimizing a merit function

of up to six variables
Optimization points are nodes of a polygon
Optimal solution is determined by:
Reflection
Expansion
Contraction
Слайд 14

Simplex Steps Reflection Expansion Contraction

Simplex Steps

Reflection

Expansion

Contraction

Слайд 15

Previous Work Investigated Operation of the MGFLO Code Designed Simple Optimization

Previous Work

Investigated Operation of the MGFLO Code
Designed Simple Optimization Routine in

Matlab
Established Algorithms to Optimize Complex Fluid-Thermal Systems
Слайд 16

Code Overview Developed Matlab Routines to Analyze MGFLO Output. Matlab Can

Code Overview

Developed Matlab Routines to Analyze MGFLO Output.
Matlab Can Compute Quantities

of Interest:
Vorticity, Divergence
Gradient, Laplacian
0th, 1st, 2nd Order Derivatives Normal to Walls
Average Quantities in Large Datasets
Слайд 17

Code Functions Initializes the solution Calls MGFLO for each simplex step

Code Functions

Initializes the solution
Calls MGFLO for each simplex step
Checks that user-specified

constraints are satisfied
Calculates the user-specified merit function
Allows user to monitor solution progression
Слайд 18

Слайд 19

Debugging & Validation Attempt to find answer to a known problem

Debugging & Validation
Attempt to find answer to a known problem
Position heat

source on top surface to maximize heat flux out of the bottom
Run on the 16-node Beowulf cluster in the CFDLab
Слайд 20

Слайд 21

Слайд 22

Optimization Path

Optimization Path

Слайд 23

Limitations Merit function dependence for pathological problems Not successful at maximizing

Limitations
Merit function dependence for pathological problems
Not successful at maximizing vorticity in

previous case
Non-smooth merit functions (too many local maxima)
Слайд 24

Applications Solve more complicated problem whose answer is not known a-priori

Applications

Solve more complicated problem whose answer is not known a-priori
System exposed

to external environment via Newton’s law of cooling (mixed boundary condition)
Use particle tracing as a visualization technique
Слайд 25

Слайд 26

Case 1: Tdesired=310K

Case 1: Tdesired=310K

Слайд 27

Particle Tracing Algorithm Heun predictor-corrector method Second-order accurate in time Allows visualization/quantification of mixing

Particle Tracing Algorithm

Heun predictor-corrector method
Second-order accurate in time
Allows visualization/quantification of

mixing
Слайд 28

Слайд 29

Слайд 30

Convergence History

Convergence History

Слайд 31

Case 2: Tdesired=340K

Case 2: Tdesired=340K

Слайд 32

Слайд 33

Слайд 34

Convergence History

Convergence History

Слайд 35

Conclusions We became familiar with the CFDLab and the MGFLO code

Conclusions

We became familiar with the CFDLab and the MGFLO code
Successfully developed

a method to optimize nonlinear fluid-thermal systems
Implemented a particle tracing algorithm in Matlab to visualize fluid mixing
Слайд 36

Recommendations Use particle tracing algorithm to optimize system mixing (currently takes

Recommendations

Use particle tracing algorithm to optimize system mixing (currently takes a

long time!)
Implement feedback control for time-varying systems
Calculate merit function interior to MGFLO
Faster
More accurate
Support unstructured grids