1 / 13

Lecture Objectives:

Lecture Objectives:. Review discretization methods for advection diffusion equation Accuracy Numerical Stability Unsteady-state CFD Explicit vs. Implicit method HW2 Turbulence http://network.bepress.com/physical-sciences-and-mathematics/physics/fluid-dynamics/

finian
Download Presentation

Lecture Objectives:

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Lecture Objectives: • Review discretization methods for advection diffusion equation • Accuracy • Numerical Stability • Unsteady-state CFD • Explicit vs. Implicit method • HW2 • Turbulence • http://network.bepress.com/physical-sciences-and-mathematics/physics/fluid-dynamics/ • http://www.transportation.anl.gov/engines/multi_dim_model_les.html

  2. Steady–state 1D example I) X direction • If Vx > 0, • If Vx < 0, Convection term - Upwind-scheme: P E W dxe dxw and a) and Dx e w Diffusion term: b) When mesh is uniform: DX = dxe = dxw Assumption: Source is constant over the control volume Source term: c)

  3. Advection diffusion equation 1-D, steady-state Dx Dx N N+1 N-1 Different notation: Dx General equation

  4. 1D example multiple (N) volumes N unknowns i N 1 2 N-1 3 Equation for volume 1 N equations Equation for volume 2 …………………………… Equation matrix: For 1D problem 3-diagonal matrix

  5. 3D problem Equation in the general format: H N W E P S L Wright this equation for each discretization volume of your discretization domain A F 60,000 elements 60,000 cells (nodes) N=60,000 = x 60,000 elements 7-diagonal matrix This is the system for only one variable ( ) When we need to solve p, u, v, w, T, k, e, C system of equation is very larger

  6. Boundary conditionsfor CFD application - indoor airflow Real geometry Model geometry Where are the boundary Conditions?

  7. CFD ACCURACY Depends on airflow in the vicinity of Boundary conditions 1) At air supply device 2) In the vicinity of occupant 3) At room surfaces • Detailed modeling • limited by • computer power

  8. Surface boundaries thickness 0.01-20 mm for forced convection Wall surface W use wall functions to model the flow in the vicinity of surface Using relatively large mesh (cell) size.

  9. momentum sources Airflow at air supply devices Complex geometry - Δ~10-4m We can spend all our computing power for one small detail

  10. Diffuser jet properties High Aspiration diffuser D D L L How small cells do you need? We need simplified models for diffusers

  11. Simulation of airflow in In the vicinity of occupants How detailed should we make the geometry?

  12. General Transport EquationUnsteady-state H N Equation in the algebraic format: W E P S L We have to solve the system matrix for each time step ! Transient term: Are these values for step  or + ? If: -  - explicit method - + - implicit method Unsteady-state 1-D

  13. General Transport Equationunsteady-state 1-D Fully explicit method: Implicit method: Value form previous time step (known value) • Make the difference between • - Calculation for different time step • - Calculation in iteration step

More Related