Do not copy paste: AF

want view factors from all the surfaces $.$ combine my two codes so I can have view factors from all surfaces

Using your Monte Carlo code, determine the view factor between a control surface on the left side of

the V$-$groove with every control surface on the opposite side of the V$-$groove. Do this for every control

surface on the left side of the V$-$groove. This means you should know the view factor between any two

control surfaces on opposing sides of the V$-$groove.

combine the codes below

import random

import numpy as np

import math

import matplotlib.pyplot as plt

# Numerical solution

def rotate$\_$vector$($vector$,$ angle$)$:

x$,$ y $=$ vector

theta $=$ math.radians$($angle$)$

x$\_$rotated $=$ math.cos$($theta$)*$ x $-$ math.sin$($theta$)*$ y

y$\_$rotated $=$ math.sin$($theta$)*$ x $+$ math.cos$($theta$)*$ y

return $[$x$\_$rotated, y$\_$rotated$]$

def find$\_$intersection$($x$,$ y$,$ xhat, yhat, surface$1,$ surface$2)$:

xi$\_$s $=$ surface$1[\text{'}$X$0\text{'}]+$ surface$1[\text{'}$XHAT$\text{'}]*$ x

yi$\_$s $=$ surface$1[\text{'}$Y$0\text{'}]+$ surface$1[\text{'}$YHAT$\text{'}]*$ y

xi$\_$r $=$ surface$2[\text{'}$X$0\text{'}]+$ surface$2[\text{'}$XHAT$\text{'}]*$ x

yi$\_$r $=$ surface$2[\text{'}$Y$0\text{'}]+$ surface$2[\text{'}$YHAT$\text{'}]*$ y

return xi$\_$s$,$ yi$\_$s$,$ xi$\_$r$,$ yi$\_$r

def monte$\_$carlo$\_$ray$\_$tracing$($num$\_$rays$)$:

# Step $1$: Define surfaces

surface$\_$length $=1$ # Length of the surfaces

surface$\_$distance $=1$ # Distance between the surfaces

surface$1=\{\text{'}$X$0\text{'}$: $0,\text{'}$Y$0\text{'}$: $0,$ 'XHAT': surface$\_$length, 'YHAT': $0\}$

surface$2=\{\text{'}$X$0\text{'}$: $0,\text{'}$Y$0\text{'}$: surface$\_$distance, 'XHAT': surface$\_$length, 'YHAT': $0\}$

rays$\_$surface$1=[]$

for $\_$ in range$($num$\_$rays$)$:

# Step $3$: Find a random spot for each ray

ts $=$ random.random$()$

tr $=$ random.random$()$

# Step $4$: Assign emission locations

x$\_$ray$\_$s$1,$ y$\_$ray$\_$s$1=$ surface$1[\text{'}$X$0\text{'}]+$ surface$1[\text{'}$XHAT$\text{'}]*$ ts$,$ surface$1[\text{'}$Y$0\text{'}]+$ surface$1[\text{'}$YHAT$\text{'}]*$ ts

surface$1\_$normal $=$ rotate$\_$vector$([$surface$1[\text{'}$XHAT$\text{'}],$ surface$1[\text{'}$YHAT$\text{'}]],90)$

R$\_$theta$\_$s$1,$ R$\_$phi$\_$s$1=$ random.random$(),$ random.random$()$

theta$\_$s$1=$ math.asin$($R$\_$theta$\_$s$1)$

if R$\_$phi$\_$s$1<0\mathrm{.}5$:

theta$\_$s$1=-$theta$\_$s$1$

xhat$\_$ray$\_$s$1=$ math.cos$($theta$\_$s$1)*$ surface$1\_$normal$[0]-$ math.sin$($theta$\_$s$1)*$ surface$1\_$normal$[1]$

yhat$\_$ray$\_$s$1=$ math.sin$($theta$\_$s$1)*$ surface$1\_$normal$[0]+$ math.cos$($theta$\_$s$1)*$ surface$1\_$normal$[1]$

rays$\_$surface$1.$append$(\{\text{'}$X$\text{'}$: x$\_$ray$\_$s$1,\text{'}$Y$\text{'}$: y$\_$ray$\_$s$1,$ 'XHAT': xhat$\_$ray$\_$s$1,$ 'YHAT': yhat$\_$ray$\_$s$1,$ 'Reflections': $0,\text{'}$ts$\text{'}$: None, $\text{'}$tr$\text{'}$: None$\})$

# Step $6$: Find intersection points and save ts and tr for each surface

for ray in rays$\_$surface$1$:

xi$\_$s$,$ yi$\_$s$,$ xi$\_$r$,$ yi$\_$r $=$ find$\_$intersection$($ray$[\text{'}$X$\text{'}],$ ray$[\text{'}$Y$\text{'}],$ ray$[\text{'}$XHAT$\text{'}],$ ray$[\text{'}$YHAT$\text{'}],$ surface$1,$ surface$2)$

denominator $=$ ray$[\text{'}$YHAT$\text{'}]*$ surface$2[\text{'}$XHAT$\text{'}]-$ ray$[\text{'}$XHAT$\text{'}]*$ surface$2[\text{'}$YHAT$\text{'}]$

if denominator $!=0$:

tr $=($xi$\_$s $*$ surface$2[\text{'}$YHAT$\text{'}]-$ yi$\_$s $*$ surface$2[\text{'}$XHAT$\text{'}]+$ surface$2[\text{'}$X$0\text{'}]*$ surface$2[\text{'}$YHAT$\text{'}]-$ surface$2[\text{'}$Y$0\text{'}]*$ surface$2[\text{'}$XHAT$\text{'}])/$ denominator

ts $=($ray$[\text{'}$X$\text{'}]+$ ray$[\text{'}$XHAT$\text{'}]*$ tr $-$ surface$1[\text{'}$X$0\text{'}])/$ surface$1[\text{'}$XHAT$\text{'}]$

ray$[\text{'}$ts$\text{'}]=$ ts

ray$[\text{'}$tr$\text{'}]=$ tr

# Step $7$: Determine impact surface

for ray in rays$\_$surface$1$:

ts $=$ ray$[\text{'}$ts$\text{'}]$

tr $=$ ray$[\text{'}$tr$\text{'}]$

if $0<=$ ts $<=1$:

if ray$[\text{'}$Reflections$\text{'}]==0$ or ts $<$ rays$\_$surface$1[$ray$[\text{'}$Reflections$\text{'}]-1][\text{'}$ts$\text{'}]$:

ray$[\text{'}$ImpactedSurface$\text{'}]=1$

else:

ray$[\text{'}$ImpactedSurface$\text{'}]=2$

else:

ray$[\text{'}$ImpactedSurface$\text{'}]=0$

# Print the results

hit$\_$count $=$ sum$($ray$[\text{'}$ImpactedSurface$\text{'}]==1$ for ray in rays$\_$surface$1)$

miss$\_$count $=$ num$\_$rays $-$ hit$\_$count

#print$($f$"$Number of rays that hit the surface: $\{$hit$\_$count$\}")$

#print$($f$"$Number of rays that did not hit the surface: $\{$miss$\_$count$\}")$

# view factor

view$\_$factor $=$ hit$\_$count $/$ num$\_$rays

print$($f$"$View Factor: $\{$view$\_$factor$\}")$

# Example usage

monte$\_$carlo$\_$ray$\_$tracing$(100000)$

second code

import numpy as np

import matplotlib.pyplot as plt

# Problem Definitions

phi $=120$ # phi $($deg$)$

phi$\_$radians $=$ phi $*$ np$.$pi $/180$ # phi $($rad$)$

Nsurf $=10$ # number of surfaces PER SIDE

length $=0\mathrm{.}1$ # total length per side $($m$)$

dL $=$ length $/$ Nsurf # length of one surface

# Set Initial Conditions

X$0=0$

Y$0=0$

# Calculate Surfaces for First Side

for i in range$($Nsurf$)$:

XHAT $=$ X$0+$ dL $*$ np$.$sin$($phi$\_$radians $/2)$

YHAT $=$ Y$0-$ dL $*$ np$.$cos$($phi$\_$radians $/2)$ # MINUS for $-$y

temp $=$ "surface" $+$ str$($i$)+"=\{\text{'}$X$0\text{'}$:X$0,\text{'}$Y$0\text{'}$:Y$0,$'XHAT':XHAT,'YHAT':YHAT$\}"$

exec$($temp$)$

X$0=$ XHAT

Y$0=$ YHAT

# Calculate Surfaces for Second Side

for i in range$($Nsurf$)$:

XHAT $=$ X$0+$ dL $*$ np$.$sin$($phi$\_$radians $/2)$

YHAT $=$ Y$0+$ dL $*$ np$.$cos$($phi$\_$radians $/2)$ # PLUS for $+$Y

temp $=$ "surface" $+$ str$($i $+$ Nsurf$)+"=\{\text{'}$X$0\text{'}$:X$0,\text{'}$Y$0\text{'}$:Y$0,$'XHAT':XHAT,'YHAT':YHAT$\}"$

exec$($temp$)$

X$0=$ XHAT

Y$0=$ YHAT

# Plot Surfaces

for i in range$(2*$ Nsurf$)$:

temp $="$plt$.$plot$([$surface$"+$ str$($i$)+"[\text{'}$X$0\text{'}],$ surface" $+$ str$($i$)+"[\text{'}$XHAT$\text{'}]],[$surface$"+$ str$($i$)+"[\text{'}$Y$0\text{'}],$ surface" $+$ str$($i$)+"[\text{'}$YHAT$\text{'}]],$ marker$=\text{'}$o$\text{'})"$

exec$($temp$)$

plt$.$show$()$

the output should have similar values or I'll down vote

F $0->10=0\mathrm{.}000862$

F $0->11=0 There are 3 Steps involved in it Step: 1 Get Instant Access to Expert-Tailored Solutions See step-by-step solutions with expert insights and AI powered tools for academic success Step: 2 Step: 3 Ace Your Homework with AI Get the answers you need in no time with our AI-driven, step-by-step assistance Get Started Recommended Textbook for Database Driven Web Sites Authors: Joline Morrison, Mike Morrison 2nd Edition ? 061906448X, 978-0619064488 More Books Students also viewed these Databases questions Question ★★★★★ An integrated test facility (ITF) is a method used by both internal and external auditors for testing IT-based system controls. Discuss the advantages and disadvantages of implementing an ITF. Answered: 1 week ago Question ★★★★★ Use MIPS not x86 4. Show how the following four consecutive instructions move through each stage of the five stage pipeline, similar to the example on slide 11 of lecture 17.This pipeline does not... Answered: 1 week ago Question ★★★★★ 123. A circular target has radius 1 foot. Assume that you hit the target (we shall ignore misses) and that the probability of hitting any region of the target is proportional to the regions area. If... Answered: 1 week ago Question ★★★★★ Presented below are incomplete manufacturing cost data. Determine the missing amounts for three differentsituations. Direct Materials Used Direct Labor Used Factory Overhead Total Manufacturing Costs... Answered: 1 week ago Question ★★★★★ Henderson Office Supply is considering a more liberal credit policy to increase sales but expects that 9 percent of the new accounts will be uncollectible. Collection costs are 6 percent of new... Answered: 1 week ago Question ★★★★★ . Safari File Edit View History Bookmarks Window Help G Q 8 Sun Jul 16 12:59 PM GO v2.cengagenow.com Problem Set: Module 4 Book Show Me How ? 1. PR.06.02A.ALGO Contribution Margin 2. EX.06.01 Molly... Answered: 1 week ago Question ★★★★★ Provide real-world examples of companies that prioritize attitude when recruiting people instead for their skills? Connect it with this phrase: "We recruit for attitude and train for skill," Atul... Answered: 1 week ago Question ★★★★★ 1. How long should a clean status last for? 2.What should happen if a clearance check is required when no QC personnel are on duty? 3. How frequently should staff wash hands or sanitise their gloves Answered: 1 week ago Question ★★★★★ Identify a product that you believe would benefit from being redesigned. What are the current problems with the product? What specific design improvements do you recommend? Why and how will your... Answered: 1 week ago Question ★★★★★ Find your own example of greenwashing. Explain why you think it constitutes greenwashing. If you were the CEO, would you have authorized the ad you chose? Why or why not? Answered: 1 week ago Question ★★★★★ Q22 (a) has been address You receive a task to create a new subnet with two servers. This subnet will be connected to an existing router. The new subnet allocated with network 192.168.20.24/29. The... Answered: 1 week ago Question ★★★★★ 9-16 How did WIM manage the change from a local MRP system to Epicors ERP system? Welding Industries Malaysia (WIM) is the only local Malasian company that manufactures welding machines for... Answered: 1 week ago Question ★★★★★ 9-15 How has Epicors ERP system helped WIM to make its production process more efficient? Welding Industries Malaysia (WIM) is the only local Malasian company that manufactures welding machines for... Answered: 1 week ago Question ★★★★★ 9-12 In MyMISLab, you will find a Collaboration and Teamwork Project dealing with the concepts in this chapter. You will be able to use Google Drive, Google Docs, Google Sites, Google+, or other... Answered: 1 week ago Previous Question Next Question View Answer in SolutionInn App Services Sitemap Fun Definitions Become Tutor Study Help Categories Recent Questions Expert Questions Company Info Security Copyrights Privacy Policy Terms & Conditions SolutionInn Fee Scholarship Get In Touch About Us Contact Us Career Jobs FAQ Campus Ambassador Secure Payment Download Our App © 2024 SolutionInn. All Rights Reserved $