{"id":1874,"date":"2017-10-25T16:11:05","date_gmt":"2017-10-25T05:11:05","guid":{"rendered":"http:\/\/casestudyhelp.com\/sample-questions\/?p=1874"},"modified":"2019-09-20T21:12:12","modified_gmt":"2019-09-20T11:12:12","slug":"matlab-programming-assignment-on-engineering-simulations-and-computations","status":"publish","type":"post","link":"https:\/\/casestudyhelp.com\/sample-questions\/matlab-programming-assignment-on-engineering-simulations-and-computations\/","title":{"rendered":"MATLAB Programming Assignment on Engineering Simulations and Computations"},"content":{"rendered":"

Get Questions in PDF Format ENG3104_assignment4<\/a><\/strong><\/p>\n

\"Hire<\/a><\/p>\n

Q 1\u00a0\u00a0 \u00a0(worth 115 marks)<\/strong><\/p>\n

1.1\u00a0\u00a0 \u00a0Introduction<\/strong><\/p>\n

Further analysis of the viaduct from Assignment 3 Question 3 is to be performed. An important reason for the construction of the Range Crossing is to reduce the gradient of the road for trucks. You will analyse the data in ass4q1in.csv to calculate gradients of the road with and without the viaduct. For this question, we are assuming that the viaduct follows a straight line between its end-points and the x<\/em>-coordinate has been reoriented to act along the length of the bridge.<\/p>\n

1.2\u00a0\u00a0 \u00a0Requirements<\/strong><\/p>\n

For this assessment item, you must perform hand calculations:<\/p>\n

1.<\/strong> Assuming that there is no viaduct, estimate the value of the gradient \u2202z\/\u2202x <\/em>at the second measurement point using backward, forward and central differences.<\/p>\n

2.<\/strong> Discuss which is the best result out of the backward, forward and central differences and why.<\/p>\n

3.<\/strong> Estimate the overall gradient of the viaduct. You must also produce MATLAB code which:<\/p>\n

4.<\/strong> Calculates backward, forward and central differences for \u2202z\/\u2202x <\/em>for the length of the road, assuming there is no viaduct. The calculations should be for each possible section (you will not calculate differences for every value of x<\/em>).<\/p>\n

5.<\/strong> Verifies the results from Requirement 4 using the results of Requirement 1.<\/p>\n

6.<\/strong> **Estimates the numerical error in the calculations of Requirement 4. Perform this using numerical derivatives for as many data points as possible; at the ends of the data, use the value closest to the end (e.g. for the 1st datum, use the second derivative calculated at the 2nd datum, and the third derivative calculated at the 3rd datum). \u00a0For data spaced evenly in x<\/em>, the necessary formulae are [the methods are O<\/em>(\u2206x<\/em>2) accurate]:<\/p>\n

\"equation\"<\/a><\/p>\n

For data unevenly \u00a0spaced in x<\/em>, the necessary formulae are [the methods are O<\/em>(\u2206x<\/em>) accu- rate]:<\/p>\n

\"1\"<\/a><\/p>\n

7.<\/strong> Plots the results from Requirements 4 and 6. Also plots \u2202z\/\u2202x <\/em>(correct Requirement 4 us- ing Requirement 6) together with the central difference from Requirement 4. On the plots of the gradient, also show the viaduct average gradient. Discusses the results, including which is the best method.<\/p>\n

8.<\/strong> Has appropriate comments throughout.<\/p>\n

1.3\u00a0\u00a0 \u00a0Assessment Criteria<\/strong><\/p>\n

Your code will be assessed using \u00a0the following scheme. Note that you are marked based on how well you perform for each category, \u00a0so the correct answer determined in a basic way will receive half marks and the correct answer determined using an excellent method\/code will receive full marks.<\/p>\n\n\n\n\n\n\n\n\n\n
Quality of hand calculations<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of difference calculations<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of verification<\/td>\n5 marks<\/td>\n<\/tr>\n
Quality of error calculations<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of the plots<\/td>\n10 marks<\/td>\n<\/tr>\n
Quality of the discussion<\/td>\n5 marks<\/td>\n<\/tr>\n
Quality of header(s) and comments<\/td>\n5 marks<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

<\/h2>\n

Q 2 \u00a0\u00a0\u00a0 (worth 90 marks)<\/strong><\/p>\n

2.1 \u00a0\u00a0\u00a0\u00a0 Introduction<\/strong><\/p>\n

The resistance of a resistor which has a constant cross-sectional area A <\/em>can be calculated using the formula<\/p>\n

R=\u03c1L\/A ,<\/em>\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0\u00a0<\/em>(5)<\/p>\n

where \u03c1 <\/em>is the resistivity and L <\/em>is the length of the resistor. If the geometric and\/or material properties are not constant, then (because each x<\/em>-location can be considered to be a resistor in series) the formula is:<\/p>\n

R=\u222bL \u03c1(x)\/\u00a0A(x)\u00a0dx\u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0 \u00a0\u00a0<\/em>(6)<\/p>\n

Thermoelectricity is a phenomenon that \u00a0relies on changing material properties: \u00a0electric currents produced by an applied voltage difference generate heat (which could replace chemical means), or temperature differences generate potential differences (most commonly applied in thermocouples). It is, however a very inefficient means of converting between heat and electricity (especially compared to heat engines using chemical reactions to drive an electric generator), which is fundamentally due to the relevant material properties.<\/p>\n

For your assignment, the following values are to be used for resistivities, radii and lengths. The properties for resistor 1 are:<\/p>\n

\u03c1<\/em>1\u00a0\u00a0\u00a0 = \u00a01.<\/em>68 \u00d7 <\/em>10\u2212<\/em>8\u00a0\u00a0\u2126m<\/p>\n

r<\/em>1\u00a0\u00a0\u00a0\u00a0 =\u00a0 \u00a01 mm<\/p>\n

L<\/em>1 \u00a0\u00a0\u00a0 = \u00a01 cm ,<\/em><\/p>\n

\u00a0resistor 2:<\/p>\n

\u03c1<\/em>2\u00a0\u00a0\u00a0 = \u00a0(1 + 0.<\/em>95179 \u00d7 <\/em>10) \u00d7 <\/em>10\u2212<\/em>8\u00a0\u00a0\u2126m<\/p>\n

r<\/em>2\u00a0\u00a0\u00a0\u00a0 =\u00a0 \u00a01 + 8.<\/em>8395 \u00f7 <\/em>10 mm<\/p>\n

L<\/em>2 \u00a0\u00a0\u00a0\u00a0 = \u00a01 + 6.<\/em>9848 \u00f7 <\/em>10 cm ,<\/em><\/p>\n

\u00a0and resistor 3:<\/p>\n

\u03c1<\/em>3\u00a0\u00a0\u00a0 = \u00a0(1 + 9.<\/em>198 \u00d7 <\/em>10) \u00d7 <\/em>10\u2212<\/em>8\u00a0\u00a0\u2126m<\/p>\n

r<\/em>3\u00a0\u00a0\u00a0\u00a0 =\u00a0 \u00a01 + 7.<\/em>2579 \u00f7 <\/em>10 mm<\/p>\n

L<\/em>3 \u00a0\u00a0\u00a0\u00a0 = \u00a01 + 2.<\/em>0632 \u00f7 <\/em>10 cm .<\/em><\/p>\n

In this question, you will calculate the resistances of a number of different resistor configu- rations:<\/p>\n

1.<\/strong> Resistor 1 with uniform cross-section.<\/p>\n

2.<\/strong> Resistor 1 with radius r<\/em>1 at the start and radius 2r<\/em>1at the end (the radius changes at a constant rate between the two ends).<\/p>\n

3.<\/strong> A resistor constructed by joining resistor 2 end-to-end with resistor 3 (the component resistors each have a constant cross-section).<\/p>\n

4.<\/strong> *A resistor constructed by joining resistor 2 end-to-end with resistor 3 (the radius changes at a constant rate along each of the two component resistors so that the radius is r<\/em>1 where they join).<\/p>\n

2.2\u00a0\u00a0 \u00a0Requirements<\/strong><\/p>\n

For this assessment item, you must perform hand calculations:<\/p>\n

1.<\/strong> Find the resistances of configurations 1\u20134 analytically.<\/p>\n

2.<\/strong> Use the trapezoidal method to determine the resistance of configurations 2 and 3 using four intervals along the total length of each configuration. \u00a0Validate the results with Requirement 1.<\/p>\n

You must also produce MATLAB code which:<\/p>\n

3.<\/strong> Repeats Requirement 2 and verifies the code.<\/p>\n

4.<\/strong> Uses the trapezoidal method to solve configurations 1\u20134 using appropriate intervals.<\/p>\n

5.<\/strong> Reports the values from Requirements 1 and 4.<\/p>\n

6.<\/strong> *Estimates, reports and discusses the overall numerical error in the calculations of each configuration in Requirement 4.<\/p>\n

7.<\/strong> Plots R<\/em>(x<\/em>) from Requirements 1 and 4.<\/p>\n

8.<\/strong> Has appropriate comments throughout.<\/p>\n

2.3\u00a0\u00a0 \u00a0Assessment Criteria<\/strong><\/p>\n

Your code will be assessed using \u00a0the following scheme. Note that you are marked based on how well you perform for each category, \u00a0so the correct answer determined in a basic way will receive half marks and the correct answer determined using an excellent method\/code will receive full marks.<\/p>\n\n\n\n\n\n\n\n
Quality of hand calculations<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of trapezoidal calculations and discussion<\/td>\n35 marks<\/td>\n<\/tr>\n
Quality of error calculations<\/td>\n10 marks<\/td>\n<\/tr>\n
Quality of the plots<\/td>\n10 marks<\/td>\n<\/tr>\n
Quality of header(s) and comments<\/td>\n5 marks<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

<\/h2>\n

3\u00a0\u00a0 \u00a0(worth 145 marks)<\/strong><\/p>\n

3.1\u00a0\u00a0 \u00a0Introduction<\/strong><\/p>\n

We return to Assignment 1. Your task is to simulate the droplet acceleration using the formulae and parameters supplied in Assignment 1, and starting with the initial \u00a0condition supplied in ass1in.csv.<\/p>\n

To obtain consistent results from your random number generation, you should initialise the seed to a fixed value using rng(seed). \u00a0For your assignment, the value is:<\/p>\n

seed = 8.2756<\/em><\/p>\n

\u00a03.2\u00a0\u00a0 \u00a0Requirements<\/strong><\/p>\n

For this assessment item, you must perform hand calculations:<\/p>\n

1.<\/strong> Use Euler\u2019s method to estimate the velocity for the 2nd data value in ass1in.csv using 5 equal steps; validate against that data value.<\/p>\n

You must also produce MATLAB code which:<\/p>\n

2.<\/strong> Repeats Requirement 1 to verify the method.<\/p>\n

3.<\/strong> Simulates the entire duration of ass1in.csv using appropriate parameters for the following methods:<\/p>\n

(a) Euler\u2019s method<\/p>\n

(b) \u00a0ode23 (c) ode45<\/p>\n

Reports the timestep used for Euler\u2019s method; if you use the default settings for the MATLAB ode solvers, report this, otherwise report what values you used for any changed settings. \u00a0Graphically validates the simulation results with \u00a0the data in ass1in.csv and reports the outcome of the validation.<\/p>\n

4.<\/strong> Uses Simulink to:<\/p>\n

(a) Simulate up to the 5th data value.<\/p>\n

(b) \u00a0***Simulate the entire duration of ass1in.csv.<\/p>\n

Plots the results and reports the final velocity within Simulink. Uses MATLAB to validate the results with the data in ass1in.csv.<\/p>\n

5.<\/strong> ***Simulates droplet collisions and accelerations.<\/p>\n

(a) Loads the data in ass2in.csv and multiplies the velocities by 50.<\/p>\n

(b) \u00a0Assigns each droplet an \u201cage\u201d. The age grows \u00a0so that, on average, for each 0.5 mm the droplet travels, it collides with another droplet. Initialise the age of each droplet to be a uniformly-distributed random number between 0 mm and 1 mm.<\/p>\n

(c) At each timestep, accelerates every droplet according to Euler\u2019s method as in Re- quirement 3.<\/p>\n

(d) \u00a0For each \u00a0timestep, increments every droplet\u2019s age. \u00a0If \u00a0the age is greater than a uniformly-distributed \u00a0random number between 0 and 1 mm, then the droplet is to collide after the new velocity has been calculated \u00a0for the timestep.<\/p>\n

(e) Pairs up droplets that are to collide and collide them according to the algorithm of Assignment 2. Randomly assign droplets to be paired with each other. \u00a0If there is an odd number of droplets marked for collision, choose one to not collide. Reset the age of the new, combined droplet according to the initialisation in Requirement 5b.<\/p>\n

(f ) Simulates up to the entire duration of ass1in.csv.<\/p>\n

(g) Plots the mean velocity, standard deviation of velocity and total number of droplets as functions of time.<\/p>\n

(h) \u00a0Reports to the Command Window the total number of droplets remaining at the end of the simulation. If there is only a single droplet remaining, reports when the last droplet collision occurred.<\/p>\n

6.<\/strong> Has appropriate comments throughout.<\/p>\n

3.3\u00a0\u00a0 \u00a0Assessment Criteria<\/strong><\/p>\n

Your code will be assessed using \u00a0the following scheme. Note that you are marked based on how well you perform for each category, \u00a0so the correct answer determined in a basic way will receive half marks and the correct answer determined using an excellent method\/code will receive full marks.<\/p>\n\n\n\n\n\n\n\n\n\n
Quality of hand calculations<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of Requirement 2<\/td>\n20 marks<\/td>\n<\/tr>\n
Quality of Requirement 3<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of using Simulink<\/td>\n30 marks<\/td>\n<\/tr>\n
Quality of simulating droplet collisions and acceleration<\/td>\n15 marks<\/td>\n<\/tr>\n
Quality of the plots<\/td>\n15 marks<\/td>\n<\/tr>\n
Quality of header(s) and comments<\/td>\n5 marks<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

 <\/p>\n