University of Wollongong

Faculty of Engineering and Information Sciences

ENGG952 Engineering Computing

Spring Session – 2017

Assignment 2

Rules:

1. The assignment may be completed individually or by a group up to 3 students. The group formation is your own responsibility. Members may be from the same or different tutorial groups.

2. No collaboration between groups permitted. Any case of plagiarism will be penalized and students should make themselves aware of the university policies regarding plagiarism (see subject outline under University and Faculty Policies)

3. The assignment is due in week 13 (Friday 27 October by 4 pm) and submitted to Central of EIS using then barcoded sheet (see subject outline under assignment submission). Late submission will incur penalty as described in the subject outline.

4. If the assignment is completed by group, statement indicating the effort or contribution to the assignment by each member and signed by all members must be included in the beginning of the report, or all students agree that they have contributed equally to the report –add a statement at the front of the report, signed by all members.

5. The assignment must be submitted as a formal report. You need to include the analysis of the problem, the procedure for solution, and discussion of the results. It is NOT sufficient to only provide the program and results. All matlab script files, function files and excel files must be included as attachments or main contents in the hard copy report.

6. In addition, all matlab codes and function files must be burned in a CD or USB and submitted with the report for checking.

Question 1: Motion of a Spring-Mass System (50%)

The motion of a damped spring-mass system is descirbed by the following ordinary differnetial equation:

d2x dx

m 2 ?c dt ? kx? 0 dt

where x is displacement from equlibrium position (meter), t is time (second), m is the mass and equal 20 kg, c is the damping coefficient (N.sec/meter) . The dampping coefficient, c, takes on two values of 5 (under damped), 40 (critically damped. The spring constant k = 20 N/meter. The initial velocity is zero, and the intial displacement x = 1 meter.

Figure 1. Damped spring-mass system

(a) Transform the problem into a system of two first order initial value ODEs. The report must clearly provide the detailed derivation of the technique.

(b) Solve for motion of a spring-mass system using the 2nd order RK Huen method over the time period 0?t ?5 sec with a step size ?t ?1.0 . You can calcualte the results manually or by matlab or by excel.

(c) Plot the displacment verus time for two values of the damping coefficient on the same figure, and discuss the results.

Question 2: Transient Heat Conduction (50%)

The non-dimensional form for the transient heat conduction in an insulated rod is

?2u ?u

2? t

?x ?

x is the nondimensional length, t is the nondimensional time, u is the nondimensional

temperature. This makes for the following boundary and initial conditions:

Boundary conditions u(0, t ) = 0.5

u(1, t ) = 2.0

Initial conditions u( x , 0) = 0.5 0? x 1

Figure 2. Heat conduction problem in an insulated rod.

x t T?To , in which L = the rod length, k = thermal Note: x?L , t ?( CL2/k) , u?TL ?To

?

conductivity of the rod material, ? = density, C = specific heat, To = temperature at x = 0, and TL = temperature at x = L.

Solve this nondimensional equation for the temperature distribution using explicit finitedifference method and implicit Crank-Nicholson method:

a) Write the finite-difference equation of the differential equation. The report must clearly provide the detailed derivation of the technique.

b) Programing to obtain the solution for time duration 0?t ?1. (For explicit finite difference method, you can use EXCEL or matlab. For implicit Crank-Nicholson method, you must use the matlab).

Please demenstrate that appropriate ?x and ?t are needed to solve u until steadystate solution is reached.

c) Plot the nondimensional temperature versus nondimensional length for a few typical values of nondimensional times, which can demonstrate the evloution of the tempeature at different time.

Dr. Hongtao Zhu

October 2017

Faculty of Engineering and Information Sciences

ENGG952 Engineering Computing

Spring Session – 2017

Assignment 2

Rules:

1. The assignment may be completed individually or by a group up to 3 students. The group formation is your own responsibility. Members may be from the same or different tutorial groups.

2. No collaboration between groups permitted. Any case of plagiarism will be penalized and students should make themselves aware of the university policies regarding plagiarism (see subject outline under University and Faculty Policies)

3. The assignment is due in week 13 (Friday 27 October by 4 pm) and submitted to Central of EIS using then barcoded sheet (see subject outline under assignment submission). Late submission will incur penalty as described in the subject outline.

4. If the assignment is completed by group, statement indicating the effort or contribution to the assignment by each member and signed by all members must be included in the beginning of the report, or all students agree that they have contributed equally to the report –add a statement at the front of the report, signed by all members.

5. The assignment must be submitted as a formal report. You need to include the analysis of the problem, the procedure for solution, and discussion of the results. It is NOT sufficient to only provide the program and results. All matlab script files, function files and excel files must be included as attachments or main contents in the hard copy report.

6. In addition, all matlab codes and function files must be burned in a CD or USB and submitted with the report for checking.

Question 1: Motion of a Spring-Mass System (50%)

The motion of a damped spring-mass system is descirbed by the following ordinary differnetial equation:

d2x dx

m 2 ?c dt ? kx? 0 dt

where x is displacement from equlibrium position (meter), t is time (second), m is the mass and equal 20 kg, c is the damping coefficient (N.sec/meter) . The dampping coefficient, c, takes on two values of 5 (under damped), 40 (critically damped. The spring constant k = 20 N/meter. The initial velocity is zero, and the intial displacement x = 1 meter.

Figure 1. Damped spring-mass system

(a) Transform the problem into a system of two first order initial value ODEs. The report must clearly provide the detailed derivation of the technique.

(b) Solve for motion of a spring-mass system using the 2nd order RK Huen method over the time period 0?t ?5 sec with a step size ?t ?1.0 . You can calcualte the results manually or by matlab or by excel.

(c) Plot the displacment verus time for two values of the damping coefficient on the same figure, and discuss the results.

Question 2: Transient Heat Conduction (50%)

The non-dimensional form for the transient heat conduction in an insulated rod is

?2u ?u

2? t

?x ?

x is the nondimensional length, t is the nondimensional time, u is the nondimensional

temperature. This makes for the following boundary and initial conditions:

Boundary conditions u(0, t ) = 0.5

u(1, t ) = 2.0

Initial conditions u( x , 0) = 0.5 0? x 1

Figure 2. Heat conduction problem in an insulated rod.

x t T?To , in which L = the rod length, k = thermal Note: x?L , t ?( CL2/k) , u?TL ?To

?

conductivity of the rod material, ? = density, C = specific heat, To = temperature at x = 0, and TL = temperature at x = L.

Solve this nondimensional equation for the temperature distribution using explicit finitedifference method and implicit Crank-Nicholson method:

a) Write the finite-difference equation of the differential equation. The report must clearly provide the detailed derivation of the technique.

b) Programing to obtain the solution for time duration 0?t ?1. (For explicit finite difference method, you can use EXCEL or matlab. For implicit Crank-Nicholson method, you must use the matlab).

Please demenstrate that appropriate ?x and ?t are needed to solve u until steadystate solution is reached.

c) Plot the nondimensional temperature versus nondimensional length for a few typical values of nondimensional times, which can demonstrate the evloution of the tempeature at different time.

Dr. Hongtao Zhu

October 2017

CASE STUDYOn 5 January 2013, Patient A, who was 81 years old, presented to a GP clinic complaining of two nights of breathlessness when lying flat and shortness of breath. On examination, Patient A was...HOLMES INSTITUTEFACULTY OFHIGHER EDUCATIONAssessment Details and Submission GuidelinesTrimester T1 2019Unit Code HI5019Unit Title Strategic Information Systems for Business and EnterpriseAssessment Type...ASSESSMENT BRIEFSubject Code and Title MKG203: Digital Marketing CommunicationsAssessment Individual DMC PlanIndividual/Group IndividualLength 2000 words (+/- 10%)Learning Outcomes a) Define digital marketing...COIT12209 - DATA SCIENCEAssessment 1 SpecificationDue date: Week 6 Monday (26 August 2019) 11:45 pm AEST ASSESSMENT1Weighting: 40%Length: No fixed lengthObjectivesAssessment 1 relates to unit learning...https://www.youtube.com/watch?v=_11xLlwKgWcSubject Learning Outcomes 1. Define health promotion and health education for the person individually, in community and within population across the lifespan;2....MGMT20132 INNOVATION AND SUSTAINABLE BUSINESS DEVELOPMENTAssessment 1 - Individual report: Value proposition proposal and reflectionAssessment 1 - Individual report: Value proposition proposal and reflectionDue...Assessment item 1 – Prioritisation of complex patient health issuesObjective(s): Learning objective: A, B, C, D and E Graduate Attribute: 2.0, 4.0 and 5.0Weight: 50%Due: Monday 26th August 2019 before...**Show All Questions**