代做MECH2700: Computational Engineering & Data Analysis (S2, 2024) Assignment 2: Convection-Diffus

MECH2700:  Computational  Engineering  &  Data  Analysis   (S2, 2024)

>  Seine  River  Contaminant  Transport

Assignment  2:  Convection-Diffusion  Modelling

Aim

The aim of this assignment is to synthesise your understanding of parabolic and hyperbolic partial differential equations by developing a numerical scheme that models a physical process where both diffusion and convection are present. Your solutions are to be implemented in Python and formally documented in a report.

The Problem

During the Paris Olympic Games, multiple events were swum in the Seine River, shown in Figure 1. Despite a massive effort to clean up pollution in the river, there were still concerns for athlete safety, in part due to fresh contamination introduced during the games. The fresh contamination generally entered the river via inflows from one of the banks, and is then transported through a combination  of water flow and turbulent diffusion.  Health risks arise if the contaminant concentration exceeds  the safe threshold in the area of the river where athletes compete.

Figure  1:    Seine  river,   Paris   (France).      Taken  3  August   2013.      Source:

http://www.panoramio.com/photo/97777268. Author: Tommie Hansen.

With a limited number of assumptions, the distribution of the contaminant in the river can be modelled as a combination of diffusive and convective processes. The unforced, two-dimensional convection- diffusion equation when convection only occurs in the x-direction is,

in which c is the contaminant concentration (kg per cubic meter) and is a function of x, y and t, Dturb  = 0.2 m2 /s is a turbulent diffusion coefficient, and v x is a x-velocity that is only dependent on y that is approximated as,

where V = 1.2 m/sand W1/2  = W/2.

Consider a two-dimensional stretch of the Seine River, shown in Figure 2, which is W = 150 m wide and L = 750 m long. Initially, the river is free of contamination, c(x,y, 0) = 0. At t = 0, a steady stream of contamination with a peak concentration of c max  = 1 kg/m3 at they = 0 bank flows into the stretch of river with a half Gaussian distribution through the upstream boundary (x = 0). The standard deviation of the distribution is σ = 10 m. The contamination is free to leave the domain from downstream boundary atx = L, while there is no convection or diffusion through the river banks at y = 0 andy = W. To model this, the appropriate boundary conditions are a Diriclet condition on the inflow boundary, x = 0 m,

zero gradient at the downstream boundary, x = L m,

and zero flux on the side boundaries

Your task is to predict the propagation of contamination in the river and report on the concentration in the competition area of the river, 0.25W < y < 0.75W, in the time following contamination release.

Figure 2: Problem domain.

The first step in the numerical solution method is to discretise the domain along its length and width into a grid of equally spaced nodes where spacing between nodes is △x = △y = L/(nx - 1). nx is the number of nodes in the x-direction. The solution at a node located atx = i△x, y = j△yand t = m△t is denoted c .

Tasks

1. Determine the discrete form. of the convection-diffusion equation, at solution node (i, j). To do this, approximate the second order space derivatives with the central difference expression (using c values at time step m) e.g.,

the first order space derivative with the upwind difference expression (using c values at time step m),

and the time derivative with a forward difference expression (equivalent to the forward Euler method for ODEs),

Present the discrete form. of the governing equation for an arbitrary interior node i,j.  Given the known errors in the finite difference approximations used (derived in an earlier section of MECH2700), what is the order of accuracy of the discrete governing equation in space and time?

2. Rearrange the discrete form. of the governing equation from Task 1 into an update equation for an

arbitrary internal node i, j (1 ≤ i ≤ nx  - 2, 1 ≤ j ≤ ny  - 2) of the form,

Next, determine the update equations for on the boundaries with i = 0, j = 0, j = ny  - 1 = jmax and i = nx  - 1 = imax. To do this, you will need to utilize the discrete form. of the boundary conditions,

Note that the corner nodes do not influence the rest of the solution for the numerical approach we are taking.

3. Write a Python function to determine a stable timestep for a numerical solution to this problem for arbitrary grid spacing △x = △y. This is to betaken to the following conservative combination of the hyperbolic and parabolic timestep expressions:

Here, conservatively take the Courant-Friedrichs-Lewy number to be CFL = 0.45. Note that when you find a solution with different resolution, the stable timestep will change.

4. Write a Python function to advance the solution for by one time step that takes the current array of c, t, △t and △x as inputs.  Utilize this function within a Python program that advances the solution in the entire domain from the initial condition at t  =  0 up to a maximum simulation time of t  = tmax. Record the full solution array every ≈200 s (to within △t).  Also record the concentration at the solution node nearest the edge of the competition zone at the end of the river segment x edge→ = (x, y) = (L, W/4) at every time step.

5. Use your numerical scheme to predict the contaminant concentration at x edge→ for t max   =  1000 seconds afer contaminant release.  Be sure to use a grid resolution that reasonably represents the contamination inflow profile. Generate contour plots of the concentration distribution in the domain at 200 s, 400 s, 600 s, 800 s and 1000 safer release. Also create a graph of the concentration atx edge→ with time. Does the concentration exceed a threshold value of c safe  = 0.01 in this period? Suggest and justify how the boundary of the competition zone could be adjusted so that c(x,y, t) ≤ c safe for all time. To assess uncertainty in the turbulent diffusion coefficient, how do your results change if the diffusion coefficient is doubled to 0.4 m2 /s?

6. Repeat the above calculation for at least one other choice of grid spacing. Compare your results and comment on the grid-dependence of your numerical predictions. Given the trend in your solution with grid resolution, is the safe competition boundary you have established conservative (e.g. the errors in the numerical method lead to the boundary being further away from the contaminant source that it needs to be in reality)? Hints: Be reasonable, here.  If you try to run with a grid spacing of 1 mm you’re most likely going to have a bad time. Also, include something in your code that tells you how far into your simulation you are.

The Report

Document your work in a formal report that includes, but is not necessarily limited to, the following:

i. Introduction: A brief description of the problem you have been asked to solve;

ii. Methodology: The definition of your approach to solving this problem, including all working and relevant assumptions;

iii. Results: The appropriate presentation of results, which might include figures, graphs and tables;

iv. Concluding Remarks: A critical discussion of your approach to the problem and your ndings. At a minimum this should include comment on the stability, convergence, accuracy and compu-tational efficiency of your scheme.

Submission

Work individually or in groups of two for this assignment. Submit a single report for the group, de- tailing your analysis, code and results for Tasks 1 to 7. This report should be brief but well structured, so that the marker can easily interpret your results. Include your Python code and combine it all together into a single pdf le with a cover page that clearly tells us who you are.  Scanned, hand- written documents with appropriate cut-and-paste items are perfectly acceptable, but plots should be computer-generated.

For assessment of your work, the markers will be looking for answers to the following questions:

• Have you done what the tasks require? Note that your report should be well structured and tidy so that the marker can easily nd your answers.  Having your plots correctly titled and axes labeled is essential.

• Are your answers accurate?

• Is your code well-structured and clearly written, with suitable documentation comments?

• Does your code correspond to the results shown?


热门主题

课程名

mktg2509 csci 2600 38170 lng302 csse3010 phas3226 77938 arch1162 engn4536/engn6536 acx5903 comp151101 phl245 cse12 comp9312 stat3016/6016 phas0038 comp2140 6qqmb312 xjco3011 rest0005 ematm0051 5qqmn219 lubs5062m eee8155 cege0100 eap033 artd1109 mat246 etc3430 ecmm462 mis102 inft6800 ddes9903 comp6521 comp9517 comp3331/9331 comp4337 comp6008 comp9414 bu.231.790.81 man00150m csb352h math1041 eengm4100 isys1002 08 6057cem mktg3504 mthm036 mtrx1701 mth3241 eeee3086 cmp-7038b cmp-7000a ints4010 econ2151 infs5710 fins5516 fin3309 fins5510 gsoe9340 math2007 math2036 soee5010 mark3088 infs3605 elec9714 comp2271 ma214 comp2211 infs3604 600426 sit254 acct3091 bbt405 msin0116 com107/com113 mark5826 sit120 comp9021 eco2101 eeen40700 cs253 ece3114 ecmm447 chns3000 math377 itd102 comp9444 comp(2041|9044) econ0060 econ7230 mgt001371 ecs-323 cs6250 mgdi60012 mdia2012 comm221001 comm5000 ma1008 engl642 econ241 com333 math367 mis201 nbs-7041x meek16104 econ2003 comm1190 mbas902 comp-1027 dpst1091 comp7315 eppd1033 m06 ee3025 msci231 bb113/bbs1063 fc709 comp3425 comp9417 econ42915 cb9101 math1102e chme0017 fc307 mkt60104 5522usst litr1-uc6201.200 ee1102 cosc2803 math39512 omp9727 int2067/int5051 bsb151 mgt253 fc021 babs2202 mis2002s phya21 18-213 cege0012 mdia1002 math38032 mech5125 07 cisc102 mgx3110 cs240 11175 fin3020s eco3420 ictten622 comp9727 cpt111 de114102d mgm320h5s bafi1019 math21112 efim20036 mn-3503 fins5568 110.807 bcpm000028 info6030 bma0092 bcpm0054 math20212 ce335 cs365 cenv6141 ftec5580 math2010 ec3450 comm1170 ecmt1010 csci-ua.0480-003 econ12-200 ib3960 ectb60h3f cs247—assignment tk3163 ics3u ib3j80 comp20008 comp9334 eppd1063 acct2343 cct109 isys1055/3412 math350-real math2014 eec180 stat141b econ2101 msinm014/msing014/msing014b fit2004 comp643 bu1002 cm2030
联系我们
EMail: 99515681@qq.com
QQ: 99515681
留学生作业帮-留学生的知心伴侣!
工作时间:08:00-21:00
python代写
微信客服:codinghelp
站长地图