ME 4943 Course Project - Part I
|
|
- Samantha Anna Green
- 6 years ago
- Views:
Transcription
1 ME 4943 Course Project - Part I Numerical Implementation of the Source Panel Method by Jeffrey A. Kornuta Mechanical Engineering Louisiana State University October 9, 27 1
2 Contents 1 Verification of Algorithm N = N = Method Applied to Arbitrary Shape The Square Appendices 8 A MATLAB Code 8 A.1 Part I A.2 Part II List of Figures 1 Visualizations of (a) N = 3 and (b) N = 5 for N = 3, R = 1 m, and V = 2 m/s Visualizations of (a) v(θ) and (b) C p (θ) for N = 3, R = 1 m, and V = 2 m/s Visualizations of (a) v(θ) and (b) C p (θ) for N = 5, R = 1 m, and V = 2 m/s Visualization of square with N = 8 panels Visualizations of (a) v(θ) and (b) C p (θ) for N = 8, (2 2) m, and V = 2 m/s
3 1 Verification of Algorithm Two cases are considered: 1. N = 3, R = 1 m, and V = 2 m/s 2. N = 5, R = 1 m, and V = 2 m/s. 1 Geometry, N = 3 1 Geometry, N = !.2.4!.4.6!.6.8! (a) N = 3!1!1! (b) N = 5 Figure 1: Visualizations of (a) N = 3 and (b) N = 5 for N = 3, R = 1 m, and V = 2 m/s. The following governing equations are utilized by the source panel method to generate the velocity and pressure coefficient distribution at the surface of a shape: V cos β i + λ N i 2 + λ j 2π I ij =, where (1) j=1,j i I ij = (ln r ij ) ds j. (2) n i Solving for λ = {λ 1, λ 2,..., λ N } using Eq.1 and Eq.2, one can solve for the velocities at the surface of the circle using the following equation: V i = V sin β i + N j=1,j i 3 j λ j 2π s j s j (ln r ij ) ds j. (3)
4 Finally, the pressure coefficients can be found using the following: C p,i = 1 ( Vi V ) 2. (4) Refer to pages in Anderson s Fundamentals of Aerodynamics, Fourth edition for details regarding notation. 1.1 N = 3 4 Velocity v(θ) 1 Pressure coefficient C p (θ) Velocity (m/s) 1 C p θ (degrees) (a) v(θ) θ (degrees) (b) C p (θ) Figure 2: Visualizations of (a) v(θ) and (b) C p (θ) for N = 3, R = 1 m, and V = 2 m/s. Note that with only three panels, the geometry forms a triangle [Figure 1a]. Although when the analytical and numerical values of v(θ) and C p (θ) are plotted against each other, the marked points lie on the curve [Figure 2], lines connecting the numerically calculated points would reveal an extremely inaccurate representation of the analytical solutions. Thus, more panels would be needed to give a good approximation of the velocity and pressure coefficient distributions. 4
5 1.2 N = 5 4 Velocity v(θ) 1 Pressure coefficient C p (θ) Velocity (m/s) 1 C p θ (degrees) (a) v(θ) θ (degrees) (b) C p (θ) Figure 3: Visualizations of (a) v(θ) and (b) C p (θ) for N = 5, R = 1 m, and V = 2 m/s. Note that with 5 panels, the geometry closely represents a circle [Figure 1b]. When the analytical and numerical values of v(θ) and C p (θ) are plotted against each other, the marked points not only lie on the curve [Figure 3], but lines that would connect the numerically calculated points would reveal an accurate representation of the analytical solutions. Thus, this number of panels gives a good approximation of the velocity and pressure coefficient distributions. In fact, reducing the number of panels to as little as 2 would yield a fair representation. These results suggest that the algorithm used for this scenario is valid. 5
6 2 Method Applied to Arbitrary Shape 2.1 The Square 1.5 Geometry, N = x (m) y (m) Figure 4: Visualization of square with N = 8 panels. The algorithm previously defined [Appendix A.1] will now be applied to an arbitrary shape in this case, a (2 2) m square [Figure 4, Appendix A.2]. For this square, the number of panels N = 8, and from the plots generated, this number seems to indicate that the surface is well-resolved [Figure 5]. As one may notice, the overall characteristics of these curves coincide with the results produced by the circle [Figure 3]. The stagnation points are located in the same locations, and the average maxima and minima values remain similar. These similarities suggest that the results for the square are sound. The only noticeable peculiarities lie with the values obtained for v(θ) and C p (θ) around the corners of the square. At these corners, lim R R, where R is the radius of curvature at the edges; thus, values of v (and correspondingly C p ) yield singularities at these points. These points are clearly seen in Figure 5 where there exist sharp spikes. 6
7 6 Velocity v(θ) 1 Pressure coefficient C p (θ) Velocity (m/s) C p θ (degrees) (a) v(θ) θ (degrees) (b) C p (θ) Figure 5: Visualizations of (a) v(θ) and (b) C p (θ) for N = 8, (2 2) m, and V = 2 m/s. In reality, this type of behavior would not occur. Since singularities of this kind are rarely seen in nature, these results suggest that a deficiency in the algorithm exists. In brief, this deficiency is characterized by two things: the absence of infinitely small radii of curvature, and the presence of viscous effects in the fluid. In actuality, edges cannot physically have an infinitely small radius of curvature, so these singularities would not exist in this scenario since R >. Also, the boundary layer of the surrounding fluid would separate at these sharp edges if they were indeed that sharp. Viscosity would cause the boundary layer to peel off at the corners, forming a wake. This scenario is void of singularities and would prove to be a much more reliable model to represent this shape, as compared to the inviscid model described by the source panel method. 7
8 Appendices A MATLAB Code A.1 Part I % Jeff Kornuta % ME Project Part I-1 clear all; clf; clc; % ask user for required information N = input( Enter number of panels: ); fprintf( \n ); U = input( Enter freestream velocity [m/s]: ); fprintf( \n ); R = 1; % plot a plane-jane circle for t = 1:361 theta = (t-1)*pi/18; u(t) = cos(theta); v(t) = sin(theta); %hold off %plot(u,v, r ) %axis square % now plot polygon (endpoints) hold on for i = 1:(N+1) dtheta = 2*pi/N; beta(i) = pi-(i-1)*dtheta; Phi(i) = beta(i)-pi/2; theta(i) = beta(i)+dtheta/2; tempx(i) = R*cos(theta(i)); tempy(i) = R*sin(theta(i)); %plot(tempx,tempy) 8
9 %pause(.25) % resize X, Y, and Beta to have one less element X(i) = tempx(i); Y(i) = tempy(i); if (beta(i) < ) beta(i) = beta(i)+2*pi; Beta(i) = beta(i)*18/pi; % plot midpoints of panels if (i ~= N) x(i) = (X(i+1)+X(i))/2; y(i) = (Y(i+1)+Y(i))/2; else x(i) = (X(1)+X(N))/2; y(i) = (Y(1)+Y(N))/2; % find length of segments S = sqrt((x(1)-x(n))^2+(y(1)-y(n))^2); %plot(x(i),y(i), o ) %pause(.25) % that was fun, now let s find I(i,j) for j = 1:N if (i ~= j) A(i,j) = -(x(i)-x(j))*cos(phi(j)) - (y(i)-y(j))*sin(phi(j)); B(i,j) = (x(i)-x(j))^2 + (y(i)-y(j))^2; % don t divide by zero if (B(i,j) == ) B(i,j) = 1e-12; C(i,j) = sin(phi(i)-phi(j)); D(i,j) = (y(i)-y(j))*cos(phi(i)) - (x(i)-x(j))*sin(phi(i)); E(i,j) = (x(i)-x(j))*sin(phi(j)) - (y(i)-y(j))*cos(phi(j)); 9
10 % please, stop dividing by zero if (E(i,j) == ) E(i,j) = 1e-12; I(i,j) = C(i,j)/2*log((S^2+2*A(i,j)*S+B(i,j))/B(i,j)) +... (D(i,j)-A(i,j)*C(i,j))/E(i,j)*(atan((S+A(i,j))/E(i,j))-... atan(a(i,j)/e(i,j))); else I(i,j) = pi; % set up system of equations for j = 1:N M(i,j) = I(i,j)/(2*pi); b(i,1) = -U*cos(beta(i)); % get those lambdas & check em lambda = inv(m)*b Sum = sum(lambda) % find velocities & pressure coefficients Vs = ; for j = 1:N if (i ~= j) M2(i,j) = ((D(i,j)-A(i,j)*C(i,j))/(2*E(i,j)))... *log((s^2+2*a(i,j)*s+b(i,j))/b(i,j))-c(i,j)... *(atan((s+a(i,j))/e(i,j))-atan(a(i,j)/e(i,j))); Vs = Vs+lambda(j)/(2*pi)*M2(i,j); V(i) = -(U*sin(beta(i))+Vs); Cp(i) = 1-(V(i)/U)^2; 1
11 % plot velocities & pressure coefficients for i = 1:361 t = (i-1)*pi/18; Cpreal(i) = 1-4*(sin(t))^2; Vreal(i) = -2*U*sin(t); clear t; t=[:36]; hold off %plot(t,vreal) plot(t,cpreal) axis square hold on %plot(beta,v, + ) plot(beta,cp, + ) xlabel( \theta (degrees) ); %ylabel( Velocity (m/s) ); %title( Velocity v(\theta) ); ylabel( C_p ); title( Pressure coefficient C_p(\theta) ); 11
12 A.2 Part II % Jeff Kornuta % ME Project Part I-2 clear all; clf; clc; % set required information N = 8; U = 2; S =.1; % plot a square hold off t = [-1:.1:1]; one = ones(1,21); negone = -1*ones(1,21); % plot(negone,t) % hold on % plot(t,one) % plot(one,t) % plot(t,negone) % axis([ ]) % axis square % find endpoints, midpoints, and angles for i = 1:11 X(i) = -1; Y(i) = (i-1)*s; if (i ~= 11) x(i) = -1; y(i) = Y(i)+S/2; beta(i) = pi; Phi(i) = pi/2; for i = 12:31 X(i) = -1+S*(i-11); Y(i) = 1; 12
13 x(i-1) = X(i)-S/2; y(i-1) = Y(i); beta(i-1) = pi/2; Phi(i-1) = ; for i = 32:51 X(i) = 1; Y(i) = 1-S*(i-31); x(i-1) = 1; y(i-1) = Y(i)+S/2; beta(i-1) = ; Phi(i-1) = 3*pi/2; for i = 52:71 X(i) = 1-S*(i-51); Y(i) = -1; x(i-1) = X(i)+S/2; y(i-1) = Y(i); beta(i-1) = 3*pi/2; Phi(i-1) = pi; for i = 72:8 X(i) = -1; Y(i) = -1+S*(i-71); x(i-1) = X(i); y(i-1) = Y(i)-S/2; beta(i-1) = pi; Phi(i-1) = pi/2; % get those last ones x(n) = -1; y(n) = -S/2; beta(n) = pi; Phi(N) = pi/2; % that was fun, now find theta values to keep things tidy for i = 1:2 theta(i) = pi-atan(abs(y(i))/abs(x(i))); for i = 21:4 13
14 theta(i) = atan(y(i)/x(i)); for i = 41:6 theta(i) = 2*pi-atan(abs(y(i))/abs(x(i))); for i = 61:8 theta(i) = pi+atan(abs(y(i))/abs(x(i))); % plot endpoints & midpoints % plot(x,y, r+,x,y, g. ) % axis square % axis([ ]) % title( Geometry, N = 8 ) % xlabel( y (m) ) % ylabel( x (m) ) %that was fun, now let s find I(i,j) for j = 1:N if (i ~= j) A(i,j) = -(x(i)-x(j))*cos(phi(j)) - (y(i)-y(j))*sin(phi(j)); B(i,j) = (x(i)-x(j))^2 + (y(i)-y(j))^2; %don t divide by zero if (B(i,j) == ) B(i,j) = 1e-12; C(i,j) = sin(phi(i)-phi(j)); D(i,j) = (y(i)-y(j))*cos(phi(i)) - (x(i)-x(j))*sin(phi(i)); E(i,j) = (x(i)-x(j))*sin(phi(j)) - (y(i)-y(j))*cos(phi(j)); %please, stop dividing by zero if (E(i,j) == ) E(i,j) = 1e-12; I(i,j) = C(i,j)/2*log((S^2+2*A(i,j)*S+B(i,j))/B(i,j)) +... (D(i,j)-A(i,j)*C(i,j))/E(i,j)*(atan((S+A(i,j))/E(i,j))-... atan(a(i,j)/e(i,j))); else I(i,j) = pi; 14
15 %set up system of equations for j = 1:N M(i,j) = I(i,j)/(2*pi); b(i,1) = -U*cos(beta(i)); %get those lambdas & check em lambda = inv(m)*b Sum = sum(lambda) %find velocities & pressure coefficients Vs = ; for j = 1:N if (i ~= j) M2(i,j) = ((D(i,j)-A(i,j)*C(i,j))/(2*E(i,j)))... *log((s^2+2*a(i,j)*s+b(i,j))/b(i,j))-c(i,j)... *(atan((s+a(i,j))/e(i,j))-atan(a(i,j)/e(i,j))); Vs = Vs+lambda(j)/(2*pi)*M2(i,j); V(i) = -(U*sin(beta(i))+Vs); Cp(i) = 1-(V(i)/U)^2; axis square subplot(1,2,1),plot(theta*18/pi,v, o ) xlabel( \theta (degrees) ); ylabel( Velocity (m/s) ); title( Velocity v(\theta) ); subplot(1,2,2),plot(theta*18/pi,cp, sq ) xlabel( \theta (degrees) ); ylabel( C_p ); title( Pressure coefficient C_p(\theta) ); 15
THE EFFECT OF REPLACING THE JOUKOWSKI MAP WITH THE GENERALIZED KARMAN-TREFFTZ MAP IN THE METHOD OF ZEDAN
GSJ: VOLUME 6, ISSUE 2, FEBRUARY 2018 1 GSJ: Volume 6, Issue 2, February 2018, Online: ISSN 2320-9186 THE EFFECT OF REPLACING THE JOUKOWSKI MAP WITH THE GENERALIZED KARMAN-TREFFTZ MAP IN THE METHOD OF
More informationOutput Primitives Lecture: 4. Lecture 4
Lecture 4 Circle Generating Algorithms Since the circle is a frequently used component in pictures and graphs, a procedure for generating either full circles or circular arcs is included in most graphics
More informationMAT 343 Laboratory 4 Plotting and computer animation in MATLAB
MAT 4 Laboratory 4 Plotting and computer animation in MATLAB In this laboratory session we will learn how to. Plot in MATLAB. The geometric properties of special types of matrices (rotations, dilations,
More informationNumerical Methods in Engineering Sciences
Numerical Methods in Engineering Sciences Lecture 1: Brief introduction to MATLAB Pablo Antolin pablo.antolinsanchez@unipv.it October 29th 2013 How many of you have used MATLAB before? How many of you
More informationMAT 275 Laboratory 1 Introduction to MATLAB
MATLAB sessions: Laboratory 1 1 MAT 275 Laboratory 1 Introduction to MATLAB MATLAB is a computer software commonly used in both education and industry to solve a wide range of problems. This Laboratory
More information1 Trigonometry. Copyright Cengage Learning. All rights reserved.
1 Trigonometry Copyright Cengage Learning. All rights reserved. 1.1 Radian and Degree Measure Copyright Cengage Learning. All rights reserved. Objectives Describe angles. Use radian measure. Use degree
More informationAOE 5204: Homework Assignment 4 Due: Wednesday, 9/21 in class. Extended to Friday 9/23 Online Students:
AOE 5204: Homework Assignment 4 Due: Wednesday, 9/21 in class. Extended to Friday 9/23 Online Students: Email (cdhall@vt.edu) by 5 PM Suppose F b is initially aligned with F i and at t = 0 begins to rotate
More informationMCG 4345 Aerodynamics Computational Assignment I. Report Presented to Dr. Stavros Tavoularis. Prepared By
MCG 4345 Aerodynamics Computational Assignment I Report Presented to Dr. Stavros Tavoularis Prepared By University of Ottawa November 21 st 2014 Table of Contents Table of Contents... ii List of Figures...
More informationMeasuring Lengths The First Fundamental Form
Differential Geometry Lia Vas Measuring Lengths The First Fundamental Form Patching up the Coordinate Patches. Recall that a proper coordinate patch of a surface is given by parametric equations x = (x(u,
More informationChapter 8 Complex Numbers & 3-D Plots
EGR115 Introduction to Computing for Engineers Complex Numbers & 3-D Plots from: S.J. Chapman, MATLAB Programming for Engineers, 5 th Ed. 2016 Cengage Learning Topics Introduction: Complex Numbers & 3-D
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 27 / 45
: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ) Chapter 10: Parametric Equations
More informationStrömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationE0005E - Industrial Image Analysis
E0005E - Industrial Image Analysis The Hough Transform Matthew Thurley slides by Johan Carlson 1 This Lecture The Hough transform Detection of lines Detection of other shapes (the generalized Hough transform)
More informationStudy Guide for Test 2
Study Guide for Test Math 6: Calculus October, 7. Overview Non-graphing calculators will be allowed. You will need to know the following:. Set Pieces 9 4.. Trigonometric Substitutions (Section 7.).. Partial
More informationMath 56 Homework 1. Matthew Jin. April 3, e n 10+ne. = O(n 1 ) n > n 0, where n 0 = 0 be- 10+ne
Math 56 Homework 1 Matthew Jin April 3, 2014 1a) e n 10+ne is indeed equal to big O of n 1 e as n approaches infinity. Let n n 10+ne n C n 1 for some constant C. Then ne n 10+ne C n Observe that Choose
More informationEfficacy of Numerically Approximating Pi with an N-sided Polygon
Peter Vu Brewer MAT66 Honors Topic Efficacy of umerically Approximating Pi with an -sided Polygon The quest for precisely finding the irrational number pi has been an endeavor since early human history.
More informationSolution 2. ((3)(1) (2)(1), (4 3), (4)(2) (3)(3)) = (1, 1, 1) D u (f) = (6x + 2yz, 2y + 2xz, 2xy) (0,1,1) = = 4 14
Vector and Multivariable Calculus L Marizza A Bailey Practice Trimester Final Exam Name: Problem 1. To prepare for true/false and multiple choice: Compute the following (a) (4, 3) ( 3, 2) Solution 1. (4)(
More informationEach point P in the xy-plane corresponds to an ordered pair (x, y) of real numbers called the coordinates of P.
Lecture 7, Part I: Section 1.1 Rectangular Coordinates Rectangular or Cartesian coordinate system Pythagorean theorem Distance formula Midpoint formula Lecture 7, Part II: Section 1.2 Graph of Equations
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationBasic Graphs. Dmitry Adamskiy 16 November 2011
Basic Graphs Dmitry Adamskiy adamskiy@cs.rhul.ac.uk 16 November 211 1 Plot Function plot(x,y): plots vector Y versus vector X X and Y must have the same size: X = [x1, x2 xn] and Y = [y1, y2,, yn] Broken
More informationAppendix E. Plane Geometry
Appendix E Plane Geometry A. Circle A circle is defined as a closed plane curve every point of which is equidistant from a fixed point within the curve. Figure E-1. Circle components. 1. Pi In mathematics,
More informationPractice problems from old exams for math 233
Practice problems from old exams for math 233 William H. Meeks III October 26, 2012 Disclaimer: Your instructor covers far more materials that we can possibly fit into a four/five questions exams. These
More informationModeling the Fluid Flow around Airfoils Using. Conformal Mapping
Modeling the Fluid Flow around Airfoils Using Conformal Mapping Nitin R. Kapania, Katherine Terracciano, Shannon Taylor August 29, 2008 Abstract The modeling of fluid interactions around airfoils is difficult
More informationInviscid Flows. Introduction. T. J. Craft George Begg Building, C41. The Euler Equations. 3rd Year Fluid Mechanics
Contents: Navier-Stokes equations Inviscid flows Boundary layers Transition, Reynolds averaging Mixing-length models of turbulence Turbulent kinetic energy equation One- and Two-equation models Flow management
More informationMath 7 Elementary Linear Algebra PLOTS and ROTATIONS
Spring 2007 PLOTTING LINE SEGMENTS Math 7 Elementary Linear Algebra PLOTS and ROTATIONS Example 1: Suppose you wish to use MatLab to plot a line segment connecting two points in the xy-plane. Recall that
More informationCONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION USED AS A TEST CASE FOR THE INVISCID PART OF RANS SOLVERS
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONFORMAL MAPPING POTENTIAL FLOW AROUND A WINGSECTION
More informationThe Humble Tetrahedron
The Humble Tetrahedron C. Godsalve email:seagods@hotmail.com November 4, 010 In this article, it is assumed that the reader understands Cartesian coordinates, basic vectors, trigonometry, and a bit of
More informationLinear algebra deals with matrixes: two-dimensional arrays of values. Here s a matrix: [ x + 5y + 7z 9x + 3y + 11z
Basic Linear Algebra Linear algebra deals with matrixes: two-dimensional arrays of values. Here s a matrix: [ 1 5 ] 7 9 3 11 Often matrices are used to describe in a simpler way a series of linear equations.
More informationA Brief Introduction to MATLAB
A Brief Introduction to MATLAB MATLAB (Matrix Laboratory) is an interactive software system for numerical computations and graphics. As the name suggests, MATLAB was first designed for matrix computations:
More informationReview of Trigonometry
Worksheet 8 Properties of Trigonometric Functions Section Review of Trigonometry This section reviews some of the material covered in Worksheets 8, and The reader should be familiar with the trig ratios,
More informationSupplementary Information. Design of Hierarchical Structures for Synchronized Deformations
Supplementary Information Design of Hierarchical Structures for Synchronized Deformations Hamed Seifi 1, Anooshe Rezaee Javan 1, Arash Ghaedizadeh 1, Jianhu Shen 1, Shanqing Xu 1, and Yi Min Xie 1,2,*
More information1.8 Coordinate Geometry. Copyright Cengage Learning. All rights reserved.
1.8 Coordinate Geometry Copyright Cengage Learning. All rights reserved. Objectives The Coordinate Plane The Distance and Midpoint Formulas Graphs of Equations in Two Variables Intercepts Circles Symmetry
More informationTopics in Analytic Geometry Part II
Name Chapter 9 Topics in Analytic Geometry Part II Section 9.4 Parametric Equations Objective: In this lesson you learned how to evaluate sets of parametric equations for given values of the parameter
More informationParametric Curves, Polar Plots and 2D Graphics
Parametric Curves, Polar Plots and 2D Graphics Fall 2016 In[213]:= Clear "Global`*" 2 2450notes2_fall2016.nb Parametric Equations In chapter 9, we introduced parametric equations so that we could easily
More information13.1. Functions of Several Variables. Introduction to Functions of Several Variables. Functions of Several Variables. Objectives. Example 1 Solution
13 Functions of Several Variables 13.1 Introduction to Functions of Several Variables Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. Objectives Understand
More informationAIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN ABSTRACT
AIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN Adeel Khalid *, Daniel P. Schrage + School of Aerospace Engineering, Georgia Institute of Technology
More informationPARAMETRIC EQUATIONS AND POLAR COORDINATES
10 PARAMETRIC EQUATIONS AND POLAR COORDINATES PARAMETRIC EQUATIONS & POLAR COORDINATES A coordinate system represents a point in the plane by an ordered pair of numbers called coordinates. PARAMETRIC EQUATIONS
More informationMATH 2023 Multivariable Calculus
MATH 2023 Multivariable Calculus Problem Sets Note: Problems with asterisks represent supplementary informations. You may want to read their solutions if you like, but you don t need to work on them. Set
More informationTrig Practice 09 & Nov The diagram below shows a curve with equation y = 1 + k sin x, defined for 0 x 3π.
IB Math High Level Year : Trig: Practice 09 & 0N Trig Practice 09 & Nov 0. The diagram below shows a curve with equation y = + k sin x, defined for 0 x. The point A, lies on the curve and B(a, b) is the
More informationKinematics of the Stewart Platform (Reality Check 1: page 67)
MATH 5: Computer Project # - Due on September 7, Kinematics of the Stewart Platform (Reality Check : page 7) A Stewart platform consists of six variable length struts, or prismatic joints, supporting a
More informationPractice problems from old exams for math 233 William H. Meeks III December 21, 2009
Practice problems from old exams for math 233 William H. Meeks III December 21, 2009 Disclaimer: Your instructor covers far more materials that we can possibly fit into a four/five questions exams. These
More information1 Vector Functions and Space Curves
ontents 1 Vector Functions and pace urves 2 1.1 Limits, Derivatives, and Integrals of Vector Functions...................... 2 1.2 Arc Length and urvature..................................... 2 1.3 Motion
More informationPART 1 PROGRAMMING WITH MATHLAB
PART 1 PROGRAMMING WITH MATHLAB Presenter: Dr. Zalilah Sharer 2018 School of Chemical and Energy Engineering Universiti Teknologi Malaysia 23 September 2018 Programming with MATHLAB MATLAB Environment
More informationSection 10.1 Polar Coordinates
Section 10.1 Polar Coordinates Up until now, we have always graphed using the rectangular coordinate system (also called the Cartesian coordinate system). In this section we will learn about another system,
More informationSolution Notes. COMP 151: Terms Test
Family Name:.............................. Other Names:............................. ID Number:............................... Signature.................................. Solution Notes COMP 151: Terms
More informationEXTRA-CREDIT PROBLEMS ON SURFACES, MULTIVARIABLE FUNCTIONS AND PARTIAL DERIVATIVES
EXTRA-CREDIT PROBLEMS ON SURFACES, MULTIVARIABLE FUNCTIONS AND PARTIAL DERIVATIVES A. HAVENS These problems are for extra-credit, which is counted against lost points on quizzes or WebAssign. You do not
More informationJim Lambers MAT 169 Fall Semester Lecture 33 Notes
Jim Lambers MAT 169 Fall Semester 2009-10 Lecture 33 Notes These notes correspond to Section 9.3 in the text. Polar Coordinates Throughout this course, we have denoted a point in the plane by an ordered
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationGoals: Course Unit: Describing Moving Objects Different Ways of Representing Functions Vector-valued Functions, or Parametric Curves
Block #1: Vector-Valued Functions Goals: Course Unit: Describing Moving Objects Different Ways of Representing Functions Vector-valued Functions, or Parametric Curves 1 The Calculus of Moving Objects Problem.
More informationSpline Curves. Spline Curves. Prof. Dr. Hans Hagen Algorithmic Geometry WS 2013/2014 1
Spline Curves Prof. Dr. Hans Hagen Algorithmic Geometry WS 2013/2014 1 Problem: In the previous chapter, we have seen that interpolating polynomials, especially those of high degree, tend to produce strong
More informationPhysics I : Oscillations and Waves Prof. S Bharadwaj Department of Physics & Meteorology Indian Institute of Technology, Kharagpur
Physics I : Oscillations and Waves Prof. S Bharadwaj Department of Physics & Meteorology Indian Institute of Technology, Kharagpur Lecture - 20 Diffraction - I We have been discussing interference, the
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 28 / 46
Polar Coordinates Polar Coordinates: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ)
More informationPoints Lines Connected points X-Y Scatter. X-Y Matrix Star Plot Histogram Box Plot. Bar Group Bar Stacked H-Bar Grouped H-Bar Stacked
Plotting Menu: QCExpert Plotting Module graphs offers various tools for visualization of uni- and multivariate data. Settings and options in different types of graphs allow for modifications and customizations
More informationOdysseas Skartsis ENGN2340. Project LEFM Convergence Study 12/11/13
Odysseas Skartsis ENGN2340 Project LEFM Convergence Study 2//3 Thanks to Tianyang Zhang, Davide Grazioli, and Ruike Zhao for their helpful tips on using ABAQUS. Structure of this report Main body: Appendix:
More informationSingularity Loci of Planar Parallel Manipulators with Revolute Joints
Singularity Loci of Planar Parallel Manipulators with Revolute Joints ILIAN A. BONEV AND CLÉMENT M. GOSSELIN Département de Génie Mécanique Université Laval Québec, Québec, Canada, G1K 7P4 Tel: (418) 656-3474,
More informationFind the specific function values. Complete parts (a) through (d) below. f (x,y,z) = x y y 2 + z = (Simplify your answer.) ID: 14.1.
. Find the specific function values. Complete parts (a) through (d) below. f (x,y,z) = x y y 2 + z 2 (a) f(2, 4,5) = (b) f 2,, 3 9 = (c) f 0,,0 2 (d) f(4,4,00) = = ID: 4..3 2. Given the function f(x,y)
More informationCHAPTER 2 REVIEW COORDINATE GEOMETRY MATH Warm-Up: See Solved Homework questions. 2.2 Cartesian coordinate system
CHAPTER 2 REVIEW COORDINATE GEOMETRY MATH6 2.1 Warm-Up: See Solved Homework questions 2.2 Cartesian coordinate system Coordinate axes: Two perpendicular lines that intersect at the origin O on each line.
More information12 - THREE DIMENSIONAL GEOMETRY Page 1 ( Answers at the end of all questions ) = 2. ( d ) - 3. ^i - 2. ^j c 3. ( d )
- THREE DIMENSIONAL GEOMETRY Page ( ) If the angle θ between the line x - y + x + y - z - and the plane λ x + 4 0 is such that sin θ, then the value of λ is - 4-4 [ AIEEE 00 ] ( ) If the plane ax - ay
More informationUse Parametric notation. Interpret the effect that T has on the graph as motion.
Learning Objectives Parametric Functions Lesson 3: Go Speed Racer! Level: Algebra 2 Time required: 90 minutes One of the main ideas of the previous lesson is that the control variable t does not appear
More informationGeneralized barycentric coordinates
Generalized barycentric coordinates Michael S. Floater August 20, 2012 In this lecture, we review the definitions and properties of barycentric coordinates on triangles, and study generalizations to convex,
More information6 Appendix B: Quick Guide to MATLAB R
6 Appendix B: Quick Guide to MATLAB R 6.1 Introduction In this course we will be using the software package MATLAB R. Version 17.12.0, Release R2011a has been installed in Foster 100, the labs on the third
More informationComplex Numbers, Polar Equations, and Parametric Equations. Copyright 2017, 2013, 2009 Pearson Education, Inc.
8 Complex Numbers, Polar Equations, and Parametric Equations Copyright 2017, 2013, 2009 Pearson Education, Inc. 1 8.5 Polar Equations and Graphs Polar Coordinate System Graphs of Polar Equations Conversion
More informationMATH 31A HOMEWORK 9 (DUE 12/6) PARTS (A) AND (B) SECTION 5.4. f(x) = x + 1 x 2 + 9, F (7) = 0
FROM ROGAWSKI S CALCULUS (2ND ED.) SECTION 5.4 18.) Express the antiderivative F (x) of f(x) satisfying the given initial condition as an integral. f(x) = x + 1 x 2 + 9, F (7) = 28.) Find G (1), where
More informationMATH 261 EXAM I PRACTICE PROBLEMS
MATH 261 EXAM I PRACTICE PROBLEMS These practice problems are pulled from actual midterms in previous semesters. Exam 1 typically has 6 problems on it, with no more than one problem of any given type (e.g.,
More informationLIGHT: Two-slit Interference
LIGHT: Two-slit Interference Objective: To study interference of light waves and verify the wave nature of light. Apparatus: Two red lasers (wavelength, λ = 633 nm); two orange lasers (λ = 612 nm); two
More informationCollege Pre Calculus A Period. Weekly Review Sheet # 1 Assigned: Monday, 9/9/2013 Due: Friday, 9/13/2013
College Pre Calculus A Name Period Weekly Review Sheet # 1 Assigned: Monday, 9/9/013 Due: Friday, 9/13/013 YOU MUST SHOW ALL WORK FOR EVERY QUESTION IN THE BOX BELOW AND THEN RECORD YOUR ANSWERS ON THE
More information14.5 Directional Derivatives and the Gradient Vector
14.5 Directional Derivatives and the Gradient Vector 1. Directional Derivatives. Recall z = f (x, y) and the partial derivatives f x and f y are defined as f (x 0 + h, y 0 ) f (x 0, y 0 ) f x (x 0, y 0
More informationSPC 307 Aerodynamics. Lecture 1. February 10, 2018
SPC 307 Aerodynamics Lecture 1 February 10, 2018 Sep. 18, 2016 1 Course Materials drahmednagib.com 2 COURSE OUTLINE Introduction to Aerodynamics Review on the Fundamentals of Fluid Mechanics Euler and
More informationPrime Time (Factors and Multiples)
CONFIDENCE LEVEL: Prime Time Knowledge Map for 6 th Grade Math Prime Time (Factors and Multiples). A factor is a whole numbers that is multiplied by another whole number to get a product. (Ex: x 5 = ;
More informationSUPPLEMENTARY FILE S1: 3D AIRWAY TUBE RECONSTRUCTION AND CELL-BASED MECHANICAL MODEL. RELATED TO FIGURE 1, FIGURE 7, AND STAR METHODS.
SUPPLEMENTARY FILE S1: 3D AIRWAY TUBE RECONSTRUCTION AND CELL-BASED MECHANICAL MODEL. RELATED TO FIGURE 1, FIGURE 7, AND STAR METHODS. 1. 3D AIRWAY TUBE RECONSTRUCTION. RELATED TO FIGURE 1 AND STAR METHODS
More informationLONG-TIME SIMULATIONS OF VORTEX SHEET DYNAMICS
LONG-TIME SIMULATIONS OF VORTEX SHEET DYNAMICS BY DANIEL SOLANO 1, LING XU 2, AND SILAS ALBEN 2 1 Rutgers University, New Brunswick, NJ 891, USA. 2 University of Michigan, Ann Arbor, MI 4819, USA. Abstract.
More informationMAC Learning Objectives. Module 12 Polar and Parametric Equations. Polar and Parametric Equations. There are two major topics in this module:
MAC 4 Module 2 Polar and Parametric Equations Learning Objectives Upon completing this module, you should be able to:. Use the polar coordinate system. 2. Graph polar equations. 3. Solve polar equations.
More informationAPPM 2460 PLOTTING IN MATLAB
APPM 2460 PLOTTING IN MATLAB. Introduction Matlab is great at crunching numbers, and one of the fundamental ways that we understand the output of this number-crunching is through visualization, or plots.
More informationRyian Hunter MAE 598
Setup: The initial geometry was produced using the engineering schematics provided in the project assignment document using the ANSYS DesignModeler application taking advantage of system symmetry. Fig.
More informationCalculators ARE NOT Permitted On This Portion Of The Exam 28 Questions - 55 Minutes
1 of 11 1) Give f(g(1)), given that Calculators ARE NOT Permitted On This Portion Of The Exam 28 Questions - 55 Minutes 2) Find the slope of the tangent line to the graph of f at x = 4, given that 3) Determine
More information3.1 Maxima/Minima Values
3.1 Maxima/Minima Values Ex 1: Find all critical points for the curve given by f (x)=x 5 25 3 x3 +20x 1 on the interval [-3, 2]. Identify the min and max values. We're guaranteed max and min points if
More informationClasses 7-8 (4 hours). Graphics in Matlab.
Classes 7-8 (4 hours). Graphics in Matlab. Graphics objects are displayed in a special window that opens with the command figure. At the same time, multiple windows can be opened, each one assigned a number.
More informationTransforms. COMP 575/770 Spring 2013
Transforms COMP 575/770 Spring 2013 Transforming Geometry Given any set of points S Could be a 2D shape, a 3D object A transform is a function T that modifies all points in S: T S S T v v S Different transforms
More informationTopic 4: Vectors Markscheme 4.6 Intersection of Lines and Planes Paper 2
Topic : Vectors Markscheme. Intersection of Lines and Planes Paper. Using an elimination method, x y + z x y z x y x + y 8 y Therefore x, y, z Using matrices, x y z x y z 5 (using a graphic display calculator)
More informationChapter 10 Similarity
Chapter 10 Similarity Def: The ratio of the number a to the number b is the number. A proportion is an equality between ratios. a, b, c, and d are called the first, second, third, and fourth terms. The
More information-SOLUTION- ME / ECE 739: Advanced Robotics Homework #2
ME / ECE 739: Advanced Robotics Homework #2 Due: March 5 th (Thursday) -SOLUTION- Please submit your answers to the questions and all supporting work including your Matlab scripts, and, where appropriate,
More informationThis theorem concerned a spherical gradient index in which the refractive density of a medium
APPENDIX Theorem 461 This theorem concerned a spherical gradient index in which the refractive density of a medium varied with a given power, q, of the distance, x, from a point, P. Young was more interested
More informationHOUGH TRANSFORM CS 6350 C V
HOUGH TRANSFORM CS 6350 C V HOUGH TRANSFORM The problem: Given a set of points in 2-D, find if a sub-set of these points, fall on a LINE. Hough Transform One powerful global method for detecting edges
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationBellman s Escape Problem for Convex Polygons
Bellman s Escape Problem for Convex Polygons Philip Gibbs philegibbs@gmail.com Abstract: Bellman s challenge to find the shortest path to escape from a forest of known shape is notoriously difficult. Apart
More informationAppendix E Calculating Normal Vectors
OpenGL Programming Guide (Addison-Wesley Publishing Company) Appendix E Calculating Normal Vectors This appendix describes how to calculate normal vectors for surfaces. You need to define normals to use
More informationf xx (x, y) = 6 + 6x f xy (x, y) = 0 f yy (x, y) = y In general, the quantity that we re interested in is
1. Let f(x, y) = 5 + 3x 2 + 3y 2 + 2y 3 + x 3. (a) Final all critical points of f. (b) Use the second derivatives test to classify the critical points you found in (a) as a local maximum, local minimum,
More information3D Geometry and Camera Calibration
3D Geometry and Camera Calibration 3D Coordinate Systems Right-handed vs. left-handed x x y z z y 2D Coordinate Systems 3D Geometry Basics y axis up vs. y axis down Origin at center vs. corner Will often
More information18.02 Multivariable Calculus Fall 2007
MIT OpenCourseWare http://ocw.mit.edu 18.02 Multivariable Calculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 18.02 Problem Set 4 Due Thursday
More informationECE251DN: Homework #3 Solutions
ECE251DN: Homework #3 Solutions Problem 3.7.2 (a) In this problem, we only have one null constraint. So N 1 j C = V(ψ ) = [e 2 ψ,..., 1,..., e j N 1 2 ψ ] T The weights of the least squares approximation
More information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More information0 Graphical Analysis Use of Excel
Lab 0 Graphical Analysis Use of Excel What You Need To Know: This lab is to familiarize you with the graphing ability of excels. You will be plotting data set, curve fitting and using error bars on the
More information3. The three points (2, 4, 1), (1, 2, 2) and (5, 2, 2) determine a plane. Which of the following points is in that plane?
Math 4 Practice Problems for Midterm. A unit vector that is perpendicular to both V =, 3, and W = 4,, is (a) V W (b) V W (c) 5 6 V W (d) 3 6 V W (e) 7 6 V W. In three dimensions, the graph of the equation
More information4.2 and 4.6 filled in notes.notebook. December 08, Integration. Copyright Cengage Learning. All rights reserved.
4 Integration Copyright Cengage Learning. All rights reserved. 1 4.2 Area Copyright Cengage Learning. All rights reserved. 2 Objectives Use sigma notation to write and evaluate a sum. Understand the concept
More informationPC-MATLAB PRIMER. This is intended as a guided tour through PCMATLAB. Type as you go and watch what happens.
PC-MATLAB PRIMER This is intended as a guided tour through PCMATLAB. Type as you go and watch what happens. >> 2*3 ans = 6 PCMATLAB uses several lines for the answer, but I ve edited this to save space.
More information2. A straightedge can create straight line, but can't measure. A ruler can create straight lines and measure distances.
5.1 Copies of Line Segments and Angles Answers 1. A drawing is a rough sketch and a construction is a process to create an exact and accurate geometric figure. 2. A straightedge can create straight line,
More informationTricurve Basics January 10, Tricurve Basics Timothy C. Lexen Copyright 2018
Tricurve Basics Timothy C. Lexen Copyright 2018 Figure 1 Tricurve Introduction When we think of substituting arcs for the straight sides of a triangle, we usually think of circular triangles based on an
More informationHOMEWORK ASSIGNMENT #4, MATH 253
HOMEWORK ASSIGNMENT #4, MATH 253. Prove that the following differential equations are satisfied by the given functions: (a) 2 u 2 + 2 u y 2 + 2 u z 2 =0,whereu =(x2 + y 2 + z 2 ) /2. (b) x w + y w y +
More informationPi at School. Arindama Singh Department of Mathematics Indian Institute of Technology Madras Chennai , India
Pi at School rindama Singh epartment of Mathematics Indian Institute of Technology Madras Chennai-600036, India Email: asingh@iitm.ac.in bstract: In this paper, an attempt has been made to define π by
More informationMath 265 Exam 3 Solutions
C Roettger, Fall 16 Math 265 Exam 3 Solutions Problem 1 Let D be the region inside the circle r 5 sin θ but outside the cardioid r 2 + sin θ. Find the area of D. Note that r and θ denote polar coordinates.
More information