Project 1

Unit Number: 301023

Unit Name: Advanced Thermal and Fluid Engineering

Instruction

• Students will develop CFD programs to solve simple problems in this project.

• There are two parts in this project, each student needs to complete the work assigned in each part and submit a written report.

• The report must be in the format of Miscrosft Word file and submitted through vUWS. Reports in hard copy will not be accepted.

Note:

Plagiarism

Plagiarism is absolutely unacceptable. If plagiarism is suspected to be involved in a report, it will be reported to the school for further investigation.

Late submission

- Submitting the report late due to your own internet connection problem will be treated as late submission without a

valid reason. To avoid delay, you need to submit the report well before the due date.

Part one: Simple initial value problem

Assignment (please read the problem description before doing the assignment)

(1) Find out the equation of motion of the cylinder. The derivation procedure of the equation must be included in the report.

(2) Develop a numerical method for predicting the vibration of the cylinder. The procedure of deriving the FDM formula must be included in the report.

(3) Develop a MATLAB program for solving the problem. The MATLAB program code must be included in the report.

(4) Keeping the damping coefficient of c = 100 N·m/s, conduct numerical simulations for KC=2 to 20 with an interval of 2 and find out the variation of the power with the KC number.

(5) Keeping the KC number to be constant of KC=10, conduct numerical simulations for c=0 to 2000 N·m/s with an interval of 50 and find the variation of the power with c.

(6) In the report, show the time histories of the vibration for all the calculated cases.

(7) Discuss how KC and c affects the power generation.

Problem description

An energy extraction device includes a cylinder elastically mounted on a spring and the electricity generator. The cylinder vibrates only in a direction parallel to the flow and drives an electricity generator. The electricity generator extracts the energy from the cylinder in the same way as a damper. It is modelled as a damper with a damping constant c when the vibration is studied. The flow velocity is a combined steady and oscillatory flow u(t) ?U0 ?Um sin(2?t /T)

where U0 is the steady flow velocity, Um is the amplitude of the oscillatory flow velocity, T is the period of the oscillatory flow and t is the time.

The KC number is defined as

UmT

KC ?

D

Flow velocity u=U0+Umsin(?t)

The electricity generator can modelled by a damper that extract energy from the motion of the cylinder with a constant damping constant of c. The structural damping is considered to be negligibly small. The energy extracted from the motion by the electricity generator (damper) can be calculated by

1 T 2 P ? T ?0 cV dt

where P is the power of the system, F is the hydrodynamic force on the cylinder, V is the vibration speed of the cylinder. In the numerical simulation, the power can be calculated by

1 N ? n?2

P? N ?n?1cV

where N is the total time step in one stable oscillatory flow period, Fn and Vn represent n-th step within one oscillatory flow period. The power must be calculated by stable numerical solution. The hydrodynamic force on the cylinder can be predicted by the Morison equation

dVr 1

Fwater ? CAmd ? ?CDApVr Vr

dt 2

? ?

Finertia Fdrag

where Vr=u-V is the velocity of the flow relative to the cylinder, CA and CD and inertia and drag coefficients, respectively, md is the displaced fluid mass and Ap is the projected area. The calculating parameters are listed in the table below.

Table 1 Parameter table

Note: In the table, the number j stands for the last digit of your student ID. For example, if your student ID is 2345673, 10+j = 13.

Amplitude of the oscillatory flow velocity, Um j

1? m/s

10

Steady flow velocity U0 0.1×Um m/s

Diameter of the cylinder, D 10+j cm

Length of the cylinder, L 1 m

Mass of the cylinder, m 50 kg

Density of the water, ? 1024 kg/m3

Stiffness of the spring (total), K 200 N/m

Damping coefficient (total), c 10 N?m/s

Added mass coefficient (assume to be constant) 1

Drag coefficient (assume to be constant) 1.8

Note: Stable vibration can be achieved only if the simulated time is long enough.

Part two: One-dimensional convection diffusion problem

Assignment (please read the problem description before doing the assignment)

• Develop a finite difference method (FDM) formula for solving the equation. The procedure of deriving the FDM formula must be included in the report.

• Develop a MATLAB program and use this program to do the simulation and answer the following questions. Divide the soil depth into 200 cells when solving the problem numerically. You must attach the MATLAB program in your report.

• If a=0.0002 m2/s, find out when the temperature at the depth of 1.5 m reaches 25°C.

• If the temperature at the depth of 1.5 m reaches 25°C is defined as T1.5, do the simulations for a in the range of 0.0002 m2/s to 0.001 m2/s with an interval of 0.0001 m2/s and discuss how a on T1.5 and why.

Problem description

Note: In the paragraph below, the number j stands for the last digit of your student ID. For example, if your student ID is 2345673, 10+j = 13.

The water on the ground flows into the soil through the pore space of the soil (the space between soil particles). The flow velocity is very small of 0.002×(1+j) m/s and it does not change along the vertical direction. Initially, the temperature is (10+j) °C in the whole soil volume (including soil and the water in the soil). The temperature of the water on top of the ground level is 30°C and remains constant. The warm water moves into the soil and increases the temperature in the soil.

A coordinate x is defined and its origin is on the ground level and pointing downwards. The governing equation is the convection diffusion equation

?T ?T ?2T

??u ?? 2 ?t ?x ?x

where

T is the temperature in the soil t is the time u is the water flow velocity a is the heat diffusion coefficient (the influence of the soil has been considered).

?T

Boundary condition is T=30°C at the top boundary and ?0 at the bottom boundary (at x= 2 m).

?x

Unit Number: 301023

Unit Name: Advanced Thermal and Fluid Engineering

Instruction

• Students will develop CFD programs to solve simple problems in this project.

• There are two parts in this project, each student needs to complete the work assigned in each part and submit a written report.

• The report must be in the format of Miscrosft Word file and submitted through vUWS. Reports in hard copy will not be accepted.

Note:

Plagiarism

Plagiarism is absolutely unacceptable. If plagiarism is suspected to be involved in a report, it will be reported to the school for further investigation.

Late submission

- Submitting the report late due to your own internet connection problem will be treated as late submission without a

valid reason. To avoid delay, you need to submit the report well before the due date.

Part one: Simple initial value problem

Assignment (please read the problem description before doing the assignment)

(1) Find out the equation of motion of the cylinder. The derivation procedure of the equation must be included in the report.

(2) Develop a numerical method for predicting the vibration of the cylinder. The procedure of deriving the FDM formula must be included in the report.

(3) Develop a MATLAB program for solving the problem. The MATLAB program code must be included in the report.

(4) Keeping the damping coefficient of c = 100 N·m/s, conduct numerical simulations for KC=2 to 20 with an interval of 2 and find out the variation of the power with the KC number.

(5) Keeping the KC number to be constant of KC=10, conduct numerical simulations for c=0 to 2000 N·m/s with an interval of 50 and find the variation of the power with c.

(6) In the report, show the time histories of the vibration for all the calculated cases.

(7) Discuss how KC and c affects the power generation.

Problem description

An energy extraction device includes a cylinder elastically mounted on a spring and the electricity generator. The cylinder vibrates only in a direction parallel to the flow and drives an electricity generator. The electricity generator extracts the energy from the cylinder in the same way as a damper. It is modelled as a damper with a damping constant c when the vibration is studied. The flow velocity is a combined steady and oscillatory flow u(t) ?U0 ?Um sin(2?t /T)

where U0 is the steady flow velocity, Um is the amplitude of the oscillatory flow velocity, T is the period of the oscillatory flow and t is the time.

The KC number is defined as

UmT

KC ?

D

Flow velocity u=U0+Umsin(?t)

The electricity generator can modelled by a damper that extract energy from the motion of the cylinder with a constant damping constant of c. The structural damping is considered to be negligibly small. The energy extracted from the motion by the electricity generator (damper) can be calculated by

1 T 2 P ? T ?0 cV dt

where P is the power of the system, F is the hydrodynamic force on the cylinder, V is the vibration speed of the cylinder. In the numerical simulation, the power can be calculated by

1 N ? n?2

P? N ?n?1cV

where N is the total time step in one stable oscillatory flow period, Fn and Vn represent n-th step within one oscillatory flow period. The power must be calculated by stable numerical solution. The hydrodynamic force on the cylinder can be predicted by the Morison equation

dVr 1

Fwater ? CAmd ? ?CDApVr Vr

dt 2

? ?

Finertia Fdrag

where Vr=u-V is the velocity of the flow relative to the cylinder, CA and CD and inertia and drag coefficients, respectively, md is the displaced fluid mass and Ap is the projected area. The calculating parameters are listed in the table below.

Table 1 Parameter table

Note: In the table, the number j stands for the last digit of your student ID. For example, if your student ID is 2345673, 10+j = 13.

Amplitude of the oscillatory flow velocity, Um j

1? m/s

10

Steady flow velocity U0 0.1×Um m/s

Diameter of the cylinder, D 10+j cm

Length of the cylinder, L 1 m

Mass of the cylinder, m 50 kg

Density of the water, ? 1024 kg/m3

Stiffness of the spring (total), K 200 N/m

Damping coefficient (total), c 10 N?m/s

Added mass coefficient (assume to be constant) 1

Drag coefficient (assume to be constant) 1.8

Note: Stable vibration can be achieved only if the simulated time is long enough.

Part two: One-dimensional convection diffusion problem

Assignment (please read the problem description before doing the assignment)

• Develop a finite difference method (FDM) formula for solving the equation. The procedure of deriving the FDM formula must be included in the report.

• Develop a MATLAB program and use this program to do the simulation and answer the following questions. Divide the soil depth into 200 cells when solving the problem numerically. You must attach the MATLAB program in your report.

• If a=0.0002 m2/s, find out when the temperature at the depth of 1.5 m reaches 25°C.

• If the temperature at the depth of 1.5 m reaches 25°C is defined as T1.5, do the simulations for a in the range of 0.0002 m2/s to 0.001 m2/s with an interval of 0.0001 m2/s and discuss how a on T1.5 and why.

Problem description

Note: In the paragraph below, the number j stands for the last digit of your student ID. For example, if your student ID is 2345673, 10+j = 13.

The water on the ground flows into the soil through the pore space of the soil (the space between soil particles). The flow velocity is very small of 0.002×(1+j) m/s and it does not change along the vertical direction. Initially, the temperature is (10+j) °C in the whole soil volume (including soil and the water in the soil). The temperature of the water on top of the ground level is 30°C and remains constant. The warm water moves into the soil and increases the temperature in the soil.

A coordinate x is defined and its origin is on the ground level and pointing downwards. The governing equation is the convection diffusion equation

?T ?T ?2T

??u ?? 2 ?t ?x ?x

where

T is the temperature in the soil t is the time u is the water flow velocity a is the heat diffusion coefficient (the influence of the soil has been considered).

?T

Boundary condition is T=30°C at the top boundary and ?0 at the bottom boundary (at x= 2 m).

?x

HNN215 Quality use of MedicinesAssessment Task 1 Case Study1500 words and 40% of overall gradePATIENT CASEPurpose of assessment taskNurses have an active role in medication management, including safe administration...Associations between Quantitative and Qualitative Job Insecurity and Well-beingResearch questionsDe Witte et al. (2010) investigated the association of employee’s perception of quantitative and qualitative...PACC6002 Individual Question A Semester 1 2019Individual Assignment ADue Date: 5.00 pm 21 August 2019Weighting: 5% to total assessmentsBackgroundChapter 1 of your text identifies a range of Australian...SER 301 ProjectIn line with the PODBL concepts, this unit is made up of a major project that is to be conducted individually while the project planning is a group exercise.MATLAB is a powerful programming...ASSESSMENT BRIEFSubject Code and Title ACCT2000: FINANCIAL ACCOUNTINGAssessment AssignmentIndividual/Group IndividualLength 1000 words (Part B only)Learning Outcomes 4. Understand methods of presenting...Assessment task 3: Response PaperObjective(s): This assessment task addresses the following subject learning objectives (SLOB):1, 2 and 3This assessment task contributes to the development of the following...Competitive Strategy H16006 —S3 / T2 - 2019 Dr Philip Booth - Holmes SydneyHandout 3b: HI6006 ASSIGNMENTS OVERVIEW - VERSION 2Assessment Tasks1. Individual Assignment: Key Strategy Development ToolsDue:...**Show All Questions**