https://community.wolfram.com RSS Feed for Wolfram Community showing any discussions tagged with Aerospace Engineering sorted by active. https://community.wolfram.com/groups/-/m/t/3502450 ![Visualizing orbital decay due to Earth's oblate shape][1] &[Wolfram Notebook][2] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=Screenshot2025-07-10at5.28.11%E2%80%AFPM.png&userId=3500732 [2]: https://www.wolframcloud.com/obj/bba46b7e-4418-4806-abf2-b10d8b8594e9 Alex Ning 2025-07-10T21:33:29Z https://community.wolfram.com/groups/-/m/t/3135874 ![enter image description here][1] &[Wolfram Notebook][2] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=Image20240306202014.jpg&userId=20103 [2]: https://www.wolframcloud.com/obj/80cb0022-4f48-468f-bf63-9d20b3313f28 Brian Beckman 2024-03-06T00:56:49Z https://community.wolfram.com/groups/-/m/t/2534173 Hi All, I'm trying to formulate and solve a pde which models thermal radiation between two bodies using a Neumann boundary condition. I've created a mesh and want to capture the thermal radiation heat transfer between element 106 and element 19 in the boundary mesh wireframe plot below. ![enter image description here][1] I've modelled convective BCs as follows: GconvPistComb = NeumannValue[hPistComb*(TPistHeat - temp[z, r]), ElementMarker == 106]; GconvBlkCool = NeumannValue[hCoolWater1*(TWater1 - temp[z, r]), ElementMarker == 75] The pde with just these included looks like: pde = {1/r D[-kMesh r D[temp[z, r], r], r] + D[-kMesh D[temp[z, r], z], z] == GconvPistComb + GconvBlkCool If I were modelling the thermal radiation from ambient gas to element 19 I would use the following Neumann BC in which TAmb is a pre-defined fixed value: GradAmbTest = NeumannValue[Epsilon Sigma (TAmb - temp[z, r])^4, ElementMarker == 19]; How can I replace TAmb with an expression that represents the solved temperature field for Element 106? [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=MathematicaThermalRadiationQuestion.png&userId=1933607 Archie Watts-Farmer 2022-05-18T18:11:52Z https://community.wolfram.com/groups/-/m/t/2484078 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/8e00e8e9-5eba-4f6b-8137-051ffacc8a96 Maher Kuzbari 2022-03-03T15:56:19Z https://community.wolfram.com/groups/-/m/t/2034161 ![Wolfram Notebook][1] &[Wolfram Notebook][2] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=IMG-2156.PNG&userId=2033168 [2]: https://www.wolframcloud.com/obj/17bd27a8-c2ca-426d-b5b9-af563c71c6d5 Anna Krzyzanska 2020-07-15T20:05:36Z https://community.wolfram.com/groups/-/m/t/2026590 ![Cassini-Huygens Probe (Artist Interpretation)][1] &[Wolfram Notebook][2] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=27159.jpg&userId=2026549 [2]: https://www.wolframcloud.com/obj/8fd0b919-bd56-4356-9059-d3c7f01e88ee Pranav Sai 2020-07-14T01:30:43Z https://community.wolfram.com/groups/-/m/t/2407750 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/685f303e-62bd-40a8-a0bd-5ed0786e5ba1 Saulo da Paz Almeida 2021-11-16T03:03:40Z https://community.wolfram.com/groups/-/m/t/2415569 Hello, I am relatively new to Wolfram. After searching appropriate solutions I want to focus here for a project of gravitational force testing. There are so many opportunities and diversity of application here. Defined by [Skylark Micro][1] rocket technology of testing facilities in reality. The purpose is to create pre-test conditions, which will be able to input materials, 3D construction and put it to a certain simulation load. For any guidance, from where to start thank you in advance. [1]: https://www.skyrora.com/skylark-micro Lariliss Liss 2021-11-30T10:31:44Z https://community.wolfram.com/groups/-/m/t/2406638 &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/3f8cce09-cf45-450e-a27f-41eab038dac5 [Original NB]: https://www.wolframcloud.com/obj/177f4cd5-e401-4b1f-84ba-905244d26dae Maher Kuzbari 2021-11-14T12:38:52Z https://community.wolfram.com/groups/-/m/t/2123140 Dear Wolfram community, a question on the definition of the quantity ‘ThermalConductivityCoefficient’ in the StandardAtmosphere Package: Provided that ‘ThermalConductivityCoefficient’ stands for the thermal conductivity lam acc. to Fourier’s law of heat conduction q = -lam grad(T), then lam should be in the order of 0.025 W/(m K) for ambient air. However, the command in the attached file yields a value of 0.000025362 J / (m s K) for condition at sea level (0 Meters). This seems to be too low by an order of 3 magnitudes. Is this quantity based on a different definition, other than assumed above ? I am using Version 11.3.0.0. ludwig.eineder 2020-11-25T12:26:27Z https://community.wolfram.com/groups/-/m/t/2011420 Has anyone in Wolfram-Land ever modeled and/or calculated the necessary parameters for spacecraft rendezvous and proximity operations? Paul Huter 2020-06-23T16:49:01Z https://community.wolfram.com/groups/-/m/t/1920848 I would like to calculate the vertical and diagonal asymptote f (x) = (x ^ 2 + x + 3) / (x-1). Will somebody help? Mateusz Janicki 2020-04-04T13:48:47Z https://community.wolfram.com/groups/-/m/t/1757394 I want to write a 'Gauss' function in my main program like the way its written in Matlab. In the following example I want to pass the variables ('aar' & 'es' are matrices and other are scalars) "aar, es, x1value, x2value, x3value, x4value, y1value, y2value, y3value, y4value" like given below. I want the output to store in detjacobs & Invdetjacobs. Kindly suggest me the correct way. [detjacobs, Invdetjacobs] = Gauss[aar, es, x1value, x2value, x3value, x4value, y1value, y2value, y3value, y4value] Do[ Do[ r = aar[[i]]; s = es[[j]]; shape1 = ((1. - r) (1. - s))/4.; shape2 = ((1. + r) (1. - s))/4.; shape3 = ((1. + r) (1. + s))/4.; shape4 = ((1. - r) (1. + s))/4.; dhdr1 = 1./4. (-1. + s); dhdr2 = (1. - s)/4.; dhdr3 = 1./4. (1. + s); dhdr4 = -(1. + s)/4.; dhds1 = 1./4. (-1. + r); dhds2 = 1./4. (-1. - r); dhds3 = (1. + r)/4.; dhds4 = (1. - r)/4.; detjacobs = {{dhdr1 x1value + dhdr2 x2value + dhdr3 x3value + dhdr4 x4value, dhdr1 y1value + dhdr2 y2value + dhdr3 y3value + dhdr4 y4value}, {dhds1 x1value + dhds2 x2value + dhds3 x3value + dhds4 x4value, dhds1 y1value + dhds2 y2value + dhds3 y3value + dhds4 y4value}}; Invdetjacobs = 1/detjacobs; ,{j,1,4}]; ,{i,1,4}]; khaja moinuddin mohammed 2019-08-09T11:59:47Z https://community.wolfram.com/groups/-/m/t/1206896 I want to rotate earth about z axis but it gives me error that object is not 3d Graphics ![ ][1] ![enter image description here][2] [1]: http://community.wolfram.com//c/portal/getImageAttachment?filename=2.PNG&userId=1199219 [2]: http://community.wolfram.com//c/portal/getImageAttachment?filename=8488Capture.PNG&userId=1199219 mohameddonia2222 2017-10-23T02:25:54Z https://community.wolfram.com/groups/-/m/t/1144481 ![enter image description here][1] Every aircraft has an elevation ceiling, which is the maximum height at which it can fly. Although most aircraft have an elevation ceiling high enough that their travel paths are not affected by geographical features, some, such as drones, fly low enough that they need to look out for mountains and other elevations. This project developed a microsite which takes two geographical locations and an elevation ceiling as input, and then outputs a map with an approximation of the shortest route around any mountains that may be in the way. Suppose we want the shortest path from LA to NYC with a ceiling of 7000 ft: loc1 = Entity["City", {"LosAngeles", "California", "UnitedStates"}]["Coordinates"]; loc2 = Entity["City", {"NewYork", "NewYork", "UnitedStates"}]["Coordinates"]; ceiling = 7000; The program first creates a "box" with the two locations at opposite corners. We then create an array of the flyable area using [`GeoElevationData`](http://reference.wolfram.com/language/ref/GeoElevationData.html?q=GeoElevationData) and our ceiling as a threshold: elevation[x1_, y1_, x2_, y2_, ceiling_] := UnitStep[QuantityMagnitude[ GeoElevationData[{{x1, y1}, {x2, y2}}, GeoZoomLevel -> 3, UnitSystem -> "Imperial"]] - ceiling] m = elevation[loc1[[1]], loc1[[2]], loc2[[1]], loc2[[2]], ceiling]; {w, l} = Dimensions[m]; From this, a network is created by connecting 'neighboring' points in our flyable area. We do this by starting out with a [`GridGraph`](http://reference.wolfram.com/language/ref/GridGraph.html), connect diagonally adjacent vertices, and then delete the invalid ones: coordTrans[w_, l_, m_] := Sort[(#2 - 1) w + w - #1 + 1 & @@@ Position[m, 1]] remove = coordTrans[w, l, m]; networkCreate[w_, l_, g_] := EdgeAdd[g, Join[Table[ If[Mod[n, w] == 1, Nothing, UndirectedEdge[n, n + w - 1]], {n, 1, w (l - 1)}], DeleteCases[ Table[UndirectedEdge[n, n + w + 1], {n, 1, w (l - 1) - 1}], Alternatives @@ Table[UndirectedEdge[n, n + w + 1], {n, w , w (l - 2), w}]]]] network = networkCreate[w, l, GridGraph[{w, l}]]; flyable = Subgraph[network, Complement[Range[l w], remove], VertexCoordinates -> AbsoluteOptions[network, VertexCoordinates][[1, 2, Complement[Range[l w], remove]]], ImageSize -> 1200] ![enter image description here][2] We then add edge weights and find the shortest path from the bottom left corner to the upper right corner: networkWeighting[UndirectedEdge[v1_, v2_]] := If[Abs[v2 - v1] == 1 || Abs[v2 - v1] == w, 1, Sqrt[2]] weights = networkWeighting /@ EdgeList[flyable]; flyable = Graph[EdgeList[flyable], EdgeWeight -> weights]; p = FindShortestPath[flyable, 1, l w]; Viola! Here's our path: coordTransBack[p_, w_, l_] := {Ceiling[#/w], l - ((Ceiling[#/w] - 1) w + l - #)} & /@ p pathCoords = coordTransBack[p, w, l]; Graphics[{Red, Line[pathCoords]}] ![enter image description here][3] Let's look at this atop a map: geoPathCoords[loc1_, loc2_, w_, l_, pathCoords_] := {(#2 - 1) (loc2[[1]] - loc1[[1]])/w + loc1[[1]], (#1 - 1) (loc2[[2]] - loc1[[2]])/l + loc1[[2]]} & @@@ pathCoords geoPath = geoPathCoords[loc1, loc2, w, l, pathCoords]; GeoGraphics[{Thick, Red, GeoPath[geoPath]}] ![enter image description here][4] There's a clear problem here. Once we're past the Rocky Mountains, we should be taking the greater circle path to our destination. Clearly we're not doing this here - this is because we were working with a discrete set of points over lat-lon space. To fix this, we can take a greedy approach and merge 2 consecutive line segments when the merge avoids the invalid regions. To do this, we express our invalid area and path as regions and use [`RegionDisjoint`](http://reference.wolfram.com/language/ref/RegionDisjoint.html): bad = ImageMesh[Image[m], Method -> "Exact", DataRange -> Reverse[Sort /@ Transpose[{loc1, loc2}]]]; linePoints[p1_, p2_] := Table[(1 - t) p1 + t p2, {t, 0, 1, 0.1`}] checkDisjoint[bad_, geoPath_, i_, j_] := RegionDisjoint[bad, Line[Reverse /@ linePoints[geoPath[[i]], geoPath[[j]]]]] i = 0; j = 1; pathFinal = {}; While[i < Length[geoPath], While[checkDisjoint[bad, geoPath, i, j] && j < Length[geoPath], j++]; AppendTo[pathFinal, geoPath[[j]]]; i = j; j = i + 1; ] The new path: GeoGraphics[{Thick, Red, GeoPath[pathFinal]}] ![enter image description here][5] Much better! A potential issue is that our solution 'hugs' the mountains too much. In the real world, this would be too dangerous. A way around this is to pad our invalid region array, say using [`Dilation`](http://reference.wolfram.com/language/ref/Dilation.html), before we build our connectivity network. This algorithm produces a good approximation of the shortest route. However, it fails when our network is disconnected, even though there could be a potential path that goes outside the box. A further continuation of this project could be done to create larger boxes to look through. [1]: http://community.wolfram.com//c/portal/getImageAttachment?filename=path.png&userId=883509 [2]: http://community.wolfram.com//c/portal/getImageAttachment?filename=network.png&userId=883509 [3]: http://community.wolfram.com//c/portal/getImageAttachment?filename=4234AdamCommunityPost.png&userId=883509 [4]: http://community.wolfram.com//c/portal/getImageAttachment?filename=6960AdamCommunityPost.png&userId=883509 [5]: http://community.wolfram.com//c/portal/getImageAttachment?filename=9170AdamCommunityPost.png&userId=883509 Adam Czarnecki 2017-07-11T15:27:41Z https://community.wolfram.com/groups/-/m/t/996238 A group of us at the International [Space Elevator Consortium (ISEC)][1] are designing a comprehensive simulation suite for the space elevator. We are trying to select a software base, and Mathematica is one of the choices. The spreadsheet shows several others. I have been using Mathematica for a few months and have filled in my scores on a scale of 1 to 10 against a list of features that we need. 1 means a lot of work is needed; 10 means it’s all there already. We would really like other users or Mathematica staff to put their own opinions in a copy of the spreadsheet. It’s OK to write in brief verbal comments or references rather than just numbers if that helps. [1]: http://www.isec.org/ John Knapman 2017-01-16T12:54:26Z https://community.wolfram.com/groups/-/m/t/935929 NASA’s Science website is a premier location for space and science cutting-edge information. The National Aeronautics and Space Administration recently chose Wolfram|Alpha to utilize the ever-growing knowledgebase for use on their Science website. NASA will use the Wolfram|Alpha API to apply this knowledge engine and offer more real-time information and data. NASA made this selection to more ably apply W|A’s natural language search abilities to enhance the user experience. Brittany Banks 2016-10-07T15:54:58Z https://community.wolfram.com/groups/-/m/t/820884 We've done a major update to the SystemModeler Examples that you can find on [www.wolfram.com/system-modeler/examples][1]. We now have two major categories: Industry Examples and Education Examples. ![enter image description here][2] When you click on one of the topics, say [Mechanical Engineering][3], you get to a page dedicated to that topic: ![Mechanical Engineering Domain][4] For pages with educational materials there are difficulty legends to guide you. There are of course examples that are suitable for more than one level. When you click through to an example, you get to the individual example page, where you can read about the example, see videos, interact with CDF dynamic content, or download the full model. The full model can also contain a Mathematica notebook with analysis and design material. ![Individual example page][5] When you find an example that you want to share on social media, there are new share buttons to make that easy. I hope you like the new pages and take a minute to browse through the examples! [1]: http://www.wolfram.com/system-modeler/examples/ [2]: http://community.wolfram.com//c/portal/getImageAttachment?filename=8188WC_ExamplesUpdate.png&userId=18621 [3]: http://www.wolfram.com/system-modeler/examples/mechanical-engineering/ [4]: http://community.wolfram.com//c/portal/getImageAttachment?filename=WC_ExamplesUpdate.png&userId=18621 [5]: http://community.wolfram.com//c/portal/getImageAttachment?filename=img3.png&userId=18621 Johan Rhodin 2016-03-10T23:40:19Z https://community.wolfram.com/groups/-/m/t/747008 I am trying to visualise contours for a simple pipe flow, I can't seem to find a notebook online. I want to write a little about the capabilities of mathematica for computational fluid dynamics as part of my masters thesis. I tried importing the code in the following documentation for the simple fluid flow at the bottom of the page, however it seems to bring errors. http://reference.wolfram.com/language/FEMDocumentation/tutorial/SolvingPDEwithFEM.html the pipe only needs to be a simple rectangle with an inlet on the left and outlet on the right. Can anyone point me in the right direction? Thank you for your time. Ahmed Hagi 2015-11-29T10:27:19Z https://community.wolfram.com/groups/-/m/t/746344 Hi guys , I have a Mathematica code that i wrote to find the natural frequencies of a triangular plate. The problem that i am facing is when i increase the number of terms above 9, i get an underdetermined system of linear equations. Any idea how to solve this problem I have attached my code below Thank you very much in advance Peter peter azizi 2015-11-27T21:48:53Z